Search NASA⌕ Search

SEARCH · Search NASA

Results for “Rainfall Runoff”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Significant Rural Contributions to Severe Coastal Urban Flooding Through Flow Connectivity

Coastal cities are increasingly threatened by hurricane-induced compound flooding. However, the relative contributions of rainfall-surge interactions and other flood drivers, such as runoff from surrounding rural and local impervious areas, in shaping urban flooding remain unclear. Here, we use a multi-scale Earth system modeling framework to disentangle the contributions of diverse flooding drivers. The results show that urban flooding is primarily affected by topographic factors (low elevation, flat slopes, and localized flow constrictions), rather than increased runoff from their impervious surfaces. Importantly, runoff from rural areas contributes more to severely flooded urban areas with inundation depth exceeding 1.7 m than local runoff. Although rainfall–surge interactions and sea-level rise exacerbate inundation along the coastline, over 90% of the compounding effects are absorbed by wetlands. By resolving flooding processes across scales and components, this study emphasizes the role of runoff from surrounding rural areas in shaping zones of severe urban flooding.

58 GEOSCIENCES↗
DOE OSTI2025

Exploring Flood Predictability in Taiwan through Coupled Atmospheric–Hydrological and High-Performance Hydrodynamic Models

Effective flood simulation capabilities can tremendously support early warning and disaster prevention. To examine the applicability of a fully physics-based and high-performance flood simulation and forecasting modeling framework for a flood-prone region in Taiwan, we conduct a numerical experiment that couples the Weather Research and Forecasting (WRF) Model, WRF-Hydrological modeling system (WRF-Hydro), and the Two-Dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) to perform integrated rainfall, streamflow, and flood simulations. Furthermore, we first use the coupled WRF and WRF-Hydro (WWH) to predict rainfall and streamflow and then drive TRITON with the predicted streamflow hydrographs to simulate flood depth and inundation area. With the refined spatial resolution and parameterization, this framework can better predict rainfall with reasonable spatial patterns. Although WWH could overestimate the amount of rainfall in some areas, the uncertain rainfall–streamflow predictions produce reasonable flood maps able to pinpoint regions at risk of flooding. In terms of model efficiency, the graphics processing unit–based computation can yield a speed-up factor as high as ∼13 compared to the central processing unit–based computation, promoting the efficacy of the coupled modeling framework in practical real-time flood forecasting.

Coupled models↗
DOE OSTI2025

Divergent responses of historic rain-on-snow flood extremes to a warmer climate

Global warming is altering flood risks induced by rain-on-snow events. However, decision-makers lack guidance on how rain-on-snow induced extreme floods could be altered with warming. Here, storyline analyses using a kilometer-scale land surface model reveal diverse responses of four historically-impactful, decision-relevant rain-on-snow induced extreme flood events over the contiguous U.S. to warming, due to alterations in their water budgets. For the 2017-Feb California floods, runoff first increases and then decreases with warming, peaking under the +3 K scenario, while runoff of the 2017-Jan California floods increases monotonically by ~53%/K. Contrastingly, runoff of the 1996-Jan Mid-Atlantic floods decreases gradually with warming. Despite these differences, warming generally shifts flood-generating regimes along elevation profiles. High elevations could experience notably increased runoff, while low elevations encounter a shift from rain-on-snow-driven to rainfall-dominated runoff. These findings underscore the need for flood control planning to quantify region- and elevation-specific changes in rain-on-snow events in a warmer climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Next-Generation Intensity-Duration-Frequency Curves for Diverse Land across the Continental United States

The current methods for designing hydrological infrastructure rely on precipitation-based intensity-duration-frequency curves. However, they cannot accurately predict flooding caused by snowmelt or rain-on-snow events, potentially leading to underdesigned infrastructure and property damage. To address these issues, next-generation intensity-duration-frequency (NG-IDF) curves have been developed for the open condition, characterizing water available for runoff from rainfall, snowmelt, and rain-on-snow. However, they lack consideration of land use land cover (LULC) factors, which can significantly affect runoff processes. We address this limitation by expanding open area NG-IDF dataset to include eight vegetated LULCs over the continental United States, including forest (deciduous, evergreen, mixed), shrub, grass, pasture, crop, and wetland. This NG-IDF 2.0 dataset offers a comprehensive analysis of hydrological extreme events and their associated drivers under different LULCs at a continental scale. It will serve as a useful resource for improving standard design practices and aiding in the assessment of infrastructure design risks. Additionally, it provides useful insights into how changes in LULC impact flooding magnitude, mechanisms, timing, and snow water supply.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Temporal persistence of postfire flood hazards under present and future climate conditions in southern Arizona, USA

