Search NASA⌕ Search

SEARCH · Search NASA

Results for “streamflow simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Can artificial intelligence and data-driven machine learning models match or even replace process-driven hydrologic models for streamflow simulation?: A case study of four watersheds with different hydro-climatic regions across the CONUS

With recent developments in computational techniques, Data-driven Machine Learning Models (DMLs) have shown great potential in simulating streamflow and capturing the rainfall-runoff relationship in given watersheds, which are traditionally fulfilled by Process-based Hydrologic Models (PHMs). There are debates on whether the DMLs can outperform and possibly replace the classical PHMs for streamflow simulation and river forecasting, but no clear conclusions have been made. This study aims to investigate whether the newer DMLs have any potential in further improving the simulation accuracy of classical PHMs, and vice versa. To do this, we compared a few popular PHMs and DMLs over four watersheds across the Continental US (CONUS) that are associated with different input, climate, and regional conditions. A total of five hydrologic models were chosen, including (1) two classical lumped models, i.e., the Sacramento Soil Moisture Accounting (SAC-SMA) and Xinanjiang (XAJ); (2) one modern distributed model, termed Coupled Routing and Excess Storage (CREST); (3) and two DMLs including an Artificial Neural Networks (ANN) and a deep learning model, termed Long Short Term Memory (LSTM). Our results demonstrated that the DMLs still significantly biased when using the baseline input scenario with the PHMs. However, the DMLs fed with delayed input scenarios had great potential and can reach high simulation accuracy. The DMLs, especially the ANN, outperformed other employed models under the rainfall-runoff relationship in which rainfall dominantly drives. Furthermore, the DMLs also showed better performance in the high-flow regime, while the PHMs had a better performance for the low-flow regime, implying both PHMs and DMLs have their own merits and are worthy of joint development. In general, our study indicated a great potential of using DMLs to simulate streamflow, but further studies are still needed to verify the transferability and scalability of DMLs in large-scale experiments, such as the Distributed Model Intercomparison Projects 1&2 conducted by National Weather Services but to compare modern DMLs and PHMs.

58 GEOSCIENCES↗

A Mass Conservation Relaxed (MCR) LSTM Model for Streamflow Simulation Across CONUS

The recent development of the physics-aware Mass-Conserving Long Short-Term Memory network (MC-LSTM) provides an alternative to other data-driven Deep Learning (DL) models in hydrology. Mass-Conserving Long Short-Term Memory incorporates mass conservation directly into the LSTM architecture. Despite the theoretical advancements, studies have reported a surprisingly limited performance of the MC-LSTM in streamflow simulation. We hypothesize that such a limitation is due to the unrealistic mass conservation scheme in MC-LSTM, which overlooks unobserved incoming water fluxes beyond precipitation. As an attempt to verify this hypothesis, we propose a Mass Conservation Relaxed LSTM (MCR-LSTM), which incorporates a bi-directional mass relaxation (MR) component to account for potential incoming water fluxes beyond precipitation. We train and test the proposed MCR-LSTM model across 531 watersheds in the contiguous United States (CONUS) against three baseline models: the Sacramento Soil Moisture Accounting, LSTM, and MC-LSTM. Our results show that MCR-LSTM outperforms MC-LSTM despite its underperformance compared to LSTM. Specifically, MCR-LSTM's advantage over MC-LSTM is mainly seen in the Plains and Western U.S., where the newly incorporated MR component better simulates water loss and suggests the likely existence of additional incoming water fluxes beyond precipitation, respectively. The novelty and contribution of this study are twofold: firstly, it introduces an alternative physics-aware DL tool (i.e., MCR-LSTM) in hydrology with higher accuracy in specific regions compared to MC-LSTM. Secondly, it provides a diagnosis of regions where strict, precipitation-based mass conservation constraints may be unrealistic in streamflow simulation.

deep learning↗

Comparing quantile regression forest and mixture density long short-term memory models for probabilistic post-processing of satellite precipitation-driven streamflow simulations

