A comparison between retrieved properties of the rain drop size distributions (DSDs) from multifrequency cloud radar observations and WRF Model results using either the Morrison or the Thompson bulk microphysics scheme is performed in order to evaluate the model’s ability to predict the rain microphysics. This comparison reveals discrepancies in the vertical profile of the rain DSDs for the stratiform region of the squall-line system observed on 12 June 2011 over Oklahoma. Based on numerical sensitivity analyses, this study addresses the bias at the top of the rain layer and the vertical evolution of the DSD properties (i.e., of Dm and N0*). In this way, the Thompson scheme is used to explore the sensitivity to the melting process. Moreover, using the Thompson and Morrison schemes, the sensitivity of the DSD vertical evolution to different breakup and self-collection parameterizations is studied. Results show that the DSDs are strongly dependent on the representation of the melting process in the Thompson scheme. In the Morrison scheme, the simulations with more efficient breakup reproduce the DSD properties with better fidelity. This study highlights how the inaccuracies in simulated Dm and N0* for both microphysics schemes can impact the evaporation rate, which is systematically underestimated in the model.

Abstract. For the last 25 years, CO-PDD (Cézeaux-Aulnat-Opme-puy de Dôme) has evolved to become a full instrumented platform for atmospheric research. It has received credentials as a national observing platform in France and is internationally recognized as a global station in the GAW (Global Atmosphere Watch) network. It is a reference site of European and national research infrastructures ACTRIS (Aerosol Cloud and Trace gases Research Infrastructure) and ICOS (Integrated Carbon Observing System). The site located on top of the puy de Dôme mountain (1465 m a.s.l.) is completed by additional sites located at lower altitudes and adding the vertical dimension to the atmospheric observations: Opme (660 m a.s.l.), Cézeaux (410 m), and Aulnat (330 m). The integration of different sites offers a unique combination of in situ and remote sensing measurements capturing and documenting the variability of particulate and gaseous atmospheric composition, but also the optical, biochemical, and physical properties of aerosol particles, clouds, and precipitations. Given its location far away from any major emission sources, its altitude, and the mountain orography, the puy de Dôme station is ideally located to sample different air masses in the boundary layer or in the free troposphere depending on time of day and seasons. It is also an ideal place to study cloud properties with frequent presence of clouds at the top in fall and winter. As a result of the natural conditions prevailing at the site and of the very exhaustive instrumental deployment, scientific studies at the puy de Dôme strongly contribute to improving knowledge in atmospheric sciences, including the characterization of trends and variability, the understanding of complex and interconnected processes (microphysical, chemical, biological, chemical and dynamical), and the provision of reference information for climate/chemistry models. In this context, CO-PDD is a pilot site to conduct instrumental development inside its wind tunnel for testing liquid and ice cloud probes in natural conditions, or in situ systems to collect aerosol and cloud. This paper reviews 25 years (1995–2020) of atmospheric observation at the station and related scientific research contributing to atmospheric and climate science.

This study investigates how multifrequency cloud radar observations can be used to evaluate the representation of rain microphysics in the WRF Model using two bulk microphysics schemes. A squall line observed over Oklahoma on 12 June 2011 is used as a case study. A recently developed retrieval technique combining observations of two vertically pointing cloud radars provides quantitative description of the drop size distribution (DSD) properties of the transition and stratiform regions of the squall-line system. For the first time, the results of this multifrequency cloud radar retrieval are compared to more conventional retrievals from a nearby polarimetric radar, and a supplementary result of this work is that this new methodology provides a much more detailed description of the DSD vertical and temporal variations. While the extent and evolution of the squall line is well reproduced by the model, the 1-h low-reflectivity transition region is not. In the stratiform region, simulations with both schemes are able to reproduce the observed downdraft and the associated significative subsaturation below the melting level, but with a slight overestimation of the relative humidity. Under this subsaturated air, the simulated rain mixing ratio continuously decreases toward the ground, in agreement with the observations. Conversely, the profiles of the mean volume diameter and the concentration parameter of the DSDs are not well reproduced. These discrepancies pinpoint at an issue in the representation of rain microphysics. The companion paper, investigates the sources of the biases in the microphysics processes in the rain layer by performing numerical sensitivity studies.