GMD Publications for 2019

Click on page icon to view abstract.

Andrews, Elisabeth, Patrick J. Sheridan, John A. Ogren, Derek Hageman, Anne Jefferson, Jim Wendell, Andrés Alástuey, Lucas Alados-Arboledas, Michael Bergin, Marina Ealo, A. Gannet Hallar, András Hoffer, Ivo Kalapov, Melita Keywood, Jeongeun Kim, Sang-Woo Kim, Felicia Kolonjari, Casper Labuschagne, Neng-Huei Lin, AnneMarie Macdonald, Olga L. Mayol-Bracero, Ian B. McCubbin, Marco Pandolfi, Fabienne Reisen, Sangeeta Sharma, James P. Sherman, Mar Sorribas and Junying Sun, (2019), Overview of the NOAA/ESRL Federated Aerosol Network, Bulletin of the American Meteorological Society, 100, 1, 123-135, 10.1175/BAMS-D-17-0175.1


To estimate global aerosol radiative forcing, measurements of aerosol optical properties are made by the National Oceanic and Atmospheric Administration (NOAA) Earth System Research Laboratory (ESRL)’s Global Monitoring Division (GMD) and their collaborators at 30 monitoring locations around the world. Many of the sites are located in regions influenced by specific aerosol types (Asian and Saharan desert dust, Asian pollution, biomass burning, etc.). This network of monitoring stations is a shared endeavor of NOAA and many collaborating organizations, including the World Meteorological Organization (WMO)’s Global Atmosphere Watch (GAW) program, the U.S. Department of Energy (DOE), several U.S. and foreign universities, and foreign science organizations. The result is a long-term cooperative program making atmospheric measurements that are directly comparable with those from all the other network stations and with shared data access. The protocols and software developed to support the program facilitate participation in GAW’s atmospheric observation strategy, and the sites in the NOAA/ESRL network make up a substantial subset of the GAW aerosol observations. This paper describes the history of the NOAA/ESRL Federated Aerosol Network, details about measurements and operations, and some recent findings from the network measurements.

Flores, Edgar, Joële Viallon, Tiphaine Choteau, Philippe Moussay, Faraz Idrees, Robert I Wielgosz, Jeongsoon Lee, Ewelina Zalewska, Gerard Nieuwenkamp, Adriaan van der Veen, L A Konopelko, Y A Kustikov, A V Kolobova, Y K Chubchenko, O V Efremova, BI Zhe, Zeyi Zhou, Zeyi Zhou, George C Rhoderick, Joseph T Hodge, Takuya Shimosaka, Nobuyuki Aoki, Brad Hall, Paul Brewer, Dariusz Cieciora, Michela Sega, Tatiana Macé, Judit Fükő, Zsófia Nagyné Szilágyi, Tamás Büki, Mudalo I Jozela, Napo G Ntsasa, Nompumelelo Leshabane, James Tshilongo, Prabha Johri and Tanil Tarhan, (2019), CCQM-K120 (Carbon dioxide at background and urban level), Metrologia, 56, 1A, 08001-08001, 10.1088/0026-1394/56/1A/08001


CCQM-K120.a comparison involves preparing standards of carbon dioxide in air which are fit for purpose for the atmospheric monitoring community, with stringent requirements on matrix composition and measurement uncertainty of the CO2 mole fraction. This represents an analytical challenge and is therefore considered as a Track C comparison. The comparison will underpin CMC claims for CO2 in air for standards and calibrations services for the atmospheric monitoring community, matrix matched to real air, over the mole fraction range of 250 μmol/mol to 520 μmol/mol.

CCQM-K120.b comparison tests core skills and competencies required in gravimetric preparation, analytical certification and purity analysis. It is considered as a Track A comparison. It will underpin CO2 in air and nitrogen claims in a mole fraction range starting at the smallest participant's reported expanded uncertainty and ending at 500 mmol/mol. Participants successful in this comparison may use their result in the flexible scheme and underpin claims for all core mixtures

This study has involved a comparison at the BIPM of a suite of 44 gas standards prepared by each of the participating laboratories. Fourteen laboratories took part in both comparisons (CCQM-K120.a, CCQM-K120.b) and just one solely in the CCQM-K120.b comparison.

The standards were sent to the BIPM where the comparison measurements were performed. Two measurement methods were used to compare the standards, to ensure no measurement method dependant bias: GC-FID and FTIR spectroscopic analysis corrected for isotopic variation in the CO2 gases, measured at the BIPM using absorption laser spectroscopy. Following the advice of the CCQM Gas Analysis Working Group, results from the FTIR method were used to calculate the key comparison reference values.

Hall, Bradley D., Andrew M. Crotwell, Benjamin R. Miller, Michael Schibig and James W. Elkins, (2019), Gravimetrically prepared carbon dioxide standards in support of atmospheric research, Atmospheric Measurement Techniques, 12, 1, 517-524, 10.5194/amt-12-517-2019