Changes to soil hydraulic properties that reduce infiltration capacity following fire can increase flash flood risks. These risks are exacerbated by rainfall intensification associated with a warming climate. However, the potential effects of climate-change-driven rainfall intensification on postfire floods remain largely unexplored. Using rainfall and runoff observations from a 49.4 km 2 watershed in southern Arizona, USA, and a hydrologic model (KINEROS2), we examined the temporal evolution following a historic fire of three crucial hydrologic parameters: soil saturated hydraulic conductivity (K sp ), net capillary drive (G p ), and hydraulic roughness (n c ). We explored how the effect of fire on these parameters may influence peak flow under future climate scenarios derived from CMIP6, specifically the medium emissions scenario (SSP2-4.5) and high emissions scenario (SSP5-8.5). Results demonstrate an increase in K sp from 11 mm h −1 in the first postfire year to 60 mm h −1 in postfire year 3. G p similarly increased from 19 mm in the first postfire year to 30 mm in the third, while nc was relatively constant. The highest simulated Q p occurred in the first postfire year. Under the SSP2-4.5 scenario, the likelihood of a 100-year flood is projected to be twice as large by the middle of the century relative to its historical magnitude. Simulations further indicate that the maximum expected discharge associated with a postfire flood, as derived from historical data, could be triggered by a 10-year rainstorm under the SSP5-8.5 scenario by the late century. Simulations also demonstrate that rainfall intensification will lead to greater persistence of elevated flood hazards following fire by the late century under both the SSP2-4.5 and the SSP5-8.5 scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Next-Generation Intensity-Duration-Frequency Curves for Climate-Resilient Infrastructure Design: Advances and Opportunities

National and international security communities (e.g., U.S. Department of Defense) have shown increasing attention for innovating critical infrastructure and installations due to recurring high-profile flooding events in recent years. The standard infrastructure design approach relies on local precipitation-based intensity-duration-frequency (PREC-IDF) curves that do not account for snow process and assume stationary climate, leading to high failure risk and increased maintenance costs. This paper reviews the recently developed next-generation IDF (NG-IDF) curves that explicitly account for the mechanisms of extreme water available for runoff including rainfall, snowmelt, and rain-on-snow under nonstationary climate. The NG-IDF curve is an enhancement to the PREC-IDF curve and provides a consistent design approach across rain- to snow-dominated regions, which can benefit engineers and planners responsible for designing climate-resilient facilities, federal emergency agencies responsible for the flood insurance program, and local jurisdictions responsible for developing design manuals and approving subsequent infrastructure designs. Further, we discuss the recent advances in climate and hydrologic science communities that have not been translated into actional information in the engineering community. To bridge the gap, we advocate that building climate-resilient infrastructure goes beyond the traditional local design scale where engineers rely on recipe-based methods only; the future hydrologic design is a multi-scale problem and requires closer collaboration between climate scientists, hydrologists, and civil engineers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Variability of Snow and Rainfall Partitioning Into Evapotranspiration and Summer Runoff Across Nine Mountainous Catchments

Abstract Understanding the partitioning of snow and rain contributing to either catchment streamflow or evapotranspiration (ET) is of critical relevance for water management in response to climate change. To investigate this partitioning, we use endmember splitting and mixing analyses based on stable isotope ( 18 O) data from nine headwater catchments in the East River, Colorado. Our results show that one third of the snow partitions to ET and 13% of the snowmelt sustains summer streamflow. Only 8% of the rainfall contributes to the summer streamflow, because most of the rain (67%) partitions to ET. The spatial variability of precipitation partitioning is mainly driven by aspect and tree cover across the sub‐catchments. Catchments with higher tree cover have a higher share of snow becoming ET, resulting in less snow in summer streamflow. Summer streamflow did not contain more rain with higher rainfall sums, but more rain was taken up in ET.

