Quantifying methane point sources from fine-scale satellite observations of atmospheric methane plumes

Anthropogenic methane emissions originate from a large number of relatively small point sources. The planned GHGSat satellite fleet aims to quantify emissions from individual point sources by measuring methane column plumes over selected similar to 10 x 10 km(2) domains with <= 50 x 50 m(2) pixel resolution and 1 %-5 % measurement precision. Here we develop algorithms for retrieving point source rates from such measurements. We simulate a large ensemble of instantaneous methane column plumes at 50 x 50 m(2) pixel resolution for a range of atmospheric conditions using the Weather Research and Forecasting model (WRF) in large eddy simulation (LES) mode and adding instrument noise. We show that standard methods to infer source rates by Gaussian plume inversion or source pixel mass balance are prone to large errors because the turbulence cannot be properly parameterized on the small scale of instantaneous methane plumes. The integrated mass enhancement (IME) method, which relates total plume mass to source rate, and the cross-sectional flux method, which infers source rate from fluxes across plume transects, are better adapted to the problem. We show that the IME method with local measurements of the 10 m wind speed can infer source rates with an error of 0.07-0.17 th(-1) + 5 %-12 % depending on instrument precision (1 %-5 %). The cross-sectional flux method has slightly larger errors (0.07-0.26 th(-1) + 8 %-12 %) but a simpler physical basis. For comparison, point sources larger than 0.3 th(-1) contribute more than 75 % of methane emissions reported to the US Greenhouse Gas Reporting Program. Additional error applies if local wind speed measurements are not available and may dominate the overall error at low wind speeds. Low winds are beneficial for source detection but detrimental for source quantification.