Experimental Site For Atmospheric Research Research Articles

Measurement report: Evaluation of the TOF-ACSM-CV for PM1.0 and PM2.5 measurements during the RITA-2021 field campaign

Abstract. The recently developed time-of-flight aerosol chemical speciation monitor with a capture vaporizer and a PM2.5 aerodynamic lens (TOF-ACSM-CV-PM2.5) aims to improve the collection efficiency and chemical characterization of aerosol particles with a diameter smaller than 2.5 µm. In this study, comprehensive cross-comparisons were performed between real-time online measurements and offline filter analysis with 24 h collection time. The goal was to evaluate the capabilities of the TOF-ACSM-CV-PM2.5 lens, as well as the accuracy of the TOF-ACSM-CV-PM2.5. The experiments were conducted at Cabauw Experimental Site for Atmospheric Research (CESAR) during the RITA-2021 campaign. The non-refractory fine particulate matter (PM1.0 and PM2.5) was measured by two collocated TOF-ACSM-CV-PM2.5 instruments by placing them behind a PM2.5 and PM1.0 inlet, respectively. A comparison between the ACSMs and PM2.5 and PM1.0 filter samples showed a much better accuracy than ±30 % less given in the previous reports, with average differences less than ±10 % for all inorganic chemical species. In addition, the ACSMs were compared to the Monitoring Instrument for Aerosol and Gas (MARGA) (slope between 0.78 and 0.97 for inorganic compounds, R2≥ 0.93) and a mobility particle size spectrometer (MPSS), measuring the particle size distribution from around 10 to 800 nm (slope was around 1.00, R2= 0.91). The intercomparison of the online measurements and the comparison between the online and offline measurements indicated a low bias (< 10 % for inorganic compounds) and demonstrated the high accuracy and stability of the TOF-ACSM-CV-PM2.5 lens for the atmospheric observations of particle matter. The two ACSMs exhibited an excellent agreement, with differences less than 7 %, which allowed a quantitative estimate of PM1.0 vs. PM2.5 chemical composition. The result showed that the PM1.0 accounted for about 70 %–80 % of the PM2.5 on average. The NO3 mass fraction increased, but the organic carbon (OC) mass fraction decreased from PM1.0 to PM2.5, indicating the size dependence on chemical composition.

Tropospheric and stratospheric ozone profiles during the 2019 TROpomi vaLIdation eXperiment (TROLIX-19)

Abstract. A TROPOspheric Monitoring Instrument (TROPOMI) validation campaign was held in the Netherlands based at the CESAR (Cabauw Experimental Site for Atmospheric Research) observatory during September 2019. The TROpomi vaLIdation eXperiment (TROLIX-19) consisted of active and passive remote sensing platforms in conjunction with several balloon-borne and surface chemical (e.g., ozone and nitrogen dioxide) measurements. The goal of this joint NASA-KNMI geophysical validation campaign was to make intensive observations in the TROPOMI domain in order to be able to establish the quality of the L2 satellite data products under realistic conditions, such as non-idealized conditions with varying cloud cover and a range of atmospheric conditions at a rural site. The research presented here focuses on using ozone lidars from NASA's Goddard Space Flight Center to better evaluate the characterization of ozone throughout TROLIX-19. Results of comparisons to the lidar systems with balloon, space-borne and ground-based passive measurements are shown. In addition, results are compared to a global coupled chemistry meteorology model to illustrate the vertical variability and columnar amounts of both tropospheric and stratospheric ozone during the campaign period.

Fifty Years of Atmospheric Boundary-Layer Research at Cabauw Serving Weather, Air Quality and Climate

