Methane fluxes from Zambian tropical wetlands

Author(s):  
Jacob Shaw ◽  
Grant Allen ◽  
Patrick Barker ◽  
Joseph Pitt ◽  
James Lee ◽  
...  

<p>Airborne measurements of methane (CH<sub>4</sub>) were recorded over three major wetland areas in Zambia in February 2019 during the MOYA (Methane Observations and Yearly Assessments) ZWAMPS field campaign. Enhancements of up to 600 ppb CH<sub>4</sub> were measured over the Bangweulu (11°36’ S, 30°05’ E), Kafue (15°43’ S, 27°17’ E), and Lukanga (14°29’ S, 27°47’ E) wetlands. Three independent methods were used to quantify methane emission fluxes; aircraft mass balance, aircraft eddy covariance, and atmospheric inversion modelling. Results yielded methane emission fluxes of 10-20 mg CH<sub>4</sub> m<sup>-2</sup> hr<sup>-1</sup>, which were up to an order of magnitude greater than the emission fluxes simulated by various wetland process models (WetCHARTs ensemble and LPX-Bern). Independent column CH<sub>4</sub> observations from the TROPOMI instrument were used to verify these fluxes, and investigate their applicability for timescales longer than the duration of the MOYA flight campaign.</p>

2020 ◽  
Author(s):  
Bruno Neininger ◽  
Jorg M. Hacker ◽  
Wolfgang Lieff

<p>Last year we described the campaign and the first results (Kelly et al., 2019; Neininger et al., 2019).</p><p>This year we will give an update on methods applied for estimating the regional methane emissions on a scale of about 10'000 km<sup>2</sup>, and sub-regions of about 2'500 km<sup>2</sup>.</p><p>Two approaches were applied:</p><ol><li>The classical mass balance, where the inflow and the outflow of an imaginary box was calculated, based on almost perfect Lagrangian cross-sections (following the air mass).</li> <li>A mass balance for the part of the boundary layer, where flight tracks were available (below 300 m above ground), supplemented by vertical turbulent fluxes to above this height.</li> </ol><p>In the best case, the two methods are leading to similar emission rates. The advantage of method (2) is, that the long flight legs can be limited to the lower boundary layer, which is especially useful when a convective boundary layer is reaching up to typically 2 km or higher above the surface.</p><p>The method worked quite well for water vapour, CO<sub>2</sub> and sensible heat, where fully resolved turbulent fluxes could be calculated based on 10 Hz measurements along the flight legs. Since CH<sub>4</sub> could only be measured with a temporal resolution of about two seconds (0.5 Hz), these a-priori results of the turbulent vertical fluxes are less consistent. However, by applying factors of turbulent versus advective fluxes from the other species, the agreement between the two methods was improved. The turbulent transport to above the 300-metre-layer during the convective conditions was about equal to the accumulation in this layer.</p><p>Since estimating the height of the convective boundary layer and the assumption that the mixing is perfect for approach (1) has many limitations, using method (2) has the advantage that less assumptions on homogeneity of the atmosphere above the densely observed layer has to be made. Even when the concentration profiles and the wind are known from vertical soundings (excursions to above the convective boundary layer), the horizontal inhomogeneity remains unknown. When using the vertical turbulent fluxes into this unknown volume above the lower layer, inhomogeneous mixing is not a problem.</p><p>The challenge of method (2) is to measure fast and precise enough for the quantification of the vertical fluxes. When concentrating on this, one could save time by omitting high soundings, improving the horizontal coverage, and therefore the statistics for the vertical fluxes.</p><p><strong>References</strong></p><p>Kelly et al.: Direct Measurement of Coal Seam Gas and Agricultural Methane Emissions in the Surat Basin, Australia. EGU 2019.</p><p>Neininger, B., J. M. Hacker and W. Lieff: Airborne Measurements for estimating Methane Emissions in the Surat Basin, Australia. EGU 2019.</p>


