An interpretable and scalable model for rapid flood extent forecasting using satellite imagery and machine learning with rotated EOF analysis
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
This study focuses on the feasibility of using SAR interferometry to support 2D hydrodynamic model calibration and provide water storage change in the floodplain. Two-dimensional (2D) flood inundation modeling has been widely studied using storage cell approaches with the availability of high resolution, remotely sensed floodplain topography. The development of coupled 1D/2D flood modeling has shown improved calculation of 2D floodplain inundation as well as channel water elevation. Most floodplain model results have been validated using remote sensing methods for inundation extent. However, few studies show the quantitative validation of spatial variations in floodplain water elevations in the 2D modeling since most of the gauges are located along main river channels and traditional single track satellite altimetry over the floodplain are limited. Synthetic Aperture Radar (SAR) interferometry recently has been proven to be useful for measuring centimeter-scale water elevation changes over the floodplain. In the current study, we apply the LISFLOOD hydrodynamic model to the central Atchafalaya River Basin, Louisiana, during a 62 day period from 1 April to 1 June 2008 using two different calibration schemes for Manning's n. First, the model is calibrated in terms of water elevations from a single in situ gauge that represents a more traditional approach. Due to the gauge location in the channel, the calibration shows more sensitivity to channel roughness relative to floodplain roughness. Second, the model is calibrated in terms of water elevation changes calculated from ALOS PALSAR interferometry during 46 days of the image acquisition interval from 16 April 2008 to 1 June 2009. Since SAR interferometry receives strongly scatters in floodplain due to double bounce effect as compared to specular scattering of open water, the calibration shows more dependency to floodplain roughness. An iterative approach is used to determine the best-fit Manning's n for the two different calibration approaches. Results suggest similar floodplain roughness but slightly different channel roughness. However, application of SAR interferometry provides a unique view of the floodplain flow gradients, not possible with a single gauge calibration. These gradients, allow improved computation of water storage change over the 46-day simulation period. Overall, the results suggest that the use of 2D SAR water elevation changes in the Atchafalaya basin offers improved understanding and modeling of floodplain hydrodynamics.
Using the 2017 Hurricane Harvey flood event as a test case, this study set up a series of sensitivity analyses to highlight three challenges associated with large-scale flood inundation modeling, including (a) model parameterization, (b) errors in digital elevation models, and (c) effects of reservoir retention. Driven by radar-based hourly rainfall data, a series of hydrologic-hydraulic models including the VIC hydrologic model, RAPID routing model, and Flood2D-GPU hydrodynamic model are set up over Harris County, Texas, to simulate flood inundation and hazards. The results demonstrate the importance of hydrologic parameters in improving flood modeling. For a large flood event such as Hurricane Harvey, the effect of the initial water depths is insignificant. The Manning's n values may increase the peak water depth by ~1%, the flood extents by 65km 2 , and the high danger zone by ~6%. On the contrary, the bathymetry correction factors may reduce the flood extent by ~1.4% and the high-danger zone by ~4%. Reducing the reservoir storage capacity to 1% may increase the flood extent by ~4% and the high-danger zone by ~17%. This study may provide supporting information to guide and prioritize the development of future high-performance computing hydrodynamic large-scale flood simulations.
Polymer rheology can have either a positive or a negative effect on polymer flooding performance under varied circumstances. Many researchers have studied the effect of polymer rheology in a vertical well, but no field scale studies have been conducted to investigate whether polymer rheology is beneficial to polymer flooding in heavy oil reservoirs developed by horizontal wells. In this paper, we conducted a numerical simulation study to examine the effect of HPAM polymer rheology on a polymer flooding pilot, which is the first-ever project conducted on a heavy oil reservoir from Alaska North Slope (ANS) developed by horizontal wells. Three rheology types were considered in the study including the apparent viscosity measured during coreflooding of using a HPAM polymer, the bulk viscosity measured with a viscometer, and a Newtonian flow model. The results suggest that using the bulk viscosity in simulation underestimates the conformance control and the water-oil-ratio reduction capability of the HPAM polymer solution. When the apparent viscosity is used, the incremental oil and sweep were largely increased, and the optimal recovery period of polymer flooding was extended greatly, especially for the heterogeneous formations. Therefore, the rheology type of polymer plays a significant role in the incremental oil recovery and injection profile of the horizontal well system given the pilot test conditions. This study has provided practical guidance to field operators for the ongoing polymer flooding pilot on ANS and will also provide valuable information for other polymer projects conducted in similar conditions.
River flooding can have a detrimental impact on a community by causing loss of life, loss or damage to property, and damage to infrastructure. Having the capability to forecast flooding events can prevent the loss of life and mitigate damage to property. A programmatic, data-driven approach using deep learning to forecast a stream’s gauge height every 6 hours has been developed by NASA’s Short-term Prediction Research and Transition Center (SPoRT) and is currently in operation. Based on our end user engagement and feedback, this medium-term product has been successfully adopted by several National Weather Service (NWS) Forecast Offices and River Forecast Centers (RFC). SPoRT is currently researching and developing a short-term hourly deep learning model which will be useful in forecasting flash flood events.
Waikīkī, Hawai‘i serves as a valuable case study for flood vulnerability due to its low elevation, present-day exposure to multiple flood sources, and role as the state’s economic engine. This study applies the high-resolution, open-source WRF-Hydro-CUFA (Coastal Urban Flood Application) model to simulate compound flood interactions among pluvial, fluvial, coastal, storm drain-driven, and subsurface processes, and to evaluate how flood characteristics in Waikīkī evolve under projected sea-level rise (SLR). Simulated water levels in an estuarine drainage canal, serving as the primary outlet for gravity-flow drainage from Waikīkī, show strong agreement with observations (r ≥ 0.935, R 2 ≥ 0.874), accurately reproducing flooding observed during recent 5-year and 50-year storm events. Simulations incorporating SLR demonstrate a shift in dominant flood drivers as sea levels rise. In addition, results identify critical thresholds at which elevated estuarine tailwater conditions cause widespread drainage backflow, eventually occurring in the absence of rainfall. These findings suggest that tidally driven flooding will progressively dominate over rainfall-driven events, increasing public exposure to highly contaminated canal and groundwater sources. The results underscore the urgent need to integrate water quality issues and drainage failure into coastal flood management strategies for low-lying urban environments, particularly those reliant on polluted estuarine waterways for stormwater management.
We evaluate performance of the Catchment Land Surface Model (CLSM) under flood conditions after the assimilation of observations of the terrestrial water storage anomaly (TWSA) from NASA's Gravity Recovery and Climate Experiment (GRACE). Assimilation offers three key benefits for the viability of GRACE observations to operational applications: (1) near-real time analysis; (2) a downscaling of GRACE's coarse spatial resolution; and (3) state disaggregation of the vertically-integrated TWSA. We select the 2011 flood event in the Missouri river basin as a case study, and find that assimilation generally made the model wetter in the months preceding flood. We compare model outputs with observations from 14 USGS groundwater wells to assess improvements after assimilation. Finally, we examine disaggregated water storage information to improve the mechanistic understanding of event generation. Validation establishes that assimilation improved the model skill substantially, increasing regional groundwater anomaly correlation from 0.58 to 0.86. For the 2011 flood event in the Missouri river basin, results show that groundwater and snow water equivalent were contributors to pre-event flood potential, providing spatially-distributed early warning information.
Evaluating extreme precipitation events (EPEs) is essential for building climate-resilient water management strategies, but it remains a major challenge in ungauged basins. Using the 26,070 km 2 Wadi al-Rummah basin in central Saudi Arabia as a case study, we developed an alternative, reliable, cost-effective satellite-based framework that combines empirically derived EPE thresholds, imagery-calibrated 2D hydrodynamic modeling, GRACE water-storage diagnostics, and bias-corrected CMIP6 projections to assess flood hazards and recharge potential under current and future climate scenarios in ungauged basins. The integrated approach and the resulting findings followed four key steps: (1) Identified a 22.5 mm EPE threshold, the 80th percentile of 3-day GPM/IMERG rainfall (2000–2024), aligned with flood-triggering events (Nov 2018: 23–28 mm; Apr, 2023: 42 mm); (2) Developed and calibrated a RiverFlow2D model using Sentinel-2 and PlanetScope imagery for the November 2018 flood, accurately reproducing flood depth and extent (RMSE ≤0.31 m; fuzzy-Dice ≥0.91), and estimating runoff (41 %), infiltration (25 %), and evaporation (34 %); (3) Independently validated the model with the April 2023 event (RMSE ≤0.35 m; fuzzy-Dice ≥0.86); (4) Conducted climate projections (2025–2100) from five bias-corrected NEX-GDDP CMIP6 models that revealed a 34 % increase in EPE intensity under SSP2-4.5 and 48 % under SSP5-8.5 scenarios, relative to 20th-century baselines. Our findings indicate that while intensifying extremes in the 21st century increase flood risk, the results highlight the potential for episodic recharge if effective retention strategies are employed, and offer a transferable model for climate-informed planning in data-scarce arid regions.
For this study, a high-resolution, summit-to-sea unstructured-grid model was used to simulate two compound flooding events in different geomorphic settings in Japan: the July 2012 flood in Kumamoto city and July 2018 flood in Okayama City, both caused by a torrential rainfall during seasonal Meiyu front but otherwise exhibiting quite different compound flood characteristics. The model shows good performance in simulating flooding extent; e.g., the Hit Rate for the inundation event in 2018 exceeds 0.9. Sensitivity tests were conducted to determine the roles/significance of levee structures along major rivers in coastal regions. Our results indicate that riverine and oceanic factors as well as levees significantly contribute to the observed flooding extents, and we are able to quantify the contributions from each factor at different locations. The model has been implemented operationally as the first nationwide 3-day compound flooding forecast system for Japan.
The 'flooded agglomerate' model of the Teflon-bonded gas diffusion electrode is discussed. A mathematical treatment of the 'flooded agglomerate' model is given; it can be used to predict the performance of the electrode as a function of measurable physical parameters.
Detailed land cover classifications were performed on the Thematic Mapper Simulator (TMS) and MSS data of the Clinton River Basin (acquired on August 19, 1981, and June 28, 1980, respectively) using supervised classification techniques. Differences in interclass separability were compared to select several promising TMS band combinations, selected from the 27 covering the Clinton River Basin. The TMS data produced a more accurate and spatially contiguous classification than MSS for this study site. While the accuracy of the 4-band TM data set was as good as the 7-band, the 3-band TMS data sets were also better than the MSS. These results indicate that both the increased spectral discrimination and spatial resolution contribute to improved classification accuracy. The possibility of reducing the data analysis burden associated with large TM data volumes through effective band selection therefore appears promising. The implications of the improved classification accuracy of TMS data are important for hydrologic and economic modeling. In particular, the higher accuracies for the developed categories (residential and commercial) should improve the predictions of runoff in flood forecasting models and of flood damage for damage calculation models appreciably.
Here, we report on methods and results for a model-based flood hazard assessment we have conducted for the Hudson River from New York City to Troy/Albany at the head of tide. Our recent work showed that neglecting freshwater flows leads to underestimation of peak water levels at up-river sites and neglecting stratification (typical with two-dimensional modeling) leads to underestimation all along the Hudson. As a result, we use a three-dimensional hydrodynamic model and merge streamflows and storm tides from tropical and extratropical cyclones (TCs, ETCs), as well as wet extratropical cyclone (WETC) floods (e.g. freshets, rain-on-snow events). We validate the modeled flood levels and quantify error with comparisons to 76 historical events. A Bayesian statistical method is developed for tropical cyclone streamflows using historical data and consisting in the evaluation of (1) the peak discharge and its pdf as a function of TC characteristics, and (2) the temporal trend of the hydrograph as a function of temporal evolution of the cyclone track, its intensity and the response characteristics of the specific basin. A k-nearest-neighbors method is employed to determine the hydrograph shape. Out of sample validation tests demonstrate the effectiveness of the method. Thus, the combined effects of storm surge and runoff produced by tropical cyclones hitting the New York area can be included in flood hazard assessment. Results for the upper Hudson (Albany) suggest a dominance of WETCs, for the lower Hudson (at New York Harbor) a case where ETCs are dominant for shorter return periods and TCs are more important for longer return periods (over 150 years), and for the middle-Hudson (Poughkeepsie) a mix of all three flood events types is important. However, a possible low-bias for TC flood levels is inferred from a lower importance in the assessment results, versus historical event top-20 lists, and this will be further evaluated as these preliminary methods and results are finalized. Future funded work will quantify the influences of sea level rise and flood adaptation plans (e.g. surge barriers). It would also be valuable to examine how streamflows from tropical cyclones and wet cool-season storms will change, as this factor will dominate at upriver locations.
The Bracknell (U.K.) Meteorological Office runs a global weather model twice a day, providing the following data: surface and radiosonde; aircraft reports; and satellite soundings and wind. A human forecast is made every six hours. The model runs on a 150 km grid with 15 levels, and takes about four minutes on a Cyber-205. The standard output from the global products are wind, temperature, height, tropopause, and maximum wind. Various experiments have been conducted to see if short-range forecasters could improve on the upper-wind forecasts over the numerical model; the numerical model remains of paramount importance. Small-scale models are being run in the U.S. and the U.K. A fine-mesh model covers Europe and the Atlantic. A mesoscale model is under development. A great deal of verification work is done to see how good the models are.
Increased intensity and frequency of floods raise concerns about the release and transport of contaminated soil and sediment to and from rivers and streams. To model these processes during flooding events, we developed an External Coupler in Python to link the Hydrologic Engineering Center-River Analysis System (HEC-RAS) 2D hydrodynamic model to the Water Quality Analysis Simulation Program (WASP). Accurate data transfer from a hydrodynamic model to a water quality model is critical. Our test results showed the External Coupler successfully linked HEC-RAS and WASP and addressed technical challenges in aggregating flow data and conserving mass during the flood event. We ran the coupled models for a 100-year flood event to calculate flood-induced transport of sediment-associated arsenic in Woodbridge Creek, NJ. Change in surface sediment and arsenic at the end of 48-h flood simulation ranged from a net loss of 13.5 cm to a net gain of 11.6 cm, and 16.2 to 2.9 mg/kg, respectively, per model segment, which demonstrates the capability of the coupled model for simulating sediment and contaminant transport in flood.
Central Texas is subject to infrequent, unusually intense rainstorms which cause extremely rapid runoff from drainage basins developed on the deeply dissected limestone and marl bedrock of the Edwards Plateau. One approach to flood hazard evaluation in this area is a parametric model relating flood hydrograph characteristics to quantitative geomorphic properties of the drainage basins. The preliminary model uses multiple regression techniques to predict potential peak flood discharge from basin magnitude, drainage density, and ruggedness number. After mapping small catchment networks from remote sensing imagery, input data for the model are generated by network digitization and analysis by a computer assisted routine of watershed analysis. The study evaluated the network resolution capabilities of the following data formats: (1) large-scale (1:24,000) topographic maps, employing Strahler's "method of v's," (2) standard low altitude black and white aerial photography (1:13,000 and 1:20,000 scales), (3) NASA - generated aerial infrared photography at scales ranging from 1:48,000 to 1:123,000, and (4) Skylab Earth Resources Experiment Package S-190A and S-190B sensors (1:750,000 and 1:500,000 respectively).
This study showcases a global, heterogeneously coupled total water level system wherein salinity and temperature outputs from a coarser-resolution (~12 km) ocean general circulation model are used to calculate density-driven terms within a global, higher-resolution (~2.5 km) depth-averaged total water level model. We demonstrate that the inclusion of baroclinic forcing in the barotropic model requires modification of the internal wave drag term to prevent excess degradation of tidal results compared to the barotropic model. By scaling the internal tide dissipation by an easy to calculate dissipation ratio, the resulting heterogeneously coupled model has complex root mean square errors (RMSE) of 2.27 cm in the deep ocean and 12.16 cm in shallow waters for the M 2 tidal constituent. While this represents a 10%–20% deterioration as compared to the barotropic model, the improvements in total water level prediction more than offset this degradation. Global median RMSE compared to observations of total water levels, 30-day sea levels, and non-tidal residuals improve by 1.86 (18.5%), 2.55 (42.5%), and 0.36 (5.3%) cm respectively. The drastic improvement in model performance highlights the importance of including density-driven effects within global hydrodynamic models and will help to improve the results of both hindcasts and forecasts in modeling extreme and nuisance flooding. With only an 11% increase in model run time compared to the fully barotropic total water level model, this approach paves the way for high resolution coastal water level and flood models to be used alongside climate models, improving operational forecasting of total water levels.
Rising sea levels due to climate change increasingly threaten medical infrastructure through flooding. This study develops machine learning models to predict flood exposure for 11,508 medical facilities in the southeastern coastal regions of the United States by integrating datasets including meteorological, hydrological, topographic, and geological data, the Natural Risk Index, and historical flood records from NASA, HIFLD, and FEMA. Six regression models, namely Linear Regression, Support Vector Regression, Random Forest, k-Nearest Neighbors, XGBoost, and Artificial Neural Networks, are trained using 16 explanatory variables identified through literature review and correlation analysis. Data preprocessing employs the SMOGN for class imbalance and Winsorization for outliers. Model performance is evaluated using MAE, MSE, and RMSE, with Random Forest and XGBoost models achieving the highest performance (MSE of 2.58e-5 and 3.69e-5, respectively). This multifactorial approach allows the models to capture complex flood-influencing relationships, enhancing adaptability and performance across geographic regions. Future work focuses on expanding across the U.S. and developing a near real-time flood monitoring system.
The 1993 US Midwest Flood produced record levels of flooding, in the Mississippi River Basin. This flooding resulted from repeated frontal passages and mesoscale convective complexes (MCCs) during the months of June and July. A better understanding of processes that influenced MCC development during the 1993 Flood may lead to improved forecasts of heavy precipitation and flooding. Here, we consider the impact of soil moisture on MCC development during a two-day period (June 23-24) of the 1993 US Midwest Flood. The purpose of this study is to assess the importance of soil moisture distribution on the timing, intensity, and location of heavy precipitation. In this study, the MM5-PLACE Atmosphere/Land-Surface Model is utilized. The atmospheric component consists of the Penn State/NCAR MM5 mesoscale model, and the land-surface component consists of the Goddard Parameterization for Land Atmosphere-Cloud Exchange (PLACE). Initial soil moisture is provided from two sources: 1) NCEP reanalysis, and 2) Antecedent Precipitation Index (API) using NOAA rain gauge measurements as a proxy for soil moisture. NCEP reanalysis provides coarse resolution initial soil moisture (2.5 degree), while API provides high resolution initial soil moisture (10-200 km depending on NOAA rain gauge spacing). Initial results indicate that the distribution of soil moisture has a significant impact on the timing and location of heavy precipitation during this two-day flood event. Precipitation in simulations with high resolution initial soil moisture agrees more closely with observed precipitation. These results suggest that high resolution soil moisture observations are necessary to accurately predict severe storm development, heavy precipitation, and subsequent flooding.