HIAPER Pole-to-Pole Observations (HIPPO): fine-grained, global-scale measurements of climatically important atmospheric gases and aerosols

S. C. Wofsy

Abstract

The HIAPER Pole-to-Pole Observations (HIPPO) programme has completed three of five planned aircraft transects spanning the Pacific from 85° N to 67° S, with vertical profiles every approximately 2.2° of latitude. Measurements include greenhouse gases, long-lived tracers, reactive species, O2/N2 ratio, black carbon (BC), aerosols and CO2 isotopes. Our goals are to address the problem of determining surface emissions, transport strength and patterns, and removal rates of atmospheric trace gases and aerosols at global scales and to provide strong tests of satellite data and global models. HIPPO data show dense pollution and BC at high altitudes over the Arctic, imprints of large N2O sources from tropical lands and convective storms, sources of pollution and biogenic CH4 in the Arctic, and summertime uptake of CO2 and sources for O2 at high southern latitudes. Global chemical signatures of atmospheric transport are imaged, showing remarkably sharp horizontal gradients at air mass boundaries, weak vertical gradients and inverted profiles (maxima aloft) in both hemispheres. These features challenge satellite algorithms, global models and inversion analyses to derive surface fluxes. HIPPO data can play a crucial role in identifying and resolving questions of global sources, sinks and transport of atmospheric gases and aerosols.

1. Introduction

Future action on climate change requires us to be able to distinguish surface emissions of greenhouse gases (GHGs) directly associated with human activities from indirect effects of climate change and from natural variations. Measurements of concentrations of GHGs (e.g. CO2, CH4 and N2O) and tracers (e.g. CO, SF6, O2/N2) have been used to derive surface fluxes, and to attribute emissions to human activity or natural processes, in a variety of conceptual frameworks [16]. A common approach is to exercise atmospheric chemistry transport models (ACTMs) to simulate concentrations of pollutants in time and space over the globe, typically combining a priori surface emission flux fields, rates for chemical reactions in the atmosphere and assimilated meteorological fields. The surface flux field is then adjusted to obtain the best match between model output and data. The model is effectively ‘inverted’ to define optimal sources and sinks from atmospheric concentrations.

Global-scale ACTMs have difficulty resolving sharp chemical gradients, such as exist at air mass boundaries, at the tropopause [7,8] or in plumes of tracers emanating from strong source regions, a limitation inherent in the numerical methods used in the Eulerian model framework [9]. ACTM outputs are commonly compared with datasets that likewise do not resolve sharp gradients or define fine-scale variance and covariance, e.g. data from sparse ground stations or satellites. These data are in some respects ideal for comparison with global simulations, since they may be presumed to represent large air masses and global-scale processes. But, because accurate knowledge of vertical and horizontal gradients is available neither from global data nor from global models, the fine-grained structure of atmospheric tracer distributions is unknown and we have no way to assess how important this structure might be for inverse models.

Air mass boundaries move in space and time. Spatially aggregated data, or measurements from widely spaced locations (e.g. remote island stations), give horizontal patterns with diffuse boundaries. This fuzzy picture might be correct if data are averaged over time (e.g. for a month), but the pattern is never actually present in the atmosphere. A model may reproduce these smoothed gradients with simulated rates of transport across air mass boundaries that are significantly in error, leading to incorrect surface fluxes inferred from an inverse analysis. Similarly, if the vertical gradient in a model is systematically in error, the atmospheric mass burden (or column integral) is improperly simulated, and significant bias can result when the model is inverted to obtain optimal surface fluxes (e.g. [10]).

Surface fluxes obtained by inverting global models are also subject to systematic errors derived from other elements of the model framework, in addition to transport errors. For example, misidentification of source processes and locations, or incorrectly specified temporal variation of a priori fluxes, introduces errors into ‘optimal’ flux fields associated directly with model constructs. Since generally the only adjustable parameters in an inverse model are surface fluxes, strong, independent tests of all elements of the model framework are needed, using atmospheric measurements, in order to reduce the impact of these diverse systematic errors and biases.

This paper introduces the data obtained by the HIAPER Pole-to-Pole Observations (HIPPO) project, a sequence of five global aircraft measurement programmes that sample the atmosphere from (almost) the North Pole to the coastal waters of Antarctica, from the surface to 14 km, spanning the seasons. Three have been completed to date—January and November 2009, April 2010—to be followed by June and August/September 2011 (denoted HIPPO-1, -2, etc.). This paper will focus primarily on results obtained in HIPPO-1, for which data have been finalized.

