Discharge and bathymetry estimations of rivers from altimetry and datasets by hybrid computational methods

2020 ◽  
Author(s):  
Jerome Monnier ◽  
Kevin Larnier ◽  
Pierre-André Garambois
Keyword(s):  
Data Assimilation ◽  
River Discharge ◽  
Hydrological Cycle ◽  
Low Complexity ◽  
Space Mission ◽  
Free Surface Flows ◽  
Learning Approaches ◽  
Flow Models ◽  
Hybrid Data Assimilation ◽  
River Models

<p>We present the Hierarchical Variational Discharge Inference (HiVDI) algorithm [1,2] and its capabilities to estimate the discharge and bathymetry of rivers from altimetry measurement, more particularly from the forthcoming SWOT space mission. The last version algorithm is based on hierarchical flow models and hybrid computational approaches : 1) a dedicated satellite-scale low-complexity model relating the discharge Q(x,t), the bathymetry b(x) and the friction parameter K [2]; 2) an advanced Variational Data Assimilation (VDA) formulation based on a relatively complete physics (Saint-Venant’s equations) [2,4] ; 3) deep neural networks based estimations obtained from recently enriched databases [1]. The resulting algorithm turns out to be robust and relatively accurate. Passed the assimilation of a hydrological cycle (~ 1 year variations, considered as a “learning period) the identified parameters (b(x), K) are identified; next given newly acquired satellite measurements, the low complexity model enables to estimate Q(x,t) in real-time [1,2].</p><p>Numerical results on numerous river datasets are analyzed in detail including for relatively complex flows and multi-satellite datasets [1,2,3].</p><p>References</p><p>[1] K. Larnier, J. Monnier. "Hybrid data assimilation - deep learning approaches to estimate rivers discharges from altimetry". Submitted.</p><p>[2] K. Larnier, J. Monnier, P.-A. Garambois, J. Verley. "River discharge and bathymetry estimations from SWOT altimetry measurements". Revised (nov. 2019).</p><p>[3] P.-A. Garambois, K. Larnier, J. Monnier, P. Finaud-Guyot, J. Verley, A. Montazem, S. Calmant. "Variational inference of effective channel and ungauged anabranching river discharge from multi-satellite water heights of different spatial sparsity". J. of Hydrology 2019.</p><p>[4] P. Brisset, J. Monnier, P.-A. Garambois, H. Roux. "On the assimilation of altimetry data in 1D Saint-Venant river models". Adv. Water Ress. 2018. </p><p>[5] "DassFlow: Data Assimilation for Free Surface Flows", open-source computational software. INSA - IMT, CNRS, CNES, CS group. http://www.math.univ-toulouse.fr/DassFlow</p>

scholarly journals A GSI-Based Coupled EnSRF–En3DVar Hybrid Data Assimilation System for the Operational Rapid Refresh Model: Tests at a Reduced Resolution

2014 ◽  
Vol 142 (10) ◽  
pp. 3756-3780 ◽  
Author(s):  
Yujie Pan ◽  
Kefeng Zhu ◽  
Ming Xue ◽  
Xuguang Wang ◽  
Ming Hu ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Hybrid System ◽  
Control Variable ◽  
Three Dimensional ◽  
Horizontal Resolution ◽  
Precipitation Forecast ◽  
Background Error ◽  
Hybrid Data Assimilation ◽  
Noticeable Improvement ◽  
Forecasting System

Abstract A coupled ensemble square root filter–three-dimensional ensemble-variational hybrid (EnSRF–En3DVar) data assimilation (DA) system is developed for the operational Rapid Refresh (RAP) forecasting system. The En3DVar hybrid system employs the extended control variable method, and is built on the NCEP operational gridpoint statistical interpolation (GSI) three-dimensional variational data assimilation (3DVar) framework. It is coupled with an EnSRF system for RAP, which provides ensemble perturbations. Recursive filters (RF) are used to localize ensemble covariance in both horizontal and vertical within the En3DVar. The coupled En3DVar hybrid system is evaluated with 3-h cycles over a 9-day period with active convection. All conventional observations used by operational RAP are included. The En3DVar hybrid system is run at ⅓ of the operational RAP horizontal resolution or about 40-km grid spacing, and its performance is compared to parallel GSI 3DVar and EnSRF runs using the same datasets and resolution. Short-term forecasts initialized from the 3-hourly analyses are verified against sounding and surface observations. When using equally weighted static and ensemble background error covariances and 40 ensemble members, the En3DVar hybrid system outperforms the corresponding GSI 3DVar and EnSRF. When the recursive filter coefficients are tuned to achieve a similar height-dependent localization as in the EnSRF, the En3DVar results using pure ensemble covariance are close to EnSRF. Two-way coupling between EnSRF and En3DVar did not produce noticeable improvement over one-way coupling. Downscaled precipitation forecast skill on the 13-km RAP grid from the En3DVar hybrid is better than those from GSI 3DVar analyses.

