Articles | Volume 21, issue 11
https://doi.org/10.5194/acp-21-8791-2021
https://doi.org/10.5194/acp-21-8791-2021
Research article
 | 
10 Jun 2021
Research article |  | 10 Jun 2021

Estimating Upper Silesian coal mine methane emissions from airborne in situ observations and dispersion modeling

Julian Kostinek, Anke Roiger, Maximilian Eckl, Alina Fiehn, Andreas Luther, Norman Wildmann, Theresa Klausner, Andreas Fix, Christoph Knote, Andreas Stohl, and André Butz
Abstract

Abundant mining and industrial activities located in the Upper Silesian Coal Basin (USCB) lead to large emissions of the potent greenhouse gas (GHG) methane (CH4). The strong localization of CH4 emitters (mostly confined to known coal mine ventilation shafts) and the large emissions of 448 and 720 kt CH4 yr−1 reported in the European Pollutant Release and Transfer Register (E-PRTR 2017) and the Emissions Database for Global Atmospheric Research (EDGAR v4.3.2), respectively, make the USCB a prime research target for validating and improving CH4 flux estimation techniques. High-precision observations of this GHG were made downwind of local (e.g., single facilities) to regional-scale (e.g., agglomerations) sources in the context of the CoMet 1.0 campaign in early summer 2018. A quantum cascade–interband cascade laser (QCL–ICL)-based spectrometer adapted for airborne research was deployed aboard the German Aerospace Center (DLR) Cessna 208B to sample the planetary boundary layer (PBL) in situ. Regional CH4 emission estimates for the USCB are derived using a model approach including assimilated wind soundings from three ground-based Doppler lidars. Although retrieving estimates for individual emitters is difficult using only single flights due to sparse data availability, the combination of two flights allows for exploiting different meteorological conditions (analogous to a sparse tomography algorithm) to establish confidence on facility-level estimates. Emission rates from individual sources not only are needed for unambiguous comparisons between bottom-up and top-down inventories but also become indispensable if (independently verifiable) sanctions are to be imposed on individual companies emitting GHGs. An uncertainty analysis is presented for both the regional-scale and facility-level emission estimates. We find instantaneous coal mine emission estimates of 451/423 ± 77/79 kt CH4 yr−1 for the morning/afternoon flight of 6 June 2018. The derived fuel-exploitation emission rates coincide (±6 %) with annual-average inventorial data from E-PRTR 2017 although they are distinctly lower (−28%/−32%) than values reported in EDGAR v4.3.2. Discrepancies in available emission inventories could potentially be narrowed down with sufficient observations using the method described herein to bridge the gap between instantaneous emission estimates and yearly averaged inventories.

Dates
1 Introduction

The growth in population and economy since the pre-industrial era has been going hand in hand with rising anthropogenic emissions, causing a strong increase in atmospheric greenhouse gas (GHG) concentrations. This is in general attributed to anthropogenic emissions from a large variety of sources omnipresent in modern life (Pachauri et al.2014), e.g., extraction, processing, and transport of fossil fuels. Despite the evident anthropogenic influence on the climate in general, large uncertainties remain in the magnitude of human-induced radiative forcing and relative contributions from different sectors (Nisbet et al.2014; Kirschke et al.2013). Approximately 20 % of the global CH4 source is estimated to arise from the sector of fossil fuel industry (Schwietzke et al.2016), which also includes activities like coal mining – an industry for which the Upper Silesian Coal Basin (USCB), located in southern Poland and the Czech Republic, is well known.