HIPPO is intended to obtain global-scale, fine-grained data for the first time, for a large number of atmospheric constituents. HIPPO data provide high-resolution, systematic, pole-to-pole curtains with deeper vertical extent than possible in the past, acquired at high temporal resolution to preserve correlations among constituents. The goal is to provide new perspectives on how to use data to wring out transport issues in ACTMs, to identify and address the errors in the spatial/temporal representation of surface fluxes and reaction rates, to discover new features of trace gas and aerosol emissions and global-scale transport and to rigorously test satellite algorithms. A related goal is to elucidate the underlying processes and controls on surface emission fluxes, in order to be able to predict emissions and atmospheric concentrations in the future.

Fine-grained data have been obtained previously from aircraft platforms, but only a few datasets have long transects covering extended time periods, principally from commercial airliners that acquire data exclusively at cruise altitude plus vertical profiles entering/leaving large urban complexes [1113]. High-resolution datasets spanning the depth of the atmosphere over large areas are available from a few surveys [1417] that covered much smaller domains and seasonal intervals than HIPPO. Thus, HIPPO data provide hitherto unavailable information for testing and refining understanding of global distributions and surface fluxes of trace gases and aerosols.

2. HIAPER Pole-to-Pole Observations payload

HIPPO measurements are summarized in table 1. The platform is the National Science Foundation’s Gulfstream V (GV or HIAPER) aircraft operated by the National Center for Atmospheric Research (NCAR). Major GHGs (CO2, CH4, N2O) and other important trace species (CO, SF6, O2/N2, H2, etc.) were measured at high frequency, many at 1 Hz, with two (or more) independent measurements for each to provide redundancy (five for CO2), check calibration and assess sensor drift. The quantum cascade laser spectrometer (QCLS), which measures CO2, CO, CH4 and N2O, is a mid-infrared (IR) sensor developed by Harvard University and Aerodyne Corp., owned by NCAR, and operated during HIPPO by the Harvard team; the Observations of the Middle Stratosphere (OMS) is a Harvard CO2 sensor using an IR gas analyser (IRGA), which has logged more than 300 flights on airborne platforms [18]. The Research Aviation Facility CO monitor (RAF-CO) refers to NCAR’s AeroLaser AL5002 vacuum–UV sensor. Data at 1 Hz were obtained in two independent measurements for O3 and for H2O (by the Vertical Cavity Surface Emitting Laser (VCSEL), a new open path near-IR multi-pass spectrometer [19], and by the Unmanned Aircraft Systems (UAS) Chromatograph for Atmospheric Trace Species (UCATS)-H2O, a National Oceanic and Atmospheric Administration (NOAA) mid-IR spectrometer with a sampling inlet). AO2 is the NCAR Airborne Oxygen Instrument, which measures O2 using a vacuum–UV absorption technique [20] and CO2 using a single-cell IRGA (Licor 850). UCATS and the PAN (PeroyAcylNitrate) and other Trace Hydrohalocarbon ExpeRiment (PANTHER) are on-board gas chromatographs with extensive flight histories, measuring at 1–3 min intervals.

View this table:
Table 1.

Measurements and species on the GV in HIPPO. Bold font indicates species with three or more measurements, sampling rates in parentheses.

A comprehensive suite of GHGs and a diverse ensemble of halocarbons, hydrocarbons and sulphur species were measured by the whole-air sampler system (WAS; [8]), using three different collection systems: stainless steel electropolished cylinders (the Advanced Whole Air Sampler (AWAS) from the University of Miami), glass vessels from the NOAA Programmable Flask Package (PFP) programme (NWAS from NOAA) and special glass vessels for measuring O2/N2 and Ar/N2 ratios and the isotopic composition of CO2 (Multiple Enclosure Device for Unfractioned Sampling of Air (MEDUSA) from the University of California at San Diego).