2010 ◽  
Vol 10 (1) ◽  
pp. 169-199 ◽  
Author(s):  
C. N. Hewitt ◽  
J. D. Lee ◽  
A. R. MacKenzie ◽  
M. P. Barkley ◽  
N. Carslaw ◽  
...  

Abstract. In April–July 2008, intensive measurements were made of atmospheric composition and chemistry in Sabah, Malaysia, as part of the "Oxidant and particle photochemical processes above a South-East Asian tropical rainforest" (OP3) project. Fluxes and concentrations of trace gases and particles were made from and above the rainforest canopy at the Bukit Atur Global Atmosphere Watch station and at the nearby Sabahmas oil palm plantation, using both ground-based and airborne measurements. Here, the measurement and modelling strategies used, the characteristics of the sites and an overview of data obtained are described. Composition measurements show that the rainforest site was not significantly impacted by anthropogenic pollution, and this is confirmed by satellite retrievals of NO2 and HCHO. The dominant modulators of atmospheric chemistry at the rainforest site were therefore emissions of BVOCs and soil emissions of reactive nitrogen oxides. At the observed BVOC:NOx volume mixing ratio (~100 pptv/pptv), current chemical models suggest that daytime maximum OH concentrations should be ca. 105 radicals cm−3, but observed OH concentrations were an order of magnitude greater than this. We confirm, therefore, previous measurements that suggest that an unexplained source of OH must exist above tropical rainforest and we continue to interrogate the data to find explanations for this.


2005 ◽  
Vol 18 (22) ◽  
pp. 4715-4730 ◽  
Author(s):  
P. Räisänen ◽  
H. W. Barker ◽  
J. N. S. Cole

Abstract The Monte Carlo Independent Column Approximation (McICA) method for computing domain-average radiative fluxes is unbiased with respect to the full ICA, but its flux estimates contain conditional random noise. Results for five experiments are used to assess the impact of McICA-related noise on simulations of global climate made by the NCAR Community Atmosphere Model (CAM). The experiment with the least noise (an order of magnitude below that of basic McICA) is taken as the reference. Two additional experiments help demonstrate how the impact of noise depends on the time interval between calls to the radiation code. Each experiment is an ensemble of seven 15-month simulations. Experiments with very high noise levels feature significant reductions to cloudiness in the lowermost model layer over tropical oceans as well as changes in highly related quantities. This bias appears immediately, stabilizes after a couple of model days, and appears to stem from nonlinear interactions between clouds and radiative heating. Outside the Tropics, insignificant differences prevail. When McICA sampling is confined to cloudy subcolumns and when, on average, 50% more samples, relative to basic McICA, are drawn for selected spectral intervals, McICA noise is much reduced and the results of the simulation are almost statistically indistinguishable from the reference. This is true both for mean fields and for the nature of fluctuations on scales ranging from 1 day to at least 30 days. While calling the radiation code once every 3 h instead of every hour allows the CAM additional time to incorporate McICA-related noise, the impact of noise is enhanced only slightly. In contrast, changing the radiative time step by itself produces effects that generally exceed the impact of McICA’s noise.


2013 ◽  
Vol 13 (2) ◽  
pp. 4331-4389 ◽  
Author(s):  
C. Jouan ◽  
J. Pelon ◽  
E. Girard ◽  
G. Ancellet ◽  
J. P. Blanchet ◽  
...  