An overview is given of 50-year Cabauw observations and research on the structure and dynamics of the atmospheric boundary layer. It is shown that over time this research site with its 200-m meteorological tower has grown into an atmospheric observatory with a comprehensive observational program encompassing almost all aspects of the atmospheric column including its boundary conditions. This is accomplished by the Cabauw Experimental Site for Atmospheric Research (CESAR) a consortium of research institutes. CESAR plays an important role in the educational programs of the CESAR universities. The current boundary-layer observational program is described in detail, and other parts of the CESAR observational program discussed more briefly. Due to an open data policy the CESAR datasets are used by researchers all over the world. Examples are given of the use of the long time series for model evaluation, satellite validation, and process studies. The role of tall towers is discussed in relation to the development of more and better ground-based remote sensing techniques. CESAR is now incorporated into the Ruisdael observatory, the large-scale atmospheric research infrastructure in the Netherlands. With Ruisdael the embedding of the Dutch atmospheric community in national policy landscape, and in the European atmospheric research infrastructures is assured for the coming decade.

Observing and Modelling the Surface Radiative Budget and Cloud Radiative Forcing at the Cabauw Experimental Site for Atmospheric Research (CESAR), the Netherlands, 2009–17

AbstractA dataset of 9 years in duration (2009–17) of clouds and radiation was obtained at the Cabauw Experimental Site for Atmospheric Research (CESAR) in the Netherlands. Cloud radiative forcings (CRF) were derived from the dataset and related to cloud cover and temperature. Also, the data were compared with RCM output. Results indicate that there is a seasonal cycle (i.e., winter, spring, summer, and autumn) in longwave (CRF-LW: 48.3, 34.4, 30.8, and 38.7 W m−2) and shortwave (CRF-SW: −23.6, −60.9, −67.8, and −32.9 W m−2) forcings at CESAR. Total CRF is positive in winter and negative in summer. The RCM has a cold bias with respect to the observations, but the model CRF-LW corresponds well to the observed CRF-LW as a result of compensating errors in the difference function that makes up the CRF-LW. The absolute value of model CRF-SW is smaller than the observed CRF-SW in summer, mostly because of albedo differences. The majority of clouds from above 2 km are present at the same time as low clouds, so the higher clouds have only a small impact on CRF whereas low clouds dominate their values. CRF-LW is a function of fractional cloudiness. CRF-SW is also a function of fractional cloudiness, if the values are normalized by the cosine of solar zenith angle. Expressions for CRF-LW and CRF-SW were derived as functions of temperature, fractional cloudiness, and solar zenith angle, indicating that CRF is the largest when fractional cloudiness is the highest but is also large for low temperature and high sun angle.

A Three‐Dimensional Array for the Study of Infrasound Propagation Through the Atmospheric Boundary Layer

AbstractThe Royal Netherlands Meteorological Institute (KNMI) operates a three‐dimensional microbarometer array at the Cabauw Experimental Site for Atmospheric Research observatory. The array consists of five microbarometers on a meteorological tower up to an altitude of 200 m. Ten ground‐based microbarometers surround the tower with an array aperture of 800 m. This unique setup allows for the study of infrasound propagation in three dimensions. The added value of the vertical dimension is the sensitivity to wind and temperature in the atmospheric boundary layer over multiple altitudes. In this study, we analyze infrasound generated by an accidental chemical explosion at the Moerdijk petrochemical plant on 3 June 2014. The recordings of the tower microbarometers show two sequential arrivals, whereas the recordings on the ground show one wavefront. This arrival structure is interpreted to be the upgoing and downgoing wavefronts. The observations are compared with propagation modeling results using global‐scale and mesoscale atmospheric models. Independent temperature and wind measurements, which are available at the Cabauw Experimental Site for Atmospheric Research, are used for comparison with model output. The modeling results explain the signal arrival times; however, the tower wavefront arrivals are not explained. This study is important for understanding the influence of the atmospheric boundary layer on infrasound detections and propagation.

Missed Fog?