BC was measured at 1 Hz by the NOAA single-particle soot photometer (SP2) [21] on air drawn into the cabin. The SP2 uses a laser beam to heat the BC component of individual fine-mode aerosol particles to vaporization, resulting in emission of thermal radiation that is interpreted to give the BC mass concentration and size distribution. Other particle probes, including a cloud droplet probe (CDP), the NCAR Ultra High Sensitivity Aerosol Sampler (UHSAS), the 2D-C particle imager and the PMS liquid water sensor (King) (Precipitable Liquid Water Content (PLWC)), resided on the wing, along with the microwave temperature profiler (MTP), which measures the temperature profile above and below the plane, and a digital camera that acquired an image every second.

3. Results

Figure 1 shows the flight tracks for HIPPO-1 and -2, plus vertical sampling and observed isentropic structure for HIPPO-1. Each deployment obtained a southbound cross section of the atmosphere near the International Date Line, covering the Arctic Ocean, Alaska and the entire Pacific to just short of Antarctica. Northbound, HIPPO-1 traversed the South Pacific to sample in the eastern tropical Pacific upwelling zone, and HIPPO-2 sampled the western Pacific warm pool. We obtained 138, 151 and 136 vertical profiles in HIPPO-1, -2 and -3, respectively, spanning 67° S–80° N. On average, consecutive samples in the mid-troposphere are separated by 2.2° of latitude, with 4.4° between consecutive near-surface or high-altitude samples. Most profiles extended from approximately 300 to 8500 m altitude, constrained by air traffic, but significant profiling extended above approximately 14 km (figure 1c).

Figure 1.

(a,b) Locations of flight tracks and vertical profiles (red points) for HIPPO deployments 1 and 2. (c) Vertical profiles in the HIPPO-1 cross section of the Pacific southbound near the date line (flights 2–7), tropopause heights (pressure altitude (km); cyan) from MTP, and stratospheric flight segments in red. (d) Cross section of potential temperature (Θ) in HIPPO-1, on the southbound leg near the date line. The white dotted lines mark the flight path of the GV, and grey lines show contours of potential temperature. Altitudes are given in metres above sea level (m.a.s.l.) (GPS). (Online version in colour.)

Figure 1d shows the cross section of potential temperature (Θ) versus latitude and GPS altitude (m.a.s.l.) during HIPPO-1. A dome of very cold, stable air, with a strong latitude gradient of Θ, covered the Arctic, with a sharp transition at the northern edge of the polar jet, near 60° N. There was a narrow zone of stratospheric influx just equatorward of 60° N, and then a belt of much weaker vertical stratification from 60° to 40° N, the region of the ‘warm conveyor belt’, where vertical transport occurs in association with jet stream dynamics and mid-latitude storms [22,23]. Equatorward of 40° N, there was a broad Θ bulge in the middle and upper troposphere, reflecting the influence of deep convection and the Hadley circulation. The Southern Hemisphere (SH) had similar structure, but the features were weaker, reflecting the summer season.

Figure 2.

Cross sections of (a) CH4 (ppb), (b) CO (ppb), (c) CO2 (ppm), (d) SF6 (ppt), (e) N2O (ppb) and (f) H2O (log10 (ppm)) on HIPPO-1, southbound along the date line, January 2009. White dotted lines show flight tracks, and grey contours show potential temperature. SF6 data represent a composite of PANTHER and UCATS data, and H2O data are from VCSEL. Altitudes are given in m.a.s.l. (GPS). (Online version in colour.)

Figure 2 shows cross sections along the date line for key species in January 2009 (HIPPO-1). Atmospheric tracers displayed sharp transitions in the horizontal, demarcating the boundaries between polar, middle latitude, subtropical and tropical air masses, and at the intertropical convergence zone (ITCZ; the persistent line of convective storms that circles the globe near the equator, demarcating the boundary between the hemispheres). Gradients were much weaker in the vertical than in the horizontal, except in the northern polar region where fresh inputs of CO, CO2, CH4 and SF6 were evidenced by high values near the surface and the presence of very short-lived pollutants (e.g. benzene). According to model simulations (discussed further below), these pollutants originated from Asia and Europe. In the Arctic, tracer isopleths broadly aligned with isentropes, with sporadic anomalies owing to proximate emissions or to stratospheric influence. Convective transport across isentropes was evident elsewhere, in mid-latitudes, tropics and subtropics, particularly for the longest lived tracers (e.g. SF6 and N2O). The peak value of SF6 over the North Pacific was observed in the middle troposphere. Concentrations of N2O showed a completely different pattern from other GHGs, with a prominent bulge in the middle and upper tropical and subtropical troposphere, plus sporadic enhancements between 10 and 14 km. Similar features were seen for N2O in all three missions flown so far.