Abstract. We have explored a one-step method for gravimetric preparation of CO2-in-air standards in aluminum cylinders. We consider both adsorption to stainless steel surfaces used in the transfer of highly pure CO2 and adsorption of CO2 to cylinder walls. We demonstrate that CO2-in-air standards can be prepared with relatively low uncertainty (∼0.04%, ∼95% confidence level) by introducing aliquots whose masses are known to high precision and by using well-characterized cylinders. Five gravimetric standards, prepared over the nominal range of 350 to 490µmol mol−1 (parts per million,ppm), showed excellent internal consistency, with residuals from a linear fit equal to 0.05ppm. This work compliments efforts to maintain the World Meteorological Organization, Global Atmosphere Watch, mole fraction scale for carbon dioxide in air, widely used for atmospheric monitoring. This gravimetric technique could be extended to other atmospheric trace gases, depending on the vapor pressure of the gas.

Hossaini, Ryan, Elliot Atlas, Sandip S. Dhomse, Martyn P. Chipperfield, Peter F. Bernath, Anton M. Fernando, Jens Mühle, Amber A. Leeson, Stephen A. Montzka, Wuhu Feng, Jeremy J. Harrison, Paul Krummel, Martin K. Vollmer, Stefan Reimann, Simon O'Doherty, Dickon Young, Michela Maione, Jgor Arduini and Chris R. Lunder, (2019), Recent Trends in Stratospheric Chlorine From Very Short‐Lived Substances, Journal of Geophysical Research: Atmospheres, 124, 4, 2318-2335, 10.1029/2018JD029400


Very short‐lived substances (VSLS), including dichloromethane (CH2Cl2), chloroform (CHCl3), perchloroethylene (C2Cl4), and 1,2‐dichloroethane (C2H4Cl2), are a stratospheric chlorine source and therefore contribute to ozone depletion. We quantify stratospheric chlorine trends from these VSLS (VSLCltot) using a chemical transport model and atmospheric measurements, including novel high‐altitude aircraft data from the NASA VIRGAS (2015) and POSIDON (2016) missions. We estimate VSLCltot increased from 69 (±14) parts per trillion (ppt) Cl in 2000 to 111 (±22) ppt Cl in 2017, with >80% delivered to the stratosphere through source gas injection, and the remainder from product gases. The modeled evolution of chlorine source gas injection agrees well with historical aircraft data, which corroborate reported surface CH2Cl2 increases since the mid‐2000s. The relative contribution of VSLS to total stratospheric chlorine increased from ~2% in 2000 to ~3.4% in 2017, reflecting both VSLS growth and decreases in long‐lived halocarbons. We derive a mean VSLCltot growth rate of 3.8 (±0.3) ppt Cl/year between 2004 and 2017, though year‐to‐year growth rates are variable and were small or negative in the period 2015–2017. Whether this is a transient effect, or longer‐term stabilization, requires monitoring. In the upper stratosphere, the modeled rate of HCl decline (2004–2017) is −5.2% per decade with VSLS included, in good agreement to ACE satellite data (−4.8% per decade), and 15% slower than a model simulation without VSLS. Thus, VSLS have offset a portion of stratospheric chlorine reductions since the mid‐2000s.

Li, Pingyang, Jens Mühle, Stephen A. Montzka, David E. Oram, Benjamin R. Miller, Ray F. Weiss, Paul J. Fraser and Toste Tanhua, (2019), Atmospheric histories, growth rates and solubilities in seawater and other natural waters of the potential transient tracers HCFC-22, HCFC-141b, HCFC-142b, HFC-134a, HFC-125, HFC-23, PFC-14 and PFC-116, Ocean Science, 15, 1, 33-60, 10.5194/os-15-33-2019


Abstract. We present consistent annual mean atmospheric histories and growth rates for the mainly anthropogenic halogenated compounds HCFC-22, HCFC-141b, HCFC-142b, HFC-134a, HFC-125, HFC-23, PFC-14 and PFC-116, which are all potentially useful oceanic transient tracers (tracers of water transport within the ocean), for the Northern and Southern Hemisphere with the aim of providing input histories of these compounds for the equilibrium between the atmosphere and surface ocean. We use observations of these halogenated compounds made by the Advanced Global Atmospheric Gases Experiment (AGAGE), the Scripps Institution of Oceanography (SIO), the Commonwealth Scientific and Industrial Research Organization (CSIRO), the National Oceanic and Atmospheric Administration (NOAA) and the University of East Anglia (UEA). Prior to the direct observational record, we use archived air measurements, firn air measurements and published model calculations to estimate the atmospheric mole fraction histories. The results show that the atmospheric mole fractions for each species, except HCFC-141b and HCFC-142b, have been increasing since they were initially produced. Recently, the atmospheric growth rates have been decreasing for the HCFCs (HCFC-22, HCFC-141b and HCFC-142b), increasing for the HFCs (HFC-134a, HFC-125, HFC-23) and stable with little fluctuation for the PFCs (PFC-14 and PFC-116) investigated here. The atmospheric histories (source functions) and natural background mole fractions show that HCFC-22, HCFC-141b, HCFC-142b, HFC-134a, HFC-125 and HFC-23 have the potential to be oceanic transient tracers for the next few decades only because of the recently imposed bans on production and consumption. When the atmospheric histories of the compounds are not monotonically changing, the equilibrium atmospheric mole fraction (and ultimately the age associated with that mole fraction) calculated from their concentration in the ocean is not unique, reducing their potential as transient tracers. Moreover, HFCs have potential to be oceanic transient tracers for a longer period in the future than HCFCs as the growth rates of HFCs are increasing and those of HCFCs are decreasing in the background atmosphere. PFC-14 and PFC-116, however, have the potential to be tracers for longer periods into the future due to their extremely long lifetimes, steady atmospheric growth rates and no explicit ban on their emissions. In this work, we also derive solubility functions for HCFC-22, HCFC-141b, HCFC-142b, HFC-134a, HFC-125, HFC-23, PFC-14 and PFC-116 in water and seawater to facilitate their use as oceanic transient tracers. These functions are based on the Clark–Glew–Weiss (CGW) water solubility function fit and salting-out coefficients estimated by the poly-parameter linear free-energy relationships (pp-LFERs). Here we also provide three methods of seawater solubility estimation for more compounds. Even though our intention is for application in oceanic research, the work described in this paper is potentially useful for tracer studies in a wide range of natural waters, including freshwater and saline lakes, and, for the more stable compounds, groundwaters.

