Sentinel-1 snow depth assimilation improves river discharge simulations in the western Alps

Author(s):  
Isis Brangers ◽  
Hans Lievens ◽  
Augusto Getirana ◽  
Sujay Kumar ◽  
Gabrielle De Lannoy

<p>In many of the world’s mountainous regions, river discharge is largely influenced by the seasonal melt of snow. Therefore, accurate information on the amount of water stored as snow is essential for water management and flood forecasting. However, there are large uncertainties in model simulations of snow depth, partly due to uncertain precipitation estimates in mountain regions with complex topography. A study by Lievens et al. (2019) showed the potential of Sentinel-1 (S1) satellite observations to provide snow depth estimates at 1 km spatial and ~weekly temporal resolution in mountain regions. In this study, we assimilated these retrievals into the Noah Multiparameterization (Noah-MP) v3.6 land surface model for the western Alps using an ensemble Kalman filter. The land surface model was coupled to the Hydrological Modeling and Analysis Platform (HyMAP) routing scheme to also provide estimates of river discharge. With S1 data assimilation, the snow depth estimates improved, reducing the bias from 0.23 m to 0.05 m compared to in situ measurements. Preliminary results also show improved discharge simulations mainly in mountain catchments at high elevations that are less prone to regulations (e.g., by dams). This study demonstrates the capability of the S1 snow depth retrievals to improve not only snow depth estimates, but also the estimation of snow melt water contributions to river discharge.</p>

2018 ◽  
Vol 22 (7) ◽  
pp. 3863-3882 ◽  
Author(s):  
Fuxing Wang ◽  
Jan Polcher ◽  
Philippe Peylin ◽  
Vladislav Bastrikov

Abstract. River discharge plays an important role in earth's water cycle, but it is difficult to estimate due to un-gauged rivers, human activities and measurement errors. One approach is based on the observed flux and a simple annual water balance model (ignoring human processes) for un-gauged rivers, but it only provides annual mean values which is insufficient for oceanic modelings. Another way is by forcing a land surface model (LSM) with atmospheric conditions. It provides daily values but with uncertainties associated with the models. We use data assimilation techniques by merging the modeled river discharges by the ORCHIDEE (without human processes currently) LSM and the observations from the Global Runoff Data Centre (GRDC) to obtain optimized discharges over the entire basin. The “model systematic errors” and “human impacts” (dam operation, irrigation, etc.) are taken into account by an optimization parameter x (with annual variation), which is applied to correct model intermediate variable runoff and drainage over each sub-watershed. The method is illustrated over the Iberian Peninsula with 27 GRDC stations over the period 1979–1989. ORCHIDEE represents a realistic discharge over the north of the Iberian Peninsula with small model systematic errors, while the model overestimates discharges by 30–150 % over the south and northeast regions where the blue water footprint is large. The normalized bias has been significantly reduced to less than 30 % after assimilation, and the assimilation result is not sensitive to assimilation strategies. This method also corrects the discharge bias for the basins without observations assimilated by extrapolating the correction from adjacent basins. The “correction” increases the interannual variability in river discharge because of the fluctuation of water usage. The E (P−E) of GLEAM (Global Land Evaporation Amsterdam Model, v3.1a) is lower (higher) than the bias-corrected value, which could be due to the different P forcing and probably the missing processes in the GLEAM model.


2021 ◽  
Author(s):  
Gabriele Arduini ◽  
Ervin Zsoter ◽  
Hannah Cloke ◽  
Elisabeth Stephens ◽  
Christel Prudhomme

<p>Snow processes, with the water stored in the snowpack and released as snowmelt, are very important components of the water balance, in particular in high latitude and mountain regions. The evolution of the snow cover and the timing of the snow melt can have major impact on river discharge. Land surface models are used in Earth System models to compute exchanges of water, energy and momentum between the atmosphere and the surface underneath, and also to compute other components of the hydrological cycle. In order to improve the snow representation, a new multi-layer snow scheme is under development in the HTESSEL land surface model of the European Centre for Medium‐Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS), to replace the current single-layer snow scheme used in HTESSEL. The new scheme has already been shown to improve snow and 2‐metre temperature, while in this study, the wider hydrological impact is evaluated and documented.</p><p>The analysis is done in the reanalysis context by comparing two ERA5-forced offline HTESSEL experiments. The runoff output of HTESSEL is coupled to the CaMa-Flood hydrodynamic model in order to derive river discharge. The analysis is done globally for the period between 1980-2018. The evaluation was carried out using over 1000 discharge observation time-series with varying catchment size. The hydrological response of the multi-layer snow scheme is generally positive, but in some areas the improvement is not clear and can even be negative with deteriorated signal in river discharge. Further investigation is needed to understand the complex hydrological impact of the new snow scheme, making sure it contributes to an improved description of all hydrological components of the Earth System.</p>