Water vapour showed the expected increases with temperature and tropical origin. Input of very wet air was notable in the South Pacific convergence zone (SPCZ; 10–20° S). The SPCZ is a persistent frontal formation of widespread cloud cover and precipitation extending in a southeast direction from New Guinea into the SH mid-latitudes, with strong convection typically extending over a large area. It is specially prominent in boreal winter [24]. The SPCZ is spatially larger, and contains more intense convection, than similar convergence zones elsewhere, such as the South Indian Convergence Zone (SICZ) and the South Atlantic Convergence Zone (SACZ), although it is not clear why [24]. Analogous, but less distinct, anomalies have been observed for other chemical tracers near the SPCZ [25,26]. In each of the three HIPPO missions so far, the water vapour plume from the SPCZ has been larger in extent, and penetrated deeper into the atmosphere, than in the ITCZ. These results suggest that the SPCZ may be an important global source of water vapour to the tropical tropopause layer over the Pacific, at least in the boreal winter season, supporting inferences from observations of cirrus clouds over the tropical western Pacific during that season [27].

Figure 3.

(a,b) Cross sections of CO2 (ppm) and CO (ppb), on HIPPO-2, southbound along the date line, November 2009. Grey contours show potential temperature. (c) Photo of dense layer of dark aerosols, looking north at 80° N, 8 km altitude, 2 November 2009 (photo: E. Kort). (d) Vertical profiles of black carbon (BC) and CO at 77.2° N, on HIPPO-2, 2 November 2009. Black lines, BC; red lines, CO. (e) Vertical profiles of CO2, CH4 and N2O at 77.2° N, on HIPPO-2, 2 November 2009. Black lines, CH4; green lines, CO2; blue lines, N2O. (Online version in colour.)

Figure 3 shows cross sections for HIPPO-2 in November 2009. Inverted vertical profiles (i.e. with higher concentrations aloft) were surprisingly strong. Pollutant levels typical of urban corridors were found in a thick layer at 6–8 km over an extensive region of the Arctic, overlying clean air. The layers preserve coherence (figure 3d,e) among tracers despite sharp gradients and transport over thousands of kilometres, with strong vertical uplift. Very high levels of BC aerosol were present, and absorption of solar radiation was visible (figure 3c). High aerosol concentrations indicate the absence of efficient precipitation scavenging, suggesting that polluted air was lofted by isentropic advection over the cold dome that was in the process of developing over the Arctic. Notable enhancements of N2O were observed in some of the plumes. In many descents, enhanced CH4 concentrations were observed near the surface of the Arctic Ocean, sometimes in otherwise pristine air, possibly signifying emission from biogenic sources or from CH4 hydrates (e.g. figure 3c, below 2000 m). Horizontal gradients were again very sharp at air mass boundaries.

Figure 4 illustrates the meridional gradients of CO2 and O2/N2. Concentrations of CO2 were enhanced above the surface (but not at the surface) from 45° S to 70° S in January (figures 2c and 4a), which might confirm a signal of elevated CO2 in this region similar to that reported by the Atmospheric Infrared Sounder (AIRS) satellite [28]; the structure of the CO2 distribution is rather complex, with probable contributions from multiple processes including fires in Australia and uptake of CO2 at the ocean surface. There is excess O2 near the surface all the way from 20° S to 60° S (figure 4b), reflecting inputs to the atmosphere from both the ocean and the land, including thermal degassing of O2 from the ocean and biological production of O2 from land and ocean in summer.

Figure 4.

(a) Meridional gradients of CO2 averaged between 0.3 and 1.5 km and between 3.5 and 6.5 km. (b) Variation of the O2:N2 ratio, as in (a). Green plus symbols, 0.3–1.5 km; open circles, 3.5–6.5 km. (Online version in colour.)