Conventional in situ observations of meteorological variables are restricted to a limited number of levels near the surface, with the lowest observation often made around 1-m height. This can result in missed observations of both shallow fog, and the initial growth stage of thicker fog layers. At the same time, numerical experiments have demonstrated the need for high vertical grid resolution in the near-surface layer to accurately simulate the onset of fog; this requires correspondingly high-resolution observational data for validation. A two-week field campaign was conducted in November 2017 at the Cabauw Experimental Site for Atmospheric Research (CESAR) in the Netherlands. The aim was to observe the growth of shallow fog layers and assess the possibility of obtaining very high-resolution observations near the surface during fog events. Temperature and relative humidity were measured at centimetre resolution in the lowest 7 m using distributed temperature sensing. Further, a novel approach was employed to estimate visibility in the lowest 2.5 m using a camera and an extended light source. These observations were supplemented by the existing conventional sensors at the site, including those along a 200-m tall tower. Comparison between the increased-resolution observations and their conventional counterparts show the errors to be small, giving confidence in the reliability of the techniques. The increased resolution of the observations subsequently allows for detailed investigations of fog growth and evolution. This includes the observation of large temperature inversions in the lowest metre (up to 5 K) and corresponding regions of (super)saturation where the fog formed. Throughout the two-week observation period, fog was observed twice at the conventional sensor height of 2.0 m. Two additional low-visibility events were observed in the lowest 0–0.5 m using the camera-based observations, but were missed by the conventional sensors. The camera observations also showed the growth of shallow radiation fog, forming in the lowest 0.5 m as early as two hours before it was observed at the conventional height of 2 m.

The transition of gravity wave or mesoscale flow to boundary layer turbulence and implications for infrasound propagation conditions

Knowledge of the atmospheric boundary layer (ABL) is important for weather and climate forecasting as well as infrasound detection. In spite of many advances, the ABL is still not represented realistically in weather models, which leads to uncertainties in weather predictions. At the same time, the influence of the ABL on infrasound detectability is not well understood, which potentially limits the capability of the technique. In this work, a network of anemometers, microbarometers, and ceilometers is used to characterize gravity wave activity and turbulence in the ABL. The sensors are co-located at the Cabauw Experimental Site for Atmospheric Research (CESAR). The sensor network allows for the estimation of velocity and pressure spectra, as well as boundary-layer structure and mixing height. The goal of this work is to study the transition from gravity wave or mesoscale flow to boundary layer turbulence and to consider implications for infrasound propagation conditions. Here, an analysis of spectra under convective and stable conditions is presented with respect to both gravity wave and boundary-layer turbulence theory. Observed distinctions in the velocity and pressure spectral transitions will be discussed, including large-scale turbulence and the mesoscale spectral gap.Knowledge of the atmospheric boundary layer (ABL) is important for weather and climate forecasting as well as infrasound detection. In spite of many advances, the ABL is still not represented realistically in weather models, which leads to uncertainties in weather predictions. At the same time, the influence of the ABL on infrasound detectability is not well understood, which potentially limits the capability of the technique. In this work, a network of anemometers, microbarometers, and ceilometers is used to characterize gravity wave activity and turbulence in the ABL. The sensors are co-located at the Cabauw Experimental Site for Atmospheric Research (CESAR). The sensor network allows for the estimation of velocity and pressure spectra, as well as boundary-layer structure and mixing height. The goal of this work is to study the transition from gravity wave or mesoscale flow to boundary layer turbulence and to consider implications for infrasound propagation conditions. Here, an analysis of spectra under...

Understanding and Reducing False Alarms in Observational Fog Prediction