Abstract. Recently, two Types of Ice Clouds (TICs) properties have been characterized using ISDAC airborne measurements (Alaska, April 2008). TIC-2B were characterized by fewer (<10 L−1) and larger (>110 μm) ice crystals, a larger ice supersaturation (>15%) and a fewer ice nuclei (IN) concentration (<2 order of magnitude) when compared to TIC-1/2A. It has been hypothesized that emissions of SO2 may reduce the ice nucleating properties of IN through acidification, resulting to a smaller concentration of larger ice crystals and leading to precipitation (e.g. cloud regime TIC-2B) because of the reduced competition for the same available moisture. Here, the origin of air masses forming the ISDAC TIC-1/2A (1 April 2008) and TIC-2B (15 April 2008) is investigated using trajectory tools and satellite data. Results show that the synoptic conditions favor air masses transport from the three potentials SO2 emission areas to Alaska: eastern China and Siberia where anthropogenic and biomass burning emission respectively are produced and the volcanic region from the Kamchatka/Aleutians. Weather conditions allow the accumulation of pollutants from eastern China/Siberia over Alaska, most probably with the contribution of acid volcanic aerosol during the TIC-2B period. OMI observations reveal that SO2 concentrations in air masses forming the TIC-2B were larger than in air masses forming the TIC-1/2A. Airborne measurements show high acidity near the TIC-2B flight where humidity was low. These results strongly support the hypothesis that acidic coating on IN are at the origin of the formation of TIC-2B.


2021 ◽  
Vol 18 (2) ◽  
pp. 557-572
Author(s):  
Sudhanshu Pandey ◽  
Sander Houweling ◽  
Alba Lorente ◽  
Tobias Borsdorff ◽  
Maria Tsivlidou ◽  
...  

Abstract. The TROPOspheric Monitoring Instrument (TROPOMI) provides observations of atmospheric methane (CH4) at an unprecedented combination of high spatial resolution and daily global coverage. Hu et al. (2018) reported unexpectedly large methane enhancements over South Sudan in these observations. Here we assess methane emissions from the wetlands of South Sudan using 2 years (December 2017–November 2019) of TROPOMI total column methane observations. We estimate annual wetland emissions of 7.4 ± 3.2 Tg yr−1, which agrees with the multiyear GOSAT inversions of Lunt et al. (2019) but is an order of magnitude larger than estimates from wetland process models. This disagreement may be explained by the underestimation (by up to 4 times) of inundation extent by the hydrological schemes used in those models. We investigate the seasonal cycle of the emissions and find the lowest emissions during the June–August season when the process models show the largest emissions. Using satellite-altimetry-based river water height measurements, we infer that this seasonal mismatch is likely due to a seasonal mismatch in inundation extent. In models, inundation extent is controlled by regional precipitation scaled to static wetland extent maps, whereas the actual inundation extent is driven by water inflow from rivers like the White Nile and the Sobat. We find the lowest emissions in the highest precipitation and lowest temperature season (June–August, JJA) when models estimate large emissions. In general, our emission estimates show better agreement in terms of both seasonal cycle and annual mean with model estimates that use a stronger temperature dependence. This suggests that temperature might be a stronger control for the South Sudan wetlands emissions than currently assumed by models. Our findings demonstrate the use of satellite instruments for quantifying emissions from inaccessible and uncertain tropical wetlands, providing clues for the improvement of process models and thereby improving our understanding of the currently uncertain contribution of wetlands to the global methane budget.


2020 ◽  
Vol 14 (1) ◽  
pp. 349-366
Author(s):  
Matthias O. Willen ◽  
Martin Horwath ◽  
Ludwig Schröder ◽  
Andreas Groh ◽  
Stefan R. M. Ligtenberg ◽  
...  

