ALBERT

Description: Ozone is a key constituent of the troposphere, where it drives photochemical processes, impacts air quality, and acts as a climate forcer. Large-scale in situ observations of ozone commensurate with the grid resolution of current Earth system models are necessary to validate model outputs and satellite retrievals. In this paper, we examine measurements from the Atmospheric Tomography (ATom; four deployments in 2016–2018) and the HIAPER Pole-to-Pole Observations (HIPPO; five deployments in 2009–2011) experiments, two global-scale airborne campaigns covering the Pacific and Atlantic basins. ATom and HIPPO represent the first global-scale, vertically resolved measurements of O3 distributions throughout the troposphere, with HIPPO sampling the atmosphere over the Pacific and ATom sampling both the Pacific and Atlantic. Given the relatively limited temporal resolution of these two campaigns, we first compare ATom and HIPPO ozone data to longer-term observational records to establish the representativeness of our dataset. We show that these two airborne campaigns captured on average 53 %, 54 %, and 38 % of the ozone variability in the marine boundary layer, free troposphere, and upper troposphere–lower stratosphere (UTLS), respectively, at nine well-established ozonesonde sites. Additionally, ATom captured the most frequent ozone concentrations measured by regular commercial aircraft flights in the northern Atlantic UTLS. We then use the repeated vertical profiles from these two campaigns to confirm and extend the existing knowledge of tropospheric ozone spatial and vertical distributions throughout the remote troposphere. We highlight a clear hemispheric gradient, with greater ozone in the Northern Hemisphere, consistent with greater precursor emissions and consistent with previous modeling and satellite studies. We also show that the ozone distribution below 8 km was similar in the extra-tropics of the Atlantic and Pacific basins, likely due to zonal circulation patterns. However, twice as much ozone was found in the tropical Atlantic as in the tropical Pacific, due to well-documented dynamical patterns transporting continental air masses over the Atlantic. Finally, we show that the seasonal variability of tropospheric ozone over the Pacific and the Atlantic basins is driven year-round by transported continental plumes and photochemistry, and the vertical distribution is driven by photochemistry and mixing with stratospheric air. This new dataset provides additional constraints for global climate and chemistry models to improve our understanding of both ozone production and loss processes in remote regions, as well as the influence of anthropogenic emissions on baseline ozone.

Description: Wetland forests around the world have been reduced to a small proportion of their original expanse due to changing climatic conditions and intensification of human land use activities. As a case in point, the Columbia bottomland hardwood forests along the Brazos–Colorado coastal basin on the Gulf coast of Texas are currently threatened by an increasingly erratic hydroclimate in the form of both extreme floods and droughts and by urban expansion. In this study, we use dendrochronology and tree ring carbon isotopes to understand the effect of changing hydroclimatic conditions on the functional attributes of these forests. We examined the tree rings of Quercus nigra at four sites within the Columbia bottomlands, of which one site experiences frequent and prolonged flooding, while the other three are less flood prone. The objectives of this study were to (i) understand the impact of hydroclimatic variation on radial growth, using tree ring width analysis, (ii) assess the magnitude of physiological stress inflicted by extreme hydroclimatic conditions, using tree ring Δ13C measurements as a proxy, and (iii) evaluate the relationship between tree ring width and Δ13C values. Radial growth across the landscape was influenced most strongly by the midgrowing season climate, while the early growing season climate had the strongest effect on Δ13C. Growth inhibition was minimal, and tree ring Δ13C values were not affected in trees at the wetter site under extreme hydrological conditions such as droughts or floods. In addition, trees at the wet site were less sensitive to precipitation and showed no response to higher temperatures. In contrast, trees at the three drier sites experienced growth inhibition and had lower tree ring Δ13C values during dry periods. Our results indicate more favorable growing conditions and lower stress in trees growing under wetter hydrological conditions. Management and conservation strategies dependent on site-specific conditions are critical for the health of these wetland forests under a rapidly changing hydroclimate. This study provides the first dendrochronological baseline for this region and a better understanding of favorable conditions for the growth and health of these forests, which can assist in management decisions such as streamflow regulation and conservation plans.