Multi-resolution Hybrid Data Assimilation Core on a Cubed-sphere Grid (HybDA)

2018 ◽  
Vol 54 (S1) ◽  
pp. 337-350 ◽  
Author(s):  
Hyo-Jong Song ◽  
Ji-Hyun Ha ◽  
In-Hyuk Kwon ◽  
Junghan Kim ◽  
Jihye Kwun
Keyword(s):  
Data Assimilation ◽  
Hybrid Data Assimilation ◽  
Hybrid Data ◽  
Cubed Sphere

Hydrological impact of the new ECMWF multi-layer snow scheme

2021 ◽  
Author(s):  
Gabriele Arduini ◽  
Ervin Zsoter ◽  
Hannah Cloke ◽  
Elisabeth Stephens ◽  
Christel Prudhomme
Keyword(s):  
Land Surface ◽  
River Discharge ◽  
Hydrological Cycle ◽  
Land Surface Model ◽  
Single Layer ◽  
Observation Time ◽  
Surface Model ◽  
Earth System ◽  
Hydrological Impact ◽  
Energy And Momentum

<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>

Abstract TP458: High Accuracy of Predictive Models for SAH Using Different Machine Learning Approaches

Stroke ◽  
2020 ◽  
Vol 51 (Suppl_1) ◽  
Author(s):  
Paul Litvak ◽  
Jeevan Medikonda ◽  
Girish Menon ◽  
Pitchaiah Mandava
Keyword(s):  
Machine Learning ◽  
Predictive Models ◽  
Prediction Accuracy ◽  
Support Vector ◽  
World Federation ◽  
Learning Approaches ◽  
Flow Models ◽  
Multi Stage ◽  
Stage 1 ◽  
Categorical Scale

Background: Patients suffering from subarachnoid hemorrhage (SAH) have poor long-term outcomes. There are predictive models for ischemic and hemorrhagic stroke. However, there is paucity of models for SAH. Machine learning concepts were applied to build multi-stage Neural Networks (NN), Support Vector Machines (SVM) and Keras/Tensor Flow models to predict SAH outcomes. Methods: A database of ~800 aneurysmal SAH patients from Kasturba Medical College was utilized. Baseline variables of World Federation of Neurosurgeons 5-point scale (WFNS 1-5), age, gender, and presence/absence of hypertension and diabetes were considered in Stage 1. Stage 2 included all Stage 1 variables along with presence/absence of radiologic signs vasospasm and ischemia. Stage 3 includes earlier 2 stages and discharge Glasgow Outcome Scale (GOS 1-5). GOS at 3 months was predicted using 2-layer NN/SVM/Keras-TensorFlow models on the five point categorical scale as well as dichotomized to dead/alive and favorable (GOS 4-5) or unfavorable (GOS 1-3). Prediction accuracy of models was compared to the recorded GOS. Results: Prediction accuracy shown as percentages (See Table) for all three stages was similar for SVM, NN and Keras/TensorFlow models. Accuracy was remarkably higher with dichotomization compared to the complete five point GOS categorical scale. Conclusions: SVM, NN, and Keras-TensorFlow based machine learning models can be used to predict SAH outcomes to a high degree of accuracy. These powerful predictive models can be used to prognosticate and select patients into trials.

A Data Assimilation Framework for Generating Space‐Time Continuous Daily SWOT River Discharge Data Products

2020 ◽  
Vol 56 (6) ◽  
Author(s):  
Dongyue Li ◽  
Konstantinos M. Andreadis ◽  
Steven A. Margulis ◽  
Dennis P. Lettenmaier
Keyword(s):  
Data Assimilation ◽  
River Discharge ◽  
Space Time ◽  
Discharge Data ◽  
Data Products

scholarly journals The effect of satellite-derived surface soil moisture and leaf area index land data assimilation on streamflow simulations over France

2017 ◽  
Vol 21 (4) ◽  
pp. 2015-2033 ◽  
Author(s):  
David Fairbairn ◽  
Alina Lavinia Barbu ◽  
Adrien Napoly ◽  
Clément Albergel ◽  
Jean-François Mahfouf ◽  
...  
Keyword(s):  
Soil Moisture ◽  
Data Assimilation ◽  
Leaf Area Index ◽  
Land Surface ◽  
River Discharge ◽  
Surface Soil ◽  
Land Surface Model ◽  
Surface Model ◽  
Surface Soil Moisture ◽  
Area Index

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.

E4DVar: Coupling an Ensemble Kalman Filter with Four-Dimensional Variational Data Assimilation in a Limited-Area Weather Prediction Model