Minschwaner, Kenneth, Anthony T. Giljum, Gloria L. Manney, Irina Petropavlovskikh, Bryan J. Johnson and Allen F. Jordan, (2019), Detection and classification of laminae in balloon-borne ozonesonde profiles: application to the long-term record from Boulder, Colorado, Atmospheric Chemistry and Physics, 19, 3, 1853-1865, 10.5194/acp-19-1853-2019


Abstract. We quantify ozone variability in the upper troposphere and lower stratosphere (UTLS) by investigating lamination features in balloon measurements of ozone mixing ratio and potential temperature. Laminae are defined as stratified variations in ozone that meet or exceed a 10 % threshold for deviations from a basic state vertical profile of ozone. The basic state profiles are derived for each sounding using smoothing methods applied within a vertical coordinate system relative to the World Meteorological Organization (WMO) tropopause. We present results of this analysis for the 25-year record of ozonesonde measurements from Boulder, Colorado. The mean number of ozone laminae identified per sounding is about 9±2 (1σ). The root-mean-square relative amplitude is 20 %, and laminae with much larger amplitudes (>40 %) are seen in ∼ 2 % of the profiles. The vertical scale of detected ozone laminae typically ranges between 0.5 and 1.2 km. The lamina occurrence frequency varies significantly with altitude and is largest within ∼2 km of the tropopause. Overall, ozone laminae identified in our analysis account for more than one-third of the total intra-seasonal variability in ozone. A correlation technique between ozone and potential temperature is used to classify the subset of ozone laminae that are associated with gravity wave (GW) phenomena, which accounts for 28 % of all laminar ozone features. The remaining 72 % of laminae arise from non-gravity wave (NGW) phenomena. There are differences in both the vertical distribution and seasonality of GW versus NGW ozone laminae that are linked to the contrast in main generating mechanisms for each laminae type.

Moeini, Omid, Zahra Vaziri Zanjani, C. Thomas McElroy, David W. Tarasick, Robert D. Evans, Irina Petropavlovskikh and Keh-Harng Feng, (2019), The effect of instrumental stray light on Brewer and Dobson total ozone measurements, Atmospheric Measurement Techniques, 12, 1, 327-343, 10.5194/amt-12-327-2019


Abstract. Dobson and Brewer spectrophotometers are the primary, standard instruments for ground-based ozone measurements under the World Meteorological Organization's (WMO) Global Atmosphere Watch program. The accuracy of the data retrieval for both instruments depends on a knowledge of the ozone absorption coefficients and some assumptions underlying the data analysis. Instrumental stray light causes nonlinearity in the response of both the Brewer and Dobson to ozone at large ozone slant paths. In addition, it affects the effective ozone absorption coefficients and extraterrestrial constants that are both instrument-dependent. This effect has not been taken into account in the calculation of ozone absorption coefficients that are currently recommended by WMO for the Dobson network. The ozone absorption coefficients are calculated for each Brewer instrument individually, but in the current procedure the effect of stray light is not considered. This study documents the error caused by the effect of stray light in the Brewer and Dobson total ozone measurements using a physical model for each instrument. For the first time, new ozone absorption coefficients are calculated for the Brewer and Dobson instruments, taking into account the stray light effect. The analyses show that the differences detected between the total ozone amounts deduced from Dobson AD and CD pair wavelengths are related to the level of stray light within the instrument. The discrepancy introduced by the assumption of a fixed height for the ozone layer for ozone measurements at high latitude sites is also evaluated. The ozone data collected by two Dobson instruments during the period of December 2008 to December 2014 are compared with ozone data from a collocated double monochromator Brewer spectrophotometer (Mark III). The results illustrate the dependence of Dobson AD and CD pair measurements on stray light.

Montzka, S.A., G.J.M. Velders, P.B. Krummel, J. Muhle, V.L. Orkin, S. Park, N. Shah and H. Walters-Terrinoni, (2019), Scientific Assessment of Ozone Depletion: 2018, Chapter 2: Hydrofluorocarbons (HFCs), World Meteorological Organization, Report, 58, 144-204,


