Inferring surface energy fluxes using drone data assimilation in large eddy simulations

• Main Authors: Pirk, Norbert, Aalstad, Kristoffer, Westermann, Sebastian, Vatne, Astrid, van Hove, Alouette, Tallaksen, Lena Merete, Cassiani, Massimo, Katul, Gabriel G
• Format: Article
• Language: English
• ISSN: 7293-7314

Spatially representative estimates of surface energy exchange from field measurements are required for improving and validating Earth system models and satellite remote sensing algorithms. The scarcity of flux measurements can limit understanding of ecohydrological responses to climate warming, especially in remote regions with limited infrastructure. Direct field measurements often apply the eddy covariance method on stationary towers, but recently, drone-based measurements of temperature, humidity, and wind speed have been suggested as a viable alternative to quantify the turbulent fluxes of sensible (H) and latent heat (LE). A data assimilation framework to infer uncertainty-aware surface flux estimates from sparse and noisy drone-based observations is developed and tested using a turbulence-resolving large eddy simulation (LES) as a forward model to connect surface fluxes to drone observations. The proposed framework explicitly represents the sequential collection of drone data, accounts for sensor noise, includes uncertainty in boundary and initial conditions, and jointly estimates the posterior distribution of a multivariate parameter space. Assuming typical flight times and observational errors of light-weight, multi-rotor drone systems, we first evaluate the information gain and performance of different ensemble-based data assimilation schemes in experiments with synthetically generated observations. It is shown that an iterative ensemble smoother outperforms both the non-iterative ensemble smoother and the particle batch smoother in the given problem, yielding well-calibrated posterior uncertainty with continuous ranked probability scores of 12 W m−2 for both H and LE, with standard deviations of 37 W m−2 (H) and 46 W m−2 (LE) for a 12 min vertical step profile by a single drone. Increasing flight times, using observations from multiple drones, and further narrowing the prior distributions of the initial conditions are viable for reducing the posterior spread. Sampling strategies prioritizing space–time exploration without temporal averaging, instead of hovering at fixed locations while averaging, enhance the non-linearities in the forward model and can lead to biased flux results with ensemble-based assimilation schemes. In a set of 18 real-world field experiments at two wetland sites in Norway, drone data assimilation estimates agree with independent eddy covariance estimates, with root mean square error values of 37 W m−2 (H), 52 W m−2 (LE), an