Abstract. Deep learning (DL) and machine learning (ML) are widely used in hydrological modelling, which plays a critical role in improving the accuracy of hydrological predictions. However, the trade-off between model performance and computational cost has always been a challenge for hydrologists when selecting a suitable model, particularly for probabilistic post-processing with large ensemble members. This study aims to systematically compare the quantile regression forest (QRF) model and countable mixtures of asymmetric Laplacians long short-term memory (CMAL-LSTM) model as hydrological probabilistic post-processors. Specifically, we evaluate their ability in dealing with biased streamflow simulations driven by three satellite precipitation products across 522 nested sub-basins of the Yalong River basin in China. Model performance is comprehensively assessed using a series of scoring metrics from both probabilistic and deterministic perspectives. Our results show that the QRF model and the CMAL-LSTM model are comparable in terms of probabilistic prediction, and their performances are closely related to the flow accumulation area (FAA) of the sub-basin. The QRF model outperforms the CMAL-LSTM model in most sub-basins with smaller FAA, while the CMAL-LSTM model has an undebatable advantage in sub-basins with FAA larger than 60 000 km2 in the Yalong River basin. In terms of deterministic predictions, the CMAL-LSTM model is preferred, especially when the raw streamflow is poorly simulated and used as input. However, setting aside the differences in model performance, the QRF model with 100-member quantiles demonstrates a noteworthy advantage by exhibiting a 50 % reduction in computation time compared to the CMAL-LSTM model with the same ensemble members in all experiments. As a result, this study provides insights into model selection in hydrological post-processing and the trade-offs between model performance and computational efficiency. The findings highlight the importance of considering the specific application scenario, such as the catchment size and the required accuracy level, when selecting a suitable model for hydrological post-processing.

Geology↗

Machine learning assisted hybrid models can improve streamflow simulation in diverse catchments across the conterminous US

Incomplete representations of physical processes often lead to structural errors in process-based (PB) hydrologic models. Machine learning (ML) algorithms can reduce streamflow modeling errors but do not enforce physical consistency. As a result, ML algorithms may be unreliable if used to provide future hydroclimate projections where climates and land use patterns are outside the range of training data. Here we test hybrid models built by integrating PB model outputs with a ML algorithm known as Long Short-Term Memory (LSTM) network on their ability to simulate streamflow in 531 catchments representing diverse conditions across the Conterminous United States. Model performance of hybrid models as measured by Nash-Sutcliffe efficiency (NSE) improved relative to standalone PB and LSTM models. More importantly, hybrid models provide highest improvement in catchments where PB models fail completely (i.e., NSE < 0). However, all models performed poorly in catchments with extended low flow periods, suggesting need for additional research.

54 ENVIRONMENTAL SCIENCES↗

Streamflow simulation in data-scarce basins using Bayesian and physics-informed machine learning models

Hydrologic predictions at rural watersheds are important but also challenging due to data shortage. Long short-term memory (LSTM) networks are a promising machine learning approach and have demonstrated good performance in streamflow predictions. However, due to its data-hungry nature, most LSTM applications focus on well-monitored catchments with abundant and high-quality observations. In this work, we investigate predictive capabilities of LSTM in poorly monitored watersheds with short observation records. To address three main challenges of LSTM applications in data-scarce locations, i.e., overfitting, uncertainty quantification (UQ), and out-of-distribution prediction, we evaluate different regularization techniques to prevent overfitting, apply a Bayesian LSTM for UQ, and introduce a physics-informed hybrid LSTM to enhance out-of-distribution prediction. Through case studies in two diverse sets of catchments with and without snow influence, we demonstrate that 1) when hydrologic variability in the prediction period is similar to the calibration period, LSTM models can reasonably predict daily streamflow with Nash–Sutcliffe efficiency above 0.8, even with only 2 years of calibration data; 2) when the hydrologic variability in the prediction and calibration periods is dramatically different, LSTM alone does not predict well, but the hybrid model can improve the out-of-distribution prediction with acceptable generalization accuracy; 3) L2 norm penalty and dropout can mitigate overfitting, and Bayesian and hybrid LSTM have no overfitting; and 4) Bayesian LSTM provides useful uncertainty information to improve prediction understanding and credibility. In conclusion, these insights have vital implications for streamflow simulation in watersheds where data quality and availability are a critical issue.