The Montreal Protocol is an international agreement designed to heal the ozone layer. It outlines schedules for the phase-out of ozone-depleting substances (ODSs) such as chlorofluorocarbons (CFCs), hydrochlorofluorocarbons (HCFCs), chlorinated solvents, halons, and methyl bromide. As a result of this phase-out, alternative chemicals and procedures were developed by industry for use in many applications including refrigeration, air-conditioning, foam-blowing, electronics, medicine, agriculture, and fire protection. Hydrofluorocarbons (HFCs) were used as ODS alternatives in many of these applications because they were suitable substitutes and they do not contain ozone-depleting chlorine or bromine; in addition, most HFCs have smaller climate impacts per molecule than the most widely used ODSs they replaced. Long-lived HFCs, CFCs, and HCFCs, however, are all potent greenhouse gases, and concerns were raised that uncontrolled future use of HFCs would lead to substantial climate warming.  As a result of these concerns, HFCs were included as one group of greenhouse gases for which emissions controls were
adopted by the 1997 Kyoto Protocol under the 1992 United Nations Framework Convention on Climate Change (UNFCCC). Consequently, developed countries (those listed in Annex I to this Convention, or “Annex I” Parties) supply annual emission estimates of HFCs to the UNFCCC. Since the Kyoto Protocol only specified limits on the sum of all controlled greenhouse gases, emissions of HFCs were not explicitly controlled. However, following the Kyoto Protocol, some countries enacted additional controls specifically limiting HFC use based on their global warming potentials (GWPs). Ultimately the Kigali Amendment to the
Montreal Protocol was agreed upon in 2016, and this Amendment supplies schedules for limiting the production and consumption of specific HFCs. Although the radiative forcing supplied by HFCs is currently small, this Amendment was designed to ensure that the radiative forcing from HFCs will not grow uncontrollably in the future. The Kigali Amendment will come into force at the start of 2019. HFC concentrations are currently monitored through atmospheric measurements. All HFCs with large abundances are monitored, as are most with small abundances. Most HFCs that are emitted to the atmosphere are intentionally produced for use in a variety of applications that were once dependent on ODSs. An exception is HFC-23, which is emitted to the atmosphere primarily as a by-product of HCFC-22 production. HFC-23 is also unique in that it has a substantially longer atmospheric lifetime and higher GWP than nearly all other HFCs. As a result, the Kigali Amendment includes different control schedules for HFC-23 production than for other HFCs. To date, HFC-23 emissions have been partially abated in developed countries through regulations or voluntary measures and in developing countries with assistance from the UNFCCC’s Clean Development Mechanism (CDM).

Naus, Stijn, Stephen A. Montzka, Sudhanshu Pandey, Sourish Basu, Ed J. Dlugokencky and Maarten Krol, (2019), Constraints and biases in a tropospheric two-box model of OH, Atmospheric Chemistry and Physics, 19, 1, 407-424, 10.5194/acp-19-407-2019


Abstract. The hydroxyl radical (OH) is the main atmospheric oxidant and the primary sink of the greenhouse gas CH4. In an attempt to constrain atmospheric levels of OH, two recent studies combined a tropospheric two-box model with hemispheric-mean observations of methyl chloroform (MCF) and CH4. These studies reached different conclusions concerning the most likely explanation of the renewed CH4 growth rate, which reflects the uncertain and underdetermined nature of the problem. Here, we investigated how the use of a tropospheric two-box model can affect the derived constraints on OH due to simplifying assumptions inherent to a two-box model. To this end, we derived species- and time-dependent quantities from a full 3-D transport model to drive two-box model simulations. Furthermore, we quantified differences between the 3-D simulated tropospheric burden and the burden seen by the surface measurement network of the National Oceanic and Atmospheric Administration (NOAA). Compared to commonly used parameters in two-box models, we found significant deviations in the magnitude and time-dependence of the interhemispheric exchange rate, exposure to OH, and stratospheric loss rate. For MCF these deviations can be large due to changes in the balance of its sources and sinks over time. We also found that changes in the yearly averaged tropospheric burden of CH4 and MCF can be obtained within 0.96ppbyr−1 and 0.14%yr−1 by the NOAA surface network, but that substantial systematic biases exist in the interhemispheric mixing ratio gradients that are input to two-box model inversions.

To investigate the impact of the identified biases on constraints on OH, we accounted for these biases in a two-box model inversion of MCF and CH4. We found that the sensitivity of interannual OH anomalies to the biases is modest (1%–2%), relative to the uncertainties on derived OH (3%–4%). However, in an inversion where we implemented all four bias corrections simultaneously, we found a shift to a positive trend in OH concentrations over the 1994–2015 period, compared to the standard inversion. Moreover, the absolute magnitude of derived global mean OH, and by extent, that of global CH4 emissions, was affected much more strongly by the bias corrections than their anomalies (∼10%). Through our analysis, we identified and quantified limitations in the two-box model approach as well as an opportunity for full 3-D simulations to address these limitations. However, we also found that this derivation is an extensive and species-dependent exercise and that the biases were not always entirely resolvable. In future attempts to improve constraints on the atmospheric oxidative capacity through the use of simple models, a crucial first step is to consider and account for biases similar to those we have identified for the two-box model.