The reduction in visibility that accompanies fog events presents a hazard to human safety and navigation. However, accurate fog prediction remains elusive, with numerical methods often unable to capture the conditions of fog formation, and observational methods having high false-alarm rates in order to obtain high hit rates of prediction. In this work, 5 years of observations from the Cabauw Experimental Site for Atmospheric Research are used to further investigate how false alarms may be reduced using the statistical method for diagnosing radiation-fog events from observations developed by Menut et al. (Boundary-Layer Meteorol 150:277–297, 2014). The method is assessed for forecast lead times of 1–6 h and implementing four optimization schemes to tune the prediction for different needs, compromising between confidence and risk. Prediction scores improve significantly with decreased lead time, with the possibility of achieving a hit rate of over 90% and a false-alarm rate of just 13%. In total, a further 31 combinations of predictive variables beyond the original combination are explored (including mostly, e.g., variables related to moisture and static stability of the boundary layer). Little change to the prediction scores indicates any appropriate combination of variables that measure saturation, turbulence, and near-surface cooling can be used. The remaining false-alarm periods are manually assessed, identifying the lack of spatio–temporal information (such as the temporal evolution of the local conditions and the advective history of the airmass) as the ultimate limiting factor in the methodology’s predictive capabilities. Future observational studies are recommended that investigate the near-surface evolution of fog and the role of non-local heterogeneity on fog formation.

Analysis and validation of Weather Research and Forecasting model tendencies for meso-to-microscale modelling of the atmospheric boundary layer

The different terms of momentum and energy budget simulated by the Weather Research and Forecasting model (WRF) are examined at two sites that depict diverse climate and orographic environments. The analysis focuses both on particular weather episodes of one-to-few days duration and in one-year statistics useful for the assessment of the annual wind climate of a particular location. In order to estimate their accuracy, a validation exercise is performed in the Cabauw Experimental Site for Atmospheric Research (CESAR) in the Netherlands, where a fair agreement is found when comparing the mean surface pressure gradient term extracted from WRF with the pressure-derived geostrophic winds observed by the CESAR project. This study is carried out under the New European Wind Atlas (NEWA) project as part of the effort to develop a model-chain system that connects mesoscale with microscale wind flow models with different fidelity ranges for wind resource assessment applications. Therefore the analysis deepens in the horizontal and vertical variability of the tendencies in each of the sites with the ultimate goal of providing some guidelines of their best usage as large-scale forcing of the microscale models.

Critical scales to explain urban hydrological response: an application in Cranbrook, London

Abstract. Rainfall variability in space and time, in relation to catchment characteristics and model complexity, plays an important role in explaining the sensitivity of hydrological response in urban areas. In this work we present a new approach to classify rainfall variability in space and time and we use this classification to investigate rainfall aggregation effects on urban hydrological response. Nine rainfall events, measured with a dual polarimetric X-Band radar instrument at the CAESAR site (Cabauw Experimental Site for Atmospheric Research, NL), were aggregated in time and space in order to obtain different resolution combinations. The aim of this work was to investigate the influence that rainfall and catchment scales have on hydrological response in urban areas. Three dimensionless scaling factors were introduced to investigate the interactions between rainfall and catchment scale and rainfall input resolution in relation to the performance of the model. Results showed that (1) rainfall classification based on cluster identification well represents the storm core, (2) aggregation effects are stronger for rainfall than flow, (3) model complexity does not have a strong influence compared to catchment and rainfall scales for this case study, and (4) scaling factors allow the adequate rainfall resolution to be selected to obtain a given level of accuracy in the calculation of hydrological response.

Impacts of afternoon and evening sea‐breeze fronts on local turbulence, and on CO2 and radon‐222 transport