Preliminary simulations of HIPPO-1 data have been made with global models including Goddard Earth Observing System (GEOS) Chemical Tracer Model (CHEM) [1], ACTM [6] and Global and regional Earth-system Monitor using Satellite and in situ data/Monitoring and Forecasting Atmospheric Composition (GEMS/MACC) [29]. Examples are shown in figure 5ad. The models extracted computed values along the flight track for detailed comparison with HIPPO. Model simulations performed reasonably well in many respects, as can be seen by comparing with figure 2. Simulations of CO2 from the GEOS-CHEM model ([30]; figure 5a) captured many features of the observations (figure 2), including the relatively well-mixed profiles in subpolar and middle latitudes. Vertical contrasts occurred mostly over the Arctic and Antarctic/Southern Ocean, reflecting inputs from fossil fuels and marine biological uptake, respectively. But horizontal gradients were much sharper in the atmosphere, and vertical gradients weaker, than in the models. Most satellite data likewise do not resolve the observed gradients. Enhanced CO2 poleward of 40° S seen by HIPPO is not reproduced in the model; the major pollution signature appears at 60° N instead of 80° N; enhanced surface CO2 in the SH deep tropics, probably from equatorial upwelling, is not simulated in the model.

Figure 5.

Model simulations for HIPPO-1, southbound along the date line, January 2009. (a,b) Cross sections of CO2 (ppm) and SF6 (ppt) from the GEOS-CHEM and ACTM models, respectively. (c) Vertical profiles of N2O in the subtropics (ACTM and observed in HIPPO-1). Blue filled squares, HIPPO-1; black open circles, ACTM. (d) Cross section for N2O from the ACTM for the HIPPO-1 flight track (flights 2–7). (Online version in colour.)

Both the ACTM and GEOS-CHEM models appear to have vertical transport rates that are too weak at high latitudes (compare CO2, SF6 and CO in figures 25), both along and across isentropes. The GEMS/MACC model gives a somewhat better appearance for vertical profiles of CO and CH4 (not shown), but has anomalies elsewhere and misses global-scale North–South gradients significantly: GEMS/MACC predicts 70 ppb for ΔCO between the Antarctic and the Arctic, compared with approximately 95 ppb observed.

HIPPO data for N2O differ radically from model results. Simulations completely miss the major feature of the global distribution, the tropical/sub-tropical maxima at altitude (figure 5c,d), with greater enhancement north of the ITCZ. This pattern has been observed in all six HIPPO cross sections flown to date, although the magnitude of the enhancement aloft was largest in HIPPO-1 southbound. HYSPLIT [31] trajectories showed that air near the surface originated from the east, whereas the air with enhanced N2O concentrations, above 2 km, came from the west, suggesting sources in regions of convection in the western Pacific and/or tropical South Asia or Indonesia. Samples collected on the CARIBIC flights of commercial aircraft from Germany to India [32] also showed diffuse maxima of N2O in the upper troposphere (8.5–12 km) during June–August, peaking in the same general latitude range as in HIPPO flights (15–25° N). It appears that much of the excess N2O seen by HIPPO was generated in the tropics, but it is difficult to distinguish land-based sources from possible production by corona discharge or lightning in convective storms [33] or by other mechanisms.

Previous studies [34, 35] using sparse station data have already suggested that unaccounted tropical sources were needed to close the N2O budget. HIPPO data strongly support that view. Elevated concentrations high in the atmosphere found in HIPPO indicate a notably stronger tropical source than could be inferred from surface data only. Even in the north polar region, HIPPO data often showed maximum concentrations of N2O aloft (cf. figures 2, 3 and 5), emphasizing the difficulty of using surface measurements for modelling of this important species.

4. Summary and conclusions

HIPPO flights provide a unique dataset for atmospheric research: simultaneous, fine-grained, global, high-frequency measurements of major GHGs, tracers with a wide range of chemical lifetimes and diverse source processes, O2/N2 and Ar/N2 ratios, BC and aerosols, plus the isotopic composition of CO2. Concentration data from HIPPO are linked to world standards, and multiple measurements provide checks on data quality. The data from all missions will be publicly available within 12–18 months of collection. HIPPO data have already played a major role in calibrating total column measurements of CO2, CO and CH4 from the Total Carbon Column Observing Network (TCCON) of Fourier transform spectrometers [36].

First impressions already highlight new phenomena: dense pollution high over the Arctic in late autumn/early winter, with a notable component of BC; the imprint of large N2O sources in the tropical or subtropical areas of Asia or the western Pacific and sources of CH4 in the Arctic both from fossil-fuel extraction and from non-industrial sources. Although not discussed here, we also observed short-lived gases emanating from various marine environments across the Pacific (e.g. methyl nitrate, haloforms) and from major industrial areas.