2012 ◽  
Vol 140 (2) ◽  
pp. 587-600 ◽  
Author(s):  
Meng Zhang ◽  
Fuqing Zhang
Keyword(s):  
Kalman Filter ◽  
Data Assimilation ◽  
Prediction Model ◽  
Ensemble Kalman Filter ◽  
Weather Prediction ◽  
Coupled System ◽  
Limited Area ◽  
Ensemble Forecasts ◽  
Background Error ◽  
Hybrid Data Assimilation

A hybrid data assimilation approach that couples the ensemble Kalman filter (EnKF) and four-dimensional variational (4DVar) methods is implemented for the first time in a limited-area weather prediction model. In this coupled system, denoted E4DVar, the EnKF and 4DVar systems run in parallel while feeding into each other. The multivariate, flow-dependent background error covariance estimated from the EnKF ensemble is used in the 4DVar minimization and the ensemble mean in the EnKF analysis is replaced by the 4DVar analysis, while updating the analysis perturbations for the next cycle of ensemble forecasts with the EnKF. Therefore, the E4DVar can obtain flow-dependent information from both the explicit covariance matrix derived from ensemble forecasts, as well as implicitly from the 4DVar trajectory. The performance of an E4DVar system is compared with the uncoupled 4DVar and EnKF for a limited-area model by assimilating various conventional observations over the contiguous United States for June 2003. After verifying the forecasts from each analysis against standard sounding observations, it is found that the E4DVar substantially outperforms both the EnKF and 4DVar during this active summer month, which featured several episodes of severe convective weather. On average, the forecasts produced from E4DVar analyses have considerably smaller errors than both of the stand-alone EnKF and 4DVar systems for forecast lead times up to 60 h.

Do high resolution GCMs overestimate precipitation over land?

2020 ◽  
Author(s):  
Omar Müller ◽  
Pier Luigi Vidale ◽  
Patrick McGuire ◽  
Benoît Vannière ◽  
Reinhard Schiemann ◽  
...  
Keyword(s):  
High Resolution ◽  
Land Surface ◽  
River Discharge ◽  
Hydrological Cycle ◽  
Transport Model ◽  
Meteorological Data ◽  
Grid Point ◽  
Climate Simulation ◽  
Surface Model ◽  
Orographic Precipitation

<p>Previous studies showed that high resolution GCMs overestimate land precipitation when compared against gridded observations or reanalysis (Demory et al. 2014, Vannière et al. 2019). In particular, grid point models (eg. HadGEM3) show a significant increase of precipitation on regions dominated by complex orography, where the scarcity of gauge stations increase the uncertainty of gridded observations. The goal of this work is to assess the effect of such differences in precipitation on river discharge, considering it as an integrator of the water balance at catchment scale. A set of JULES and CLM simulations have been conducted turning rivers on with Total Runoff Integrating Pathways (TRIP) and the River Transport Model (RTM) respectively. The simulations form three ensembles for each land surface model (LSM) which main difference is given by the forcing dataset. The forcings are WFDEI (reanalysis), LR (~1° resolution in meteorological data from GCMs) and HR (~0.25° resolution in meteorological data from GCMs). These ensembles are evaluated in a set of 280 catchments distributed around the world.</p><p>In terms of correlation between simulated and observed river discharge observations, the results show that LSMs forced by reanalysis have higher performance than LSMs forced by GCMs as expected. In terms of biases, the river discharge is underestimated in eight out of eleven major basins when LSMs are forced by reanalysis. On those basins, the extra precipitation estimated by GCMs help to simulate an amount of river discharge closer to observations (Eg. Yenisey and Lena). Moreover, 37 small basins with a strong component of orographic precipitation over the Andes, the Rocky Mountains, the Alps and in the Maritime Continent were evaluated. In most cases HR offers notably better results than LR and WFDEI, suggesting that high resolution models produce orographic precipitation in the correct place and time.</p><p>In future works offline TRIP simulations will be carried out directly forced by runoff and subsurface runoff from GCMs. It will allow to discard errors in evapotranspiration produced by JULES or CLM when they are used to simulate river discharge. This work is part of the European Process-based climate sIMulation: AdVances in high resolution modelling and European climate Risk Assessment (PRIMAVERA) Project. PRIMAVERA is a collaboration between 19 funded by the European Union’s Horizon 2020 Research & Innovation Programme.</p><p>Demory, M. E., Vidale, P. L., Roberts, M. J., Berrisford, P., Strachan, J., Schiemann, R., & Mizielinski, M. S. (2014). The role of horizontal resolution in simulating drivers of the global hydrological cycle. CLIM DYNAM, 42(7-8), 2201-2225.</p><p>Vannière, B., Demory, M. E., Vidale, P. L., Schiemann, R., Roberts, M. J., Roberts, C. D., ... & Senan, R. (2018). Multi-model evaluation of the sensitivity of the global energy budget and hydrological cycle to resolution. CLIM DYNAM, 1-30.</p>

Sign in / Sign up

Export Citation Format

Share Document