Search NASA⌕ Search

Engineering topics

Michel Bechtold

Publications and source records attributed to Michel Bechtold.

Contributions of irrigation modeling, soil moisture and snow data assimilation to high-resolution water budget estimates over the Po basin: progress towards digital replicas

High-resolution water budget estimates benefit from modeling of human water management and satellite data assimilation (DA) in river basins with a large human footprint. Utilizing the Noah-MP land surface model with dynamic vegetation growth and river routing, in combination with an irrigation module, Sentinel-1 backscatter and snow depth retrievals, we produce a set of 0.7-km2 water budget estimates of the Po river basin (Italy) for 2015–2023. The results demonstrate that irrigation modeling improves the seasonal soil moisture variation and summer streamflow at all gauges in the valley after withdrawal of irrigation water from the streamflow in postprocessing (12% error reduction relative to observed low summer streamflow), even if the basin-wide irrigation amount is underestimated. Sentinel-1 backscatter DA for soil moisture updating strongly interacts with irrigation modeling: when both are activated, the soil moisture updates are limited, and the simulated irrigation amounts are reduced. Backscatter DA systematically reduces soil moisture in the spring, which improves downstream spring streamflow. Assimilating Sentinel-1 snow depth retrievals over the surrounding Alps and Apennines further improves spring streamflow in a complementary way (2% error reduction relative to observed high spring streamflow). Despite the seasonal improvements, irrigation modeling and Sentinel-1 backscatter DA cannot significantly improve short-term or interannual variations in soil moisture, irrigation modeling causes a systematically prolonged high vegetation productivity, and snow depth DA only impacts the deep snowpacks. This study helps advancing the design of digital water budget replicas for river basins.

Gabrielle J M De Lannoy↗

Improving the Fire Weather Index System for Peatlands Using Peat-Specific Hydrological Input Data

The Canadian Fire Weather Index (FWI) system, even though originally developed and calibrated for an upland Jack pine forest, is used globally to estimate fire danger for any fire environment. However, for some environments, such as peatlands, the applicability of the FWI in its current form, is often questioned. In this study, we replaced the original moisture codes of the FWI with hydrological estimates resulting from the assimilation of satellite-based L-band passive microwave observations into a peatland-specific land surface model. In a conservative approach that maintains the integrity of the original FWI structure, the distributions of the hydrological estimates were first matched to those of the corresponding original moisture codes before replacement. The resulting adapted FWI, hereafter called FWIpeat, was evaluated using satellite-based information on fire presence over boreal peatlands from 2010 through 2018. Adapting the FWI with model- and satellite-based hydrological information was found to be beneficial in estimating fire danger, especially when replacing the deeper moisture codes of the FWI. For late-season fires, further adaptations of the fine fuel moisture code show even more improvement due to the fact that late-season fires are more hydrologically driven. The proposed FWIpeat should enable improved monitoring of fire risk in boreal peatlands.

Canadian Fire Weather Index↗

Impact of Design Factors for ESA CCI Satellite Soil Moisture Data Assimilation over Europe

In this study, soil moisture retrievals of the combined active–passive ESA Climate Change Initiative (CCI) soil moisture product are assimilated into the Noah-MP land surface model over Europe using a one-dimensional ensemble Kalman filter and an 18-yr study period. The performance of the data assimilation (DA) system is evaluated by comparing it with a model-only experiment (at in situ sites) and by assessing statistics of innovations and increments as DA diagnostics (over the entire domain). For both assessments, we explore the impact of three design choices, resulting in the following insights. 1) The magnitude of the assumed observation errors strongly affects the skill improvements evaluated against in situ stations and internal diagnostics. 2) Choosing between climatological or monthly cumulative distribution function matching as the observation bias correction method only has a marginal effect on the in situ skill of the DA system. However, the internal diagnostics suggest a more robust system parameterization if the observations are rescaled monthly. 3) The choice of atmospheric reanalysis dataset to force the land surface model affects the model-only skill and the DA skill improvements. The model-only skill is higher with input from the MERRA-2 than with input from the ERA5 reanalysis, resulting in larger DA skill improvements for the latter. Additionally, we show that the added value of the DA strongly depends on the quality of the satellite retrievals and land cover, with the most substantial soil moisture skill improvements occurring over croplands and skill degradation occurring over densely forested areas.

Soil moisture↗

Lightning Over Central Canada: Skill Assessment for Various Land-Atmosphere Model Configurations and Lightning Indices Over a Boreal Study Area