2017 ◽  
Vol 21 (4) ◽  
pp. 2015-2033 ◽  
Author(s):  
David Fairbairn ◽  
Alina Lavinia Barbu ◽  
Adrien Napoly ◽  
Clément Albergel ◽  
Jean-François Mahfouf ◽  
...  

Abstract. This study evaluates the impact of assimilating surface soil moisture (SSM) and leaf area index (LAI) observations into a land surface model using the SAFRAN–ISBA–MODCOU (SIM) hydrological suite. SIM consists of three stages: (1) an atmospheric reanalysis (SAFRAN) over France, which forces (2) the three-layer ISBA land surface model, which then provides drainage and runoff inputs to (3) the MODCOU hydro-geological model. The drainage and runoff outputs from ISBA are validated by comparing the simulated river discharge from MODCOU with over 500 river-gauge observations over France and with a subset of stations with low-anthropogenic influence, over several years. This study makes use of the A-gs version of ISBA that allows for physiological processes. The atmospheric forcing for the ISBA-A-gs model underestimates direct shortwave and long-wave radiation by approximately 5 % averaged over France. The ISBA-A-gs model also substantially underestimates the grassland LAI compared with satellite retrievals during winter dormancy. These differences result in an underestimation (overestimation) of evapotranspiration (drainage and runoff). The excess runoff flowing into the rivers and aquifers contributes to an overestimation of the SIM river discharge. Two experiments attempted to resolve these problems: (i) a correction of the minimum LAI model parameter for grasslands and (ii) a bias-correction of the model radiative forcing. Two data assimilation experiments were also performed, which are designed to correct random errors in the initial conditions: (iii) the assimilation of LAI observations and (iv) the assimilation of SSM and LAI observations. The data assimilation for (iii) and (iv) was done with a simplified extended Kalman filter (SEKF), which uses finite differences in the observation operator Jacobians to relate the observations to the model variables. Experiments (i) and (ii) improved the median SIM Nash scores by about 9 % and 18 % respectively. Experiment (iii) reduced the LAI phase errors in ISBA-A-gs but had little impact on the discharge Nash efficiency of SIM. In contrast, experiment (iv) resulted in spurious increases in drainage and runoff, which degraded the median discharge Nash efficiency by about 7 %. The poor performance of the SEKF originates from the observation operator Jacobians. These Jacobians are dampened when the soil is saturated and when the vegetation is dormant, which leads to positive biases in drainage and/or runoff and to insufficient corrections during winter, respectively. Possible ways to improve the model are discussed, including a new multi-layer diffusion model and a more realistic response of photosynthesis to temperature in mountainous regions. The data assimilation should be advanced by accounting for model and forcing uncertainties.


2011 ◽  
Vol 12 (4) ◽  
pp. 508-530 ◽  
Author(s):  
Natacha B. Bernier ◽  
Stéphane Bélair ◽  
Bernard Bilodeau ◽  
Linying Tong

Abstract A high-resolution 2D near-surface and land surface model was developed to produce snow and temperature forecasts over the complex alpine region of the Vancouver 2010 Winter Olympic and Paralympic Games. The model is driven by downscaled operational outputs from the Meteorological Service of Canada’s regional and global forecast models. Downscaling is applied to correct forcings for elevation differences between the operational forecast models and the high-resolution surface model. The high-resolution near-surface and land surface model is then used to further refine the forecasts. The model was validated against temperature and snow depth observations. The largest improvements were found in regions where low-resolution (i.e., on the order of 10 km or more) operational models typically lack the spatial resolution to capture rapid elevation changes. The model was found to better reproduce the intermittent snow cover at low-lying stations and to reduce snow depth error by as much as 3 m at alpine stations.


2015 ◽  
Vol 9 (6) ◽  
pp. 6733-6790
Author(s):  
B. Decharme ◽  
E. Brun ◽  
A. Boone ◽  
C. Delire ◽  
P. Le Moigne ◽  
...  