HIPPO data clearly delineate the atmospheric imprint of summertime sinks for CO2 and sources for O2 in the SH, providing a new tool to infer the strength of ocean and land fluxes in the carbon cycle. Curtain plots reveal large-scale features invisible to satellites and blurred by models, revealing the internal structure of global CO2 distributions hitherto observed fuzzily by satellites and sparse remote stations.

The most novel aspect of HIPPO is the imaging of the signatures of global atmospheric transport modes with high clarity. From the SPCZ and ITCZ to the cold dome over the winter pole, HIPPO data show the influence of convection, isentropic transport and stratosphere–troposphere exchange, with fine resolution at global scale. Remarkably sharp gradients in tracer concentrations are observed at air mass boundaries. Inverted tracer gradients, with maxima aloft, were found to be common in both hemispheres over the Pacific. These features present major challenges to global models and to observation systems using sparse surface stations and satellites. If not represented in global inverse model studies, they can lead to major biases and inconsistencies. HIPPO data can play a crucial role in identifying and resolving these issues.

Acknowledgements

The HIPPO Programme was supported by NSF grants ATM-0628575, ATM-0628519 and ATM-0628388 to Harvard University, University of California (San Diego), University Corporation for Atmospheric Research, University of Colorado/CIRES and by the NCAR. The NCAR is supported by the National Science Foundation. Participation by several instruments (SP-2, Ozone, UCATS, PANTHER, NWAS flasks), and weather forecasting, were supported by offices and programmes of the National Oceanic and Atmospheric Administration: the Atmospheric Composition and Climate Programme, the Office of Oceanic and Atmospheric Research and the Environmental Research Laboratory. The AWAS flasks system was supported by NSF grants NSF ATM0849086 and AGS0959853 to the University of Miami. VCSEL was supported by NSF grant AGS-1036275 to Princeton University. We are grateful to the crew of the GV for their dedication and professional skill in making the flights possible and in taking the GV to places not previously visited by a jet aircraft, and to the NCAR Earth Observing Laboratory for support of logistics and public outreach. The authors gratefully acknowledge the NOAA Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and dispersion model.

Footnotes

  • S. C. Wofsy, B. C. Daube, R. Jimenez, E. Kort, J. V. Pittman, S. Park, R. Commane, B. Xiang, G. Santoni, D. Jacob, J. Fisher, C. Pickett-Heaps, H. Wang, K. Wecht, Q.-Q. Wang (School of Engineering and Applied Science, Harvard University, Cambridge, MA, USA); B. B. Stephens, S. Shertz, P. Romashkin, T. Campos, J. Haggerty, W. A. Cooper, D. Rogers, S. Beaton, R. Hendershot (National Center for Atmospheric Research, Boulder, CO, USA); J. W. Elkins, D. W. Fahey, R. S. Gao, F. Moore, S. A. Montzka, J. P. Schwarz, D. Hurst, B. Miller, C. Sweeney, S. Oltmans, D. Nance, E. Hintsa, G. Dutton, L. A. Watts, J. R. Spackman, K. H. Rosenlof, E. A. Ray, B. Hall (NOAA ESRL and CIRES, Boulder, CO, USA); M. A. Zondlo, M. Diao (Department of Civil and Environmental Engineering, Princeton University, Princeton, NJ, USA); R. Keeling, J. Bent (Scripps Institution of Oceanography, University of California at San Diego, CA, USA); E. L. Atlas, R. Lueb (Rosenstiel School of Marine and Atmospheric Science, University of Miami, Miami, FL, USA); M. J. Mahoney, M. Chahine, E. Olsen (Jet Propulsion Laboratory, Pasadena, CA, USA); P. Patra, K. Ishijima (Research Institute for Global Change, JAMSTEC, Yokohama, Japan); R. Engelen, J. Flemming (European Centre for Medium-Range Weather Forecasts, Reading, UK); R. Nassar, D. B. A. Jones (Department of Physics, University of Toronto, Toronto, Canada); S. E. Mikaloff Fletcher (National Institute for Water and Atmospheric Research, Wellington, New Zealand).

  • One contribution of 17 to a Discussion Meeting Issue ‘Greenhouse gases in the Earth system: setting the agenda to 2030’.

References

View Abstract