54 ENVIRONMENTAL SCIENCES↗

Hydrologic Model Data for the East Fork Poplar Creek Watershed Simulated with the Advanced Terrestrial Simulator (ATS): Streamflow and Network Expansion–Contraction Dynamics

This dataset supports hydrologic modeling and stream network expansion–contraction analysis for the East Fork Poplar Creek (EFPC) Watershed in Tennessee. It includes a Jupyter notebook for model setup, model configuration files, simulation outputs, and derived products used to evaluate model performance and investigate stream dynamics under varying hydrologic conditions. The dataset was generated using the Watershed Workflow Python package and the Advanced Terrestrial Simulator (ATS), enabling integrated surface–subsurface hydrologic simulations using a stream-aligned mesh. Outputs include high-resolution time series of streamflow, active network length, water table depth, and related hydrologic variables. Also included are spatially explicit stream persistency indices and classifications of reaches as perennial or non-perennial. These data facilitate reproducibility and support further research on stream intermittency and variability in network extent.The model data archive is organized in following directories:1) model_setup_inputsContains the Watershed Workflow Jupyter notebooks (accessed through any open source code editor), selected input datasets, and resulting ATS input files, including XML files (access through any open source code editor), computational mesh (.exo files can be viewed using Paraview), and meteorological forcing files (.h5 files can be accessed through h5py python package and HDFView open source software). 2) model_outputsIncludes ATS simulation outputs relevant to this study. Time series of spatially integrated or averaged variables (e.g., streamflow, water table depth) are provided as CSV files. Select spatial fields (e.g., ponded depth and water table depth) are saved as pickled Python objects to reduce file size, and can be accessed through pickle package in Python. Key geometry objects from Watershed Workflow—such as the surface mesh and river tree—are also included to support analysis of streamflow persistency and expansion–contraction dynamics. These files can also be accessed through Watershed Workflow Python package.3) model_evaluationProvides observed streamflow time series and field survey-based flow regime classifications used to evaluate model performance. Jupyter notebooks for processing ATS outputs and comparing model predictions with observations to build confidence in the model prior to scientific analysis are also included.4) Q_L_relationshipsContains workflows for generating time series of discharge, active network length, and related hydrologic variables used in the stream network expansion–contraction analysis. Includes routines for delineating baseflow-dominated periods. For each catchment, notebooks and processed data (as pickled DataFrames accessed through Pandas Python package) are provided. 5) figure_scriptsProvides the Jupyter notebooks used to generate the figures presented in the paper.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of distributed process-based hydrologic model performance using only a priori information to define model inputs

Fully distributed, integrated surface–subsurface hydrological models (ISSHMs) have seen renewed interest due to availability of better software, high performance computing facilities, and high-resolution, spatially extensive data products. ISSHMs are valuable as tools for advancing system understanding as they can resolve multiple processes defined on the plot scale including three-dimensional interaction of surface water and groundwater. Here, we evaluated the performance of an ISSHM, the Advanced Terrestrial Simulator (ATS), on seven diverse catchments across the continental US using widely available data products to define model inputs without calibration. We compare the ATS-simulated streamflow and evapotranspiration with gauge observations and MODIS-derived evapotranspiration, respectively. Using the Kling-Gupta Efficiency (KGE) as metric, ATS with default data products performed reasonably well at 6 of 7 catchments for streamflow. However, in one of those 6 catchments ATS had poor performance on baseflow and ATS’s overall performance was thus judged to be inadequate despite the acceptable KGE. ATS performance for evapotranspiration was good in all 7 catchments using default data products. In the two catchments where ATS streamflow performance using default data products was not acceptable, the performance was significantly improved by using local information on subsurface properties below the soil. We also compare the model-simulated streamflow and evapotranspiration with the Sacramento soil moisture accounting (SAC-SMA) model, a semi-distributed model that was calibrated on a catchment-by-catchment basis. Uncalibrated ATS performance is comparable to the calibrated SAC-SMA model in terms of streamflow while ATS performance is similar to or better (much better in certain catchments) in reproducing MODIS-derived evapotranspiration. Reasonably good performance of ATS without catchment-specific calibration provides new confidence in the ISSHM class of models and community data products as tools for advancing understanding of watershed function in a changing environment.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the effects of urbanization on floods in a changing environment to promote water security — A case study of two adjacent basins in Texas

