An Inverse Modeling Approach for Retrieving High-Resolution Surface Fluxes of Greenhouse Gases from Measurements of Their Concentrations in the Atmospheric Boundary Layer

This study explores the potential of using Unmanned Aircraft Vehicles (UAVs) as a measurement platform for estimating greenhouse gas (GHG) fluxes over complex terrain. We proposed and tested an inverse modeling approach for retrieving GHG fluxes based on two-level measurements of GHG concentrations and airflow properties over complex terrain with high spatial resolution. Our approach is based on a three-dimensional hydrodynamic model capable of determining the airflow parameters that affect the spatial distribution of GHG concentrations within the atmospheric boundary layer. The model is primarily designed to solve the forward problem of calculating the steady-state distribution of GHG concentrations and fluxes at different levels over an inhomogeneous land surface within the model domain. The inverse problem deals with determining the unknown surface GHG fluxes by minimizing the difference between measured and modeled GHG concentrations at two selected levels above the land surface. Several numerical experiments were conducted using surrogate data that mimicked UAV observations of varying accuracies and density of GHG concentration measurements to test the robustness of the approach. Our primary modeling target was a 6 km² forested area in the foothills of the Greater Caucasus Mountains in Russia, characterized by complex topography and mosaic vegetation. The numerical experiments show that the proposed inverse modeling approach can effectively solve the inverse problem, with the resulting flux distribution having the same spatial pattern as the required flux. However, the approach tends to overestimate the mean value of the required flux over the domain, with the maximum errors in flux estimation associated with areas of maximum steepness in the surface topography. The accuracy of flux estimates improves as the number of points and the accuracy of the concentration measurements increase. Therefore, the density of UAV measurements should be adjusted according to the complexity of the terrain to improve the accuracy of the modeling results.