Current lightning predictions are uncertain because they rely on empirical diagnostic relationships and often use coarse-scale climate scenario simulations in which deep convection is parameterized. Previous studies demonstrated that simulations with convection permitting resolutions improve lightning predictions compared to coarser-grid simulations using convection parameterizations for different geographical locations but not over the boreal zone. In this study, lightning simulations with the NASA Unified-Weather Research and Forecasting (NU-WRF) model are evaluated over a domain including the Great Slave Lake in Canada, for six lightning seasons. The simulations are performed at convection parameterized (9 km) and convection-permitting (3 km) resolution using the Goddard 4ICE and the Thompson microphysics schemes. Four lightning indices are evaluated against observations from the Canadian Lightning Detection Network, in terms of spatiotemporal frequency distribution, spatial pattern, daily climatology, and an event-based overall skill assessment. The Thompson scheme is, regardless of the spatial resolution, superior to the Goddard 4ICE scheme in predicting daily climatology but worse in predicting the spatial patterns of lightning occurrence. Results indicate that lightning estimation benefits from modeling at convection-permitting resolution, in particular for the ice based lightning indices. In contrast, the product of convective available potential energy and precipitation rate proved to be the most robust index that was largely invariant to varying spatial resolution. Finally, this study reveals issues of the models to reproduce the observed spatial pattern of lightning well, which might be related to an insufficient representation of land surface heterogeneity, including peatlands, in the study area.

Jonas Mortelmans↗

Perspective on Satellite-Based Land Data Assimilation to Estimate Water Cycle Components in an Era of Advanced Data Availability and Model Sophistication

The beginning of the 21st century is marked by a rapid growth of land surface satellite data and model sophistication. This offers new opportunities to estimate multiple components of the water cycle via satellite-based land data assimilation (DA) across multiple scales. By resolving more processes in land surface models and by coupling the land, the atmosphere, and other Earth system compartments, the observed information can be propagated to constrain additional unobserved variables. Furthermore, access to more satellite observations enables the direct constraint of more and more components of the water cycle that are of interest to end users. However, the finer level of detail in models and data is also often accompanied by an increase in dimensions, with more state variables, parameters, or boundary conditions to estimate, and more observations to assimilate. This requires advanced DA methods and efficient solutions. One solution is to target specific observations for assimilation based on a sensitivity study or coupling strength analysis, because not all observations are equally effective in improving subsequent forecasts of hydrological variables, weather, agricultural production, or hazards through DA. This paper offers a perspective on current and future land DA development, and suggestions to optimally exploit advances in observing and modeling systems.

Data Assimilation↗

L-Band Microwave Satellite Data and Model Simulations Over the Dry Chaco to Estimate Soil Moisture, Soil Temperature, Vegetation and Soil Salinity

The Dry Chaco in South America is a semi-arid ecoregion prone to dryland salinization. In this region, we investigated coarse-scale surface soil moisture (SM), soil temperature, soil salinity and vegetation, using L-band microwave brightness temperature (TB) observations and retrievals from the Soil Moisture Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) satellite missions, Catchment Land Surface Model (CLSM) simulations, and in situ measurements within 26 sampled satellite pixels. Across these 26 sampled pixels, the satellite-based SM outperformed CLSM SM compared to field data, and forward L-band TB simulations derived from in situ SM and temperature performed better than those derived from CLSM estimates relative to SMOS TB observations. The surface salinity for the sampled pixels was on average only 4 mg/g and only locally influenced the TB simulations, when including salinity in the dielectric mixing model of the forward radiative transfer model (RTM) simulations. To explore the potential of retrieving salinity together with other RTM parameters to optimize TB simulations over the entire Dry Chaco, the RTM was inverted using 10 years of multi-angular SMOS TB data and constraints of CLSM SM and temperature. However, the latter modeled SM was not sufficiently accurate and factors such as open surface water were missing in the background constraints, so that the salinity retrievals effectively represented a bulk correction of the dielectric constant, rather than salinity per se. However, the retrieval of vegetation, scattering albedo and surface roughness resulted in realistic values.

L-Band↗

FORWARD AND INVERSE L-BAND RADIATIVE TRANSFER MODELING OVER THE DRY CHACO, USING SMOS OBSERVATIONS, LAND SURFACE MODELING AND IN SITU DATA

Passive microwave L-band remote sensing is well known for its sensitivity to surface soil moisture over land.The signal is also affected by other dynamic variables such as vegetation and soil salinity. In this research, L-band microwave observations of the Soil Moisture Ocean Salinity (SMOS) mission are used to explore soil surface salinity, moisture and vegetation, in the Argentinean Dry Chaco, an area with possible emerging dryland salinity. A Radiative Transfer Model (RTM) with inclusion of a correction to the soil’s dielectric constant for salinity was used in forward and inverse mode, using either in situ data or Catchment Land Surface Model (CLSM) simulations as RTM input. The forward analysis pointed out shortcomings in the modeled soil moisture and soil surface temperature estimates. The impact of salinity on forward simulations over the Dry Chaco was limited.The RTM inversion using 10 years of SMOS brightness temperature observations resulted in realistic estimates of vegetation and roughness, both with and withouto ptimizing a correction term for the dielectric constant in terms of salinity equivalents. However,the retrieval of the correction term to the dielectric constant was very uncertain and not representative of soil surface salinity, but rather of open water, texture uncertainty or soil moisture bias.

L-band↗