The increased occurrence of flood events resulting from urbanization and global climate change is a great threat to water security. To systematically evaluate the impacts of urbanization on floods, here we applied a paired catchments approach to two adjacent river basins in south-central Texas — the San Antonio River Basin (SARB), with fast urbanization; and the Guadalupe River Basin (GRB), with little land cover change. A physics-based distributed hydrological model — the Distributed Hydrology Soil Vegetation Model, embedded with a multi-purpose reservoir module (DHSVM-Res) — was used to simulate streamflow and reservoir storage. The simulations were conducted under different land cover scenarios, including a newly developed continuous land cover series (CLCS). Holistic analyses were then conducted for the paired basins using three methods: analyzing the selected flood events, detecting change points (CP) of monthly floods, and testing the elasticity of long-term flood regimes. The results suggest that: (1) urbanization may reduce lag time and elevate flood peaks significantly by 3–30% in our study area; (2) when there is little land cover change, changing climate is the major driver of variations in the monthly maximum streamflow (MMS); (3) fast urbanization can amplify streamflow variability, increase MMS significantly, and thus alter the timing of CP; and (4) the mean MMS of observed streamflow in the SARB has increased by as much as 75.7% from the pre-CP to post-CP periods. This comprehensive study fills in a gap in our current understanding of the isolated impacts of urbanization on flooding and is expected to support future explorations of anthropogenic influences on floods.

54 ENVIRONMENTAL SCIENCES↗

Agro-IBIS/THMB Modeled Yield, Nitrogen Loss, and Profitability Maps for the Raccoon River Basin, Iowa, USA

This data was used to assess the potential for improving water quality in this watershed when integrating miscanthus, a bioenergy perennial grass, on to areas of low profit under current corn-soybean management and areas of high nitrogen leaching under both current and future climate conditions. This dataset includes all observations (yield, streamflow, nitrogen loads) and model simulation data used from Dr. Kelsie Ferin’s (former PhD student with Dr. Andy VanLoocke, ISU Dept. of Agronomy) dissertation and publication in Global Change Biology – Bioenergy (DOI: 10.1111/gcbb.13078). Modeled output included in this dataset includes corn, soybean, and miscanthus yields, nitrogen leaching rates, net mineralization rates, land use cover fractions with and without the inclusion of miscanthus, profitability maps, additional soil property maps, and simulated streamflow and nitrogen loads near the outlet of the Raccoon River Basin (Van Meter, Iowa) for both historical and future climate conditions. This dataset includes a mixture of netCDF, .tiff, .csv, and .xlsx files. See README.pdf for more information.

Ferin, Kelsie↗

Near-Real-Time Statistical Analysis and Visualization of Streamflow from a Deep-Learning Rainfall-Runoff Model

Near-real-time (NRT) streamflow data are critical importance for timely water resources management. Here, we developed an open-source tool, FlowStats, for NRT streamflow analysis and visualization in Germany, based on NRT meteorological data from the German Weather Service and simulated streamflow from a long short-term memory neural network (LSTM). The LSTM model achieved very good overall performance, median NSE of 0.80 for the test period across 1,479 catchments. FlowStats provides options for deriving various streamflow statistics, from normal and abnormal streamflow detection to drought and flood analyses. An example analysis from FlowStats revealed widespread below-normal to extreme low-flow conditions across Germany from March to May 2025, which weakened from June to September 2025. Drought analysis for September 2025 highlighted severe to extreme drought conditions in northwestern Germany, while flood classifications indicated that high-flow events occurred in southwestern Germany. FlowStats can be used for various hydrological assessments to support water resources management.

Hydrological modeling↗

Uncertainties in Simulating Flooding During Hurricane Harvey Using 2D Shallow Water Equations

