skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: 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.  more » « less
Award ID(s):
2053429
PAR ID:
10566883
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
61
Issue:
1
ISSN:
0043-1397
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Numerous studies have examined the reliability of various precipitation products over the Mekong River Basin (MRB) and modeled its basin hydrology. However, there is a lack of comprehensive studies on precipitation‐induced uncertainties in hydrological simulations using process‐based land surface models. This study examines the propagation of precipitation uncertainty into hydrological simulations over the entire MRB using the Community Land Model version 5 (CLM5) at a high spatial resolution of 0.05° (∼5 km) and without any parameter calibration. Simulations conducted using different precipitation datasets are compared to investigate the discrepancies in streamflow, terrestrial water storage (TWS), soil moisture, and evapotranspiration (ET) caused by precipitation uncertainty. Results indicate that precipitation is a key determinant of simulated streamflow in the MRB; peak flow and soil moisture are particularly sensitive to precipitation input. Further, precipitation data with a higher spatial resolution did not improve the simulations, contrary to the common perception that using meteorological forcing with higher spatial resolution would improve hydrological simulations. In addition, since high flow indicators are particularly influenced by precipitation data, the choice of precipitation data could directly impact flood pulse simulations in the MRB. Notable differences are also found among TWS, soil moisture, and ET simulated using different precipitation products. Moreover, TWS, soil moisture, and ET exhibit a varying degree of sensitivity to precipitation uncertainty. This study provides crucial insights on precipitation‐induced uncertainties in process‐based hydrological modeling and uncovers these uncertainties in the MRB. 
    more » « less
  2. Abstract. Two-dimensional (2D), depth-averaged shallow water equation (SWE) models are routinely used to simulate flooding in coastal areas – areas that often include vast networks of channels and flood-control topographic features and/or structures, such as barrier islands and levees. Adequately resolving these features within the confines of a 2D model can be computationally expensive, which has led to coupling 2D simulation tools to less expensive one-dimensional (1D) models. Under certain 1D–2D coupling approaches, this introduces internal constraints that must be considered in the generation of the 2D computational mesh used. In this paper, we further develop an existing automatic unstructured mesh generation tool for SWE models, ADMESH+, to sequentially (i) identify 1D constraints from the raw input data used in the mesh generation process, namely the digital elevation model (DEM) and land–water delineation data; (ii) distribute grid points along these internal constraints, according to feature curvature and user-prescribed minimum grid spacing; and (iii) integrate these internal constraints into the 2D mesh size function and mesh generation processes. The developed techniques, which include a novel approach for determining the so-called medial axis of a polygon, are described in detail and demonstrated on three test cases, including two inland watersheds with vast networks of channels and a complex estuarine system on the Texas, USA, coast. 
    more » « less
  3. Abstract Uncertainty attribution in water supply forecasting is crucial to improve forecast skill and increase confidence in seasonal water management planning. We develop a framework to quantify fractional forecast uncertainty and partition it between (1) snowpack quantification methods, (2) variability in post‐forecast precipitation, and (3) runoff model errors. We demonstrate the uncertainty framework with statistical runoff models in the upper Tuolumne and Merced River basins (California, USA) using snow observations at two endmember spatial resolutions: a simple snow pillow index and full‐catchment snow water equivalent (SWE) maps at 50 m resolution from the Airborne Snow Observatories. Bayesian forecast simulations demonstrate a nonlinear decrease in the skill of statistical water supply forecasts during warm snow droughts, when a low fraction of winter precipitation remains as SWE. Forecast skill similarly decreases during dry snow droughts, when winter precipitation is low. During a shift away from snow‐dominance, the uncertainty of forecasts using snow pillow data increases about 1.9 times faster than analogous forecasts using full‐catchment SWE maps in the study area. Replacing the snow pillow index with full‐catchment SWE data reduces statistical forecast uncertainty by 39% on average across all tested climate conditions. Attributing water supply forecast uncertainty to reducible error sources reveals opportunities to improve forecast reliability in a warmer future climate. 
    more » « less
  4. Abstract Unprecedented floods from extreme rainfall events worldwide emphasize the need for flood inundation mapping for floodplain management and risk reduction. Access to flood inundation maps and risk evaluation tools remains challenging in most parts of the world, particularly in rural regions, leading to decreased flood resilience. The use of hydraulic and hydrodynamic models in rural areas has been hindered by excessive data and computational requirements. In this study, we mapped the flood inundation in Huron Creek watershed, Michigan, USA for an extreme rainfall event (1000-year return period) that occurred in 2018 (Father’s Day Flood) using the Height Above Nearest Drainage (HAND) model and a synthetic rating curve developed from LIDAR DEM. We compared the flood inundation extent and depth modeled by the HAND with flood inundation characteristics predicted by two hydrodynamic models, viz., HEC-RAS 2D and SMS-SRH 2D. The flood discharge of the event was simulated using the HEC-HMS hydrologic model. Results suggest that, in different channel segments, the HAND model produces different degrees of concurrence in both flood inundation extent and depth when compared to the hydrodynamic models. The differences in flood inundation characteristics produced by the HAND model are primarily due to the uncertainties associated with optimal parameter estimation of the synthetic rating curve. Analyzing the differences between the HAND and hydrodynamic models also highlights the significance of terrain characteristics in model predictions. Based on the comparable predictive capability of the HAND model to map flood inundation areas during extreme rainfall events, we demonstrate the suitability of the HAND-based approach for mitigating flood risk in data-scarce, rural regions. 
    more » « less
  5. Abstract Existing stochastic rainfall generators (SRGs) are typically limited to relatively small domains due to spatial stationarity assumptions, hindering their usefulness for flood studies in large basins. This study proposes StormLab, an SRG that simulates precipitation events at 6‐hr and 0.03° resolution in the Mississippi River Basin (MRB). The model focuses on winter and spring storms caused by water vapor transport from the Gulf of Mexico—the key flood‐generating storm type in the basin. The model generates anisotropic spatiotemporal noise fields that replicate local precipitation structures from observed data. The noise is transformed into precipitation through parametric distributions conditioned on large‐scale atmospheric fields from a climate model, reflecting spatial and temporal nonstationarity. StormLab can produce multiple realizations that reflect the uncertainty in fine‐scale precipitation arising from a specific large‐scale atmospheric environment. Model parameters were fitted monthly from December–May, based on storms identified from 1979 to 2021 ERA5 reanalysis data and Analysis of Record for Calibration (AORC) precipitation. StormLab then generated 1,000 synthetic years of precipitation events based on 10 CESM2 ensemble simulations. Empirical return levels of simulated annual maxima agree well with AORC data and show an overall increase in 1‐ to 500‐year events in the future period (2022–2050). To our knowledge, this is the first SRG simulating nonstationary, anisotropic high‐resolution precipitation over continental‐scale river basins, demonstrating the value of conditioning such stochastic models on large‐scale atmospheric variables. StormLab provides a wide range of extreme precipitation scenarios for design floods in the MRB and can be further extended to other large river basins. 
    more » « less