Oltmans, S. J., L. C. Cheadle, B. J. Johnson, R. C. Schnell, D. Helmig, A. M. Thompson, P. Cullis, E. Hall, A. Jordan, C. Sterling, A. McClure-Begley, J. T. Sullivan, T. J. McGee and D. Wolfe, (2019), Boundary layer ozone in the Northern Colorado Front Range in July–August 2014 during FRAPPE and DISCOVER-AQ from vertical profile measurements, Elem Sci Anth, 7, 1, 6, 10.1525/elementa.345


The Northern Colorado Front Range (NCFR) east of the foothills of the Rocky Mountains and north of Denver and its nearby suburbs includes the moderate sized (~100,000) population centers of Boulder, Longmont, Ft. Collins, and Greeley. To the east is the Denver-Julesburg Basin (DJB), a major oil and natural gas production region that also includes significant agricultural and livestock operations. To the west, the Rocky Mountains with a number of 3000 m peaks rise abruptly and contribute to a complex geographical and meteorological regime (Losleben et al., 2000). During summer months, the Front Range and nearby foothills and mountains regularly experience high ozone (O3) episodes that exceed the National Ambient Air Quality Standard (NAAQS). The U.S. EPA classifies the region as an O3 non-attainment area.

During July–August 2014 a multi-institution campaign was carried out to investigate contributors to the high O3 episodes and possible solutions for controlling air pollution in the region. Two major components of the campaign were the Front Range Air Pollution and Photochemistry Experiment (FRAPPE), led by the Colorado Department of Public Health and Environment (CDPHE) and the National Center for Atmospheric Research (NCAR), and a Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) deployment led by NASA. An array of measurement platforms was deployed during the study period that included multiple aircraft, instrumented surface sites, mobile instrumented vans, and profiling capabilities using lidar, tethered balloons, release balloons, and a tall tower. Along with an extensive suite of chemical measurements many of the platforms and sites included meteorological parameters as well.

In this paper several vertical profile observing platforms including tethered and free release ozonesondes, a tall tower, and surface sites along a vertical elevation gradient (Figure 1) are used to relate surface O3 observations that are the basis for O3 regulatory actions to O3 dynamics and synoptic transport within the boundary layer. The focus of this study is primarily the NCFR, stretching approximately from Boulder to Ft. Collins and the area to the east. This region (Figure 1) is subject to several O3 precursor sources including urban emissions, agricultural and landfill emissions, and notably oil and natural gas related emissions (Pétron et al., 2012; Thompson et al., 2014; McDuffie et al., 2016; Halliday et al., 2016; Abeleira et al., 2017; Cheadle et al., 2017). Elevated O3 occurrences in the NCFR have been found to be influenced by both local and synoptic meteorology patterns (Toth and Johnson, 1985; Pfister et al., 2017b; Evans and Helmig, 2016).

 Original | PPT

Figure 1 

Map of the Northern Colorado Front Range. Map of the Northern Colorado Front Range showing monitoring sites (red triangles), population centers (black dots), and oil and natural gas wells (small gray dots). The inset shows sites along the Colorado Front Range in Boulder County that measured O3 during this study. Numbers next to the site name indicate the site elevation in km. DOI:

Meteorological conditions during July and August 2014 were generally typical of those seen during the summer in the NCFR (Toth and Johnson, 1985) with prevailing daytime upslope and nighttime downslope flow (Losleben et al, 2000). On selected occasions closed cyclonic circulations NE of Denver developed under conditions when synoptic frontal system passed through the region (Pfister et al., 2017b). Especially in July several periods of high pressure dominated the region leading to the thermally driven circulation systems noted above. Transport pathways computed from back trajectories reflected both these more local circulations patterns and periods when synoptic scale disturbances dominated.

In this paper profiles of O3 and accompanying meteorological variables are used to investigate the contribution of O3 mixed down into the boundary layer, boundary layer O3 dependence on transport, and the production of O3 within the boundary layer including near the surface. Cases with enhanced afternoon O3 levels (≥75 ppb) are contrasted with days with minimal afternoon O3 enhancements (O3 ≤65 ppb).

Ortega, Ivan, Rebecca R. Buchholz, Emrys G. Hall, Dale F. Hurst, Allen F. Jordan and James W. Hannigan, (2019), Tropospheric water vapor profiles obtained with FTIR: comparison with balloon-borne frost point hygrometers and influence on trace gas retrievals, Atmospheric Measurement Techniques, 12, 2, 873-890, 10.5194/amt-12-873-2019