Abstract Flooding is one of the most impactful weather‐related natural hazards. Numerical models that solve the two dimensional (2D) shallow water equations (SWE) represent the first‐principles approach to simulate all types of spatial flooding, such as pluvial, fluvial, and coastal flooding, and their compound dynamics. High spatial resolution (e.g., () m) is needed in 2D SWE simulations to capture flood dynamics accurately, resulting in formidable computational challenges. Thus, relatively coarser spatial resolutions are used for large‐scale simulations of flooding, which introduce uncertainties in the results. It is unclear how the uncertainty associated with the model resolution compares to the uncertainties in precipitation data sets and assumptions regarding boundary conditions when channelized flows interact with other water bodies. In this study, we compare these three sources of uncertainties in 2D SWE simulations for the 2017 Houston flooding event. Our results show that precipitation uncertainty and mesh resolution have more significant impacts on the simulated streamflow and inundation dynamics than the choice of the downstream boundary condition at the watershed outlet. We point out the viability to confine the uncertainty of coarsening mesh resolution by using the variable resolution mesh (VRM) which refines critical topographic features with far fewer grid cells. Specifically, in simulations with VRM, the simulated inundation depths over the refined region are comparable to that use the finest uniform mesh. This study contributes to understanding the challenges and pathways for applying 2D SWE models to improve the realism of flood simulations over large scales.

54 ENVIRONMENTAL SCIENCES↗

Investigating Mountain Watershed Headwater‐To‐Groundwater Connections, Water Sources, and Storage Selection Behavior With Dynamic‐Flux Particle Tracking

Abstract Climate change will impact mountain watershed streamflow both directly—with changing precipitation amounts and variability—and indirectly—through temperature shifts altering snowpack, melt, and evapotranspiration. To understand how these complex processes will affect ecosystem functioning and water resources, we need tools to distinguish connections between water sources (rain/snowmelt), groundwater storage, and exit fluxes (streamflow/evapotranspiration), and to determine how these connections change seasonally and as climate shifts. Here, we develop novel watershed‐scale approaches to understand water source, storage, and exit flux connections using a dynamic‐flux particle tracking model (EcoSLIM) applied in California's Cosumnes Watershed, which connects the Sierra Nevada and Central Valley. This work develops new visualizations and applications to provide mechanistic understanding that underpins the interpretation of isotopic field data at watershed scales to distinguish sources, flow paths, residence times, and storage selection. In our simulations, streamflow comes primarily from snow‐derived water while evapotranspiration generally comes from rain. Most streamflow starts above 1,000 m while evapotranspiration is sourced relatively evenly across the watershed and is generally younger than streamflow. Modeled streamflow consists primarily of water sourced from precipitation in the previous 5 years but before the current water year, while ET consists primarily of water from precipitation in the current water year. ET, and to a lesser extent streamflow, are both younger than water in groundwater storage. However, snowmelt‐derived streamflow preferentially discharges older water from snow‐derived storage. Dynamic‐flux particle tracking and new approaches presented here enable novel model‐tracer comparisons in large‐scale watersheds to better understand watershed behavior in a changing climate.

54 ENVIRONMENTAL SCIENCES↗

An Improved Zhang's Dynamic Water Balance Model Using Budyko-Based Snow Representation for Better Streamflow Predictions

Understanding the water balance of a catchment in relation to its regional climate forcings and catchment characteristics is critical for predicting current and future water resources amid changing climate and land cover. This study intends to improve Zhang's monthly water balance model (a physics-based conceptual hydrologic model) that reflects the physical partitioning process of the hydrological cycle at the basin level based on regional climate and catchment characteristics. The existing model does not include snow process and has confronted evident limitations in snow-affected areas, which is a critical aspect since snowmelt water has been a significant source of water resources for many regions, especially in the temperate and frigid zones. We introduce a snow module based on surface energy balance and Budyko-limits on melting and combine it with the existing water balance equations. Moreover, monthly parameterization is applied to the model to better explain the time-varying hydrological characteristics of a catchment. The proposed model involves five different monthly parameters, which determine the physical partitioning process of the hydrological cycle, and they are regionally calibrated and validated under Budyko-type constraints. The model is applied to 1,210 basins across the continental United States (CONUS), and the simulated streamflow is compared to the observed data. The proposed model significantly outperformed the original model, improving the median NSE by 31% (from 0.51 to 0.67) and increasing the number of catchments with an acceptable NSE by 58%. The spatial variability of the basin characteristics across the CONUS is also investigated based on the calibrated parameters.