We investigated sharp disruptions of local turbulence and scalar transport due to the arrival of sea‐breeze fronts (SBFs). To this end, we employed a comprehensive 10‐year observational database from the Cabauw Experimental Site for Atmospheric Research (CESAR, the Netherlands). Sea‐breeze (SB) days were selected using a five‐filter algorithm, which accounts for large‐scale conditions and a clear mesoscale‐frontal signal associated with the land–sea contrast. Among those days (102 in all, 8.3%), based on the value of the sensible‐heat flux at the onset of SB, we identified three atmospheric boundary‐layer (ABL) regimes: convective, transition and stable. In the convective regime, the thermally driven convective boundary layer is only slightly altered by a small enhancement of the shear when the SBF arrives. Regarding the transition regime, we found that the ABL afternoon transition is accelerated. This was quantified by estimating the contributions of shear and buoyancy to the turbulent kinetic energy. Other relevant disruptions are the sharp reduction in ABL depth (∼250 m/hr) and the sudden increase in average wind speed (> 2 m/s). In the stable regime, the arrival of the SB leads to disturbances in the wind profile at the surface layer. We observed a deviation of more than 1 m/s in the observed surface‐layer wind profile compared with the profile calculated using Monin–Obukhov Similarity Theory (MOST). Our findings furthermore reveal the determinant role of the SB direction in the transport of water vapour, CO2 and 222Rn. The return of continental air masses driven by the SB circulation generates sharp CO2 increases (up to 14 ppm in half an hour) in a few SB events. We suggest that the variability in 222Rn evolution may also be influenced by other non‐local processes such as the large‐scale footprint from more remote sources.

A Novel Radar-Based Visibility Estimator

Remote visibility (Vis) estimation by radar is of interest to aviation, road traffic, and other fields. Millimeter-wave radars are suitable candidates because of such advantages as high spatial resolution and sensitivity to small droplets in reflection and attenuation. To investigate remote Vis estimation and to develop physics-based models and algorithms, a 35 GHz cloud radar at the Cabauw Experimental Site for Atmospheric Research (CESAR) in the western part of the Netherlands has acquired data during fog periods in “fog mode.” The advantage in using millimeter-waves for remote sensing of fog is that they interact strong enough with fog for sensing it, but not too strong so that they can penetrate the fog allowing to sense fog top. Simultaneously, fog drop size distribution (DSD) and Vis are continuously automatically measured by the in situ optical sensors at CESAR. Radar reflectivity (Z) and Vis can be linked theoretically since they are related to the sixth and the second moments of an assumed Gamma-shaped DSD. However, in reality the fog DSD is not always Gamma-shaped, leading to errors in the Vis estimation. A further development of the Vis-Z model includes the attenuation factor (La), which is proportional to the liquid water content at a given radar wavelength. This improves the estimated accuracy of Vis, in the theoretical moments-based model. Finally, we were able to arrive at a higher accuracy by introducing an empirical exponential model, estimating Vis from Z and La. A test based on DSD data sets for various fog types in the literature showed robust performance of the Vis-Z-La model for large variations in DSD. The Vis-Z-La model is also validated with the actual fog data sets that were collected by the in situ and remote sensing instruments synergy at CESAR.

Ambient and laboratory observations of organic ammonium salts in PM1.

Ambient measurements of PM1 aerosol chemical composition at Cabauw, the Netherlands, implicate higher ammonium concentrations than explained by the formation of inorganic ammonium salts. This additional particulate ammonium is called excess ammonium (eNH4). Height profiles over the Cabauw Experimental Site for Atmospheric Research (CESAR) tower, of combined ground based and airborne aerosol mass spectrometric (AMS) measurements on a Zeppelin airship show higher concentrations of eNH4 at higher altitudes compared to the ground. Through flights across the Netherlands, the Zeppelin based measurements furthermore substantiate eNH4 as a regional phenomenon in the planetary boundary layer. The excess ammonium correlates with mass spectral signatures of (di-)carboxylic acids, making a heterogeneous acid-base reaction the likely process of NH3 uptake. We show that this excess ammonium was neutralized by the organic fraction forming particulate organic ammonium salts. We discuss the significance of such organic ammonium salts for atmospheric aerosols and suggest that NH3 emission control will have benefits for particulate matter control beyond the reduction of inorganic ammonium salts.

In situ observations of the isotopic composition of methane at the Cabauw tall tower site