58 GEOSCIENCES↗
DOE OSTI2025

How Flood Hazards in a Warming Climate Could Be Amplified by Changes in Spatiotemporal Patterns and Mechanisms of Water Available for Runoff

Prior research on climate change impacts on flooding has primarily focused on changes in extreme rainfall magnitudes, often neglecting snow processes and spatiotemporal storm patterns, such as hyetograph shapes and areal reduction factors (ARFs). This study examines projected changes in extreme water available for runoff (W) events in two snow-dominated basins in the western United States: the Yakima River Basin (YRB) in Washington State and the Walker River Basin (WRB) spanning the California-Nevada border. We analyze changes in W magnitudes, mechanisms, hyetograph shapes, and ARFs, and study their compounded impacts on flood hazard. Our findings suggest increased extreme W magnitudes across a large portion of the basins, with steeper or flatter hyetographs, and higher ARF values under the future climate. These changes are driven by a shift from seasonal snowmelt to more rain-on-snow events at higher elevations and by increased rainfall at lower elevations. We then use a single event-based rainfall-runoff model to estimate flood hazard changes based on extreme W magnitudes, hyetograph shapes, ARFs, and their compounded impacts. Our analysis reveals that focusing solely on the magnitude of changes in extreme W can significantly underestimate future flood hazards and uncertainties. Ignoring future changes in spatiotemporal patterns can underestimate future flood hazards by 63% and underestimate the uncertainty in future flood events by 18% in the WRB. These results underscore the necessity of incorporating spatiotemporal dynamics into future flood hazard assessments to provide a more accurate evaluation of potential impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Measuring and modeling soil moisture and runoff at solar farms using a disconnected impervious surface approach

Abstract Ground‐mounted photovoltaic sites are often treated as impervious surfaces in stormwater permits. This ignores the pervious soils beneath and between solar arrays and leads to an overestimation of runoff. Our objective was to improve solar farm stormwater hydrology models by explicitly considering the disconnected impervious nature of solar design and site characteristics. Experimental sites established on utility scale solar farms in Colorado, Georgia, Minnesota, New York, and Oregon had perennial vegetative plantings with mean precipitation ranging from 40.6 to 124.5 cm, and soil texture ranging from loamy sand to clay. Soil moisture measurements were collected beneath arrays, under drip edges, and in the vegetated area between arrays at each site. Hydrus‐3D models for soil moisture and stormwater hydrology were developed that accounted for precipitation falling on solar panels, drip edge redistribution of rainfall, infiltration, and runoff in the pervious areas between solar arrays and beneath panels. Drip edge runoff averaged 3‐ to 10‐times incident precipitation at the New York and Minnesota sites, respectively. Root mean square error values between measured sub‐hourly soil moisture and predicted moisture for large measured single storm events averaged 0.029 across all five sites. Predicted runoff depths were strongly affected by precipitation depth, soil texture, soil profile depth, and soil bulk density. Runoff depths across the five experimental sites averaged 13%, 25%, and 45% of the 2‐, 10‐, and 100‐year design storm depths, clearly showing that these solar farms do not behave like impervious surfaces, but rather as disconnected impervious surfaces with substantial infiltration of runoff in the vegetated areas between and beneath solar arrays.

Agriculture↗
DOE OSTI2020

TRITON

TRITON is a 2D hydraulic, digital elevation model (DEM)-based model that accepts a variety of flow inputs including point-based inflow hydrographs, external boundary conditions, and runoff (i.e., excess rainfall). It provides information of flood maps and temporal evolution of water depth at each grid cell point at user-specified output time. Its unique feature is that TRITON can run on single or multiple CPUs and GPUs, proving to be very flexible in terms of system architecture.

Morales-Hernandez, Mario↗
DOE OSTI2020

Understanding the distinct impacts of MCS and non-MCS rainfall on the surface water balance in the central US using a numerical water-tagging technique