Description: In the current era of rapid climate change, accurate characterization of climate-relevant gas dynamics – namely production, consumption, and net emissions – is required for all biomes, especially those ecosystems most susceptible to the impact of change. Marine environments include regions that act as net sources or sinks for numerous climate-active trace gases including methane (CH4) and nitrous oxide (N2O). The temporal and spatial distributions of CH4 and N2O are controlled by the interaction of complex biogeochemical and physical processes. To evaluate and quantify how these mechanisms affect marine CH4 and N2O cycling requires a combination of traditional scientific disciplines including oceanography, microbiology, and numerical modeling. Fundamental to these efforts is ensuring that the datasets produced by independent scientists are comparable and interoperable. Equally critical is transparent communication within the research community about the technical improvements required to increase our collective understanding of marine CH4 and N2O. A workshop sponsored by Ocean Carbon and Biogeochemistry (OCB) was organized to enhance dialogue and collaborations pertaining to marine CH4 and N2O. Here, we summarize the outcomes from the workshop to describe the challenges and opportunities for near-future CH4 and N2O research in the marine environment.

Description: The optical properties of black carbon (BC) are a major source of uncertainty in regional and global climate studies. In the past, detailed investigation of BC absorption has been hampered by systematic biases in the measurement instrumentation. We present airborne measurements of aerosol absorption and black carbon microphysical properties in highly aged biomass burning plumes measured 4–8 d from their source over the southeast Atlantic Ocean during CLARIFY-2017, using a suite of novel photoacoustic spectrometers to measure aerosol absorption at 405, 514, and 655 nm and a single-particle soot photometer to measure the BC mass concentration, size, and mixing state. These measurements are of sufficient quality and detail to provide constraint on optical schemes used in climate models for the first time in biomass burning plumes far from their source – an aerosol environment that is one of the most important climatically. The average absorption Ångström exponents (AAE) were 1.38 over the wavelength range from 405 to 514 nm and 0.88 over the range from 514 to 655 nm, suggesting that brown carbon (BrC) contributed to 11±2 % of absorption at 405 nm. The effective organic aerosol (OA) mass absorption coefficient (MAC) was 0.31±0.09 m2 g−1 at 405 nm. The BC particles were universally thickly coated, and almost no externally mixed BC particles were detected. The average MAC of BC was 20±4, 15±3, and 12±2 m2g−1 at wavelengths of 405, 514, and 655 nm respectively, with equivalent absorption enhancements of around 1.85±0.45 at all three wavelengths, suggesting that the thick coatings acted as a lens that enhanced light absorption by the BC. We compared the measured MAC and AAE values with those calculated using several optical models and absorption parameterisations that took the measured BC mass and mixing state as inputs. Homogeneous grey-sphere Mie models were only able to replicate MAC for some low (real and imaginary) values of the complex BC refractive index (mBC) at the shortest wavelength, but they would have to use unrealistically low values of mBC to accurately replicate the AAE. A core–shell Mie model was able to generate good agreement for MAC in the green–red end of the visible spectrum for most values of mBC. However, there are no possible values of mBC that produce MAC values that agree with our observations at all three wavelengths, due to a wavelength-dependent underestimation of the MAC of the underlying BC core. Four semiempirical parameterisations from the literature were also tested, linking the BC mixing state to either the MAC or absorption enhancement. Two of these schemes produced results that agreed within a few percent with the measured MAC at all three wavelengths, and the AAE agreed well when discounting the effects of BrC. Our results uniquely demonstrate the validity of absorption parameterisations, as well as the failings of Mie calculations, in this highly aged environment. We recommend that future work should conduct similar analyses in environments where BC has different properties; future studies should also investigate the impact of implementing these types of schemes within climate models as well as the impact of developing equivalent schemes for light scattering by soot particles at visible wavelengths.

