The focus of the present study is to assess and quantify the uncertainty in predicting the steady and unsteady aerodynamic performance as well as the major mechanical characteristics of a contrarotating turbofan, primarily due to geometric variations stemming from the manufacturing process. The basis of this study is the optically scanned blisk of the first rotor, for which geometric variations from blade to blade are considered. In a first step, selected profile sections of the first rotor were evaluated aerodynamically by applying the 2D coupled Euler/boundary-layer solver mises. Statistical properties of the relevant flow quantities were calculated firstly based on the results of the nine manufactured blades. In a second step, the geometric variations were decomposed into their corresponding eigenforms by means of principal component analysis (PCA). These modes were the basis for carrying out Monte Carlo (MC) simulations in order to analyze in detail the blade's aerodynamic response to the prescribed geometric variations. By means of 3D-computational fluid dynamics (CFD) simulations of the entire fan stage for all the nine scanned rotor 1 blade geometries, the variation of the overall stage performance parameters will be quantified. The impact of the instrumentation will be discussed, here partly doubling the standard deviation of the major performance indicators for the instrumented blades and also triggering a premature laminar/turbulent transition of the boundary layer. In terms of the unsteady blade row interaction, the standard deviation of the resulting blade pressure amplitude shall be discussed based on unsteady simulations, taking advantage of a novel harmonic balance approach. It will be shown that the major uncertainty in terms of the predicted blade pressure amplitude is in the aft part of the front rotor and results from upstream shock/blade interaction. Apart from the aerodynamic performance, an analysis of the mechanical properties in terms of Campbell characteristics and eigenfrequencies was carried out for each of the scanned blades of rotor 1, reflecting the frequency scattering of each eigenmode due to geometric variability.