According to the European Pollutant Release and Transfer Register (E-PRTR 2017; https://prtr.eea.europa.eu/, last access: 6 May 2020) a total of 448 kt CH4 yr−1 is emitted into the air from the USCB region, making it one of Europe's methane emission hot spots. The intensive mining activities and the heavy industry spread around the city of Katowice lead to these significant amounts of CH4 emitted into the atmosphere, where over 99 % of the CH4 emissions reported in E-PRTR 2017, listing emitters above a threshold of 0.1 kt CH4 yr−1, are attributed to mining and related industry. These large emissions are also apparent in the Emissions Database for Global Atmospheric Research (EDGAR v4.3.2) emission inventory reporting a total of 720 kt CH4 yr−1 in 2012.

The design and subsequent control of mitigation measures to slow down the increase in atmospheric GHG concentrations require reliable verification and attribution of GHG emissions now and in the future. An established method to derive GHG emissions is known as the top-down approach. This method is based on observed GHG concentrations in the atmosphere and projects their variations (both in time and space) back onto the emissions that may have caused these variations (Nisbet and Weiss2010; Chevallier et al.2005; Peters et al.2007). For CH4, concentrations downwind of local- to regional-scale emission sources can be sampled efficiently using high-precision airborne measurements within the planetary boundary layer (PBL). Onboard meteorological instrumentation allows for concurrent sensing of important atmospheric state variables like static pressure and air temperature, as well as the local wind field, which are particularly useful to estimate emissions. In situ instruments provide point measurements at high precision that can be used well for flux estimation using techniques like the mass balance approach (Karion et al.2013; Conley et al.2017; Pitt et al.2019). Previous studies used airborne in situ instrumentation to estimate regional CH4 emissions from oil and natural gas operations in the USA and Canada (Johnson et al.2017; Karion et al.2015; Barkley et al.2017). These studies find emission inventories (EDGAR) to underestimate CH4 emissions from the respective sector. Recent studies have also targeted urban CH4 emissions (Ryoo et al.2019; Plant et al.2019; Ren et al.2018) and anthropogenic CH4 emissions from agriculture and waste treatment (Yu et al.2020). Airborne in situ data have further been used to estimate emissions on a facility level by flying closed circles around individual emitters (Lavoie et al.2015; Conley et al.2017; Hajny et al.2019; Mehrotra et al.2017; Baray et al.2018).

Recently, Luther et al. (2019) reported on XCH4 flux estimates ranging from 6 ± 1 kt CH4 yr−1 for single shafts to up to 109 ± 33 kt CH4 yr−1 for a subregion of the USCB from ground-based, portable, sun-viewing Fourier transform spectrometers mounted on a truck. Fiehn et al. (2020) investigated CH4 emissions from the USCB using a mass balance approach. They report emission estimates of 436 ± 115 and 477 ± 101 kt CH4 yr−1 from two research flights along with a detailed uncertainty analysis. The present study is based on the very same research flights and aims at contributing an advanced model approach. Previous studies have used Lagrangian models to simulate the dispersion (Tuccella et al.2017; Raut et al.2017) of plumes emanating from oil and gas platforms or identification of CH4 sources (Platt et al.2018). Atmospheric transport models have been used to infer CH4 emissions from the oil and natural gas industry (Barkley et al.2017). Here, a combination of an Eulerian atmospheric transport model and a Lagrangian particle dispersion model is used in conjunction with assimilated Doppler lidar soundings to infer instantaneous CH4 emissions for Europe's largest coal extraction region, the USCB.

Section 2 provides an introduction to the USCB as the region of prime interest followed by a research flight overview in Sect. 3. Section 4 details a model-based flux estimation approach. CH4 emission estimates will be given in the form of a case study in Sect. 4.1 for two research flights on 6 June 2018 along with an estimate of the uncertainties involved. Section 5 summarizes our findings and concludes the study.

2 The Upper Silesian Coal Basin

The Upper Silesian Coal Basin is a plateau elevated between 200 and 300 m above sea level (a.s.l.) in southern Poland. To the south it is confined by the Tatra Mountains reaching up to 2655 m a.s.l. and forming a natural border between Slovakia and Poland. To the west it extends across the national border between Poland and the Czech Republic into the Ostrava region.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f01

Figure 1Flight trajectories of two flights of the German Aerospace Center (DLR) Cessna 208B sampling the USCB area during the CoMet field campaign. The plot shows flight trajectories for the morning flight (black) and afternoon flight (red) on 6 June 2018. Red triangles mark the location of three Doppler wind lidars, deployed in the USCB area during the CoMet campaign. Gray triangles mark known coal mine ventilation shafts. Colored tiles are from the EDGAR v4.3.2 CH4 emission inventory for 2012, showing typical emissions ranging up to  100 kt yr−1.

According to Gzyl et al. (2017), the USCB is well known for its abundant mining and industrial activities, including coal, zinc, and lead ore exploitation. Coal mining activities make up the largest part, with an approximate total of 10×109 t extracted since the industrial revolution, where over 70 % of this exploitation took place after 1945. To date an approximate 75×106 t of coal is extracted every year from 27 active mines. It is these figures and the large area of approximately 7400 km2 covered that make the USCB the largest coal extraction region in Europe (Dulias2016). The intensive coal mining activities and the heavy industry spread around the city of Katowice, Poland, located in the north of the USCB, lead to significant amounts of GHG emissions in this area. Fugitive CH4 emanating from the coal mine shafts reaching several hundred meters into the ground is either actively ventilated (active mines) or degasses passively from abandoned mines. Mines located in the north of the USCB are mostly abandoned and partially flooded, while intensive, active coal exploitation is located in the southern USCB, in both Poland and the Czech Republic (Gzyl et al.2017).

Global emission inventories show large sources of methane in this area as depicted by the colored tiles in Fig. 1. The figure is based on a subset of the publicly available EDGAR v4.3.2 CH4 (https://edgar.jrc.ec.europa.eu/, last access” 6 May 2020) emission inventory (Janssens-Maenhout et al.2017). It shows CH4 emissions range up to approx. 100 kt yr−1 on a 0.1× 0.1 grid with source strengths increasing towards the southern USCB. Accordingly, the strongest sources are located near the Czech border midway between the cities of Bielsko-Biała, Poland, and Ostrava, Czech Republic.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f02

Figure 2Total non-fuel-exploitation CH4 emission from the EDGAR v4.3.2 emission inventory corresponding to the following five sectors: solid waste landfills, energy for buildings, waste water handling, enteric fermentation, and oil refineries and transformation energy.

According to EDGAR v4.3.2, these CH4 sources are among the strongest in Europe. The total CH4 emissions from this inventory amount to approximately 720 kt yr−1 for the USCB region, where  620 kt yr−1 is attributed to the fuel-exploitation sector. The EDGAR v4.3.2 inventory further includes information on sectorial partitioning of the remaining non-fuel-exploitation CH4 emissions making up for approximately 14 % of total annual CH4 emissions in the USCB. From these  14 % approximately 90 % are attributed to the following five sectors: solid waste landfills, energy for buildings, waste water handling, enteric fermentation, and oil refineries and transformation energy. The spatial distribution of the total non-fuel-exploitation CH4 emissions from these five sectors is shown in Fig. 2. Most emitters are weak in comparison to fuel exploitation and uniformly distributed and will be canceled out by the later-described background subtraction. We added stronger source tiles (threshold  4 kt yr−1) emitting a total of 33 kt yr−1 to our FLEXPART-WRF (FLEXible PARTicle dispersion model–Weather Research and Forecasting) simulation. Facility-level emission data of CH4 are provided by E-PRTR 2017. The locations of 74 documented coal mine ventilation shafts (active and inactive) have been added to Fig. 1 for reference. These locations were visually identified from satellite imagery, and emission values from E-PRTR 2017 were evenly distributed among the ventilation shafts for each company (see also Nickl et al.2020). According to E-PRTR 2017, individual contributions sum to a total CH4 emission of approximately 448 kt yr−1. This value is approximately 38 % lower compared to the EDGAR v4.3.2 inventory (28 % if only considering the fuel-exploitation sector), showing the large uncertainties present in the available data.

3 Research flight overview

The CoMet mission in early summer 2018 primarily aimed at providing observations of GHG (mainly CO2 and CH4) gradients along large-scale latitudinal transects over Europe from the coordinated operation of several state-of-the-art instruments on the ground and aboard five research aircraft. Aboard the Cessna 208B, a rich dataset of simultaneous airborne observations of CH4, C2H6, CO2, CO, N2O, and H2O was collected using the quantum cascade laser spectrometer (QCLS) instrument (see Fig. 3 and Kostinek et al.2019, for details) during 30 flight hours. In the following, a subset of these data from two research flights undertaken on 6 June 2018 was used to retrieve CH4 fluxes emanating from the USCB region. Both flight tracks are shown in Fig. 1 along with the locations of three co-deployed Leosphere Windcube 200S Doppler wind lidars (Wildmann et al.2020). The morning (black line) and afternoon (red line) flights circumvent all known ventilation shafts in the area (gray triangles) and are in fact very similar (congruent) from the top-down perspective. This is intended to enhance confidence in retrieved GHG fluxes. Moderate (3–6 m s−1) winds throughout the day from northeasterly directions drive advection of the CH4 plumes towards the Czech border and into the Ostrava region.

The morning flight (black line in Fig. 1) on 6 June 2018 starts off from Katowice Airport, located to the north of the city center, at around 09:15 UTC. Following a short constant-altitude transect, a spiral-up was flown in the east to obtain a sounding out of the boundary layer. This maneuver, revealing a boundary layer depth of approximately 1150 m above mean sea level (a.m.s.l.), was followed by an upwind leg flown at a constant altitude of 900 ma.m.s.l. showing a fairly homogeneous CH4 inflow into the area of interest, thus allowing for subtracting an out-of-plume background (as described in Sect. 4.4) from the measured mole fractions downwind of the mines. Mixing ratios decreased slightly towards free tropospheric background values when climbing above the PBL. Before returning back to Katowice Airport at around 11:45 UTC, the downwind wall maneuver, consisting of five constant-altitude flight legs to the west, was performed at altitudes of approximately 800 m, 1.1 km, 950 m, 1.4 km, and 1.8 kma.m.s.l.. During the last two flight legs, the aircraft was outside of the PBL.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f03

Figure 3The DLR Cessna 208B on the taxiway at Katowice Airport. The sample air intake is mounted underneath the right wing. The meteorological data acquisition system is underneath the left wing. The QCLS instrument (lower right panel) is located inside the cabin behind the pilot seats.

Download

The afternoon flight (red line in Fig. 1) started off from Katowice Airport around 13:15 UTC. An upwind leg flown at a constant altitude of 900 ma.m.s.l. again showed a fairly homogeneous CH4 inflow into the USCB area. Mixing ratios decrease slightly towards free tropospheric background values when climbing above the PBL during a climb and descent maneuver flown parallel to the sensed mean wind direction. During this flight, we observed a latitudinally inclined PBL with an approximate depth of 1.7 kma.m.s.l. in the northern section and 1.3 kma.m.s.l. towards the south. Before returning back to Katowice Airport at around 15:30 UTC, the downwind wall maneuver, consisting of six constant-altitude flight legs, was performed over the western USCB region at altitudes of approximately 800 m, 890 m, 975 m, 950 m, 1.1 km, 1.5 km, and 1.8 kma.m.s.l..

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f04

Figure 4The FLEXPART-WRF domain resides in the nested WRF domain D2 providing the meteorological driver data. Generous spacing towards the driver domain has been included to avoid spurious boundary effects.

4 Estimating emissions

The model-based approach developed in this work employs a combination of Eulerian and Lagrangian particle dispersion models. Due to the known locations of the coal mine ventilation shafts, their emissions are modeled forward in time with constant emission rates. Modeled data are then extracted at the aircraft position in space and time and compared to actual airborne in situ observations. This comparison depends on the quality of the a priori emission data, the quality of the measurements, and the quality of the transport model simulation, which in turn depends on the quality of the meteorological data (winds, PBL heights, etc.). Validating the meteorological data is therefore important to enable regional emission estimates based on particle dispersion models. Here, meteorological driver data are generated using the Weather Research and Forecasting (WRF) v4.0 model (Powers et al.2017) with assimilated soundings from three Leosphere Windcube 200S Doppler wind lidars. Data are then fed into the Lagrangian particle dispersion model FLEXPART-WRF – a FLEXPART (Pisso et al.2019) version adapted for WRF meteorology – and used to model the exhaust plumes emanating from ventilation shafts of the emitters listed in E-PRTR 2017.

4.1 Local-scale meteorology using WRF

Figure 4 shows a satellite map of central Europe with the two domains specified for the USCB region. The outer domain D1 (light blue box in Fig. 4) with a horizontal grid resolution of ∼15km includes large parts of central Europe.

This domain is fed by NCEP GDAS/FNL Operational Global Analysis data on a 0.25×0.25 grid, available from the NCAR–UCAR Research Data Archive at a 3 h time resolution (National Centers for Environmental Prediction/National Weather Service/NOAA/U.S. Department of Commerce2015). The grid four-dimensional data assimilation (GFDDA) module is used to nudge modeled meteorology towards the analysis data at each grid point. The outer domain is intended to catch the large-scale weather situation over Europe and to provide a smooth transition between the coarse NCEP GDAS/FNL Operational Global Analysis and the region of prime interest. The inner domain D2 (yellow box in Fig. 4) has a horizontal grid resolution of ∼3km and covers the entire USCB region. The model output from D2 is the primary product required for subsequent FLEXPART runs. Both domains are driven with the original WRF v4.0 topographic data with a resolution of 30 arcsec. Vertically, the model atmosphere is divided into 33 stacked layers, with the top layer at 200 hPa (corresponding to approximately 12 km in altitude). Vertical layers are more closely spaced at lower altitudes to enable a better resolution of boundary layer processes. The modeled atmospheric state variables are output every hour for D1 and every 5 min for D2.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f05

Figure 5Ensemble WRF runs with varying radii of influence rxy in comparison to interpolated NCEP GDAS/FNL and actual lidar soundings for 6 June 2018 at 09:00 UTC. The shaded area shows the maximum variability including soundings timed 20 min before and 20 min after the observations used. Abscissa units are m s−1 and degrees.

Download

Soundings from three Doppler lidars (marked DLR85, DLR86, and DLR89) deployed in the USCB area during the CoMet mission (see Fig. 1 for respective positions) have been used to augment the model output (Wildmann et al.2020). These data are available on a regular, continuous basis throughout the campaign period at 10 min time intervals with soundings typically reaching up to  2.5 kma.m.s.l. depending on the atmospheric conditions. Domains D1 and D2 are both nudged towards the Doppler soundings using the WRF-FDDA subsystem (Deng et al.2008). Sensitivity of the model output to three key parameters of the observational data assimilation subsystem, namely the radius of influence rxy, time window Δt, and horizontal wind coefficient cuv, was analyzed through numerous runs with the goal of finding an appropriate configuration. Figure 5 shows ensemble runs with varying rxy in comparison to interpolated NCEP GDAS/FNL and the actual lidar soundings for 6 June 2018 at 09:00 UTC. The gray-shaded area beneath the orange-colored lidar soundings shows the maximum variability including soundings timed 20 min before and 20 min after the observations used. Figure 5 demonstrates that modeled data are in good agreement with observed Doppler soundings when using WRF-FDDA. It also shows discrepancies between NCEP GDAS/FNL driver data and observations in wind direction and more importantly in the wind speed in the lower troposphere and the PBL depth. To further enhance compatibility between model and observations, the WRFDA submodule (Barker et al.2012) was used in 3DVar cycling mode similarly to in Liu et al. (2013) using the NCAR CV3 background error covariance (Barker et al.2004). The required observational error covariances are taken from the measurement uncertainties.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f06

Figure 6Comparison of 1 Hz wind speed (a), wind direction (b), static pressure (c), and static air temperature (d) as measured aboard the Cessna 208B on 6 June 2018 between 10:00 and 11:20 UTC to ensemble runs with varying rxy from above. The graphs are plotted as a function of time in minutes since 06:00 UTC.

Download

To verify and validate the observational FDDA approach, non-assimilated meteorological in situ data collected aboard the Cessna 208B are compared to modeled data in Fig. 6. In particular, Fig. 6 compares the 1 Hz wind speed, wind direction, static pressure, and static air temperature as measured on 6 June 2018 between 10:00 and 11:20 UTC with the underwing boom-mounted data acquisition system to ensemble runs with varying rxy from above. These data were collected approximately 35 km (minimum distance) to 65 km (maximum distance) to the west of the nearest wind lidar during the downwind wall phase of the morning flight (see Fig. 1). Simulated data, extracted at the aircraft positions in space and time, agree with 1 Hz observations of wind speed and direction to within an RMSE of ±0.7 m s−1 (1σ) and ±5 (1σ), respectively. Here, the NCAR Command Language (NCL; Brown et al.2012) has been used to interpolate from gridded model output to the exact aircraft position in space and time. Modeled wind speed deviates from observed winds during the last 20 min of the downwind wall. A possible reason for this might be that the flight leg was in close vicinity to the PBL top height. A bias of modeled static pressure and static air temperature is evident from the lower panels in Fig. 6. Modeled pressure has a consistent offset of −5hPa compared to in situ data, and modeled temperature is biased by approximately 2.2 K towards lower values.

4.2 Plume dispersion using FLEXPART-WRF

FLEXPART-WRF version 3.3.2 (Brioude et al.2013) was used to model the exhaust plumes of known emitters forward in time using the meteorological data (including PBLH) obtained from the WRF simulations described above (see Sect. 4.1) as a driver. The model is set to release 50 000 particles and an arbitrarily chosen total mass of me=1×105kg for each release location during the total simulated time of τe= 9 h. The model output is gridded into 100×100 horizontal tiles and 24 vertically stacked layers ranging from ground level up to 3 km in altitude. This results in a horizontal resolution of approximately 1.3 km and a vertical resolution of 50 m near the ground, gradually increasing to 500 m above 2 km in altitude. The domain has been placed inside the nested WRF domain D2 with generous spacing towards the domain boundaries as indicated in Fig. 4 to avoid spurious boundary effects. The main product of the FLEXPART-WRF runs is concentration fields for each release location in units of ng m−3, which are scaled a posteriori to deduce the emission rate of each modeled release. Each coal mine ventilation shaft is modeled as a constant, continuous volume source φi with a 10 m× 10 m horizontal footprint and extending 10m in the vertical direction. The volume emitter sizes are based on the construction of typical ventilation shafts in the USCB (Swolkień2020).

Mass densities in units of ng m−3 can be extracted for the aircraft positions from the model output. The result is an m×n linear forward model matrix Kji that links scaling factors for emission rates to atmospheric mass density enhancements at the measurement instances, where m is the number of observations available and n is the number of modeled release locations; i.e., Kji is the mass density that source i contributes to observation j. A scaling coefficient xi is assigned to each of the n sources φi=me,iτe,i-1, with the total emission time τe,i in seconds and the total mass emitted me,i in kilograms for each simulated source. These last two parameters are both assigned in the FLEXPART-WRF input file.

Following a maximum a posteriori (MAP) approach, the scaling coefficients xi can be found for each of the n modeled sources φi and for each of the m observed enhancements yj making use of a priori information xa on the emissions of the individual shafts. Following Bayes' theorem the MAP solution is given by the minimum of the cost function (Tarantola2004; Jacob2007; Rodgers2000):

(1) J x = x - x a T S a - 1 x - x a + y - K x T S ϵ - 1 y - K x ,

with later-defined a priori and observational error covariance matrices Sa and Sϵ, respectively. The MAP solution can be found by solving for xJ(x)=0 and is given by

(2) x ^ = x a + G y - K x ,

with the gain matrix

(3) G = S a K T KS a K T + S ϵ - 1 .

By exploiting the averaging kernel A=GK, the number of degrees of freedom for signal ds can be computed as

(4) d s = t r A .

This number describes the reduction in the normalized error in x introduced by the available observations and hence provides a measure for the improvement in knowledge of x, relative to the a priori, due to the observations.

The total emission estimate Φ in units of kg s−1 follows from the scaled sum of the individual contributions φi:

(5) Φ = i = 1 n x i φ i .

Here, the non-negative least squares (NNLS) algorithm (Lawson and Hanson1995) has been used to minimize the MAP cost function subject to the constraint x>0. This constraint is equivalent to the absence of negative sources. The NNLS algorithm solves the constrained least squares problem by splitting into active and passive subsets, where active and passive refer to the state of the constraint. The algorithm subsequently solves the unconstrained least squares problem for the passive set.

4.3 Estimating total uncertainty

The outlined approach is based on assumptions, of which the most important ones are a constant emission rate over the timescale of transport from the source to the aircraft; an appropriate atmospheric background vector b (used to compute CH4 enhancements y=ρ-b from measured mass densities ρ); long lifetime of the species of interest, i.e., no chemical and physical removal on the timescale of a flight; and the model being able to adequately represent the meteorological state variables. To assess uncertainty on the retrieved emission rates, several variables have been selected as most influencing systematic error sources: wind speed, wind direction, PBL height, source dislocation, and an error in sensed mole fractions that is further intended to include an error due to chosen background. Individual contributions of these error sources to total uncertainty can be identified from ensemble model runs with systematically perturbed parameters.

In addition to derived systematic uncertainties, statistical errors related to the MAP fit are to be acknowledged. The statistical uncertainty ϵi in the retrieved parameters xi can be expressed in terms of the parameter covariance matrix S^ as

(6) ϵ i = S ^ i i

for individual scaling coefficients. For regional estimates we further included off-diagonal elements of S^. The parameter covariance matrix S^ is computed from the m-by-n dimensional forward model K using

(7) S ^ = K T S ϵ - 1 K + S a - 1 - 1 .
https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f07

Figure 7(a) Observational covariance augmented with a first-order autoregressive model AR(1) structure with φ=0.7. (b) A priori covariance matrix highlighting correlated clusters from individual ventilation shafts. The scaling factor results from the scaling of the FLEXPART-WRF plume to the designated flux unit kt yr−1.

Download

Due to the different timescales of model output and observations, the observational error covariance matrix Sϵ cannot simply be taken as a diagonal matrix as this would neglect any correlations. Although the influence on regional estimates is small, these correlations have a significant effect on subregional estimates. The main diagonal of Sϵ has been estimated from the squared observation uncertainties σ2 inflated with a transport model error χm obtained from ensemble runs with perturbed parameters (see Sect. 4.6). In order to compensate for tempo-spatio-autocorrelation, we simulated a puff release to check the impulse response of the simulation on a 1 s release from individual sources. We sampled simulated observations from this puff release at the aircraft location in space and time and computed the autocorrelation function (ACF) from a single dispersed puff. The exponential decay of the correlogram suggested to augment the observational covariance matrix with a first-order autoregressive model AR(1) structure with φ=0.7. Figure 7 shows a zoom on the first 400 × 400 elements of this 14 936 × 14 936 matrix to highlight the introduced off-diagonal elements. The measurement uncertainties σi were obtained via standard error propagation from the uncertainties associated with different instruments aboard the aircraft needed for the computation of the CH4 mass density observations (see Eq. 8). Static air temperature, static air pressure, and wind speed can be probed with an uncertainty of σT= 0.15 K, σp= 1 hPa, and σu= 0.3 m s−1 (1s-1σ; Mallaun et al.2015), respectively. CH4 mole fractions were sampled with a total uncertainty better than 1.85 ppb (1s-1σ). The main diagonal of the a priori error covariance matrix Sa (i.e., the a priori variances) contains the squared a priori uncertainties, estimated with 50 % of the nominal value. As several shafts cluster around individual mines at distances of not much more than a kilometer, we further introduced a +0.5 correlation into the mines belonging to the same cluster and mining company and added local information (see review comment by Jaroslaw Necki, https://doi.org/10.5194/acp-2020-962-SC1). The final a priori error covariance matrix is depicted in Fig. 7. The scaling factor results from the scaling of the FLEXPART-WRF plume to the designated flux unit kt yr−1.

4.4 Case study – 6 June 2018

Figure 8 shows a time series of the measured and modeled CH4 mass density enhancement as a function of flight time during the downwind wall phase with the atmospheric background subtracted and source coefficients xi already optimized (according to Eq. 1).

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f08

Figure 8Time series of the in situ measured (gray) and modeled CH4 mass density (solid black) with subtracted background as a function of flight time during the downwind wall phase of the morning (a) and afternoon (b) flight on 6 June 2018 with optimized source coefficients xi. The dotted light blue line corresponds to the same forward simulation using scaling coefficients deduced from E-PRTR 2017.

Download

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f09

Figure 9(a) Interpolated cross section of the model output along the downwind wall including scattered in situ observations of CH4 for the morning flight on 6 June 2018. The color bar also applies to the isolines on the right. (b) Top-down view of the model output and the downwind wall observations ρ at 750 ma.g.l. with underlaid topography. Triangles mark simulated emitters with colors corresponding to the optimized source strengths. Gray isolines correspond to non-fuel-exploitation emissions. Both panels show a snapshot of the model output at one fixed time chosen as the center time of the downwind wall phase.

The measured scalar mass density ρ has been deduced from the ideal gas law pV =mRsT (mass m, specific gas constant Rs=R/M, and molar mass of CH4 M) using the in situ measured static air temperature and static air pressure according to

(8) ρ x = m x V air = m x m air p R T air ,

where mx denotes the total mass of the species of interest. The unitless coefficient mxmair-1=cxMxMair-1 is obtained from the sensed CH4 mole fractions cx in units of mol mol−1.

Prior to conversion, the atmospheric background has to be subtracted from the observed cx. The choice of background is to some extent ambiguous because there is no clear edge between background and in-plume sampling. This contributes to total flux estimation uncertainty, as will be discussed later. Here, a piecewise linear interpolation between the outermost boundaries of each of the four flight legs (see Fig. 9) has been considered the best guess of atmospheric background. The mean value of 20 samples has been used on both edges of each flight leg. Using this approach, latitudinal and longitudinal gradients in background CH4 mole fractions are accounted for by using both edges of each flight leg. Vertical gradients in background CH4 levels are accounted for by treating each constant-altitude flight leg separately.

From Fig. 8 a good overall match between model (solid black line) and in situ observations (gray) is apparent with a mean bias of 2.5×10-9kg m−3 and a root mean square error of 1.6×10-8kg m−3. Some of the minor structure is not reproduced in detail by the model, which is to be expected due to the model's 3 km horizontal grid resolution. The reason for the discrepancy between model and the first few hundred observations in Fig. 8 becomes more obvious when looking at the 2D scene shown in Fig. 9. The left panel of Fig. 9 shows a cross section of the model output along the downwind wall including the in situ observations of ρ. The right panel of Fig. 9 depicts the top-down view of the model output and the downwind wall observations at a fixed altitude of 750 ma.g.l. It should be noted here that both panels show a snapshot of the model output at one fixed time chosen as the center time of the downwind wall.

The discrepancy between model and observation at the lowermost (first in time) flight leg, corresponding to the southernmost trajectory section in Fig. 9 (right panel), cannot be reproduced by any of the included emission sources. A possible source is urban CH4 emissions of Krakow, located to the east of the USCB region. An area source, covering the greater city area, has therefore been included in the model. Its influence can be seen at the rightmost edge of Fig. 9 (right panel). Although we identified the city of Krakow as a possible source, we omitted it in the USCB emission estimates, as it does not officially belong to the USCB area. There are other parts in the time series where the model either underestimates (e.g., times around observation numbers 2000–3000 and 4000–4500) or overestimates (e.g., around observation number 6000) emissions. This might well be related to sources not taken into account or deficiencies in wind speed, wind direction, PBL height, etc.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f10

Figure 10(a) Interpolated cross section of the model output along the downwind wall including scattered in situ observations of CH4 for the afternoon flight on 6 June 2018. The color bar also applies to the isolines on the right. (b) Top-down view of the model output and the downwind wall observations ρ at 750 ma.g.l. with underlaid topography. Gray isolines correspond to non-fuel-exploitation emissions. Triangles mark simulated emitters with colors corresponding to the optimized source strengths. Both panels show a snapshot of the model output at one fixed time chosen as the center time of the downwind wall phase.

The instantaneous fuel-exploitation emission estimate directly follows from the optimized parameters xi via Eq. (5). The emission estimate obtained for the morning flight on 6 June 2018 using the model-based approach amounts to Φ= 451 ± 77 kt yr−1. It differs from the yearly averaged inventorial emission estimates for the USCB region by approximately −28% for EDGAR v4.3.2 and +1% for the E-PRTR inventory (excluding simulated non-fuel-exploitation sources). In addition to the coal mine emissions, approximately 27 kt yr−1 of CH4 are estimated to emanate from simulated non-coal-mine sources. The retrieval for the morning flight yields 48 degrees of freedom for signal and a total of 35 out of 74 modeled sources actively emitting. Here, the large number of degrees of freedom for signal is indicative of the validity of the total emission estimate. The latter can be confirmed by the negligible impact of the a priori 1-Axa versus the observations Ax in the emission estimate. The algorithm makes use of the large number of modeled sources to enable a total emission estimate plus additional information on individual sources.

To enhance confidence in the emission estimate, an afternoon flight of the DLR Cessna 208B was carried out a few hours after the morning flight ended on 6 June 2018. Due to consistent wind directions on that day, the flight pattern was kept as close as possible to the morning flight. The flight trajectories for both flights are depicted in Fig. 1. The right panel in Fig. 8 shows the corresponding time series of the measured and modeled CH4 mass density as a function of flight time during the downwind wall phase with source coefficients xi already optimized. Like for the morning flight a good overall match between model and in situ observations can be observed with a mean bias of 3×10-10kg m−3 and a root mean square error of 1×10-8kg m−3. The sensed mixing ratios are lower compared to the morning flight due to a further developed and hence more diluted boundary layer in the afternoon. The corresponding snapshot 2D scene is depicted in Fig. 10, with the left panel showing a cross section of the model output along the downwind wall including the in situ observations of ρ. The right panel of Fig. 10 shows the top-down view and the downwind wall observations as before. Both panels show a snapshot of the model output at the center time of the downwind wall phase. It is evident from Fig. 10 that the inclined boundary layer height observed during the flight is nicely captured by the model. Boundary layer depth is generally enhanced compared to the morning flight. Plume trajectories are streamlined implying consistent winds over time.

The fuel-exploitation emission estimate obtained for the afternoon flight on 6 June 2018 using the model-based approach amounts to Φ= 423 ± 79 kt yr−1. The obtained instantaneous emission estimate differs from the yearly averaged inventorial emission values for the USCB region by approximately −32 % for EDGAR v4.3.2 and −6 % for the E-PRTR inventory (excluding simulated non-fuel-exploitation sources). The retrieval for the afternoon flight yields 31 kt yr−1 of CH4 from simulated non-coal-mine sources, 42 degrees of freedom for the signal, and a total of 38 out of 74 simulated coal mine sources actively emitting. Both flights yield similar ds values, indicating that not all information stems from observations alone. Hence, neither flight can be used alone to retrieve all modeled sources. In an effort to minimize the dependency on the a priori, both flights will be analyzed together in the next section.

4.5 Subregional emission estimates

The model-based approach provides a unique advantage over established mass balance techniques in terms of spatial information, as it enables attributing sensed CH4 mole fractions to remote sources at distances of tens to hundreds of kilometers. The achievable level of confidence for these subregional estimates however strongly depends on the observational data. The total emission estimate has been introduced in Sect. 4.2 as the sum over n sources φi that are individually scaled with a coefficient xi. The emission rate Φi corresponding to the ith source is thus given by xiφi. By including all n sources in the state vector, individual scaling coefficients can be derived for individual sources. Here, the availability of data from two research flights on 6 June 2018 was exploited to estimate subregional emission rates Φi for individual sources. As the mean wind direction differed by  10 % between the two flights, uncertainty on the shaft level remained large and observational data are too limited for a more specific estimate. With mission planning further optimized for the Bayesian inversion from airborne in situ data, as presented in this paper, these uncertainties can potentially be narrowed down in future campaigns. The retrieval yields ds=32 with 53 sources actively emitting.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f11

Figure 11Emission estimates Φi (blue bars) in kt yr−1 for 74 individual mining shafts using the morning and afternoon flight of 6 June 2018. Slim green bars are the reported yearly average values for each mining company (E-PRTR 2017) evenly distributed among the respective ventilation shafts. The orange error bars stem from the quadrature sum of the statistical uncertainties ϵi (computed from the parameter covariance matrix S^) and the uncertainties σensemble derived from a variational ensemble with systematically perturbed parameters.

Download

Figure 11 illustrates Φi in kt yr−1 for all modeled mining shafts, taking into account both research flights on 6 June 2018. The blue bars represent the estimated Φi and are to be related to the yearly average values (slim green bars) for each mining company reporting to E-PRTR 2017 for illustrative purposes. The estimated uncertainty depicted in Fig. 11 includes systematic uncertainties derived from a variational ensemble and statistical uncertainties due to the fit algorithm used (see Sect. 4.3). The variational ensemble introduced in Sect. 4.3 includes scaling coefficients xi subject to systematic variations in key sources of uncertainty. Systematic uncertainties for each source are directly obtained from this ensemble run. Differences in estimated and reported (E-PRTR 2017) Φi values are evident. This is however to be expected due to the comparison of instantaneous emission estimates and yearly averages. Figure 12 shows the a posteriori correlation matrix as deduced from the MAP covariance matrix. The matrix indicates that there remained some uncertainty regarding to which shaft the emissions had to be assigned.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f12

Figure 12A posteriori correlation matrix as deduced from the a posteriori covariance matrix showing there remained some uncertainty regarding to which shaft the emissions had to be assigned.

Download

4.6 Uncertainty analysis

The influence of several variables on the total flux estimate Φ has been computed from eight sensitivity runs with symmetrically perturbed parameters. The systematic transport model uncertainty is subsequently estimated as the standard deviation of this ensemble. Figure 13 shows the influence of an error in wind speed (σu= 0.9 m s−1), wind direction (σd= 5 ), PBL height (σpbl= 100 m), and a source dislocation (σsd= 1 km) on total uncertainty for the flights detailed in the previous section. An assumed error in sensed mole fractions (σc= 10 ppb) is intended to include an error due to wrongly chosen background. The error in wind speed σu is taken as the standard deviation of the difference between WRF modeled wind and non-assimilated in situ observations from the data depicted in Fig. 6.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f13

Figure 13Ensemble runs to assess uncertainty in the flux estimates derived using a model-based approach. All selected error sources contribute to total uncertainty on a similar level.

Download

The same holds for the wind direction. The difference between modeled data and observations should therefore reflect overall uncertainty in these variables. Two spiral-up soundings out of the PBL revealed a boundary layer height of 1150 m at 09:37 UTC and 1300 m at 11:45 UTC. Based on these two soundings, the uncertainty in boundary layer depth is estimated with σpbl= 100 m for the downwind wall phase between 10:00 and 11:00 UTC. For this sensitivity analysis, the WRF fields were perturbed systematically during the FLEXPART read phase in readwind.f90. The source dislocation was implemented in the FLEXPART configuration file. It is evident from Fig. 13 that all selected error sources contribute on a similar level to total systematic uncertainty, which is ultimately computed as the standard deviation of the ensemble.

https://acp.copernicus.org/articles/21/8791/2021/acp-21-8791-2021-f14

Figure 14Jacobian with respect to xi and the observations yj for the morning flight on 6 June 2018. All scaling coefficients xi are sensitive to variations in yj and can therefore be deduced from a MAP fit. Measurements centered around observation 2300 are not covered by the model and are thus obsolete for flux estimation using this particular flight.

Download

In addition to the derived systematic uncertainties, statistical errors related to the MAP fit have been computed following Sect. 4.3. Figure 14 depicts the Jacobian K with respect to xi and the observations of the morning flight on 6 June 2018. It describes the change in residuals introduced by a change in parameter xi. From this figure, it can be seen that all 74 modeled coal mine sources were sampled by the aircraft using the chosen flight pattern, as all scaling coefficients are represented in the Jacobian. For the fluxes emanating from the USCB area, the statistical uncertainty ϵi computed from the Jacobian following Eqs. (6) and (7) amounts to approximately 66 kt yr−1 or 15 %, respectively.

Ultimately, the total uncertainty for the morning flight (same for the afternoon flight) on 6 June 2018 is the quadrature sum of systematic (44 kt yr−1) errors from the ensemble runs and statistical uncertainty (66 kt yr−1) from the fitting algorithm adding up to approx. 17 % relative uncertainty.

5 Conclusions

A modified Aerodyne dual QCLS instrument was deployed aboard the DLR Cessna 208B in the context of the CoMet 1.0 campaign in early summer 2018 with the goal of estimating hard-coal-mine CH4 emissions emanating from the USCB area – Europe's largest coal extraction region. Intensive mining activities and the heavy industry spread around the city of Katowice lead to significant amounts of GHGs emitted into the atmosphere. The reported inventorial CH4 emission rates for the entire USCB region amount to 720 kt yr−1 (EDGAR v4.3.2) and 448 kt yr−1 (E-PRTR 2017). The latter corresponds to 12.5 MtCO2eq.yr-1 using a value of CH4 GWP100=28 from the IPCC Fifth Assessment Report (Pachauri et al.2014). Assuming an average carbon content of 75 %, a net calorific value of 29 MJ kg−1, an emission factor of 94 t CO2 (TJ)−1, and the approximate 75×106 t of coal extracted from the USCB every year results in yearly CO2 emissions of 205 Mt CO2 yr−1 from burning of the extracted coal. The CH4 emissions from mining alone therefore make up approximately 6 % in terms of GWP.

Estimates of coal mine CH4 emissions in the USCB were derived using a model approach based on the Eulerian WRF model and the Lagrangian particle dispersion model FLEXPART-WRF. Data assimilation further exploits the availability of additional data products, e.g., wind lidar soundings during the CoMet 1.0 campaign. Due to the known locations of the coal mine ventilation shafts, sources are modeled forward in time assuming a constant emission rate. Modeled data are then extracted at the aircraft positions in space and time and compared to actual airborne in situ observations. Here, meteorological driver data were generated using the WRF v4.0 model with continuous assimilated wind lidar soundings using WRF's OBS-FDDA and WRFDA subsystems. After validation with unassimilated in situ measurements, data were fed into the Lagrangian particle dispersion model FLEXPART-WRF and used to model the exhaust plumes of the ventilation shafts. Using an inverse modeling approach, a priori emission data from E-PRTR 2017 are optimized to allow a better fit to the observations. Thereby, total emission estimates for the USCB area of Φ= 451 ± 77 kt yr−1 and Φ= 423 ± 79 kt yr−1 were obtained for a morning flight and an afternoon flight on 6 June 2018, respectively. This includes non-fuel-exploitation fluxes, estimated with 27 and 31 kt yr−1 for the morning and afternoon flights, respectively. Morning and afternoon flights differ by less than 4 % corresponding to an excellent agreement well within the uncertainty range. The obtained emission estimate differs from the inventorial emission estimates by approximately −28%/−32% for the EDGAR v4.3.2 inventory (morning flight/afternoon flight) and +1%/−6% for the E-PRTR inventory (excluding non-fuel-exploitation sources), respectively. Differences in estimated and reported emission rates are however expected due to the comparison of instantaneous estimates and yearly averages. This is in line with previous studies hinting at EDGAR v4.3.2 overestimating CH4 emissions in the USCB (Luther et al.2019; Fiehn et al.2020). Uncertainty estimates include systematic contributions from ensemble runs and statistical uncertainty introduced by the fitting algorithm. Data from both research flights are further exploited to estimate individual source contributions. Differences between individual estimates and E-PRTR reported emissions are observed. This is expected due to several reasons: a limited number of measurements relative to the yearly averages provided in the inventories, wind directions do not differ by much between the two flights, and the evenly distributed emissions among the ventilation shafts for each mining company. In general, the approach described herein delivers more information compared to the conventional mass balance, albeit at increased effort: wind lidars need to be deployed during the measurement campaign, models need to be run, wind lidar data need to be assimilated, and inverse estimation techniques need to be applied. The additional possibility of remote source attribution, however, coupled with the results obtained in Sect. 4.2 for the regional USCB anthropogenic CH4 emissions makes this approach a potent alternative to the mass balance technique. Although retrieving estimates for individual emitters is not possible using only single flights, due to sparse data availability, the combination of two or more flights allows for exploiting different meteorological conditions to enhance confidence in facility-level estimates.

Code and data availability

Data are available from HALO-DB: https://halo-db.pa.op.dlr.de/mission/94 (DLR2019). WRF v4.0 can be downloaded from https://www.mmm.ucar.edu/weather-research-and-forecasting-model (University Corporation for Atmospheric Research2020). FLEXPART-WRF can be downloaded from https://www.flexpart.eu/ (Zentralanstalt für Meteorologie und Geodynamik2019). Model setups for both models employed are available upon request.

Author contributions

AR and AB developed the research question. AR, JK, ME, AF, and AL actively took part in the CoMet field campaign by operating instrumentation and collecting and analyzing data. NW deployed wind lidars and collected and analyzed data. AS contributed substantially to the inverse method. AF coordinated the CoMet field campaign operations. TK provided emission data for individual shafts. CK provided and optimized the WRF model configuration. JK wrote the paper.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “CoMet: a mission to improve our understanding and to better quantify the carbon dioxide and methane cycles”. It is not associated with a conference.

Acknowledgements

We thank DLR VO-R for funding the young investigator research group Greenhouse Gases. We also acknowledge funding from the BMBF under project AIRSPACE (grant no. FKZ01LK170). The CoMet mission has further been supported by the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG) within the DFG Priority Programme (SPP 1294) Atmospheric and Earth System Research with the “High Altitude and Long Range Research Aircraft” (HALO). We greatly appreciate continuous support from Hans Schlager and Markus Rapp from the DLR. We thank Paul Stock, Monika Scheibe, and Michael Lichtenstern from the DLR for engineering support. Furthermore we would like to thank everyone involved during the CoMet field campaign for their relentless dedication and helpful discussions, especially Jarosław Necki and Justyna Swolkien, who provided valuable information on CH4 ventilation and took an active role as local advisers during the campaign period. This work used resources of the Deutsches Klimarechenzentrum (DKRZ) granted by its Scientific Steering Committee (WLA) under project ID 1104.

Financial support

The CoMet 1.0 campaign was funded by the MPG (Max Planck Society) and BMBF (German Federal Ministry of Education and Research) through project AIRSPACE (FK 390 01LK1701C), as well as the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG Priority Program SPP 1294).

The article processing charges for this open-access publication were covered by the German Aerospace Center (DLR).

Review statement

This paper was edited by Ilse Aben and reviewed by two anonymous referees.

References

Baray, S., Darlington, A., Gordon, M., Hayden, K. L., Leithead, A., Li, S.-M., Liu, P. S. K., Mittermeier, R. L., Moussa, S. G., O'Brien, J., Staebler, R., Wolde, M., Worthy, D., and McLaren, R.: Quantification of methane sources in the Athabasca Oil Sands Region of Alberta by aircraft mass balance, Atmos. Chem. Phys., 18, 7361–7378, https://doi.org/10.5194/acp-18-7361-2018, 2018. a

Barker, D., Huang, X.-Y., Liu, Z., Auligné, T., Zhang, X., Rugg, S., Ajjaji, R., Bourgeois, A., Bray, J., Chen, Y., Demirtas, M., Guo, Y.-R., Henderson, T., Huang, W., Lin, H.-C., Michalakes, J., Rizvi, S., and Zhang, X.: The Weather Research and Forecasting Model's Community Variational/Ensemble Data Assimilation System: WRFDA, B. Am. Meteorol. Soc., 93, 831–843, https://doi.org/10.1175/BAMS-D-11-00167.1, 2012. a

Barker, D. M., Huang, W., Guo, Y.-R., Bourgeois, A. J., and Xiao, Q. N.: A Three-Dimensional Variational Data Assimilation System for MM5: Implementation and Initial Results, Mon. Weather Rev., 132, 897–914, https://doi.org/10.1175/1520-0493(2004)132<0897:ATVDAS>2.0.CO;2, 2004. a

Barkley, Z. R., Lauvaux, T., Davis, K. J., Deng, A., Miles, N. L., Richardson, S. J., Cao, Y., Sweeney, C., Karion, A., Smith, M., Kort, E. A., Schwietzke, S., Murphy, T., Cervone, G., Martins, D., and Maasakkers, J. D.: Quantifying methane emissions from natural gas production in north-eastern Pennsylvania, Atmos. Chem. Phys., 17, 13941–13966, https://doi.org/10.5194/acp-17-13941-2017, 2017. a, b

Brioude, J., Arnold, D., Stohl, A., Cassiani, M., Morton, D., Seibert, P., Angevine, W., Evan, S., Dingwell, A., Fast, J. D., Easter, R. C., Pisso, I., Burkhart, J., and Wotawa, G.: The Lagrangian particle dispersion model FLEXPART-WRF version 3.1, Geosci. Model Dev., 6, 1889–1904, https://doi.org/10.5194/gmd-6-1889-2013, 2013. a

Brown, D., Brownrigg, R., Haley, M., and Huang, W.: The NCAR Command Language (NCL)(version 6.0.0), UCAR/NCAR Computational and Information Systems Laboratory, https://doi.org/10.5065/D6WD3XH5, 2012. a

Chevallier, F., Fisher, M., Peylin, P., Serrar, S., Bousquet, P., Bréon, F.-M., Chédin, A., and Ciais, P.: Inferring CO2 sources and sinks from satellite observations: Method and application to TOVS data, J. Geophys. Res.: Atmospheres, 110, D24309, https://doi.org/10.1029/2005JD006390, 2005. a

Conley, S., Faloona, I., Mehrotra, S., Suard, M., Lenschow, D. H., Sweeney, C., Herndon, S., Schwietzke, S., Pétron, G., Pifer, J., Kort, E. A., and Schnell, R.: Application of Gauss's theorem to quantify localized surface emissions from airborne measurements of wind and trace gases, Atmos. Meas. Tech., 10, 3345–3358, https://doi.org/10.5194/amt-10-3345-2017, 2017. a, b

Deng, A., Stauffer, D., Dudhia, J., Hunter, G., and Bruyère, C.: WRF-ARW analysis nudging update and future development plan, Conference paper, available at: https://www.researchgate.net/publication/268409039_WRF-ARW_ANALYSIS_NUDGING_UPDATE_AND_FUTURE_DEVELOPMENT_PLAN (last access: 1 April 2019), 2008. a

DLR: Carbon Dioxide and Methane Mission for HALO, Mission 94, available at: https://halo-db.pa.op.dlr.de/mission/94 (last access: 20 February 2019), https://doi.org/10.17616/R39Q0T, 2019. a

Dulias, R.: A Brief History of Mining in the Upper Silesian Coal Basin, in: The Impact of Mining on the Landscape, Springer, 31–49, 2016. a

Fiehn, A., Kostinek, J., Eckl, M., Klausner, T., Gałkowski, M., Chen, J., Gerbig, C., Röckmann, T., Maazallahi, H., Schmidt, M., Korbeń, P., Neçki, J., Jagoda, P., Wildmann, N., Mallaun, C., Bun, R., Nickl, A.-L., Jöckel, P., Fix, A., and Roiger, A.: Estimating CH4, CO2 and CO emissions from coal mining and industrial activities in the Upper Silesian Coal Basin using an aircraft-based mass balance approach, Atmos. Chem. Phys., 20, 12675–12695, https://doi.org/10.5194/acp-20-12675-2020, 2020. a, b

Gzyl, G., Janson, E., and Łabaj, P.: Mine Water Discharges in Upper Silesian Coal Basin (Poland), chap. 17, in: Assessment, Restoration and Reclamation of Mining Influenced Soils, edited by: Bech, J., Bini, C., and Pashkevich, M. A., Academic Press, 463–486, https://doi.org/10.1016/B978-0-12-809588-1.00017-7, 2017. a, b

Hajny, K. D., Salmon, O. E., Rudek, J., Lyon, D. R., Stuff, A. A., Stirm, B. H., Kaeser, R., Floerchinger, C. R., Conley, S., Smith, M. L., and Shepson, P. B.: Observations of Methane Emissions from Natural Gas-Fired Power Plants, Environ. Sci. Technol., 53, 8976–8984, https://doi.org/10.1021/acs.est.9b01875, 2019. a

Jacob, D.: Inverse modeling techniques, in: Observing Systems for Atmospheric Composition, Springer, 230–237, 2007. a

Janssens-Maenhout, G., Crippa, M., Guizzardi, D., Muntean, M., Schaaf, E., Dentener, F., Bergamaschi, P., Pagliari, V., Olivier, J. G. J., Peters, J. A. H. W., van Aardenne, J. A., Monni, S., Doering, U., and Petrescu, A. M. R.: EDGAR v4.3.2 Global Atlas of the three major Greenhouse Gas Emissions for the period 1970–2012, Earth Syst. Sci. Data Discuss. [preprint], https://doi.org/10.5194/essd-2017-79, 2017. a

Johnson, M. R., Tyner, D. R., Conley, S., Schwietzke, S., and Zavala-Araiza, D.: Comparisons of Airborne Measurements and Inventory Estimates of Methane Emissions in the Alberta Upstream Oil and Gas Sector, Environ. Sci. Technol., 51, 13008–13017, https://doi.org/10.1021/acs.est.7b03525, 2017. a

Karion, A., Sweeney, C., Pétron, G., Frost, G., Hardesty, R. M., Kofler, J., Miller, B. R., Newberger, T., Wolter, S., Banta, R., Brewer, A., Dlugokencky, E., Lang, P., Montzka, S. A., Schnell, R., Tans, P., Trainer, M., Zamora, R., and Conley, S.: Methane emissions estimate from airborne measurements over a western United States natural gas field, Geophys. Res. Lett., 40, 4393–4397, https://doi.org/10.1002/grl.50811, 2013. a

Karion, A., Sweeney, C., Kort, E. A., Shepson, P. B., Brewer, A., Cambaliza, M., Conley, S. A., Davis, K., Deng, A., Hardesty, M., Herndon, S. C., Lauvaux, T., Lavoie, T., Lyon, D., Newberger, T., Pétron, G., Rella, C., Smith, M., Wolter, S., Yacovitch, T. I., and Tans, P.: Aircraft-Based Estimate of Total Methane Emissions from the Barnett Shale Region, Environ. Sci. Technol., 49, 8124–8131, https://doi.org/10.1021/acs.est.5b00217, 2015. a

Kirschke, S., Bousquet, P., Ciais, P., Saunois, M., Canadell, J., Dlugokencky, E., Bergamaschi, P., Bergmann, D., Blake, D., Bruhwiler, L., Cameron-Smith, P., Castaldi, S., Chevallier, F., Feng, L., Fraser, A., Heimann, M., Hodson, E., Houweling, S., Josse, B., and Zeng, G.: Three decades of global methane sources and sinks, Nat. Geosci., 6, 813–823, https://doi.org/10.1038/ngeo1955, 2013. a

Kostinek, J., Roiger, A., Davis, K. J., Sweeney, C., DiGangi, J. P., Choi, Y., Baier, B., Hase, F., Groß, J., Eckl, M., Klausner, T., and Butz, A.: Adaptation and performance assessment of a quantum and interband cascade laser spectrometer for simultaneous airborne in situ observation of CH4, C2H6, CO2, CO and N2O, Atmos. Meas. Tech., 12, 1767–1783, https://doi.org/10.5194/amt-12-1767-2019, 2019. a

Lavoie, T. N., Shepson, P. B., Cambaliza, M. O. L., Stirm, B. H., Karion, A., Sweeney, C., Yacovitch, T. I., Herndon, S. C., Lan, X., and Lyon, D.: Aircraft-Based Measurements of Point Source Methane Emissions in the Barnett Shale Basin, Environ. Sci. Technol., 49, 7904–7913, https://doi.org/10.1021/acs.est.5b00410, 2015. a

Lawson, C. L. and Hanson, R. J.: Solving Least Squares Problems, Society for Industrial and Applied Mathematics, 337 pp., https://doi.org/10.1137/1.9781611971217, 1995. a

Liu, J., Bray, M., and Han, D.: A study on WRF radar data assimilation for hydrological rainfall prediction, Hydrol. Earth Syst. Sci., 17, 3095–3110, https://doi.org/10.5194/hess-17-3095-2013, 2013. a

Luther, A., Kleinschek, R., Scheidweiler, L., Defratyka, S., Stanisavljevic, M., Forstmaier, A., Dandocsi, A., Wolff, S., Dubravica, D., Wildmann, N., Kostinek, J., Jöckel, P., Nickl, A.-L., Klausner, T., Hase, F., Frey, M., Chen, J., Dietrich, F., Nȩcki, J., Swolkień, J., Fix, A., Roiger, A., and Butz, A.: Quantifying CH4 emissions from hard coal mines using mobile sun-viewing Fourier transform spectrometry, Atmos. Meas. Tech., 12, 5217–5230, https://doi.org/10.5194/amt-12-5217-2019, 2019. a, b

Mallaun, C., Giez, A., and Baumann, R.: Calibration of 3-D wind measurements on a single-engine research aircraft, Atmos. Meas. Tech., 8, 3177–3196, https://doi.org/10.5194/amt-8-3177-2015, 2015. a

Mehrotra, S., Faloona, I., Suard, M., Conley, S., and Fischer, M. L.: Airborne Methane Emission Measurements for Selected Oil and Gas Facilities Across California, Environ. Sci. Technol., 51, 12981–12987, https://doi.org/10.1021/acs.est.7b03254, 2017. a

National Centers for Environmental Prediction/National Weather Service/NOAA/U.S. Department of Commerce: NCEP GDAS/FNL 0.25 Degree Global Tropospheric Analyses and Forecast Grids, Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory, https://doi.org/10.5065/D65Q4T4Z, 2015. a

Nickl, A.-L., Mertens, M., Roiger, A., Fix, A., Amediek, A., Fiehn, A., Gerbig, C., Galkowski, M., Kerkweg, A., Klausner, T., Eckl, M., and Jöckel, P.: Hindcasting and forecasting of regional methane from coal mine emissions in the Upper Silesian Coal Basin using the online nested global regional chemistry–climate model MECO(n) (MESSy v2.53), Geosci. Model Dev., 13, 1925–1943, https://doi.org/10.5194/gmd-13-1925-2020, 2020. a

Nisbet, E. and Weiss, R.: Atmospheric science. Top-down versus bottom-up, Science, 328, 1241–1243, https://doi.org/10.1126/science.1189936, 2010. a

Nisbet, E. G., Dlugokencky, E. J., and Bousquet, P.: Methane on the Rise–Again, Science, 343, 493–495, https://doi.org/10.1126/science.1247828, 2014. a

Pachauri, R. K., Allen, M. R., Barros, V. R., Broome, J., Cramer, W., Christ, R., Church, J. A., Clarke, L., Dahe, Q., Dasgupta, P., Dubash, N. K., Edenhofer, O., Elgizouli, I., Field, C. B., Forster, P., Friedlingstein, P., Fuglestvedt, J., Gomez-Echeverri, L., Hallegatte, S., Hegerl, G., Howden, M., Jiang, K., Cisneroz, B. J., Kattsov, V., Lee, H., Mach, K. J., Marotzke, J., Mastrandrea, M. D., Meyer, L., Minx, J., Mulugetta, Y., O'Brien, K., Oppenheimer, M., Pereira, J. J., Pichs-Madruga, R., Plattner, G.-K., Pörtner, H.-O., Power, S. B., Preston, B., Ravindranath, N. H., Reisinger, A., Riahi, K., Rusticucci, M., Scholes, R., Seyboth, K., Sokona, Y., Stavins, R., Stocker, T. F., Tschakert, P., van Vuuren, D., and van Ypserle, J.-P.: Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, IPCC, Geneva, Switzerland, 2014. a, b

Peters, W., Jacobson, A. R., Sweeney, C., Andrews, A. E., Conway, T. J., Masarie, K., Miller, J. B., Bruhwiler, L. M. P., Pétron, G., Hirsch, A. I., Worthy, D. E. J., van der Werf, G. R., Randerson, J. T., Wennberg, P. O., Krol, M. C., and Tans, P. P.: An atmospheric perspective on North American carbon dioxide exchange: CarbonTracker, P. Natl. Acad. Sci. USA, 104, 18925–18930, https://doi.org/10.1073/pnas.0708986104, 2007. a

Pisso, I., Sollum, E., Grythe, H., Kristiansen, N. I., Cassiani, M., Eckhardt, S., Arnold, D., Morton, D., Thompson, R. L., Groot Zwaaftink, C. D., Evangeliou, N., Sodemann, H., Haimberger, L., Henne, S., Brunner, D., Burkhart, J. F., Fouilloux, A., Brioude, J., Philipp, A., Seibert, P., and Stohl, A.: The Lagrangian particle dispersion model FLEXPART version 10.4, Geosci. Model Dev., 12, 4955–4997, https://doi.org/10.5194/gmd-12-4955-2019, 2019. a

Pitt, J. R., Allen, G., Bauguitte, S. J.-B., Gallagher, M. W., Lee, J. D., Drysdale, W., Nelson, B., Manning, A. J., and Palmer, P. I.: Assessing London CO2, CH4 and CO emissions using aircraft measurements and dispersion modelling, Atmos. Chem. Phys., 19, 8931–8945, https://doi.org/10.5194/acp-19-8931-2019, 2019. a

Plant, G., Kort, E. A., Floerchinger, C., Gvakharia, A., Vimont, I., and Sweeney, C.: Large Fugitive Methane Emissions From Urban Centers Along the U.S. East Coast, Geophys. Res. Lett., 46, 8500–8507, https://doi.org/10.1029/2019GL082635, 2019. a

Platt, S. M., Eckhardt, S., Ferré, B., Fisher, R. E., Hermansen, O., Jansson, P., Lowry, D., Nisbet, E. G., Pisso, I., Schmidbauer, N., Silyakova, A., Stohl, A., Svendby, T. M., Vadakkepuliyambatta, S., Mienert, J., and Lund Myhre, C.: Methane at Svalbard and over the European Arctic Ocean, Atmos. Chem. Phys., 18, 17207–17224, https://doi.org/10.5194/acp-18-17207-2018, 2018. a

Powers, J. G., Klemp, J. B., Skamarock, W. C., Davis, C. A., Dudhia, J., Gill, D. O., Coen, J. L., Gochis, D. J., Ahmadov, R., Peckham, S. E., Grell, G. A., Michalakes, J., Trahan, S., Benjamin, S. G., Alexander, C. R., Dimego, G. J., Wang, W., Schwartz, C. S., Romine, G. S., Liu, Z., Snyder, C., Chen, F., Barlage, M. J., Yu, W., and Duda, M. G.: The Weather Research and Forecasting Model: Overview, System Efforts, and Future Directions, B. Am. Meteorol. Soc., 98, 1717–1737, https://doi.org/10.1175/BAMS-D-15-00308.1, 2017. a

Raut, J.-C., Marelle, L., Fast, J. D., Thomas, J. L., Weinzierl, B., Law, K. S., Berg, L. K., Roiger, A., Easter, R. C., Heimerl, K., Onishi, T., Delanoë, J., and Schlager, H.: Cross-polar transport and scavenging of Siberian aerosols containing black carbon during the 2012 ACCESS summer campaign, Atmos. Chem. Phys., 17, 10969–10995, https://doi.org/10.5194/acp-17-10969-2017, 2017. a

Ren, X., Salmon, O. E., Hansford, J. R., Ahn, D., Hall, D., Benish, S. E., Stratton, P. R., He, H., Sahu, S., Grimes, C., Heimburger, A. M. F., Martin, C. R., Cohen, M. D., Stunder, B., Salawitch, R. J., Ehrman, S. H., Shepson, P. B., and Dickerson, R. R.: Methane Emissions From the Baltimore-Washington Area Based on Airborne Observations: Comparison to Emissions Inventories, J. Geophys. Res.-Atmos., 123, 8869–8882, https://doi.org/10.1029/2018JD028851, 2018. a

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding, World scientific, 256 pp., https://doi.org/10.1142/3171, 2000. a

Ryoo, J.-M., Iraci, L. T., Tanaka, T., Marrero, J. E., Yates, E. L., Fung, I., Michalak, A. M., Tadić, J., Gore, W., Bui, T. P., Dean-Day, J. M., and Chang, C. S.: Quantification of CO2 and CH4 emissions over Sacramento, California, based on divergence theorem using aircraft measurements, Atmos. Meas. Tech., 12, 2949–2966, https://doi.org/10.5194/amt-12-2949-2019, 2019.  a

Schwietzke, S., Sherwood, O., Bruhwiler, L., Miller, J., Etiope, G., Dlugokencky, E., Englund Michel, S., Arling, V., Vaughn, B., White, J., and Tans, P.: Upward revision of global fossil fuel methane emissions based on isotope database, Nature, 538, 88–91, https://doi.org/10.1038/nature19797, 2016. a

Swolkień, J.: Polish underground coal mines as point sources of methane emission to the atmosphere, Int. J. Greenh. Gas Con., 94, 102921, https://doi.org/10.1016/j.ijggc.2019.102921, 2020. a

Tarantola, A.: Inverse Problem Theory and Methods for Model Parameter Estimation, Society for Industrial and Applied Mathematics, USA, 2004. a

Tuccella, P., Thomas, J., Law, K. S., Raut, J.-C., Marelle, L., Roiger, A., Weinzierl, B., Denier Van Der Gon, H., Schlager, H., and Onishi, T.: Air pollution impacts due to petroleum extraction in the Norwegian Sea during the ACCESS aircraft campaign, Elementa, 5, 25, https://doi.org/10.1525/elementa.124.s1, 2017. a

University Corporation for Atmospheric Research: WRF v4.0, available at: https://www.mmm.ucar.edu/weather-research-and-forecasting-model, last access: 6 May 2020. a

Wildmann, N., Päschke, E., Roiger, A., and Mallaun, C.: Towards improved turbulence estimation with Doppler wind lidar velocity-azimuth display (VAD) scans, Atmos. Meas. Tech., 13, 4141–4158, https://doi.org/10.5194/amt-13-4141-2020, 2020. a, b

Yu, X., Millet, D. B., Wells, K. C., Griffis, T. J., Chen, X., Baker, J. M., Conley, S. A., Smith, M. L., Gvakharia, A., Kort, E. A., Plant, G., and Wood, J. D.: Top-Down Constraints on Methane Point Source Emissions From Animal Agriculture and Waste Based on New Airborne Measurements in the U.S. Upper Midwest, J. Geophys. Res.-Biogeo., 125, e2019JG005429, https://doi.org/10.1029/2019JG005429, 2020. a

Zentralanstalt für Meteorologie und Geodynamik: FLEXPART-WRF, available at: https://www.flexpart.eu/, last access: 1 February 2019. a

Download
Short summary
Abundant mining and industrial activities in the Upper Silesian Coal Basin lead to large emissions of the potent greenhouse gas methane. This study quantifies these emissions with continuous, high-precision airborne measurements and dispersion modeling. Our emission estimates are in line with values reported in the European Pollutant Release and Transfer Register (E-PRTR 2017) but significantly lower than values reported in the Emissions Database for Global Atmospheric Research (EDGAR v4.3.2).
Altmetrics
Final-revised paper
Preprint