Abstract. High-precision analyses of the isotopic composition of methane in ambient air can potentially be used to discriminate between different source categories. Due to the complexity of isotope ratio measurements, such analyses have generally been performed in the laboratory on air samples collected in the field. This poses a limitation on the temporal resolution at which the isotopic composition can be monitored with reasonable logistical effort. Here we present the performance of a dual isotope ratio mass spectrometric system (IRMS) and a quantum cascade laser absorption spectroscopy (QCLAS)-based technique for in situ analysis of the isotopic composition of methane under field conditions. Both systems were deployed at the Cabauw Experimental Site for Atmospheric Research (CESAR) in the Netherlands and performed in situ, high-frequency (approx. hourly) measurements for a period of more than 5 months. The IRMS and QCLAS instruments were in excellent agreement with a slight systematic offset of (+0.25 ± 0.04) ‰ for δ13C and (−4.3 ± 0.4) ‰ for δD. This was corrected for, yielding a combined dataset with more than 2500 measurements of both δ13C and δD. The high-precision and high-temporal-resolution dataset not only reveals the overwhelming contribution of isotopically depleted agricultural CH4 emissions from ruminants at the Cabauw site but also allows the identification of specific events with elevated contributions from more enriched sources such as natural gas and landfills. The final dataset was compared to model calculations using the global model TM5 and the mesoscale model FLEXPART-COSMO. The results of both models agree better with the measurements when the TNO-MACC emission inventory is used in the models than when the EDGAR inventory is used. This suggests that high-resolution isotope measurements have the potential to further constrain the methane budget when they are performed at multiple sites that are representative for the entire European domain.

Aerosol source apportionment from 1-year measurements at the CESAR tower in Cabauw, the Netherlands

Abstract. Intensive measurements of submicron aerosol particles and their chemical composition were performed with an Aerosol Chemical Speciation Monitor (ACSM) at the Cabauw Experimental Site for Atmospheric Research (CESAR) in Cabauw, the Netherlands, sampling at 5 m height above ground. The campaign lasted nearly 1 year from July 2012 to June 2013 as part of the EU-FP7-ACTRIS project (Q-ACSM Network). Including equivalent black carbon an average particulate mass concentration of 9.50 µg m−3 was obtained during the whole campaign with dominant contributions from ammonium nitrate (45 %), organic aerosol (OA, 29 %), and ammonium sulfate (19 %). There were 12 exceedances of the World Health Organization (WHO) PM2.5 daily mean limit (25 µg m−3) observed at this rural site using PM1 instrumentation only. Ammonium nitrate and OA represented the largest contributors to total particulate matter during periods of exceedance. Source apportionment of OA was performed season-wise by positive matrix factorization (PMF) using the multilinear engine 2 (ME-2) controlled via the source finder (SoFi). Primary organic aerosols were attributed mainly to traffic (8–16 % contribution to total OA, averaged season-wise) and biomass burning (0–23 %). Secondary organic aerosols (SOAs, 61–84 %) dominated the organic fraction during the whole campaign, particularly on days with high mass loadings. A SOA factor which is attributed to humic-like substances (HULIS) was identified as a highly oxidized background aerosol in Cabauw. This shows the importance of atmospheric aging processes for aerosol concentration at this rural site. Due to the large secondary fraction, the reduction of particulate mass at this rural site is challenging on a local scale.

Forecasting radiation fog at climatologically contrasting sites: evaluation of statistical methods and WRF