Description: The detection and attribution of high background ozone (O3) events in the southwestern US is challenging but relevant to the effective implementation of the lowered National Ambient Air Quality Standard (NAAQS; 70 ppbv). Here we leverage intensive field measurements from the Fires, Asian, and Stratospheric Transport−Las Vegas Ozone Study (FAST-LVOS) in May–June 2017, alongside high-resolution simulations with two global models (GFDL-AM4 and GEOS-Chem), to study the sources of O3 during high-O3 events. We show possible stratospheric influence on 4 out of the 10 events with daily maximum 8 h average (MDA8) surface O3 above 65 ppbv in the greater Las Vegas region. While O3 produced from regional anthropogenic emissions dominates pollution events in the Las Vegas Valley, stratospheric intrusions can mix with regional pollution to push surface O3 above 70 ppbv. GFDL-AM4 captures the key characteristics of deep stratospheric intrusions consistent with ozonesondes, lidar profiles, and co-located measurements of O3, CO, and water vapor at Angel Peak, whereas GEOS-Chem has difficulty simulating the observed features and underestimates observed O3 by ∼20 ppbv at the surface. On days when observed MDA8 O3 exceeds 65 ppbv and the AM4 stratospheric ozone tracer shows 20–40 ppbv enhancements, GEOS-Chem simulates ∼15 ppbv lower US background O3 than GFDL-AM4. The two models also differ substantially during a wildfire event, with GEOS-Chem estimating ∼15 ppbv greater O3, in better agreement with lidar observations. At the surface, the two models bracket the observed MDA8 O3 values during the wildfire event. Both models capture the large-scale transport of Asian pollution, but neither resolves some fine-scale pollution plumes, as evidenced by aerosol backscatter, aircraft, and satellite measurements. US background O3 estimates from the two models differ by 5 ppbv on average (greater in GFDL-AM4) and up to 15 ppbv episodically. Uncertainties remain in the quantitative attribution of each event. Nevertheless, our multi-model approach tied closely to observational analysis yields some process insights, suggesting that elevated background O3 may pose challenges to achieving a potentially lower NAAQS level (e.g., 65 ppbv) in the southwestern US.

Description: Low clouds persist in the summer Arctic with important consequences for the radiation budget. In this study, we simulate the linear relationship between liquid water content (LWC) and cloud droplet number concentration (CDNC) observed during an aircraft campaign based out of Resolute Bay, Canada, conducted as part of the Network on Climate and Aerosols: Addressing Key Uncertainties in Remote Canadian Environments study in July 2014. Using a single-column model, we find that autoconversion can explain the observed linear relationship between LWC and CDNC. Of the three autoconversion schemes we examined, the scheme using continuous drizzle (Khairoutdinov and Kogan, 2000) appears to best reproduce the observed linearity in the tenuous cloud regime (Mauritsen et al., 2011), while a scheme with a threshold for rain (Liu and Daum, 2004) best reproduces the linearity at higher CDNC. An offline version of the radiative transfer model used in the Canadian Atmospheric Model version 4.3 is used to compare the radiative effects of the modelled and observed clouds. We find that there is no significant difference in the upward longwave cloud radiative effect at the top of the atmosphere from the three autoconversion schemes (p=0.05) but that all three schemes differ at p=0.05 from the calculations based on observations. In contrast, the downward longwave and shortwave cloud radiative effect at the surface for the Wood (2005b) and Khairoutdinov and Kogan (2000) schemes do not differ significantly (p=0.05) from the observation-based radiative calculations, while the Liu and Daum (2004) scheme differs significantly from the observation-based calculation for the downward shortwave but not the downward longwave fluxes.

Description: Evaluation of photosynthetic competency in time and space is critical for better estimates and models of oceanic primary productivity. This is especially true for areas where the lack of iron (Fe) limits phytoplankton productivity, such as the Southern Ocean. Assessment of photosynthetic competency on large scales remains challenging, but phytoplankton chlorophyll a fluorescence (ChlF) is a signal that holds promise in this respect as it is affected by, and consequently provides information about, the photosynthetic efficiency of the organism. A second process affecting the ChlF signal is heat dissipation of absorbed light energy, referred to as non-photochemical quenching (NPQ). NPQ is triggered when excess energy is absorbed, i.e. when more light is absorbed than can be used directly for photosynthetic carbon fixation. The effect of NPQ on the ChlF signal complicates its interpretation in terms of photosynthetic efficiency, and therefore most approaches relating ChlF parameters to photosynthetic efficiency seek to minimize the influence of NPQ by working under conditions of sub-saturating irradiance. Here, we propose that NPQ itself holds potential as an easily acquired optical signal indicative of phytoplankton physiological state with respect to Fe limitation. We present data from a research voyage to the Subantarctic Zone south of Australia. Incubation experiments confirmed that resident phytoplankton were Fe-limited, as the maximum quantum yield of primary photochemistry, Fv∕Fm, measured with a fast repetition rate fluorometer (FRRf), increased significantly with Fe addition. The NPQ “capacity” of the phytoplankton also showed sensitivity to Fe addition, decreasing with increased Fe availability, confirming previous work. The fortuitous presence of a remnant warm-core eddy in the vicinity of the study area allowed comparison of fluorescence behaviour between two distinct water masses, with the colder water showing significantly lower Fv∕Fm than the warmer eddy waters, suggesting a difference in Fe limitation status between the two water masses. Again, NPQ capacity measured with the FRRf mirrored the behaviour observed in Fv∕Fm, decreasing as Fv∕Fm increased in the warmer water mass. We also analysed the diel quenching of underway fluorescence measured with a standard fluorometer, such as is frequently used to monitor ambient chlorophyll a concentrations, and found a significant difference in behaviour between the two water masses. This difference was quantified by defining an NPQ parameter akin to the Stern–Volmer parameterization of NPQ, exploiting the fluorescence quenching induced by diel fluctuations in incident irradiance. We propose that monitoring of this novel NPQ parameter may enable assessment of phytoplankton physiological status (related to Fe availability) based on measurements made with standard fluorometers, as ubiquitously used on moorings, ships, floats and gliders.