Abstract. Glacial isostatic adjustment (GIA) is a major source of uncertainty for ice and ocean mass balance estimates derived from satellite gravimetry. In Antarctica the gravimetric effect of cryospheric mass change and GIA are of the same order of magnitude. Inverse estimates from geodetic observations hold some promise for mass signal separation. Here, we investigate the combination of satellite gravimetry and altimetry and demonstrate that the choice of input data sets and processing methods will influence the resultant GIA inverse estimate. This includes the combination that spans the full GRACE record (April 2002–August 2016). Additionally, we show the variations that arise from combining the actual time series of the differing data sets. Using the inferred trends, we assess the spread of GIA solutions owing to (1) the choice of different degree-1 and C20 products, (2) viable candidate surface-elevation-change products derived from different altimetry missions corresponding to different time intervals, and (3) the uncertainties associated with firn process models. Decomposing the total-mass signal into the ice mass and the GIA components is strongly dependent on properly correcting for an apparent bias in regions of small signal. Here our ab initio solutions force the mean GIA and GRACE trend over the low precipitation zone of East Antarctica to be zero. Without applying this bias correction, the overall spread of total-mass change and GIA-related mass change using differing degree-1 and C20 products is 68 and 72 Gt a−1, respectively, for the same time period (March 2003–October 2009). The bias correction method collapses this spread to 6 and 5 Gt a−1, respectively. We characterize the firn process model uncertainty empirically by analysing differences between two alternative surface mass balance products. The differences propagate to a 10 Gt a−1 spread in debiased GIA-related mass change estimates. The choice of the altimetry product poses the largest uncertainty on debiased mass change estimates. The spread of debiased GIA-related mass change amounts to 15 Gt a−1 for the period from March 2003 to October 2009. We found a spread of 49 Gt a−1 comparing results for the periods April 2002–August 2016 and July 2010–August 2016. Our findings point out limitations associated with data quality, data processing, and correction for apparent biases.


2021 ◽  
Author(s):  
Christian Scharun ◽  
Roland Ruhnke ◽  
Michael Weimer ◽  
Peter Braesicke

&lt;p&gt;Methane (CH&lt;sub&gt;4&lt;/sub&gt;) is the second most important greenhouse gas after CO&lt;sub&gt;2&lt;/sub&gt; affecting global warming. Various sources (e.g. fossil fuel production, agriculture and waste, biomass burning and natural wetlands) and sinks (the reaction with the OH-radical as the main sink contributes to tropospheric ozone production) determine the methane budget. Due to its long lifetime in the atmosphere methane can be transported over long distances.&lt;/p&gt;&lt;p&gt;Disused and active offshore platforms can emit methane, the amount being difficult to quantify. In addition, explorations of the sea floor in the North Sea showed a release of methane near the boreholes of both, oil and gas producing platforms. The basis of this study is the established emission data base EDGAR (Emission Database for Global Atmospheric Research), an inventory that includes methane emission fluxes in the North Sea region. While methane emission fluxes in the EDGAR inventory and platform locations are matching for most of the oil platforms almost all of the gas platform sources are missing in the database. We develop a method for estimating the missing emission sources based on the EDGAR inventory and the known locations of gas platforms as additional point sources will be inserted in the model.&lt;/p&gt;&lt;p&gt;In this study the global model ICON-ART (ICOsahedral Nonhydrostatic model - Aerosols and Reactive Trace gases) is used. ART is an online-coupled model extension for ICON that includes chemical gases and aerosols. One aim of the model is the simulation of interactions between the trace substances and the state of the atmosphere by coupling the spatiotemporal evolution of tracers with atmospheric processes. ICON-ART sensitivity simulations are performed with inserted and adjusted sources to access their influence on the methane and OH-radical distribution on regional (North Sea) and global scales.&lt;/p&gt;


2015 ◽  
Vol 15 (16) ◽  
pp. 9237-9251 ◽  
Author(s):  
M. Sörgel ◽  
I. Trebs ◽  
D. Wu ◽  
A. Held

