Understanding the magnitude of tropical gross primary production (GPP) is critical for carbon cycle modeling and climate projections, but this quantity is poorly constrained at regional scales. ...Biospheric uptake of carbonyl sulfide (OCS) provides a promising approach to estimating regional GPP. Here, we simulate OCS concentrations driven by surface flux scenarios encompassing a wide range of GPP estimates for the Amazon basin. We compare the model output to satellite retrievals and find a regional GPP estimate of 2375 ± 914 g(C) m−2 yr−1, consistent with previous estimates, including the TRENDY model ensemble range of 1431–3812 g(C) m−2 yr−1.
Plain Language Summary
Carbonyl sulfide is an atmospheric gas that can be used to estimate how much carbon plants assimilate during photosynthesis. One of the most important regions for understanding the carbon cycle is the Amazon rainforest. In order to determine how much carbon is being taken up in the Amazon, we use models of the chemistry and movement of the atmosphere to estimate the uptake of carbonyl sulfide and compare our model output to satellite observations. Our research shows that estimates of Amazonian GPP using carbonyl sulfide uptake are consistent with ecosystem models and other types of observations. This builds confidence in our understanding of carbon uptake in the Amazon, which has major implications for climate predictions.
Key Points
Remotely sensed atmospheric carbonyl sulfide can be used to estimate surface processes, namely carbon uptake through photosynthesis
Models of carbonyl sulfide indicate that gross primary productivity in the Amazon basin closely follows lower‐productivity ecosystem models
High productivity and deep convection over the Amazon enable the surface uptake signal to be picked up in upper tropospheric observations
We present a global carbonyl sulfide (OCS) data set covering the period June 2002 to April 2012, derived from FTIR (Fourier transform infrared) limb emission spectra measured with the Michelson ...Interferometer for Passive Atmospheric Sounding (MIPAS) on the ENVISAT satellite. The vertical resolution is 4–5 km in the height region 6–15 km and 15 at 40 km altitude. The total estimated error amounts to 40–50 pptv between 10 and 20 km and to 120 pptv at 40 km altitude. MIPAS OCS data show no systematic bias with respect to balloon observations, with deviations mostly below ±50 pptv. However, they are systematically higher than the OCS volume mixing ratios of the ACE-FTS instrument on SCISAT, with maximum deviations of up to 100 pptv in the altitude region 13–16 km. The data set of MIPAS OCS exhibits only moderate interannual variations and low interhemispheric differences. Average concentrations at 10 km altitude range from 480 pptv at high latitudes to 500–510 pptv in the tropics and at northern mid-latitudes. Seasonal variations at 10 km altitude amount to up to 35 pptv in the Northern and up to 15 pptv in the Southern Hemisphere. Northern hemispheric OCS abundances at 10 km altitude peak in June in the tropics and around October at high latitudes, while the respective southern hemispheric maxima were observed in July and in November. Global OCS distributions at 250 hPa (∼ 10–11 km) show enhanced values at low latitudes, peaking during boreal summer above the western Pacific and the Indian Ocean, which indicates oceanic release. Further, a region of depleted OCS amounts extending from Brazil to central and southern Africa was detected at this altitude, which is most pronounced in austral summer. This depletion is related to seasonally varying vegetative uptake by the tropical forests. Typical signatures of biomass burning like the southern hemispheric biomass burning plume are not visible in MIPAS data, indicating that this process is only a minor source of upper tropospheric OCS. At the 150 hPa level (∼ 13–14 km) enhanced amounts of OCS were also observed inside the Asian monsoon anticyclone, but this enhancement is not especially outstanding compared to other low latitude regions at the same altitude. At the 80 hPa level (∼ 17–18 km), equatorward transport of mid-latitude air masses containing lower OCS amounts around the summertime anticyclones was observed. A significant trend could not be detected in upper tropospheric MIPAS OCS amounts, which points to globally balanced sources and sinks. Simulations with the ECHAM-MESSy model reproduce the observed latitudinal cross sections fairly well.
Peroxyacetyl nitrate (PAN) is an important chemical species in the troposphere as it aids the long-range transport of NOx and subsequent formation of O3 in relatively clean remote regions. Over the ...past few decades observations from aircraft campaigns and surface sites have been used to better understand the regional distribution of PAN. However, recent measurements made by satellites allow for a global assessment of PAN in the upper troposphere–lower stratosphere (UTLS). In this study, we investigate global PAN distributions from two independent retrieval methodologies, based on measurements from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) instrument, on board Envisat from the Institute of Meteorology and Climate Research (IMK), Karlsruhe Institute of Technology, and the Department of Physics and Astronomy, University of Leicester (UoL). Retrieving PAN from MIPAS is challenging due to the weak signal in the measurements and contamination from other species. Therefore, we compare the two MIPAS datasets with observations from the Atmospheric Chemistry Experiment Fourier transform spectrometer (ACE-FTS), in situ aircraft data and the 3-D chemical transport model TOMCAT. MIPAS shows peak UTLS PAN concentrations over the biomass burning regions (e.g. ranging from 150 to > 200 pptv at 150 hPa) and during the summertime Asian monsoon as enhanced convection aids the vertical transport of PAN from the lower atmosphere. At 150 hPa, we find significant differences between the two MIPAS datasets in the tropics, where IMK PAN concentrations are larger by 50–100 pptv. Comparisons between MIPAS and ACE-FTS show better agreement with the UoL MIPAS PAN concentrations at 200 hPa, but with mixed results above this altitude. TOMCAT generally captures the magnitude and structure of climatological aircraft PAN profiles within the observational variability allowing it to be used to investigate the MIPAS PAN differences. TOMCAT–MIPAS comparisons show that the model is both positively (UoL) and negatively (IMK) biased against the satellite products. These results indicate that satellite PAN observations are able to detect realistic spatial variations in PAN in the UTLS, but further work is needed to resolve differences in existing retrievals to allow quantitative use of the products.
Ozone loss in the lower and middle stratosphere in spring and summer, in particular over polar regions, is driven mainly by halogens and nitrogen oxides (NOx). Whereas the stratospheric chlorine ...levels are expected to decrease in the future, the role of NOx for the O3 budget in a changing climate is not well quantified. Here we combine satellite measurements and model simulations to diagnose the accumulated O3 loss during winter and spring 2002–2003 in the Arctic polar stratosphere. We show that in a winter stratosphere strongly disturbed by warmings, O3 loss processes driven by halogens and NOx can significantly overlap within the polar column and become comparable in magnitude even if a significant, halogen‐induced O3 loss has occurred. Whereas, until the beginning of March 2003, polar column O3 loss was mainly caused by the halogen chemistry within the vortex at an altitude around 18 km, the chemical O3 destruction in March and April was dominated by the NOx chemistry in O3‐rich air masses transported from the subtropics and mixed with the polar air above the region affected by the halogens. This NOx‐related O3 loss started around mid‐December 2002 in subtropical air masses above 30 km that moved poleward after the major warming in January, descended to 22 km with an increasing magnitude of O3 loss and reached surprisingly high values of up to 50% local loss around the end of April. To some extent, the NOx‐driven O3 loss was enhanced by mesospheric air trapped in the vortex at the beginning of the winter as a layer of few km in the vertical and transported downward within the vortex. The effect of NOx transported from the subtropics dominated the O3 loss processes in the polar stratosphere in spring 2003, both relative to the effect of the halogens and relative to the contribution of the mesospheric NOx sources. A comparison with the 1999/2000 Arctic winter and with the Antarctic vortex split event in 2002 shows that wave events triggered by stratospheric warmings may significantly enhance O3 loss driven by NOx when O3‐ and NOx‐rich air masses from the subtropics are transported poleward and are mixed with the vortex air.
Carbonyl sulfide (COS) has the potential to be used as a climate diagnostic due to its close coupling to the biospheric uptake of CO2 and its role in the formation of stratospheric aerosol. The ...current understanding of the COS budget, however, lacks COS sources, which have previously been allocated to the tropical ocean. This paper presents a first attempt at global inverse modelling of COS within the 4-dimensional variational data-assimilation system of the TM5 chemistry transport model (TM5-4DVAR) and a comparison of the results with various COS observations. We focus on the global COS budget, including COS production from its precursors carbon disulfide (CS2) and dimethyl sulfide (DMS). To this end, we implemented COS uptake by soil and vegetation from an updated biosphere model (Simple Biosphere Model – SiB4). In the calculation of these fluxes, a fixed atmospheric mole fraction of 500 pmol mol−1 was assumed. We also used new inventories for anthropogenic and biomass burning emissions. The model framework is capable of closing the COS budget by optimizing for missing emissions using NOAA observations in the period 2000–2012. The addition of 432 Gg a−1 (as S equivalents) of COS is required to obtain a good fit with NOAA observations. This missing source shows few year-to-year variations but considerable seasonal variations. We found that the missing sources are likely located in the tropical regions, and an overestimated biospheric sink in the tropics cannot be ruled out due to missing observations in the tropical continental boundary layer. Moreover, high latitudes in the Northern Hemisphere require extra COS uptake or reduced emissions. HIPPO (HIAPER Pole-to-Pole Observations) aircraft observations, NOAA airborne profiles from an ongoing monitoring programme and several satellite data sources are used to evaluate the optimized model results. This evaluation indicates that COS mole fractions in the free troposphere remain underestimated after optimization. Assimilation of HIPPO observations slightly improves this model bias, which implies that additional observations are urgently required to constrain sources and sinks of COS. We finally find that the biosphere flux dependency on the surface COS mole fraction (which was not accounted for in this study) may substantially lower the fluxes of the SiB4 biosphere model over strong-uptake regions. Using COS mole fractions from our inversion, the prior biosphere flux reduces from 1053 to 851 Gg a−1, which is closer to 738 Gg a−1 as was found by Berry et al. (2013). In planned further studies we will implement this biosphere dependency and additionally assimilate satellite data with the aim of better separating the role of the oceans and the biosphere in the global COS budget.
Carbonyl sulfide (COS), a trace gas in our atmosphere that leads to the formation of aerosols in the stratosphere, is largely taken up by terrestrial ecosystems. Quantifying the biosphere uptake of ...COS could provide a useful quantity to estimate gross primary productivity (GPP). Some COS sources and sinks still contain large uncertainties, and several top-down estimates of the COS budget point to an underestimation of sources, especially in the tropics. We extended the inverse model TM5-4DVAR to assimilate Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) satellite data, in addition to National Oceanic and Atmospheric Administration (NOAA) surface data as used in a previous study. To resolve possible discrepancies among the two observational data sets, a bias correction scheme is necessary and implemented. A set of inversions is presented that explores the influence of the different measurement streams and the settings of the prior fluxes. To evaluate the performance of the inverse system, the HIAPER Pole-to-Pole Observations (HIPPO) aircraft observations and NOAA airborne profiles are used. All inversions reduce the COS biosphere uptake from a prior value of 1053 GgS a−1 to much smaller values, depending on the inversion settings. These large adjustments of the biosphere uptake often turn parts of Amazonia into a COS source. Only inversions that exclusively use MIPAS observations, or strongly reduce the prior errors on the biosphere flux, maintain the Amazon as a COS sink. Inclusion of MIPAS data in the inversion leads to a better separation of land and ocean fluxes. Over the Amazon, these inversions reduce the biosphere uptake from roughly 300 to 100 GgS a−1, indicating a strongly overestimated prior uptake in this region. Although a recent study also reported reduced COS uptake over the Amazon, we emphasise that a careful construction of prior fluxes and their associated errors remains important. For instance, an inversion that gives large freedom to adjust the anthropogenic and ocean fluxes of CS2, an important COS precursor, also closes the budget satisfactorily with much smaller adjustments to the biosphere. We achieved better characterisation of biosphere prior and uncertainty, better characterisation of combined ocean and land fluxes, and better constraint of both by combining surface and satellite observations. We recommend more COS observations to characterise biosphere and ocean fluxes, especially over the data-poor tropics.
To avoid unnecessary data traffic it is sometimes desirable to apply mean averaging kernels to mean profiles of atmospheric state variables. Unfortunately, application of averaging kernels and ...averaging are not commutative in cases when averaging kernels and state variables are correlated. That is to say, the application of individual averaging kernels to individual profiles and subsequent averaging will, in general, lead to different results than averaging of the original profiles prior to the application of the mean averaging kernels, unless profiles and averaging kernels are fully independent. The resulting error, however, can be corrected by subtraction of the covariance between the averaging kernel and the vertical profile. Thus, it is recommended to calculate the covariance profile along with the mean profile and the mean averaging kernel.
Using the IMK–IAA data processor, methane and nitrous oxide distributions were retrieved from version-8 limb emission spectra recorded with the Michelson Interferometer for Passive Atmospheric ...Sounding (MIPAS). The dataset includes measurements from the nominal, upper troposphere–lower stratosphere, middle-atmosphere, upper-atmosphere and noctilucent-cloud observation modes. The processing differs from the previous version-5 data with respect to the atmospheric state variables that are jointly retrieved along with the target gases CH4 and N2O, the treatment of the radiance offset, the selection of microwindows, the regularization, the spectroscopic data used and the treatment of horizontal variability of the atmospheric state. Besides the regular data product, a coarse-grid representation of the profiles with unity averaging kernels is available, as well as a specific research product for middle-atmosphere measurements resulting from a slightly different retrieval approach. The CH4 errors are dominated by the large spectroscopic uncertainty for line intensities, which probably is too pessimistic, and estimated to be 21 %–34 % in the altitude range 6–68 km for northern midlatitude summer day conditions. The N2O errors are 7 %–17 % below 45 km. At higher altitudes they increase strongly due to nearly vanishing N2O amounts. Analysis of the horizontal averaging kernels reveals that for both gases the horizontal resolution is sampling-limited; i.e., information is not smeared over consecutive limb scans. Zonal-mean seasonal composites of both CH4 and N2O exhibit the typical distribution of source gases with strong upwelling in the tropics and subsidence above the winter poles. Comparison with the previous data version shows several improvements: first, the vertical resolution of the retrieved CH4 (N2O) profiles has generally been significantly enhanced and varies between 2.5 (2.5) and 4 (5) km at altitudes between 10 and 60 km, with the best resolution around 30 km for both species. Secondly, the number of non-converged retrievals has been clearly reduced, and thirdly, formerly strongly oscillating profiles are now considerably smoother.
Using the IMK-IAA data processor, methane and nitrous oxide distributions were retrieved from version-8 limb emission spectra recorded with the Michelson Interferometer for Passive Atmospheric ...Sounding (MIPAS). The dataset includes measurements from the nominal, upper troposphere-lower stratosphere, middle-atmosphere, upper-atmosphere and noctilucent-cloud observation modes. The processing differs from the previous version-5 data with respect to the atmospheric state variables that are jointly retrieved along with the target gases CH.sub.4 and N.sub.2 O, the treatment of the radiance offset, the selection of microwindows, the regularization, the spectroscopic data used and the treatment of horizontal variability of the atmospheric state. Besides the regular data product, a coarse-grid representation of the profiles with unity averaging kernels is available, as well as a specific research product for middle-atmosphere measurements resulting from a slightly different retrieval approach. The CH.sub.4 errors are dominated by the large spectroscopic uncertainty for line intensities, which probably is too pessimistic, and estimated to be 21 %-34 % in the altitude range 6-68 km for northern midlatitude summer day conditions. The N.sub.2 O errors are 7 %-17 % below 45 km. At higher altitudes they increase strongly due to nearly vanishing N.sub.2 O amounts. Analysis of the horizontal averaging kernels reveals that for both gases the horizontal resolution is sampling-limited; i.e., information is not smeared over consecutive limb scans. Zonal-mean seasonal composites of both CH.sub.4 and N.sub.2 O exhibit the typical distribution of source gases with strong upwelling in the tropics and subsidence above the winter poles. Comparison with the previous data version shows several improvements: first, the vertical resolution of the retrieved CH.sub.4 (N.sub.2 O) profiles has generally been significantly enhanced and varies between 2.5 (2.5) and 4 (5) km at altitudes between 10 and 60 km, with the best resolution around 30 km for both species. Secondly, the number of non-converged retrievals has been clearly reduced, and thirdly, formerly strongly oscillating profiles are now considerably smoother.
In this study, we present simultaneous airborne measurements of peroxyacetyl nitrate (PAN), ethane (C2H6), formic acid (HCOOH), methanol (CH3OH), and ethylene (C2H4) above the South Atlantic ...in September and October 2019. Observations were obtained from the Gimballed Limb Observer for Radiance Imaging of the Atmosphere (GLORIA), as two-dimensional altitude cross sections along the flight path. The flights were part of the SouthTRAC (Transport and Composition in the Southern Hemisphere Upper Troposphere/Lower Stratosphere) campaign with the German High Altitude and Long Range Research Aircraft (HALO). On two flights (8 September 2019 and 7 October 2019), large enhancements of all these substances were found between 7 and 14 km altitude with maximum volume mixing ratios (VMRs) of 1000 pptv for PAN, 1400 pptv for C2H6, 800 pptv for HCOOH, 4500 pptv for CH3OH, and 200 pptv for C2H4. One flight showed a common filamentary structure in the trace gas distributions, while the second flight is characterized by one large plume. Using backward trajectories, we show that measured pollutants likely reached upper troposphere and lower stratosphere (UTLS) altitudes above South America and central Africa, where elevated PAN VMRs are visible at the surface layer of the Copernicus Atmosphere Monitoring Service (CAMS) model during the weeks before both measurements. In comparison to results of the CAMS reanalysis interpolated onto the GLORIA measurement geolocations, we show that the model is able to reproduce the overall structure of the measured pollution trace gas distributions. For PAN, the absolute VMRs are in agreement with the GLORIA measurements. However, C2H6 and HCOOH are generally underestimated by the model, while CH3OH and C2H4, the species with the shortest atmospheric lifetimes of the pollution trace gases discussed, are overestimated by CAMS. The good agreement between model and observations for PAN suggests that the general transport pathways and emissions locations are well captured by the model. The poorer agreement for other species is therefore most likely linked to model deficiencies in the representation of loss processes and emission strength.