Abstract. In this study we analysed how an improved representation of snowpack processes and soil properties in the multi-layer snow and soil schemes of the ISBA land surface model impacts the simulation of soil temperature profiles over North-Eurasian regions. For this purpose, we refine ISBA's snow layering algorithm and propose a parameterization of snow albedo and snow compaction/densification adapted from the detailed Crocus snowpack model. We also include a dependency on soil organic carbon content for ISBA's hydraulic and thermal soil properties. First, changes in the snowpack parameterization are evaluated against snow depth, snow water equivalent, surface albedo, and soil temperature at a 10 cm depth observed at the Col de Porte field site in the French Alps. Next, the new model version including all of the changes is used over Northern-Eurasia to evaluate the model's ability to simulate the snow depth, the soil temperature profile and the permafrost characteristics. The results confirm that an adequate simulation of snow layering and snow compaction/densification significantly impacts the snowpack characteristics and the soil temperature profile during winter, while the impact of the more accurate snow albedo computation is dominant during the spring. In summer, the accounting for the effect of soil organic carbon on hydraulic and thermal soil properties improves the simulation of the soil temperature profile. Finally, the results confirm that this last process strongly influences the simulation of the permafrost active layer thickness and its spatial distribution.


2013 ◽  
Vol 10 (8) ◽  
pp. 11093-11128 ◽  
Author(s):  
N. C. MacKellar ◽  
S. J. Dadson ◽  
M. New ◽  
P. Wolski

Abstract. Land surface models (LSMs) are advanced tools which can be used to estimate energy, water and biogeochemical exchanges at regional scales. The inclusion of a river flow routing module in an LSM allows for the simulation of river discharge from a catchment and offers an approach to evaluate the response of the system to variations in climate and land-use, which can provide useful information for regional water resource management. This study offers insight into some of the pragmatic considerations of applying an LSM over a regional domain in Southern Africa. The objectives are to identify key parameter sensitivities and investigate differences between two runoff production schemes in physically contrasted catchments. The Joint UK Land Environment Simulator (JULES) LSM was configured for a domain covering Southern Africa at a 0.5° resolution. The model was forced with meteorological input from the WATCH Forcing Data for the period 1981–2001 and sensitivity to various model configurations and parameter settings were tested. Both the PDM and TOPMODEL sub-grid scale runoff generation schemes were tested for parameter sensitivities, with the evaluation focussing on simulated river discharge in sub-catchments of the Orange, Okavango and Zambezi rivers. It was found that three catchments respond differently to the model configurations and there is no single runoff parameterization scheme or parameter values that yield optimal results across all catchments. The PDM scheme performs well in the upper Orange catchment, but poorly in the Okavango and Zambezi, whereas TOPMODEL grossly underestimates discharge in the upper Orange and shows marked improvement over PDM for the Okavango and Zambezi. A major shortcoming of PDM is that it does not realistically represent subsurface runoff in the deep, porous soils typical of the Okavango and Zambezi headwaters. The dry-season discharge in these catchments is therefore not replicated by PDM. TOPMODEL, however, simulates a more realistic seasonal cycle of subsurface runoff and hence improved dry-season flow.


2021 ◽  
Author(s):  
Ebony Lee ◽  
Seon Ki Park

<p>The Noah Land Surface Model (Noah LSM) estimates snow depth using snow water equivalent and snow density. The snow density is determined by snow compaction, snowmelt water storing, and density of fresh snowfall. The Noah LSM usually underestimates snow depth compared to the ground observations in Korea, which occurs from the beginning of snowfall. We performed an optimal estimation of parameters related to the density of fresh snowfall, using micro-genetic algorithm (μ-GA) that uses the evolution process concept through natural selection and mutation mechanism. Ground observations from 36 sites of the Korea Meteorological Administration, for the recent 10 years (May 2009 – April 2019), are used for offline forcing of the Noah LSM and evaluating the fitness function in μ-GA. Optimized parameters reduced the density of fresh snowfall, and improved the simulated snow depth. The root-mean-square error of snow depth decreased from 8.1 cm to 7.1 cm.</p>


2018 ◽  
Vol 22 (6) ◽  
pp. 3515-3532 ◽  
Author(s):  
Clement Albergel ◽  
Emanuel Dutra ◽  
Simon Munier ◽  
Jean-Christophe Calvet ◽  
Joaquin Munoz-Sabater ◽  
...  