Abstract. Retrievals of vertical profiles of key atmospheric gases provide a critical long-term record from ground-based Fourier transform infrared (FTIR) solar absorption measurements. However, the characterization of the retrieved vertical profile structure can be difficult to validate, especially for gases with large vertical gradients and spatial–temporal variability such as water vapor. In this work, we evaluate the accuracy of the most common water vapor isotope (H216O, hereafter WV) FTIR retrievals in the lower and upper troposphere–lower stratosphere. Coincident high-quality vertically resolved WV profile measurements obtained from 2010 to 2016 with balloon-borne NOAA frost point hygrometers (FPHs) are used as reference to evaluate the performance of the retrieved profiles at two sites: Boulder (BLD), Colorado, and at the mountaintop observatory of Mauna Loa (MLO), Hawaii. For a meaningful comparison, the spatial–temporal variability has been investigated. We present results of comparisons among FTIR retrievals with unsmoothed and smoothed FPH profiles to assess WV vertical gradients. Additionally, we evaluate the quantitative impact of different a priori profiles in the retrieval of WV. An orthogonal linear regression analysis shows the best correlation among tropospheric layers using ERA-Interim (ERA-I) a priori profiles and biases are lower for unsmoothed comparisons. In Boulder, we found a negative bias of 0.02±1.9 % (r=0.95) for the 1.5–3 km layer. A larger negative bias of 11.1±3.5 % (r=0.97) was found in the lower free troposphere layer of 3–5 km attributed to rapid vertical change of WV, which is not always captured by the retrievals. The bias improves in the 5–7.5 km layer (1.0±5.3 %, r=0.94). The bias remains at about 13 % for layers above 7.5 km but below 13.5 km. At MLO the spatial mismatch is significantly larger due to the launch of the sonde being farther from the FTIR location. Nevertheless, we estimate a negative bias of 5.9±4.6 % (r=0.93) for the 3.5–5.5 km layer and 9.9±3.7 % (r=0.93) for the 5.5–7.5 km layer, and we measure positive biases of 6.2±3.6 % (r=0.95) for the 7.5–10 km layer and 12.6 % and greater values above 10 km. The agreement for the first layer is significantly better at BLD because the air masses are similar for both FTIR and FPH. Furthermore, for the first time we study the influence of different WV a priori profiles in the retrieval of selected gas profiles. Using NDACC standard retrievals we present results for hydrogen cyanide (HCN), carbon monoxide (CO), and ethane (C2H6) by taking NOAA FPH profiles as the ground truth and evaluating the impact of other WV profiles. We show that the effect is minor for C2H6 (bias <0.5 % for all WV sources) among all vertical layers. However, for HCN we found significant biases between 6 % for layers close to the surface and 2 % for the upper troposphere depending on the WV profile source. The best results (reduced bias and precision and r values closer to unity) are always found for pre-retrieved WV. Therefore, we recommend first retrieving WV to use in subsequent retrieval of gases.

Park, Soojin, Sang-Woo Kim, Neng-Huei Lin, Shantanu Kumar Pani, Patrick J. Sheridan and Elisabeth Andrews, (2019), Variability of Aerosol Optical Properties Observed at Polluted Marine (Gosan, Korea) and High-altitude Mountain (Lulin, Taiwan) Sites in the Asian Continental Outflow, Aerosol and Air Quality Research, 10.4209/aaqr.2018.11.0416


We investigated variability of in-situ and columnar aerosol optical properties (AOPs) at two regional background sites in the Asian continental outflow from 2012 to 2014. The monthly variability of AOPs at each site was explained using regional-scale atmospheric circulation patterns and the resulting changes in air mass sources. The median aerosol scattering and absorption coefficient values for sub-10 μm particles (σsp10 μm and σap10 μm) at the polluted marine site, Gosan (GSN; σsp10 μm: 59.94 Mm-1σap10 μm: 4.73 Mm-1), were larger than those of the high-altitude mountain site, Lulin (LLN; σsp10 μm: 17.07 Mm-1, σap10 μm: 1.72 Mm-1). Elevated σsp10 μm and σap10 μm at GSN in May and June can be explained by the accumulation of locally emitted aerosols together with long-range transported aerosols under stagnant synoptic patterns, enhanced particle formation and subsequent growth. The significant peaks of σsp10 μm and σap10 μm, with median values of 50.76 Mm-1 and 5.92 Mm-1, at LLN during March and April can be attributed to biomass burning aerosols transported from Northern Indochina and Southern China. The LLN site was mostly influenced by clean free tropospheric air and maritime air masses during the other months, showing medians of 14.65±19.81 Mm-1 and 1.46±1.86 Mm-1 in σsp10 μm and σap10 μm, respectively. Smaller median values of the sub-micron to sub-10 μm ratio of aerosol scattering (0.60) and absorption (0.81) at GSN, compared to LLN (0.81 and 0.91, respectively), are indicative of the larger proportion of coarse aerosols such as sea salt and dust at GSN. Single scattering albedo for sub-10 μm particles showed similar median values at both sites (GSN: 0.93±0.02, LLN: 0.91±0.03).

Petropavlovskikh, I., S. Godin-Beekmann, D. Hubert, R. Damadeo, B. Hassler and V. Sofieva, (2019), SPARC/IO3C/GAW Report on Long-term Ozone Trends and Uncertainties in the Stratosphere, SPARC LOTUS Activity, SPARC , i-78, 10.17874/f899e57a20b