A six‐year climatology of radiation fog has been compiled at two sites: the Research Centre for the Lower Atmosphere (CIBA, Spain) and the Cabauw Experimental Site for Atmospheric Research (CESAR, The Netherlands). These sites are contrasted in terms of geographical situation, climate zone, altitude, humidity and soil water availability. Therefore, several climatological differences in fog abundance, onset, dissipation and duration have been quantified between the two sites. The more humid site (CESAR) is characterised by relatively short radiation fog events distributed throughout the year. However, radiation fog at the drier site (CIBA) is more persistent and appears during late autumn/winter months. In general, its formation requires more time after sunset (∼2 h more), since further cooling is required to reach saturation. The forecast of these fog events has been evaluated through two different approaches. First, we extend the statistical method presented by Menut et al. (2014) (M14). This method uses statistics to define threshold values on key variables for fog formation (pre‐fog) and verifies its predictability using observations and numerical model output. We present some of the most appropriate threshold values for the forecasting of pre‐fog periods at both sites, which differ from those presented in M14 and depend on the optimisation of the hit rate or the false‐alarm rate. Additionally, we also extend M14 by suggesting other variables as potential predictors for fog formation (friction velocity and visibility tendency). Finally, we focus on fog simulation by the Weather Research and Forecasting (WRF) model in terms of liquid water content. The WRF model was able to simulate radiation fog when configured with sophisticated physical options and high resolution. However it failed in simulating the onset, dissipation and the vertical extent of fog (which was overestimated). The model results were extremely sensitive to the spin‐up time.

Role of synoptic- and meso-scales on the evolution of the boundary-layer wind profile over a coastal region: the near-coast diurnal acceleration

The contributions of synoptic- and meso-scales to the boundary layer wind profile evolution in a coastal environment are examined. The analysis is based on observations of the wind profile within the first 200 m of the atmosphere continuously recorded during a 10 year period (2001–2010) at the 213-m meteorological tower at the Cabauw Experimental Site for Atmospheric Research (CESAR, The Netherlands). The analysis is supported by a numerical experiment based on the Weather Research and Forecasting (WRF) model performed at high horizontal resolution of 2 km and spanning the complete observational period (10 years). Results indicate that WRF is able to reproduce the inter-annual wind variability but with a tendency to be too geostrophic. At seasonal scales, we find a differentiated behavior between Winter and Summer seasons with the Spring and Autumn transition periods more similar to the Summer and Winter modes, respectively. The winter momentum budget shows a weak intradiurnal variability. The synoptic scale controls the shape of the near surface wind profile that is characterized by weaker and more ageostrophic winds near the surface than at higher altitudes within the planetary boundary layer (PBL) as a result of the frictional turning. In turn, during summer, mesoscale circulations associated with the differential heating of land and sea become important. As a result, the PBL winds show a stronger intradiurnal component that is characterized by an oscillation of the near surface winds around the geostrophic direction with the maximum departure in the afternoon. Although also driven by thermal land-sea differences, this mesoscale component is not associated with the classical concept of a sea-breeze front. It originates from the thermal expansion of the boundary layer over land and primarily differs from the sea-breeze in its propagation speed resulting in a wind rotation far ahead of any coastal front. We refer to it as the near-coast diurnal acceleration (NCDA). The contribution of the NCDA depends on the specific orientation of the coast (NE-SW at CESAR). Our findings stress the importance of evaluating and understanding the performance of mesoscale models with multi-year observational/simulated data sets in order to provide a statistically robust characterization of the limitations of surface layer and boundary-layer parameterizations and thus compensate for the scarceness of upper level wind observations.

Study of aerosol hygroscopic events over the Cabauw experimental site for atmospheric research (CESAR) using the multi-wavelength Raman lidar Caeli

This article presents a study of aerosol optical and microphysical properties under different relative humidity (RH) but well mixed layer conditions using optical and microphysical aerosol properties from multi-wavelength (MW) Raman lidar and in-situ aerosol observations collected at the Cabauw Experimental Site for Atmospheric Research (CESAR). Two hygroscopic events are described through 3 backscatter (β) and 2 extinction (α) coefficients which in turn provide intensive parameters such as the backscatter-related Ångström exponent (åβ) and the lidar ratio (LR). Along with it, profiles of RH were inferred from Raman lidar observations and therefore, as a result of varying humidity conditions, a shift on the aerosol optical properties can be described. Thus, it is observed that as RH increases, aerosols uptake water vapour, augment their size and consequently the åβ diminishes whereas the LR increases. The enhancement factor based on the backscatter coefficient at 532 nm, which characterizes the aerosol from hygroscopic standpoint, is also estimated. Finally, microphysical properties that are necessary for aerosol radiative forcing estimates – such as volume, effective radii, refractive index and size distribution, all vertically resolved – are retrieved using the inversion with regularization. Using this method, two hygroscopic events are described in detail.