Abstract. The European Centre for Medium-Range Weather Forecasts (ECMWF) recently released the first 7-year segment of its latest atmospheric reanalysis: ERA-5 over the period 2010–2016. ERA-5 has important changes relative to the former ERA-Interim atmospheric reanalysis including higher spatial and temporal resolutions as well as a more recent model and data assimilation system. ERA-5 is foreseen to replace ERA-Interim reanalysis and one of the main goals of this study is to assess whether ERA-5 can enhance the simulation performances with respect to ERA-Interim when it is used to force a land surface model (LSM). To that end, both ERA-5 and ERA-Interim are used to force the ISBA (Interactions between Soil, Biosphere, and Atmosphere) LSM fully coupled with the Total Runoff Integrating Pathways (TRIP) scheme adapted for the CNRM (Centre National de Recherches Météorologiques) continental hydrological system within the SURFEX (SURFace Externalisée) modelling platform of Météo-France. Simulations cover the 2010–2016 period at half a degree spatial resolution. The ERA-5 impact on ISBA LSM relative to ERA-Interim is evaluated using remote sensing and in situ observations covering a substantial part of the land surface storage and fluxes over the continental US domain. The remote sensing observations include (i) satellite-driven model estimates of land evapotranspiration, (ii) upscaled ground-based observations of gross primary production, (iii) satellite-derived estimates of surface soil moisture and (iv) satellite-derived estimates of leaf area index (LAI). The in situ observations cover (i) soil moisture, (ii) turbulent heat fluxes, (iii) river discharges and (iv) snow depth. ERA-5 leads to a consistent improvement over ERA-Interim as verified by the use of these eight independent observations of different land status and of the model simulations forced by ERA-5 when compared with ERA-Interim. This is particularly evident for the land surface variables linked to the terrestrial hydrological cycle, while variables linked to vegetation are less impacted. Results also indicate that while precipitation provides, to a large extent, improvements in surface fields (e.g. large improvement in the representation of river discharge and snow depth), the other atmospheric variables play an important role, contributing to the overall improvements. These results highlight the importance of enhanced meteorological forcing quality provided by the new ERA-5 reanalysis, which will pave the way for a new generation of land-surface developments and applications.


2020 ◽  
Vol 12 (4) ◽  
pp. 645 ◽  
Author(s):  
Sujay Kumar ◽  
David Mocko ◽  
Carrie Vuyovich ◽  
Christa Peters-Lidard

Surface albedo has a significant impact in determining the amount of available net radiation at the surface and the evolution of surface water and energy budget components. The snow accumulation and timing of melt, in particular, are directly impacted by the changes in land surface albedo. This study presents an evaluation of the impact of assimilating Moderate Resolution Imaging Spectroradiometer (MODIS)-based surface albedo estimates in the Noah multi-parameterization (Noah-MP) land surface model, over the continental US during the time period from 2000 to 2017. The evaluation of simulated snow depth and snow cover fields show that significant improvements from data assimilation (DA) are obtained over the High Plains and parts of the Rocky Mountains. Earlier snowmelt and reduced agreements with reference snow depth measurements, primarily over the Northeast US, are also observed due to albedo DA. Most improvements from assimilation are observed over locations with moderate vegetation and lower elevation. The aggregate impact on evapotranspiration and runoff from assimilation is found to be marginal. This study also evaluates the relative and joint utility of assimilating fractional snow cover and surface albedo measurements. Relative to surface albedo assimilation, fractional snow cover assimilation is found to provide smaller improvements in the simulated snow depth fields. The configuration that jointly assimilates surface albedo and fractional snow cover measurements is found to provide the most beneficial improvements compared to the univariate DA configurations for surface albedo or fractional snow cover. Overall, the study also points to the need for improving the albedo formulations in land surface models and the incorporation of observational uncertainties within albedo DA configurations.


2013 ◽  
Vol 14 (1) ◽  
pp. 203-219 ◽  
Author(s):  
Eric Brun ◽  
Vincent Vionnet ◽  
Aaron Boone ◽  
Bertrand Decharme ◽  
Yannick Peings ◽  
...  

Abstract The Crocus snowpack model within the Interactions between Soil–Biosphere–Atmosphere (ISBA) land surface model was run over northern Eurasia from 1979 to 1993, using forcing data extracted from hydrometeorological datasets and meteorological reanalyses. Simulated snow depth, snow water equivalent, and density over open fields were compared with local observations from over 1000 monitoring sites, available either once a day or three times per month. The best performance is obtained with European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-Interim). Provided blowing snow sublimation is taken into account, the simulations show a small bias and high correlations in terms of snow depth, snow water equivalent, and density. Local snow cover durations as well as the onset and vanishing dates of continuous snow cover are also well reproduced. A major result is that the overall performance of the simulations is very similar to the performance of existing gridded snow products, which, in contrast, assimilate local snow depth observations. Soil temperature at 20-cm depth is reasonably well simulated. The methodology developed in this study is an efficient way to evaluate different meteorological datasets, especially in terms of snow precipitation. It reveals that the temporal disaggregation of monthly precipitation in the hydrometeorological dataset from Princeton University significantly impacts the rain–snow partitioning, deteriorating the simulation of the onset of snow cover as well as snow depth throughout the cold season.


Sign in / Sign up

Export Citation Format

Share Document