Warm-season rainfall associated with mesoscale convective systems (MCSs) in the central US is characterized by higher intensity and nocturnal timing compared to rainfall from non-MCS systems, suggesting their potentially different footprints on the land surface. To differentiate the impacts of MCS and non-MCS rainfall on the surface water balance, a water tracer tool embedded in the Noah land surface model with multi-parameterization options (WT-Noah-MP) is used to numerically “tag” water from MCS and non-MCS rainfall separately during April to August (1997-2018) and track their transit in the terrestrial system. From the water-tagging results, over 50% of warm-season rainfall leaves the surface-subsurface system through evapotranspiration by the end of August, but non-MCS rainfall contributes a larger fraction. However, MCS rainfall plays a more important role in generating surface runoff. These differences are mostly attributed to the rainfall intensity differences. The higher intensity MCS rainfall tends to produce more surface runoff through infiltration excess flow and drives a deeper penetration of the rainwater into the soil. Over 70% of the top 10 percentile runoff is contributed by MCS rainfall, demonstrating its important contribution to local flooding. In contrast, lighter intensity non-MCS rainfall resides mostly in the top layer and contributes more to evapotranspiration through soil evaporation. Diurnal timing of rainfall has negligible effects on the flux partitioning for both MCS and non-MCS rainfall. Differences in soil moisture profiles for MCS and non-MCS rainfall and the resultant evapotranspiration suggest differences in their roles in soil moisture-precipitation feedbacks and ecohydrology.

Hu, Huancui↗
DOE OSTI2025

Stream Temperature Responses to Summer Urban Rain Events Along the Savannah River

The hydrological urban heat island (HUHI) is a recent facet of the urban heat island (UHI), describing thermal enhancement of bodies of water in response to urbanization. Although the forefront of this work has been studied for large Metropolitan areas, this effect on developing cities is currently unknown. To locally measure and quantify HUHI effects in a developing city, we utilized Multi-Radar Multi Sensor (MRMS) radar-derived rainfall data to estimate rainfall in the Augusta Metropolitan Area (AMA) and the Savannah Metropolitan Area (SMA), two developing communities adjacent to the Savannah River. We analyzed mean temperature trends of two stream gauges measuring temperature at, and downstream of the central AMA and the SMA. Results show overall peak stream temperatures within AMA are reached within a 2hr timeframe with surge temperatures between 1-2 K, perhaps assisted by the presence of the Augusta Shoals further upstream. The SMA appears to have little HUHI effects due to urban greening and impacts to coastal meteorology. It is suggested that Moderate rainfall events (10-25mm of rainfall) have a significant correlation (~0.1 correlation coefficient with a significance level of 0.05) between stream temperature and rainfall, indicating potential warmer runoff input into the Savannah River. Although results seem promising, further research into moderate rainfall events are needed to determine the extent of the HUHI for AMA.