54 ENVIRONMENTAL SCIENCES↗

Enhanced flood synchrony and downstream severity in the Delaware River under rising temperatures

River floods threaten life and economic stability, with risks increasing globally, especially in densely populated coastal areas. In mountainous coastal watersheds like the Delaware River Basin, rising temperature is projected to reduce snowpack, reshaping upstream–downstream flood dynamics. However, the impact on flood synchronization between upland tributaries and estuarine mainstems remains poorly understood. Using multidecadal streamflow simulations from a high-resolution hydrological model, we find significant increases in the frequency and magnitude of synchronized floods under future warming scenarios, particularly severe floods (above the 75 th percentile). Under higher warming scenarios, snowpack and rain-on-snow floods in headwater subbasins nearly vanish. Surprisingly, this regime shift amplifies estuarine flooding by enhancing synchronization between historically snow-dominated subbasins and their downstream counterparts. Despite intensified synchronization and flood magnitude, cold-season flood risk declines due to fewer rain-on-snow events. Conversely, summer floods grow larger, more frequent, and synchronized, driving a seasonal reconfiguration of flood risk that should be considered in future mitigation strategies.

59 BASIC BIOLOGICAL SCIENCES↗

Thermal extremes in regulated river systems under climate change: an application to the southeastern U.S. rivers

High river temperatures, or “thermal extremes”, can cause fish mortality and thermoelectric powerplant derating. Under climate change, projected higher air temperature and stronger surface energy fluxes will lead to increased water temperatures, exacerbating thermal extremes. However, cold hypolimnetic releases from thermally stratified reservoirs can depress tailwater temperatures and therefore alleviate thermal extremes. Thermal extremes are more harmful when they coincide with low flows, which we refer to as “hydrologic hot-dry events”. To assess multisectoral impacts of climate change over large regions, we evaluate thermal events according to three impact attributes: duration (D), intensity (I), and severity (S). We apply an established model framework to simulate streamflow and stream temperature over the southeastern US regulated river system. We quantify climate change impacts (by the 2080s under RCP8.5) by comparing historical and future periods and quantify regulation impacts by comparing unregulated and regulated scenarios. We find that climate change will exacerbate thermal extremes (all three metrics) in both unregulated and regulated scenarios, albeit less in the regulated scenario. Thermal mitigation from reservoir regulation will be stronger under climate change, decreasing the three metrics compared to the unregulated case. Even so, thermal extremes in the regulated scenario will still be more severe under climate change, and only 12.2%, 19.7%, and 26.0% of D, I, and S can be mitigated by reservoirs. Despite stronger reservoir stratification, the number of regulated river segments that experience simultaneous high temperature and low flow events (hydrologic hot-dry events) will increase by 21.4% by the 2080s under RCP8.5. These events will have a median annual duration of 10.3 days/year, over 10 times the historical value.

54 ENVIRONMENTAL SCIENCES↗

Error Characteristics and Scale Dependence of Current Satellite Precipitation Estimates Products in Hydrological Modeling