Description: Satellite retrievals of column-averaged dry-air mole fractions of carbon dioxide (CO2) and methane (CH4), denoted XCO2 and XCH4, respectively, have been used in recent years to obtain information on natural and anthropogenic sources and sinks and for other applications such as comparisons with climate models. Here we present new data sets based on merging several individual satellite data products in order to generate consistent long-term climate data records (CDRs) of these two Essential Climate Variables (ECVs). These ECV CDRs, which cover the time period 2003–2018, have been generated using an ensemble of data products from the satellite sensors SCIAMACHY/ENVISAT and TANSO-FTS/GOSAT and (for XCO2) for the first time also including data from the Orbiting Carbon Observatory 2 (OCO-2) satellite. Two types of products have been generated: (i) Level 2 (L2) products generated with the latest version of the ensemble median algorithm (EMMA) and (ii) Level 3 (L3) products obtained by gridding the corresponding L2 EMMA products to obtain a monthly 5∘×5∘ data product in Obs4MIPs (Observations for Model Intercomparisons Project) format. The L2 products consist of daily NetCDF (Network Common Data Form) files, which contain in addition to the main parameters, i.e., XCO2 or XCH4, corresponding uncertainty estimates for random and potential systematic uncertainties and the averaging kernel for each single (quality-filtered) satellite observation. We describe the algorithms used to generate these data products and present quality assessment results based on comparisons with Total Carbon Column Observing Network (TCCON) ground-based retrievals. We found that the XCO2 Level 2 data set at the TCCON validation sites can be characterized by the following figures of merit (the corresponding values for the Level 3 product are listed in brackets) – single-observation random error (1σ): 1.29 ppm (monthly: 1.18 ppm); global bias: 0.20 ppm (0.18 ppm); and spatiotemporal bias or relative accuracy (1σ): 0.66 ppm (0.70 ppm). The corresponding values for the XCH4 products are single-observation random error (1σ): 17.4 ppb (monthly: 8.7 ppb); global bias: −2.0 ppb (−2.9 ppb); and spatiotemporal bias (1σ): 5.0 ppb (4.9 ppb). It has also been found that the data products exhibit very good long-term stability as no significant long-term bias trend has been identified. The new data sets have also been used to derive annual XCO2 and XCH4 growth rates, which are in reasonable to good agreement with growth rates from the National Oceanic and Atmospheric Administration (NOAA) based on marine surface observations. The presented ECV data sets are available (from early 2020 onwards) via the Climate Data Store (CDS, https://cds.climate.copernicus.eu/, last access: 10 January 2020) of the Copernicus Climate Change Service (C3S, https://climate.copernicus.eu/, last access: 10 January 2020).

Description: Nonuniform rock uplift in the form of tilting has been documented in convergent margins, postorogenic landscapes, and extensional provinces. Despite the prevalence of tilting, the transient fluvial response to tilting has not been quantified such that tectonic histories involving tilt can be extracted from river network forms. We used numerical landscape evolution models to characterize the transient erosional response of a river network initially at equilibrium to rapid tilting. We focus on the case of punctuated rigid-block tilting, though we explore longer-duration tilting events and nonuniform uplift that deviates from perfect rigid-block tilting such as that observed when bending an elastic plate or with more pronounced internal deformation of a fault-bounded block. Using a model river network composed of linked 1-D river longitudinal profile evolution models, we show that the transient response to a punctuated rigid-block tilting event creates a suite of characteristic forms or geomorphic signatures in mainstem and tributary profiles that collectively are distinct from those generated by other perturbations, such as a step change in the uniform rock uplift rate or a major truncation of the headwater drainage area, that push a river network away from equilibrium. These signatures include (1) a knickpoint in the mainstem that separates a downstream profile with uniform steepness (i.e., channel gradient normalized for drainage area) from an upstream profile with nonuniform steepness, with the mainstem above the knickpoint more out of equilibrium than the tributaries following forward tilting toward the outlet, versus the mainstem less out of equilibrium than the tributaries following back tilting toward the headwaters; (2) a pattern of mainstem incision below paleo-topography markers that increases linearly up to the mainstem knickpoint or vice versa following back tilting; and (3) tributary knickzones with nonuniform steepness that mirrors that of the mainstem upstream of the slope-break knickpoint. Immediately after a punctuated tilting event, knickpoints form at the mainstem outlet and each mainstem–tributary junction. Time since the cessation of rapid tilting is recorded by the mainstem knickpoint location relative to base level and by the upstream end of tributary knickzones relative to the mainstem–tributary junction. Tilt magnitude is recorded in the spatial gradient of mainstem incision depth and, in the forward tilting case, also by the spatial gradient in tributary knickzone drop height. Heterogeneous lithology can modulate the transient response to tilting and, post tilt, knickpoints can form anywhere in a stream network where more erodible rock occurs upstream of less erodible rock. With a full 2-D model, we show that stream segments flowing in the tilt direction have elevated channel gradient early in the transient response. Tilting is also reflected in network topologic changes via stream capture oriented in the direction of tilt. As an example of how these geomorphic signatures can be used in concert with each other to estimate the timing and magnitude of a tilting event, we show a sample of rivers from two field sites: the Sierra Nevada, California, USA, and the Sierra San Pedro Mártir, Baja California, Mexico, two ranges thought to have been tilted westward toward river outlets in the late Cenozoic.

Description: A key limitation of volcanic forcing datasets for the Pinatubo period is the large uncertainty that remains with respect to the extent of the optical depth of the Pinatubo aerosol cloud in the first year after the eruption, the saturation of the SAGE-II instrument restricting it to only be able to measure the upper part of the aerosol cloud in the tropics. Here we report the recovery of stratospheric aerosol measurements from two shipborne lidars, both of which measured the tropical reservoir of volcanic aerosol produced by the June 1991 Mount Pinatubo eruption. The lidars were on board two Soviet vessels, each ship crossing the Atlantic, their measurement datasets providing unique observational transects of the Pinatubo cloud across the tropics from Europe to the Caribbean (∼ 40 to 8∘ N) from July to September 1991 (the Professor Zubov ship) and from Europe to south of the Equator (∼ 40∘ N to 8∘ S) between January and February 1992 (the Professor Vize ship). Our philosophy with the data recovery is to follow the same algorithms and parameters that appear in the two peer-reviewed articles that presented these datasets in the same issue of GRL in 1993, and here we provide all 48 lidar soundings made from the Professor Zubov and 11 of the 20 conducted from the Professor Vize, ensuring we have reproduced the aerosol backscatter and extinction values in the figures of those two papers. These original approaches used thermodynamic properties from the CIRA-86 standard atmosphere to derive the molecular backscattering, vertically and temporally constant values applied for the aerosol backscatter-to-extinction ratio, and the correction factor of the aerosol backscatter wavelength dependence. We demonstrate this initial validation of the recovered stratospheric aerosol extinction profiles, providing full details of each dataset in this paper's Supplement S1, the original profiles of backscatter ratio, and the calculated profiles of aerosol backscatter and extinction. We anticipate these datasets will provide potentially important new observational case studies for modelling analyses, including a 1-week series of consecutive soundings (in September 1991) at the same location showing the progression of the entrainment of part of the Pinatubo plume into the upper troposphere and the formation of an associated cirrus cloud. The Zubov lidar dataset illustrates how the tropically confined Pinatubo aerosol cloud transformed from a highly heterogeneous vertical structure in August 1991, maximum aerosol extinction values around 19 km for the lower layer and 23–24 for the upper layer, to a more homogeneous and deeper reservoir of volcanic aerosol in September 1991. We encourage modelling groups to consider new analyses of the Pinatubo cloud, comparing the recovered datasets, with the potential to increase our understanding of the evolution of the Pinatubo aerosol cloud and its effects. Data described in this work are available at https://doi.org/10.1594/PANGAEA.912770 (Antuña-Marrero et al., 2020).