Wermter, Joseph E. [Savannah River National Labora↗
DOE OSTI2026

Flood Estimation under Snowmelt and Rain-on-Snow Processes in Alaska: A Military Installation Perspective

Accurate flood estimation in snow-dominated and high-latitude regions remains challenging due to complex interactions among rainfall, snowmelt, and rain-on-snow (ROS) processes, which are not captured in conventional precipitation-based intensity-duration-frequency (PREC-IDF) curves. This study evaluates the Next-Generation IDF (NG-IDF) framework, an extension of PREC-IDF that incorporates total water available for runoff (precipitation minus changes in snow water equivalent), in two contrasting Alaskan watersheds of Department of Defense (DoD) significance: the Little Chena River Basin (LCRB) in interior Alaska, a tributary of the Chena River that flows through Fort Wainwright and near Eielson Air Force Base, and the Upper Ship Creek Basin (USCB), which drains Joint Base Elmendorf-Richardson near Anchorage. Using long-term SNOTEL observations and event-based rainfall-runoff modeling, NG-IDF and PREC-IDF flood estimates were compared against observation-based flood frequency analyses. Results show that NG-IDF consistently reduces flood-estimation bias by 15–20% relative to PREC-IDF, particularly for snowmelt- and ROS-dominated events. In the interior LCRB, permafrost conditions can substantially amplify flood responses during snowmelt events, an effect not explicitly represented in standard design tools. These findings demonstrate that NG-IDF provides a more physically consistent and transferable framework for flood estimation in cold regions, with potential relevance to mission-critical DoD installation resilience. Projected increases in ROS frequency and permafrost degradation across Alaska further emphasize the need to integrate physics-based hydrologic models that explicitly represent snow and permafrost processes to enhance design resilience and operational readiness for military and other critical infrastructure.

Yan, Hongxiang (ORCID:000000022387403X)↗
DOE OSTI2022

Augmentation of WRF-Hydro to simulate overland-flow- and streamflow-generated debris flow susceptibility in burn scars

In steep wildfire-burned terrains, intense rainfall can produce large runoff that can trigger highly destructive debris flows. However, the ability to accurately characterize and forecast debris flow susceptibility in burned terrains using physics-based tools remains limited. Here, we augment the Weather Research and Forecasting Hydrological modeling system (WRF-Hydro) to simulate both overland and channelized flows and assess postfire debris flow susceptibility over a regional domain. We perform hindcast simulations using high-resolution weather-radar-derived precipitation and reanalysis data to drive non-burned baseline and burn scar sensitivity experiments. Our simulations focus on January 2021 when an atmospheric river triggered numerous debris flows within a wildfire burn scar in Big Sur – one of which destroyed California's famous Highway 1. Compared to the baseline, our burn scar simulation yields dramatic increases in total and peak discharge and shorter lags between rainfall onset and peak discharge, consistent with streamflow observations at nearby US Geological Survey (USGS) streamflow gage sites. For the 404 catchments located in the simulated burn scar area, median catchment-area-normalized peak discharge increases by ~ 450 % compared to the baseline. Catchments with anomalously high catchment-area-normalized peak discharge correspond well with post-event field-based and remotely sensed debris flow observations. We suggest that our regional postfire debris flow susceptibility analysis demonstrates WRF-Hydro as a compelling new physics-based tool whose utility could be further extended via coupling to sediment erosion and transport models and/or ensemble-based operational weather forecasts. Given the high-fidelity performance of our augmented version of WRF-Hydro, as well as its potential usage in probabilistic hazard forecasts, we argue for its continued development and application in postfire hydrologic and natural hazard assessments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Role of Snowmelt Temporal Pattern in Flood Estimation for a Small Snow‐Dominated Basin in the Sierra Nevada

Abstract Prior research confirmed the substantial bias from using precipitation‐based intensity‐duration‐frequency curves (PREC‐IDF) in design flood estimates and proposed next‐generation IDF curves (NG‐IDF) that represent both rainfall and snow processes in runoff generation. This study improves the NG‐IDF technology for a snow‐dominated test basin in the Sierra Nevada. A well‐validated physics‐based hydrologic model, the Distributed Hydrology Soil Vegetation Model (DHSVM), is used to continuously simulate snowmelt and streamflow that are used as benchmark data sets to systematically assess the NG‐IDF technology. We find that, for the studied small snow‐dominated basin, the use of standard rainfall hyetographs in the NG‐IDF technology leads to substantial underestimation of design floods. Thus, we propose probabilistic hyetographs that can represent unique patterns of events with different underlying mechanisms. For the test basin where flooding events are generated entirely by snowmelt, we develop a hyetograph that characterizes snowmelt temporal patterns, which greatly improves the performance of NG‐IDF technology in design flood estimates. In contrast to the standard rainfall hyetographs characterized by a symmetrically peaked, bell‐shaped curve, the snowmelt hyetograph displays a more rapid rise (i.e., greater intensity) and a distinct diurnal pattern influenced by solar energy input. The results also show that the uncertainty of hyetography plays an important role in design flood estimation and can have important implications for future flood projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Topographic hydro-conditioning to resolve surface depression storage and ponding in a fully distributed hydrologic model

Land surface depressions play a central role in the transformation of rainfall to ponding, infiltration and runoff, yet digital elevation models (DEMs) used by spatially distributed hydrologic models that resolve land surface processes rarely capture land surface depressions at spatial scales relevant to this transformation. Methods to generate DEMs through processing of remote sensing data, such as optical and light detection and ranging (LiDAR) have favored surfaces without depressions to avoid adverse slopes that are problematic for many hydrologic routing methods. Here, in this study, we present a new topographic conditioning workflow, Depression-Preserved DEM Processing (D2P) algorithm, which is designed to preserve physically meaningful surface depressions for depression-integrated and efficient hydrologic modeling. D2P includes several features: (1) an adaptive screening interval for delineation of depressions, (2) the ability to filter out anthropogenic land surface features (e.g., bridges), (3) the ability to blend river smoothing (e.g., a general downslope profile) and depression resolving functionality. From a case study in the Goodwin Creek Experimental Watershed, D2P successfully resolved 86% of the ponds at a DEM resolution of 10 m. Topographic conditioning was achieved with minimum impact as D2P reduced the number of modified cells from the original DEM by 51% compared to a conventional algorithm. Furthermore, hydrologic simulation using a D2P processed DEM resulted in a more robust characterization on surface water dynamics based on higher surface water storage as well as an attenuated and delayed peak streamflow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Predicting river turbidity in Pine Island Bayou using machine learning techniques coupled with variational mode decomposition

Elevated turbidity levels pose significant public health risks by facilitating the transport of harmful pollutants, including metals, organic compounds, and pathogenic microorganisms into the surface water. These conditions create serious challenges for public recreational water use and drinking water treatment, leading to economic losses and health risks. This study utilizes water monitoring data in Pine Island Bayou, Texas, and develops a Sequence-to-Sequence (S2S) model to predict turbidity using Attention-based Gated Recurrent Units with Encoder-Decoder (AT-GRU-ED) and Long Short-Term Memory (LSTM), coupled with Variational Mode Decomposition (VMD). Compared to the model without VMD, the model demonstrates satisfactory 72-hour turbidity prediction performance, achieving MAEs of 2.60 and 3.29 NTU (reductions of 53% and 58%), RMSEs of 21.08 and 31.49 NTU (reductions of 82% and 80%), and R² values of 0.96 and 0.84 on the validation and test sets, respectively. Feature importance analysis reveals that water temperature is the dominant factor influencing seasonal turbidity patterns, while real-time hourly rainfall significantly contributes to short-term variability. Turbidity typically peaks within 48 hours after rainfall events due to lagged effects from surface runoff and upstream flow. Findings suggest suspending recreational water use and water supply pumping for three days after heavy rainfall can benefit public health and improve water treatment processes. Discharges above 100 m3/s are found to accelerate sediment dilution and transport, reducing turbidity levels more quickly after the peak. In conclusion, the proposed model demonstrates reliable 72-hour turbidity prediction, supporting decision-making for water treatment plant operations and providing early warning for public recreational water use.

Deep learning↗
DOE OSTI2023

Coupled effects of future rainfall and land use on urban stormwater drainage system in Tampa, Florida (USA)

Future rainfall and land use are two important factors for flood management since they both directly and indirectly affect the functionality and performance of urban stormwater drainage infrastructures. Although some studies have already examined the coupled effects of future rainfall and land use conditions, they have not concluded how such effects on both surface runoff and flood areas would change with regional surface elevation variation and infrastructure conditions. This paper analyzed the coupled effects of future rainfall and land use on urban drainage systems in terms of surface runoff quantity and flood area changes using EPA SWMM. Future downscaled and bias-corrected precipitation projections for 2040–2060 and 2080–2099 from “Downscaled CMIP3 and CMIP5 Climate and Hydrology Projections” Archive and future land use for 2050 s and 2090 s from EPA ICLUS project were used for the City of Tampa in Florida (USA). It was shown that potential runoff volume and flood area changes ranged from –40% to 160%, and –40% to 400% due to rainfall change, while potential changes induced by land use change ranged from 0% to 3.5%, and 0% to 18% respectively. Additionally, this study examined the coupled effects of future rainfall and land use changes on flood area variations, considering the capacity of drainage infrastructure and elevation variation of study site. Results showed that flood area variations are not solely influenced by runoff quantity changes due to the two external drivers, but also by other factors. Specifically, in low-relief areas with inadequate drainage infrastructure and less elevation variation, these two drivers have no additive effects on the percentage of flood area changes. In contrast, in areas with adequate drainage infrastructure and greater surface elevation variation, there are additive effects from rainfall and land use changes on flood area expansions. These findings have important implications for stakeholders involved in city planning and flooding control, particularly for urban areas facing future precipitation and urbanization challenges.

54 ENVIRONMENTAL SCIENCES↗