Despite the apparent simplicity, it is notoriously difficult to measure rainfall accurately because of the challenging environment within which it is measured. Systematic bias caused by wind is inherent in rainfall measurement and introduces an inconvenient unknown into hydrological science that is generally ignored. This paper examines the role of rain gauge shape and mounting height on catch efficiency (CE), where CE is defined as the ratio between nonreference and reference rainfall measurements. Using a pit gauge as a reference, we have demonstrated that rainfall measurements from an exposed upland site, recorded by an adjacent conventional cylinder rain gauge mounted at 0.5 m, were underestimated by more than 23% on average. At an exposed lowland site, with lower wind speeds on average, the equivalent mean undercatch was 9.4% for an equivalent gauge pairing. An improved‐aerodynamic gauge shape enhanced CE when compared to a conventional cylinder gauge shape. For an improved‐aerodynamic gauge mounted at 0.5 m above the ground, the mean undercatch was 11.2% at the upland site and 3.4% at the lowland site. The mounting height of a rain gauge above the ground also affected CE due to the vertical wind gradient near to the ground. Identical rain gauges mounted at 0.5 and 1.5 m were compared at an upland site, resulting in a mean undercatch of 11.2% and 17.5%, respectively. By selecting three large rainfall events and splitting them into shorter‐duration intervals, a relationship explaining 81% of the variance was established between CE and wind speed.

Recent studies have used numerical models to estimate the collection efficiency of solid precipitation gauges when exposed to the wind in both shielded and unshielded configurations. The models used computational fluid dynamics (CFD) simulations of the airflow pattern generated by the aerodynamic response to the gauge-shield geometry. These are used as initial conditions to perform Lagrangian tracking of solid precipitation particles. Validation of the results against field observations yielded similarities in the overall behavior, but the model output only approximately reproduced the dependence of the experimental collection efficiency on wind speed. This paper presents an improved snowflake trajectory modeling scheme due to the inclusion of a dynamically determined drag coefficient. The drag coefficient was estimated using the local Reynolds number as derived from CFD simulations within a time-independent Reynoldsaveraged Navier-Stokes approach. The proposed dynamic model greatly improves the consistency of results with the field observations recently obtained at the Marshall Field winter precipitation test bed in Boulder, Colorado.

The use of windshields to reduce the impact of wind on snow measurements is common. This paper investigates the catching performance of shielded and unshielded gauges using numerical simulations. In Part II, the role of the windshield and gauge aerodynamics, as well as the varying flow field due to the turbulence generated by the shield-gauge configuration, in reducing the catch efficiency is investigated. This builds on the computational fluid dynamics results obtained in Part I, where the airflow patterns in the proximity of an unshielded and single Alter shielded Geonor T-200B gauge are obtained using both time-independent [Reynolds-averaged Navier-Stokes (RANS)] and time-dependent [large-eddy simulation (LES)] approaches. A Lagrangian trajectory model is used to track different types of snowflakes (wet and dry snow) and to assess the variation of the resulting gauge catching performance with the wind speed. The collection efficiency obtained with the LES approach is generally lower than the one obtained with the RANS approach. This is because of the impact of the LES-resolved turbulence above the gauge orifice rim. The comparison between the collection efficiency values obtained in case of shielded and unshielded gauge validates the choice of installing a single Alter shield in a windy environment. However, time-dependent simulations show that the propagating turbulent structures produced by the aerodynamic response of the upwind single Alter blades have an impact on the collection efficiency. Comparison with field observations provides the validation background for the model results.

The airflow surrounding any catching‐type rain gauge when impacted by wind is deformed by the presence of the gauge body, resulting in the acceleration of wind above the orifice of the gauge, which deflects raindrops and snowflakes away from the collector (the wind‐induced undercatch). The method of mounting a gauge with the collector at or below the level of the ground, or the use of windshields to mitigate this effect, is often not practicable. The physical shape of a gauge has a significant impact on its collection efficiency. In this study, we show that appropriate “aerodynamic” shapes are able to reduce the deformation of the airflow, which can reduce undercatch. We have employed computational fluid‐dynamic simulations to evaluate the time‐averaged airflow realized around “aerodynamic” rain gauge shapes when impacted by wind. Terms of comparison are provided by the results obtained for two standard “conventional” rain gauge shapes. The simulations have been run for different wind speeds and are based on a time‐averaged Reynolds‐Averaged Navier‐Stokes model. The shape of the aerodynamic gauges is shown to have a positive impact on the time‐averaged airflow patterns observed around the orifice compared to the conventional shapes. Furthermore, the turbulent air velocity fields for the aerodynamic shapes present “recirculating” structures, which may improve the particle‐catching capabilities of the gauge collector.

The aerodynamic response of snow gauges when exposed to the wind is responsible for a significant reduction of their collection performance. The modifications induced by the gauge and the windshield onto the space-time patterns of the undisturbed airflow deviate the snowflake trajectories. In Part I, the disturbed air velocity field in the vicinity of shielded and unshielded gauge configurations is investigated. In Part II, the airflow is the basis for a particle tracking model of snowflake trajectories to estimate the collection efficiency. A Geonor T-200B gauge inside a single Alter shield is simulated for wind speeds varying from 1 to 8 m s 21. Both time-averaged and time-dependent computational fluid dynamics simulations are performed, based on Reynolds-averaged Navier-Stokes (RANS) and large-eddy simulation (LES) models, respectively. A shear stress tensor k-V model (where k is the turbulent kinetic energy and V is the turbulent specific dissipation rate) is used for the RANS formulation and solved within a finite-volume method. The LES is implemented with a Smagorinsky subgrid-scale method that models the subgrid stresses as a gradient-diffusion process. The RANS simulations confirm the attenuation of the airflow velocity above the gauge when using a single Alter shield, but the generated turbulence above the orifice rim is underestimated. The intensity and spatial extension of the LES-resolved turbulent region show a dependency on the wind speed that was not detected by the RANS. The time-dependent analysis showed the propagation of turbulent structures and the impact on the turbulent kinetic energy above the gauge collecting section.