An Operational Product for Peatland Applications: The Version 7 of the SMAP Level-4 Soil Moisture Data Assimilation Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 9-km resolution, 3-hourly surface (0-5 cm) and root-zone (0-100 cm) soil moisture estimates from April 2015 to present with a mean latency of 2.5 days from the time of observation. The L4_SM estimates are derived from the assimilation of SMAP L-band (1.4 GHz) brightness temperature (Tb) observations into the NASA Catchment Land Surface Model (CLSM). The recently released Version 7 of the L4_SM product comprises a key advancement that is of interest for peatland applications. CLSM now includes the recently developed PEATCLSM hydrology module for peatlands and uses an updated global peatland distribution. In this presentation, we first give an overview of the operational product for the peatland community. We further show how the incorporation of PEATCLSM considerably improves the dynamics of water table depth, surface soil moisture and evapotranspiration in L4_SM Version 7 over Version 6 when evaluated against in situ measurements in peatlands. These improvements are also manifested in smaller Tb observation-minus-forecast residuals. Eventually, we provide three research examples in which the data assimilation product is used to gain insights into: peatland hydrological processes, the peat moisture dependency of the carbon cycle, and wildfire occurrence in peatlands.

Michel Bechtold↗

Insights Into the Hydrology of the Congo Peatlands Through Land Surface Modeling and Data Assimilation

The 16.8 million ha of peatlands in the Cuvette Centrale wetland complex in the Congo Basin is one of the largest peatland regions on Earth but still highly understudied. Understanding the hydrological functioning of these peatlands and the effects of external disturbances thereon remains a major challenge. Recent research suggested fundamental hydrological differences between the Congo peatlands and the well-studied Southeast Asian peatlands. The Congo peatlands have a doming gradient that is up to ten times smaller, and they are influenced by river hydrology to some extent. In this study, we explore the Congo peatland hydrology through land surface modeling and data assimilation. We build upon our recently developed tropical PEATCLSM module (Apers et al., 2022) that was parameterized based on data from Southeast Asian peatlands due to the lack of field data from other tropical peatland regions. In a first step, we derive Congo-specific peat hydraulic and discharge function parameters from a scalar parametrization of water level dynamics in the Congo peatlands, using observed water level data at two locations. These Congo-specific parameters differ considerably from the original literature-based parameters from Southeast Asian peatlands. In a second step, we apply our original and Congo-specific parameters in an assimilation scheme for L-band brightness temperature (Tb) data from the Soil Moisture and Ocean Salinity (SMOS) mission. The data assimilation results are used in two ways. First, the effect of these parameters on the simulated peatland hydrology and the observation-minus-forecast Tb residuals is evaluated. It is hypothesized that the new parameters reduce the previously reported modeling errors over the Congo peatlands and reduce the residuals in Tb as well. Second, we analyze the data assimilation diagnostics to learn about other model improvement possibilities. In preliminary results, we found long periods of temporally autocorrelated total water storage increments (difference of pre- and post-update) that coincided with anomalies in river stages measured upstream of the peatlands. Since PEATLCSM neglects possible river influence, this concurrence suggests that the typically used grid-based approach of land surface models should be combined with a river routing scheme over the Congo peatlands.

Sebastian Apers↗

Version 8 of the SMAP Level-4 Soil Moisture Data Assimilation Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 9-km resolution, 3-hourly surface (0-5 cm) and root-zone (0-100 cm) soil moisture from April 2015 to present with a mean latency of 2.5 days from the time of observation. The product is based on the assimilation of SMAP L-band (1.4 GHz) brightness temperature (Tb) observations into the NASA Catchment land surface model as the model is driven with observations-based precipitation forcing. This presentation discusses the improvements in the forthcoming Version 8 of L4_SM, including updates in the precipitation forcing, the Catchment model parameters, and the L-band microwave radiative transfer model (mwRTM). The precipitation observations used in L4_SM Version 8 outside of North America and the high latitudes are from the latest (Version 7) NASA Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement mission (IMERG) products. Moreover, occasionally excessive precipitation rates in earlier versions of L4_SM along certain longitudes in North America were eliminated by a bug fix in the precipitation corrections algorithm. The Catchment model in L4_SM Version 8 uses climatological snow albedo values based on observations from the Moderate Resolution Imaging Spectroradiometer, replacing the look-up table parameterization of earlier versions. Additionally, corrected soil parameters were implemented for a small region in Argentina that had erroneously been classified as peat because of an error in the ancillary soil data. Finally, the mwRTM in L4_SM Version 8 uses the Mironov soil mixing approach and updated values of the L-band scattering albedo, soil roughness, and vegetation opacity climatology obtained from the latest (Release 19) SMAP Level-2 dual-channel soil moisture retrieval product. During the development of L4_SM Version 8, the change in the mwRTM parameterization resulted in a reduced unbiased RMSE of surface soil moisture when verified against in situ measurements. It also reduced the standard deviation of the Tb observation-minus-forecast residuals by ~0.15 K, highlighting the importance of the mwRTM for successful data assimilation. The bug fixes in the precipitation corrections algorithm and the Catchment model soil parameters in the Argentina region result in locally large improvements of the simulated land surface states. In summary, the ongoing refinements of the L4_SM product continue to improve its science quality and performance for global soil moisture monitoring

Rolf Reichle↗