Impact of spatial and temporal resolution of rainfall inputs on urban hydrodynamic modelling outputs: A multi-catchment investigation

Urban catchments are typically characterised by high spatial variability and fast runoff processes resulting in short response times. Hydrological analysis of such catchments requires high resolution precipitation and catchment information to properly represent catchment response. This study investigated the impact of rainfall input resolution on the outputs of detailed hydrodynamic models of seven urban catchments in North-West Europe. The aim was to identify critical rainfall resolutions for urban catchments to properly characterise catchment response. Nine storm events measured by a dual-polarimetric X-band weather radar, located in the Cabauw Experimental Site for Atmospheric Research (CESAR) of the Netherlands, were selected for analysis. Based on the original radar estimates, at 100 m and 1 min resolutions, 15 different combinations of coarser spatial and temporal resolutions, up to 3000 m and 10 min, were generated. These estimates were then applied to the operational semi-distributed hydrodynamic models of the urban catchments, all of which have similar size (between 3 and 8 km2), but different morphological, hydrological and hydraulic characteristics. When doing so, methodologies for standardising model outputs and making results comparable were implemented. Results were analysed in the light of storm and catchment characteristics. Three main features were observed in the results: (1) the impact of rainfall input resolution decreases rapidly as catchment drainage area increases; (2) in general, variations in temporal resolution of rainfall inputs affect hydrodynamic modelling results more strongly than variations in spatial resolution; (3) there is a strong interaction between the spatial and temporal resolution of rainfall input estimates. Based upon these results, methods to quantify the impact of rainfall input resolution as a function of catchment size and spatial–temporal characteristics of storms are proposed and discussed.

Numerical simulation of the interaction between ammonium nitrate aerosol and convective boundary-layer dynamics

We investigate the interaction between the ammonium nitrate aerosol (NAO3) abundance and convective boundary-layer (CBL) dynamics by means of a large-eddy simulation (LES) framework. In our LES model the CBL dynamics is solved coupled with radiation, chemistry, and surface exchange. Concerning the aerosol coupling we assume a simplified representation that accounts for black carbon, aerosol water and inorganic aerosols, focusing on the semi-volatile ammonium nitrate aerosol within the CBL. The aerosol absorption and scattering of shortwave radiation is also taken into consideration. We use a data set of observations taken at the Cabauw Experimental Site for Atmospheric Research during the IMPACT/EUCAARI (European Integrated Project on Aerosol, Cloud, Climate, and Air Quality Interactions) campaign to successfully evaluate our LES approach. We highlight that our LES framework reproduces the observations of the ratio between gas-phase nitrate and total nitrate at the surface, with a diurnally-averaged overestimation of only ≈12%. We show that the dependence between gas-aerosol conversion of nitrate and CBL (thermo)dynamics produces highly non-linear concentration and turbulent flux vertical profiles. The flux profiles maximize at around 1/3 of the CBL. Close to the surface, we show that the outgassing of NAO3 affects the dry deposition of nitrate. This outgassing is responsible for the high deposition velocities obtained from the concentration and flux measurements during observational campaigns. To account for the influence of CBL (thermo)dynamics on gas-aerosol conversion we propose an effective turbulent exchange coefficient based on an analysis of the flux budget equation of aerosol nitrate calculated by our LES. The implementation of this effective turbulent exchange coefficient in a 1D model leads to a better agreement with the LES results and with surface observations.