Abstract. Vertical mixing ratio profiles of nitrous acid (HONO) were measured in a clearing and on the forest floor in a rural forest environment. For the forest floor, HONO was found to predominantly deposit, whereas for the clearing, net deposition dominated only during nighttime and net emissions were observed during daytime. For selected days, net fluxes of HONO were calculated from the measured profiles using the aerodynamic gradient method. The emission fluxes were in the range of 0.02 to 0.07 nmol m−2 s−1 and thus were in the lower range of previous observations. These fluxes were compared to the strengths of postulated HONO sources. Laboratory measurements of different soil samples from both sites revealed an upper limit for soil biogenic HONO emission fluxes of 0.025 nmol m−2 s−1. HONO formation by light-induced NO2 conversion was calculated to be below 0.03 nmol m−2 s−1 for the investigated days, which is comparable to the potential soil fluxes. Due to light saturation at low irradiance, this reaction pathway was largely found to be independent of light intensity, i.e. it was only dependent on ambient NO2. We used three different approaches based on measured leaf nitrate loadings for calculating HONO formation from HNO3 photolysis. While the first two approaches based on empirical HONO formation rates yielded values in the same order of magnitude as the estimated fluxes, the third approach based on available kinetic data of the postulated pathway failed to produce noticeable amounts of HONO. Estimates based on reported cross sections of adsorbed HNO3 indicate that the lifetime of adsorbed HNO3 was only about 15 min, which would imply a substantial renoxification. Although the photolysis of HNO3 was significantly enhanced at the surface, the subsequent light-induced conversion of the photolysis product NO2 did not produce considerable amounts of HONO. Consequently, this reaction might occur via an alternative mechanism. By explicitly calculating HONO formation based on available kinetic data and simple parameterizations, we showed that (a) for low NOx the light-induced conversion of NO2 on humic acids is already light saturated by the early morning, (b) HONO formation from photolysis of adsorbed HNO3 appears to proceed via an alternative mechanism and (c) estimates of HONO emissions from soil are very sensitive to mass transfer and acidic soils do not necessarily favour HONO emissions.


2020 ◽  
Author(s):  
James L. France ◽  
Anna Jones ◽  
Tom Lachlan-Cope ◽  
Alex Weiss ◽  
Marcos Andrade ◽  
...  

&lt;p&gt;Tropical wetlands have been proposed as a potential driver for the recent rise in global atmospheric methane. However, direct access and quantification of emissions is difficult. In March 2019, a pilot study was given permission to overfly the Bolivian Llanos de Moxos wetlands to measure atmospheric mixing ratios of methane and collect spot samples for isotopic analysis. Combined with this was a short ground campaign to collect isotopic samples directly above the wetland edge to compare with the integrated atmospheric signature.&lt;/p&gt;&lt;p&gt;Atmospheric mixing ratios of methane reached a maximum of 2400 ppb (500 ppb above baseline concentrations) in the well mixed boundary layer flying at 400m above the wetland. Upwind and downwind transects were a maximum of 300 km, and methane mixing ratios increased roughly linearly with distance downwind. The isotopic data from the airborne surveys and ground surveys give a bulk isotopic signature for &amp;#948;&lt;sup&gt;13&lt;/sup&gt;C&lt;sub&gt;CH4&lt;/sub&gt; of ~-59 &amp;#8240; &amp;#177; 4, which is less negative than Amazon floodplain work focusing on emission of methane through trees, but match well with bulk isotopic values from the Amazon Basin. Ground based wetland samples taken concurrently near Trinidad, Bolivia, gave a source signature of -56 &amp;#8240; &amp;#177; 4 re-enforcing the likelihood that the atmospheric enhancements measured are related to the wetland emissions. For comparison, tropical wetlands measured at ground level during a recent Ugandan and Zambian campaign gave heavier &amp;#948;&lt;sup&gt;13&lt;/sup&gt;C&lt;sub&gt;CH4&lt;/sub&gt; isotopic source signatures of -50 to -54 &amp;#8240;. Based on this snap shot study, flux estimations suggest that the Bolivian wetlands could be emitting ~10mg CH&lt;sub&gt;4&lt;/sub&gt; m&lt;sup&gt;-2 &lt;/sup&gt;h&lt;sup&gt;-1&lt;/sup&gt;. The observed mole fractions will be compared to model simulations to determine how well the Bolivian wetland methane fluxes are represented.&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document