The assessment of long-term observations by LOTUS confirms the significant decline of ozone concentrations in the upper stratosphere (at altitudes above the 10–5 hPa level) between January 1985 and December 1996. The strongest trends are observed near 2 hPa (~42 km) with values of 5.9–6.2 % per decade at mid-latitudes and 4.8 % per decade in the tropics. Trends are significant at more than 5 standard deviations in this altitude range.
Trends derived from satellite and ground-based records in the pre-1997 time period agree with climate model simulations within respective uncertainties thus confirming our understanding of ozone loss processes in the upper stratosphere during that period.
Between January 2000 and December 2016, positive trends are obtained throughout the upper stratosphere for satellite and ground-based records. The combined trends from six merged satellite records are larger in the Northern Hemisphere mid-latitudes (2–3 % per decade between ~5–1 hPa) than in the tropics (1–1.5 % per decade between ~3–1 hPa) and Southern Hemisphere mid-latitudes (~2 % per decade near 2 hPa). Statistical confidence is largest for trends in the Northern Hemisphere mid-latitudes.
For altitudes below the 4 hPa level, ozone trends in the post-2000 time period are not significant. Though not significant, negative ozone trends of 0.5–1.5 % per decade are consistently detected by multiple satellite combined records in the 50–15 hPa altitude range over the tropics. Trends derived from ground-based data and Chemistry-Climate Model Initiative (CCMI) model simulations are generally consistent but more variable in this region. The mean CCMI model trend is negative at altitudes below 30 hPa, but the range of individual model trends is large; trends in ground-based records tend to be negative at 20 hPa but increase at lower altitudes (except trends from the microwave records). At mid-latitudes, the trends are close to zero down to 50 hPa.
Larger differences in post-2000 trends from the various records are observed in the lowermost stratosphere (100–50 hPa) in all latitude bands. Non-significant negative trends are derived from merged satellite records over the tropics and the Northern Hemisphere mid-latitudes. Model simulations show positive trends in the mid-latitudes in both hemispheres in this altitude range, although the trends are not statistically significant.
LOTUS estimates of past and recent ozone trends are in fairly good agreement with results from previous studies. For the post-2000 period, the largest differences are found throughout the middle stratosphere. These differences stem primarily from extensions of and revisions to existing data records, the addition of new data records, and in some cases the use of a different trend model.
While trend values in recent studies are fairly similar, the uncertainties and hence significances of the combined trends in broad latitude bands differ substantially. The LOTUS approach, based on both error propagation and standard error of the mean, also explicitly accounts for correlation between the data sets, which results in more conservative uncertainties and thus lower, but more realistic, confidence in positive upper stratospheric trend values compared to the most recently published assessment of merged satellite data set trends.

Have ozone concentrations in the stratosphere significantly increased since the end of the 1990s when levels of ozone depleting substances (ODSs) started to decline? Finding an answer to this question is of great societal importance to ensure that the measures taken by the Montreal Protocol and subsequent amendments to reduce ODSs continue to adequately protect the ozone layer. However, the confidence with which we can assess changes in stratospheric ozone since the mid-1990s has been the subject of considerable scientific debate in recent years, as it depends on the data sets and the analysis methods used. Settling this scientific debate is one of the main objectives of the LOTUS activity, short for Long-term Ozone Trends and Uncertainties in the Stratosphere. Below, we summarise the main results obtained during the first phase of LOTUS, which was primarily targeted at providing timely input to the 2018 World Meteorological Organization (WMO) Ozone Assessment (WMO, 2018). During this phase we reevaluated the satellite and ground-based data records as well as the time series analysis methods commonly used to derive long-term trends. Using a single “LOTUS regression” model, we reassessed past and recent trends in the vertical distribution of stratospheric ozone from the updated individual data records. We then developed a new approach for combining the individual trend estimates from satellite-based records into a single best estimate of ozone profile trends with associated uncertainty estimates. Finally, we compared the satellite-based profile trends in broad latitude bands to trends from ground-based data, from the collection of CCMI-1 model simulations, and from past evaluations of satellite-based trends in peer-reviewed literature.  Some regions in the stratosphere have not been considered (e.g., polar) or have not been analysed in full detail (e.g.,
lower stratosphere) because of the timeline for the 2018 WMO Ozone Assessment (WMO, 2018).
Sorribas, M., E. Andrews, J.A. Ogren, A. del Aguila, R. Fraile, P. Sheridan and M. Yela, (2019), Climatological study for understanding the aerosol radiative effects at southwest Atlantic coast of Europe, Elsevier, 205, 52-66, 10.1016/j.atmosenv.2019.02.017