Satellite precipitation estimates (SPEs) are promising alternatives to gauge observations for hydrological applications (e.g., streamflow simulation), especially in remote areas with sparse observation networks. However, the existing SPEs products are still biased due to imperfections in retrieval algorithms, data sources and post-processing, which makes the effective use of SPEs a challenge, especially at different spatial and temporal scales. In this study, we used a distributed hydrological model to evaluate the simulated discharge from eight quasi-global SPEs at different spatial scales and explored their potential scale effects of SPEs on a cascade of basins ranging from approximately 100 to 130,000 km 2 . The results indicate that, regardless of the difference in the accuracy of various SPEs, there is indeed a scale effect in their application in discharge simulation. Specifically, when the catchment area is larger than 20,000 km2, the overall performance of discharge simulation emerges an ascending trend with the increase of catchment area due to the river routing and spatial averaging. Whereas below 20,000 km 2 , the discharge simulation capability of the SPEs is more randomized and relies heavily on local precipitation accuracy. Our study also highlights the need to evaluate SPEs or other precipitation products (e.g., merge product or reanalysis data) not only at the limited station scale, but also at a finer scale depending on the practical application requirements. Here we have verified that the existing SPEs are scale-dependent in hydrological simulation, and they are not enough to be directly used in very fine scale distributed hydrological simulations (e.g., flash flood). More advanced retrieval algorithms, data sources and bias correction methods are needed to further improve the overall quality of SPEs.

DTVGM↗

Disentangling the hydrological and hydraulic controls on streamflow variability in Energy Exascale Earth System Model (E3SM) V2 – a case study in the Pantanal region

Abstract. Streamflow variability plays a crucial role in shaping the dynamics and sustainability of Earth's ecosystems, which can be simulated and projected by a river routing model coupled with a land surface model. However, the simulation of streamflow at large scales is subject to considerable uncertainties, primarily arising from two related processes: runoff generation (hydrological process) and river routing (hydraulic process). While both processes have impacts on streamflow variability, previous studies only calibrated one of the two processes to reduce biases in the simulated streamflow. Calibration focusing only on one process can result in unrealistic parameter values to compensate for the bias resulting from the other process; thus other water-related variables remain poorly simulated. In this study, we performed several experiments with the land and river components of the Energy Exascale Earth System Model (E3SM) over the Pantanal region to disentangle the hydrological and hydraulic controls on streamflow variability in coupled land–river simulations. Our results show that the generation of subsurface runoff is the most important factor for streamflow variability contributed by the runoff generation process, while floodplain storage effect and main-channel roughness have significant impacts on streamflow variability through the river routing process. We further propose a two-step procedure to robustly calibrate the two processes together. The impacts of runoff generation and river routing on streamflow are appropriately addressed with the two-step calibration, which may be adopted by developers of land surface and earth system models to improve the modeling of streamflow.

54 ENVIRONMENTAL SCIENCES↗

The impact of forest-controlled snow variability on late-season streamflow varies by climatic region and forest structure

Previous studies have documented how forests influence snow at fine spatial scales, but none have documented the influence that existing forest-snow variability has on streamflow. Here to test how much forest-controlled snow variability influences streamflow, a tiling parameterization based on classifications from high-resolution (1–3 m) vegetation maps was incorporated into the Distributed Hydrology and Soil Vegetation Model (DHSVM). Within each grid cell (90–150 m), the tiling parameterization simulated forest-snow variability with four independently evolving snowpacks. Each tile had unique radiation conditions to represent conditions underneath the canopy, in exposed areas, and along north- and south-facing forest edges. This tiled parameterization was used to test where and when detailed forest-snow modelling should be considered further and where and when the impacts are too small to be worth the effort. To test this, tiled model simulations of streamflow were compared to non-tiled model simulations in the Sierra Nevada, CA, the Jemez Mountains, NM, and the Eastern Cascades, WA. In Tuolumne, CA, the tiled model simulated little difference in grid cell average SWE, and late-season streamflow decreased by only 3%–4% compared to the non-tiled model. In Jemez, NM, the tiled model decreased late-season streamflow by 18% due to increased sublimation. In Chiwawa, WA, the tiled model increased late-season streamflow by 15% due to high shortwave radiation attenuation and less longwave radiation enhancement from the forest. Furthermore, within the Chiwawa, a substantial silvicultural practice was synthetically implemented to increase the north-facing edge's fractional area. This silvicultural experiment, which used the same fractional forest area in all simulations increased late-season streamflow by 35% compared to tiled model simulations that did not represent forest edges. In conclusion, representing forest-SWE variability had an effect on late-season streamflow in some watersheds but not in others based on the fractional area of the forest edges, forest characteristics, and climate conditions.

54 ENVIRONMENTAL SCIENCES↗