-
The spatio-temporal variability of integrated water vapour (IWV) on small scales of less than 10 km and hours is assessed with data from the 2 months of the High Definition Clouds and Precipitation for advancing Climate Prediction (HD(CP)<sup>2</sup>) Observational Prototype Experiment (HOPE). The statistical intercomparison of the unique set of observations during HOPE (microwave radiometer (MWR), Global Positioning System (GPS), sun photometer, radiosondes, Raman lidar, infrared and near-infrared Moderate Resolution Imaging Spectroradiometer (MODIS) on the satellites Aqua and Terra) measuring close together reveals a good agreement in terms of random differences (standard deviation ≤1 kg m<sup>−2</sup>) and correlation coefficient (≥ 0.98). The exception is MODIS, which appears to suffer from insufficient cloud filtering. <br><br> For a case study during HOPE featuring a typical boundary layer development, the IWV variability in time and space on scales of less than 10 km and less than 1 h is investigated in detail. For this purpose, the measurements are complemented by simulations with the novel ICOsahedral Nonhydrostatic modelling framework (ICON), which for this study has a horizontal resolution of 156 m. These runs show that differences in space of 3–4 km or time of 10–15 min induce IWV variabilities on the order of 0.4 kg m<sup>−2</sup>. This model finding is confirmed by observed time series from two MWRs approximately 3 km apart with a comparable temporal resolution of a few seconds. <br><br> Standard deviations of IWV derived from MWR measurements reveal a high variability (> 1 kg m<sup>−2</sup>) even at very short time scales of a few minutes. These cannot be captured by the temporally lower-resolved instruments and by operational numerical weather prediction models such as COSMO-DE (an application of the Consortium for Small-scale Modelling covering Germany) of Deutscher Wetterdienst, which is included in the comparison. However, for time scales larger than 1 h, a sampling resolution of 15 min is sufficient to capture the mean standard deviation of IWV. The present study shows that instrument sampling plays a major role when climatological information, in particular the mean diurnal cycle of IWV, is determined.
-
We use a global coupled chemistry–climate–land model (CESM) to assess the integrated effect of climate, emissions and land use changes on annual surface O<sub>3</sub> and PM<sub>2.5</sub> in the United States with a focus on national parks (NPs) and wilderness areas, using the RCP4.5 and RCP8.5 projections. We show that, when stringent domestic emission controls are applied, air quality is predicted to improve across the US, except surface O<sub>3</sub> over the western and central US under RCP8.5 conditions, where rising background ozone counteracts domestic emission reductions. Under the RCP4.5 scenario, surface O<sub>3</sub> is substantially reduced (about 5 ppb), with daily maximum 8 h averages below the primary US Environmental Protection Agency (EPA) National Ambient Air Quality Standards (NAAQS) of 75 ppb (and even 65 ppb) in all the NPs. PM<sub>2.5</sub> is significantly reduced in both scenarios (4 μg m<sup>−3</sup>; ~50%), with levels below the annual US EPA NAAQS of 12 μg m<sup>−3</sup> across all the NPs; visibility is also improved (10–15 dv; >75 km in visibility range), although some western US parks with Class I status (40–74 % of total sites in the US) are still above the 2050 planned target level to reach the goal of natural visibility conditions by 2064. We estimate that climate-driven increases in fire activity may dominate summertime PM<sub>2.5</sub> over the western US, potentially offsetting the large PM<sub>2.5</sub> reductions from domestic emission controls, and keeping visibility at present-day levels in many parks. Our study indicates that anthropogenic emission patterns will be important for air quality in 2050. However, climate and land use changes alone may lead to a substantial increase in surface O<sub>3</sub> (2–3 ppb) with important consequences for O<sub>3</sub> air quality and ecosystem degradation at the US NPs. Our study illustrates the need to consider the effects of changes in climate, vegetation, and fires in future air quality management and planning and emission policy making.
-
The HadGEM2 earth system climate model was used to assess the impact of biomass burning on surface ozone concentrations over the Amazon forest and its impact on vegetation, under present-day climate conditions. Here we consider biomass burning emissions from wildfires, deforestation fires, agricultural forest burning, and residential and commercial combustion. Simulated surface ozone concentration is evaluated against observations taken at two sites in the Brazilian Amazon forest for years 2010 to 2012. The model is able to reproduce the observed diurnal cycle of surface ozone mixing ratio at the two sites, but overestimates the magnitude of the monthly averaged hourly measurements by 5–15 ppb for each available month at one of the sites. We vary biomass burning emissions over South America by ±20, 40, 60, 80 and 100% to quantify the modelled impact of biomass burning on surface ozone concentrations and ozone damage on vegetation productivity over the Amazon forest. We used the ozone damage scheme in the "high" sensitivity mode to give an upper limit for this effect. Decreasing South American biomass burning emissions by 100% (i.e. to zero) reduces surface ozone concentrations (by about 15 ppb during the biomass burning season) and suggests a 15% increase in monthly mean net primary productivity averaged over the Amazon forest, with local increases up to 60%. The simulated impact of ozone damage from present-day biomass burning on vegetation productivity is about 230 TgC yr<sup>−1</sup>. Taking into account that uncertainty in these estimates is substantial, this ozone damage impact over the Amazon forest is of the same order of magnitude as the release of carbon dioxide due to fire in South America; in effect it potentially doubles the impact of biomass burning on the carbon cycle.
-
Given significant challenges with available measurements of aerosol acidity, proxy methods are frequently used to estimate the acidity of atmospheric particles. In this study, four of the most common aerosol acidity proxies are evaluated and compared: (1) the ion balance method, (2) the molar ratio method, (3) thermodynamic equilibrium models, and (4) the phase partitioning of ammonia. All methods are evaluated against predictions of thermodynamic models and against direct observations of aerosol–gas equilibrium partitioning acquired in Mexico City during the Megacity Initiative: Local and Global Research Objectives (MILAGRO) study. The ion balance and molar ratio methods assume that any deficit in inorganic cations relative to anions is due to the presence of H<sup>+</sup> and that a higher H<sup>+</sup> loading and lower cation / anion ratio both correspond to increasingly acidic particles (i.e., lower pH). Based on the MILAGRO measurements, no correlation is observed between H<sup>+</sup> levels inferred with the ion balance and aerosol pH predicted by the thermodynamic models and NH<sub>3</sub>–NH<sub>4</sub><sup>+</sup> partitioning. Similarly, no relationship is observed between the cation / anion molar ratio and predicted aerosol pH. Using only measured aerosol chemical composition as inputs without any constraint for the gas phase, the E-AIM (Extended Aerosol Inorganics Model) and ISORROPIA-II thermodynamic equilibrium models tend to predict aerosol pH levels that are inconsistent with the observed NH<sub>3</sub>–NH<sub>4</sub><sup>+</sup> partitioning. The modeled pH values from both E-AIM and ISORROPIA-II run with gas + aerosol inputs agreed well with the aerosol pH predicted by the phase partitioning of ammonia. It appears that (1) thermodynamic models constrained by gas + aerosol measurements and (2) the phase partitioning of ammonia provide the best available predictions of aerosol pH. Furthermore, neither the ion balance nor the molar ratio can be used as surrogates for aerosol pH, and previously published studies with conclusions based on such acidity proxies may need to be reevaluated. Given the significance of acidity for chemical processes in the atmosphere, the implications of this study are important and far reaching.
-
Although ozone is an atmospheric gas with high spatial and temporal variability, mesoscale numerical weather prediction (NWP) models simplify the specification of ozone concentrations used in their shortwave schemes by using a few ozone profiles. In this paper, a two-part study is presented: (i) an evaluation of the quality of the ozone profiles provided for use with the shortwave schemes in the Advanced Research version of the Weather Research and Forecasting (WRF-ARW) model and (ii) an assessment of the impact of deficiencies in those profiles on the performance of model simulations of direct solar radiation. The first part compares simplified data sets used to specify the total ozone column in six schemes (i.e., Goddard, New Goddard, RRTMG, CAM, GFDL and Fu–Liou–Gu) with the Multi-Sensor Reanalysis data set during the period 1979–2008 examining the latitudinal, longitudinal and seasonal limitations in the ozone profile specifications of each parameterization. The results indicate that the maximum deviations are over the poles and show prominent longitudinal patterns in the departures due to the lack of representation of the patterns associated with the Brewer–Dobson circulation and the quasi-stationary features forced by the land–sea distribution, respectively. In the second part, the bias in the simulated direct solar radiation due to these deviations from the simplified spatial and temporal representation of the ozone distribution is analyzed for the New Goddard and CAM schemes using the Beer–Lambert–Bouguer law and for the GFDL using empirical equations. For radiative applications those simplifications introduce spatial and temporal biases with near-zero departures over the tropics throughout the year and increasing poleward with a maximum in the high middle latitudes during the winter of each hemisphere.
-
We measured total dissolved reactive nitrogen in precipitation samples collected at Uljin, a Korean coastal site upwind of the southern East Sea/Sea of Japan (EJS), selected as a representative study site of atmospheric deposition over the northwestern Pacific margin. NO<sub>3</sub><sup>-</sup> was found to be the most abundant nitrogen species, followed by NH<sub>4</sub><sup>+</sup> and dissolved organic nitrogen (DON). Air-mass back-trajectory (AMBT) analysis revealed that a significant fraction of the inorganic nitrogen (NO<sub>3</sub><sup>-</sup> and NH<sub>4</sub><sup>+</sup>) originated from mainland Asia, whereas the DON was primarily derived from Korea. Using varimax-rotated factor analysis in combination with major ions as tracers, agricultural activities in Korea were identified as the primary sources of DON in these samples. In addition, a positive correlation was found at Uljin between the size of organic fraction in total reactive nitrogen and nitrogen to carbon atomic ratio in organic matter. This correlation has also been observed at other locations worldwide, implying the utilization potential of atmospheric organic nitrogen might increase with its proportion in total nitrogen. Combining wet deposition measurements in this study with literature values for dry deposition observed at a remote island in the EJS, the total atmospheric depositional flux of reactive nitrogen was estimated to be 115 mmol N m<sup>−2</sup> yr<sup>−1</sup> over the southern EJS. Our study sheds new light on the potentially significant contribution to primary productivity of the northwestern Pacific Ocean by atmospheric deposition of nitrogen, especially the organic fraction.
-
Sources and propagation characteristics of high-frequency gravity waves observed in the mesosphere using airglow emissions from Gadanki (13.5° N, 79.2° E) and Hyderabad (17.5° N, 78.5° E) are investigated using reverse ray tracing. Wave amplitudes are also traced back, including both radiative and diffusive damping. The ray tracing is performed using background temperature and wind data obtained from the MSISE-90 and HWM-07 models, respectively. For the Gadanki region, the suitability of these models is tested. Further, a climatological model of the background atmosphere for the Gadanki region has been developed using nearly 30 years of observations available from a variety of ground-based (MST radar, radiosondes, MF radar) and rocket- and satellite-borne measurements. ERA-Interim products are utilized for constructing background parameters corresponding to the meteorological conditions of the observations. With the reverse ray-tracing method, the source locations for nine wave events could be identified to be in the upper troposphere, whereas for five other events the waves terminated in the mesosphere itself. Uncertainty in locating the terminal points of wave events in the horizontal direction is estimated to be within 50–100 km and 150–300 km for Gadanki and Hyderabad wave events, respectively. This uncertainty arises mainly due to non-consideration of the day-to-day variability in the tidal amplitudes. Prevailing conditions at the terminal points for each of the 14 events are provided. As no convection in and around the terminal points is noticed, convection is unlikely to be the source. Interestingly, large (~9 m s<sup>−1</sup>km<sup>−1</sup>) vertical shears in the horizontal wind are noticed near the ray terminal points (at 10–12 km altitude) and are thus identified to be the source for generating the observed high-phase-speed, high-frequency gravity waves.
-
Trends in air quality across the Northern Hemisphere over a 21-year period (1990–2010) were simulated using the Community Multiscale Air Quality (CMAQ) multiscale chemical transport model driven by meteorology from Weather Research and Forecasting (WRF) simulations and internally consistent historical emission inventories obtained from EDGAR. Thorough comparison with several ground observation networks mostly over Europe and North America was conducted to evaluate the model performance as well as the ability of CMAQ to reproduce the observed trends in air quality over the past 2 decades in three regions: eastern China, the continental United States and Europe. <br><br> The model successfully reproduced the observed decreasing trends in SO<sub>2</sub>, NO<sub>2</sub>, 8 h O<sub>3</sub> maxima, SO<sub>4</sub><sup>2−</sup> and elemental carbon (EC) in the US and Europe. However, the model fails to reproduce the decreasing trends in NO<sub>3</sub><sup>−</sup> in the US, potentially pointing to uncertainties of NH<sub>3</sub> emissions. The model failed to capture the 6-year trends of SO<sub>2</sub> and NO<sub>2</sub> in CN-API (China – Air Pollution Index) from 2005 to 2010, but reproduced the observed pattern of O<sub>3</sub> trends shown in three World Data Centre for Greenhouse Gases (WDCGG) sites over eastern Asia. Due to the coarse spatial resolution employed in these calculations, predicted SO<sub>2</sub> and NO<sub>2</sub> concentrations are underestimated relative to all urban networks, i.e., US-AQS (US – Air Quality System; normalized mean bias (NMB) = −38% and −48%), EU-AIRBASE (European Air quality data Base; NMB = −18 and −54%) and CN-API (NMB = −36 and −68%). Conversely, at the rural network EU-EMEP (European Monitoring and Evaluation Programme), SO<sub>2</sub> is overestimated (NMB from 4 to 150%) while NO<sub>2</sub> is simulated well (NMB within ±15%) in all seasons. Correlations between simulated and observed O<sub>3</sub> wintertime daily 8 h maxima (DM8) are poor compared to other seasons for all networks. Better correlation between simulated and observed SO<sub>4</sub><sup>2−</sup> was found compared to that for SO<sub>2</sub>. Underestimation of summer SO<sub>4</sub><sup>2−</sup> in the US may be associated with the uncertainty in precipitation and associated wet scavenging representation in the model. The model exhibits worse performance for NO<sub>3</sub><sup>−</sup> predictions, particularly in summer, due to high uncertainties in the gas/particle partitioning of NO<sub>3</sub><sup>−</sup> as well as seasonal variations of NH<sub>3</sub> emissions. There are high correlations (<i>R</i> > 0.5) between observed and simulated EC, although the model underestimates the EC concentration by 65% due to the coarse grid resolution as well as uncertainties in the PM speciation profile associated with EC emissions. <br><br> The almost linear response seen in the trajectory of modeled O<sub>3</sub> changes in eastern China over the past 2 decades suggests that control strategies that focus on combined control of NO<sub>x</sub> and volatile organic compound (VOC) emissions with a ratio of 0.46 may provide the most effective means for O<sub>3</sub> reductions for the region devoid of nonlinear response potentially associated with NO<sub>x</sub> or VOC limitation resulting from alternate strategies. The response of O<sub>3</sub> is more sensitive to changes in NO<sub>x</sub> emissions in the eastern US because the relative abundance of biogenic VOC emissions tends to reduce the effectiveness of VOC controls. Increasing NH<sub>3</sub> levels offset the relative effectiveness of NO<sub>x</sub> controls in reducing the relative fraction of aerosol NO<sub>3</sub><sup>−</sup> formed from declining NO<sub>x</sub> emissions in the eastern US, while the control effectiveness was assured by the simultaneous control of NH<sub>3</sub> emission in Europe.
-
The temperature measurements of the rotational Raman lidar of the University of Hohenheim (UHOH RRL) during the High Definition of Clouds and Precipitation for advancing Climate Prediction (HD(CP)<sup>2</sup>) Observation Prototype Experiment (HOPE) in April and May 2013 are discussed. The lidar consists of a frequency-tripled Nd:YAG laser at 355 nm with 10 W average power at 50 Hz, a two-mirror scanner, a 40 cm receiving telescope, and a highly efficient polychromator with cascading interference filters for separating four signals: the elastic backscatter signal, two rotational Raman signals with different temperature dependence, and the vibrational Raman signal of water vapor. The main measurement variable of the UHOH RRL is temperature. For the HOPE campaign, the lidar receiver was optimized for high and low background levels, with a novel switch for the passband of the second rotational Raman channel. The instrument delivers atmospheric profiles of water vapor mixing ratio as well as particle backscatter coefficient and particle extinction coefficient as further products. As examples for the measurement performance, measurements of the temperature gradient and water vapor mixing ratio revealing the development of the atmospheric boundary layer within 25 h are presented. As expected from simulations, a reduction of the measurement uncertainty of 70% during nighttime was achieved with the new low-background setting. A two-mirror scanner allows for measurements in different directions. When pointing the scanner to low elevation, measurements close to the ground become possible which are otherwise impossible due to the non-total overlap of laser beam and receiving telescope field of view in the near range. An example of a low-level temperature measurement is presented which resolves the temperature gradient at the top of the stable nighttime boundary layer 100 m above the ground.
-
East Asia has experienced rapid development with increasing carbon monoxide (CO) emission in the past decades. Therefore, uplifting CO from the boundary layer to the free troposphere in East Asia can have great implications on regional air quality around the world. It can also influence global climate due to the longer lifetime of CO at higher altitudes. In this study, three cases of high CO episodes in the East China Sea and the Sea of Japan from 2003 to 2005 are examined with spaceborne Measurements of Pollution in the Troposphere (MOPITT) data, in combination with aircraft measurements from the Measurement of Ozone and Water Vapor by Airbus In-Service Aircraft (MOZAIC) program. High CO abundances of 300–550 ppbv are observed in MOZAIC data in the free troposphere during these episodes. These are among the highest CO abundances documented at these altitudes. On average, such episodes with CO over 400 ppbv (in the 2003 and 2004 cases) and between 200 and 300 ppbv (in the 2005 case) may occur 2–5 and 10–20% in time, respectively, in the respective altitudes over the region. Correspondingly, elevated CO is shown in MOPITT daytime data in the middle to upper troposphere in the 2003 case, in the lower to middle troposphere in the 2004 case, and in the upper troposphere in the 2005 case. Through analyses of the simulations from a chemical transport model GEOS-Chem and a trajectory dispersion model FLEXPART, we found different CO signatures in the elevated CO and distinct transport pathways and mechanisms for these cases. In the 2003 case, emissions from large forest fires near Lake Baikal dominated the elevated CO, which had been rapidly transported upward by a frontal system from the fire plumes. In the 2004 case, anthropogenic CO from the North China Plain experienced frontal lifting and mostly reached ~ 700 hPa near the East China Sea, while CO from biomass burning over Indochina experienced orographic lifting, lee-side-trough-induced convection, and frontal lifting through two separate transport pathways, leading to two distinct CO enhancements around 700 and 300 hPa. In the 2005 case, the observed CO of ~ 300 ppbv around 300 hPa originated from anthropogenic sources over the Sichuan Basin and the North China Plain and from forest fires over Indochina. The high CO was transported to such altitudes through strong frontal lifting, interacting with convection and orographic lifting. These cases show that topography affects vertical transport of CO in East Asia via different ways, including orographic uplifting over the Hengduan Mountains, assisting frontal lifting in the North China Plain, and facilitating convection in the Sichuan Basin. In particular, topography-induced lee-side troughs over Indochina led to strong convection that assisted CO uplifting to the upper troposphere. This study shows that the new daytime MOPITT near-infrared (NIR) and thermal-infrared (TIR) data (version 5 or above) have enhanced vertical sensitivity in the free troposphere and may help qualitative diagnosis of vertical transport processes in East Asia.