In order to describe the means, variability and trends of the aerosol radiative effects on the southwest Atlantic coast of Europe, 11 years of aerosol light scatteringsp) and 4 years of aerosol light absorptionap) are analyzed. A 2006–2016 trend analysis of σsp for D < 10 μm indicates statistically significant trends for March, May–June and September–November, with a decreasing trend ranging from −1.5 to −2.8 Mm−1/year. In the 2009–2016 period, the decreasing trend is only observed for the months of June and September. For scattering Ångström exponent (SAE) there is an increasing trend during June with a rate of 0.059/year and a decreasing trend during October with −0.060/year. The trends observed may be caused by a reduction of Saharan dust aerosol or a drop in particle loading in anthropogenic influenced air masses. The relationship between SAE and absorption Ångström exponent is used to assess the aerosol typing. Based on this typing, the sub-micron particles are dominated by black carbon, mixed black and brown carbon or marine with anthropogenic influences, while the super-micrometer particles are desert dust and sea spray aerosol. The mean and standard deviation of the dry aerosol direct radiative effect at the top of the atmosphere (DRETOA) are −4.7 ± 4.2 W m−2. DRETOA for marine aerosol shows all observations more negative than −4 W m−2 and for anthropogenic aerosol type, DRETOA ranges from −5.0 to −13.0 W m−2. DRETOA of regional marine aerosol ranges from −3 to −7 W m−2, as it consists of a mixture of sea salt and anthropogenic aerosol. The variability in DRETOA is mainly dependent on AOD, given that variations in backscatter fraction and the single scattering albedo tend to counteract each other in the radiative forcing efficiency equation. The results shown here may help in interpretation of satellite retrieval products and provide context for model evaluation.

Suntharalingam, P., E. Buitenhuis, L. J. Carpenter, J. H. Butler, M. J. Messias, S. J. Andrews and S. C. Hackenberg, (2019), Evaluating Oceanic Uptake of Atmospheric CCl: A Combined Analysis of Model Simulations and Observations , Geophysical Research Letters, 46, 1, 472-482, 10.1029/2018GL080612


We provide new estimates of the air‐sea flux of CCl4 using simulations from a global ocean biogeochemistry model (NEMO‐PlankTOM) in combination with depth‐resolved CCl4 observations from global oceanic databases. Estimates of global oceanic CCl4 uptake are derived from a range of model analyses, including prescribed parameterizations using reported values on hydrolysis and degradation, and analyses optimized using the global observational databases. We evaluate the sensitivity of our results to uncertainties in air‐sea gas exchange parameterization, estimation period, and circulation processes. Our best constrained estimate of ocean CCl4 uptake for the period 1996–2000 is 20.1 Gg/year (range 16.6–22.7), corresponding to estimates of the partial atmospheric lifetime with respect to ocean uptake of 124 (110–150) years. This new oceanic lifetime implies higher emissions of CCl4 than currently estimated and therefore a larger missing atmospheric source of CCl4.

Yang, Huang, Darryn W. Waugh, Clara Orbe, Prabir K. Patra, Patrick Jöckel, Jean-Francois Lamarque, Simone Tilmes, Douglas Kinnison, James W. Elkins and Edward J. Dlugokencky, (2019), Evaluating Simulations of Interhemispheric Transport: Interhemispheric Exchange Time Versus SF Age , Geophysical Research Letters, 46, 2, 1113-1120, 10.1029/2018GL080960


Two recent studies using sulfur hexafluoride (SF6) observations to evaluate interhemispheric transport in two different ensembles of atmospheric chemistry models reached different conclusions on model performance. We show here that the different conclusions are due to the use of different metrics and not differences in the performance of the models. For both model ensembles, the multimodel mean interhemispheric exchange time τex agrees well with observations, but in nearly all models the SF6 age in the southern hemisphere is older than observed. This occurs because transport from the northern extratropics into the tropics is too slow in most models, and the SF6 age is more sensitive to this bias than τex. Thus, simulating τex correctly does not necessarily mean that transport from northern midlatitudes into the southern hemisphere is correct. It also suggests that more attention needs to be paid to evaluating transport from northern midlatitudes into the tropics.

Zhang, Taiping, Paul W. Stackhouse, Stephen J. Cox, J. Colleen Mikovitz and Charles N. Long, (2019), Clear-sky shortwave downward flux at the Earth's surface: Ground-based data vs. satellite-based data, Journal of Quantitative Spectroscopy and Radiative Transfer, 224, 247-260, 10.1016/j.jqsrt.2018.11.015


The radiative flux data and other meteorological data in the BSRN archive start in 1992, but the RadFlux data, the clear-sky radiative fluxes at the BSRN sites empirically inferred through regression analyses of actually observed clear-sky fluxes, did not come into existence until the early 2000s, and at first, they were limited to the 7 NOAA SURFRAD and 4 DOE ARM sites, a subset of the BSRN sites. Recently, the RadFlux algorithm was applied more extensively to the BSRN sites for the production of clear-sky ground-based fluxes. At the time of this writing, there are 7119 site-months of clear-sky fluxes at 42 BSRN sites spanning from 1992 to late 2017. These data provide an unprecedented opportunity to validate the satellite-based clear-sky fluxes. In this paper, the GEWEX SRB GSW(V3.0) clear-sky shortwave downward fluxes spanning 24.5 years from July 1983 to December 2007, the CERES SYN1deg(Ed4A) and EBAF(Ed4.0) clear-sky shortwave fluxes spanning March 2000 to mid-2017 are compared with their RadFlux counterparts on the hourly, 3-hourly, daily and monthly time scales. All the three datasets show reasonable agreement with their ground-based counterparts. Comparison of the satellite-based surface shortwave clear-sky radiative fluxes to the BSRN RadFlux analysis shows negative biases (satellite-based minus RadFlux). Further analysis shows that the satellite-based atmosphere contains greater aerosol loading as well as more precipitable water than RadFlux analysis estimates.