Search NASA⌕ Search

SEARCH · Search NASA

Results for “2D flood model”

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

TRITON: A Multi-GPU open source 2D hydrodynamic flood model

A new open source multi-GPU 2D flood model called TRITON is presented in this work. The model solves the 2D shallow water equations with source terms using a time-explicit first order upwind scheme based on an Augmented Roe's solver that incorporates a careful estimation of bed strengths and a local implicit formulation of friction terms. Here, the scheme is demonstrated to be first order accurate, robust and able to solve for flows under various conditions. TRITON is implemented such that the model effectively utilizes heterogeneous architectures, from single to multiple CPUs and GPUs. Different test cases are shown to illustrate the capabilities and performance of the model, showing promising runtimes for large spatial and temporal scales when leveraging the computer power of GPUs. Under this hardware configuration, communication and input/output subroutines may impact the scalability. The code is developed under an open source license and can be freely downloaded in https://code.ornl.gov/hydro/triton.

2D flood model↗

Extremely rapid, Lagrangian modeling of 2D flooding: A rivulet-based approach

Estimates of potential flood inundation areas and depths are critical to informing the preparedness, response, and investment decisions of many government agencies and private sector organizations, especially under a changing climate. The standard modeling approaches, however, are often either computationally intensive or constrained in their accuracy or applicability. A novel, rivulet-based, 2D model of pluvial flooding is described in this article that is 10,000 to 10 million times less computationally complex than the full solution of the shallow water equations, yet achieves inundation area hit rates of between 0.8 and 0.9 and relative absolute mean errors of 10%-20% across a wide range of flow depths. This combination of accuracy and efficiency will enable real-time depth estimates during flood events, detailed sensitivity analyses, and the generation of large ensembles to support broad uncertainty analyses.

54 ENVIRONMENTAL SCIENCES↗

Enhancing 2D hydrodynamic flood models through machine learning and urban drainage integration

Two-dimensional hydrodynamic flood models are commonly employed for simulating flood extent and inundation depth. However, the influence of urban drainage network (UDN) is frequently overlooked in these models, potentially compromising their accuracy. Furthermore, the expensive computational costs and longer processing times make them challenging for large-scale hydrodynamic simulation. To address these challenges, this paper develops a machine learning (ML)-driven emulator for an open-source flood model, the Two-dimensional Runoff Inundation Toolkit for Operational Needs (TRITON). A TRITON-ML Emulator (TR-Emulator) that utilizes Convolutional Long Short-Term Memory is developed to capture the spatiotemporal features of flood events based on the outputs from TRITON. We further enhance the emulator by integrating UDN parameters (TR-UDN), such as the flow capacity of drainage pipes, pipe size, and pipe length, via an ML stacking technique to improve the water surface elevation (WSE) simulation. Hurricane Harvey 2017 in Houston, TX is used as the case study. We compare WSE results from TRITON, TR-Emulator, TR-UDN, and the United States Geological Survey (USGS) observations to evaluate the performance of these models. The results indicate that the TR-Emulator effectively replicates the WSE simulated by TRITON. Additionally, TR-UDN performs well in capturing WSE patterns and peak flows, aligning more closely with USGS observations, except in areas with milder slopes where conveyance discrepancies are observed. We further test the generalizability of our ML-based models using another smaller event. This paper shows that the TR-Emulator is effective for users and engineers to emulate a 2D hydrodynamic model, and the enhanced version of the TR-Emulator, TR-UDN, can be an efficient tool for predicting WSEs during urban flooding.

54 ENVIRONMENTAL SCIENCES↗

Integrated Analysis of Interferometric SAR, Satellite Altimetry and Hydraulic Modeling to Quantify Louisiana Wetland Dynamics

Wetland loss in Louisiana has been accelerating due primarily to anthropogenic and nature processes, and is being advocated as a problem with national importance. Accurate measurement or modeling of wetland-wide water level changes, its varying extent, its storage and discharge changes resulting in part from sediment loads, erosion and subsidence are fundamental to assessment of hurricane-induced flood hazards and wetland ecology. Here, we use innovative method to integrate interferometric SAR (InSAR) and satellite radar altimetry for measuring absolute or geocentric water level changes and applied the methodology to remote areas of swamp forest in coastal Louisiana. Coherence analysis of InSAR pairs suggested that the HH polarization is preferred for this type of observation, and polarimetric analysis can help to identi:fy double-bonnce backscattering areas in the wetland. Envisat radar altimeter-measured 18- Hz (along-track sampling of 417 m) water level data processed with regional stackfile method have been used to provide vertical references for water bodies separated by levees. The high-resolution (approx.40 m) relative water changes measured from ALOS PALSAR L-band and Radarsat-l C-band InSAR are then integrated with Envisat radar altimetry to obtain absolute water level. The resulting water level time series were validated with in situ gauge observations within the swamp forest. Furthermore, we compare our water elevation changes with 2D flood modeling from LISFLOOD hydrodynamic model. Our study demonstrates that this new technique allows retrospective reconstruction and concurrent monitoring of water conditions and flow dynamics in wetlands, especially those lacking gauge networks.

Lee, Hyongki↗

Improved Calibration of Modeled Discharge and Storage Change in the Atchafalaya Floodplain Using SAR Interferometry

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.

Jung, Hahn Chul↗

A coupled hydrodynamic (HEC-RAS 2D) and water quality model (WASP) for simulating flood-induced soil, sediment, and contaminant transport

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.

54 ENVIRONMENTAL SCIENCES↗

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↗

Geological study of an outburst flood event in the upper Yangtze River and risk of similar extreme events

Here, outburst flood sediments are studied on the upper course of the Yangtze River (known as the Jinsha River) where is a key region for the development of hydroelectric power. As the history and magnitude of outburst flooding in this stretch of the Jinsha River during historic times has remained largely undocumented, we have endeavored as part of this study to conduct field investigations into, and provide numerical dating (using OSL and 14 C dating techniques) of outburst flood evidence found in the First Bend area of the Yangtze River. Carbonized nutshells and charcoal from the ancient settlement yielded 14 C ages that precisely pinpointed an outburst flood event ~1200–1300 cal yr BP, during China's Tang Dynasty (618–907 CE). However, OSL ages using the minimum age model for all outburst flood samples ranged from 9.1 to 1.8 ka, which was ~0.5–7.8 ka older than the 14 C-dated outburst flood. This discrepancy suggests a potential overestimation of quartz age when dating outburst flood sediments from the Holocene period. Furthermore, using 2D hydraulic modeling, we estimated the magnitude of the aforementioned outburst flood, revealing a reconstructed peak flow discharge of 51,000–55,300 m 3 /s. This discharge was three times greater than the 10,000-year flood for the hydroelectric power station nearest to the study area. Our study underscores the significant role of combining archaeological and geological evidence in enhancing paleoflood hydrology research. Through a multidisciplinary approach, our findings emphasize the importance of paleoflood hydrology studies in comprehending both the latent and catastrophic flood risks in alpine mountain areas amid the context of global warming.

14C dating↗

Los Alamos National Laboratory Design Basis Flood Level Modeling Report

This report describes a project conducted by Los Alamos National Laboratory (LANL) and 3AE Green to determine the base flood elevations for the return periods associated with Flood Design Categories (FDC) 1-3 as identified in DOE-STD-1020-2016. As part of this work, 3AE Green performed two-dimensional (2D) hydraulic modeling of the Design Basis Flood Level (DBFL) for the primary watersheds within the LANL boundary and site-specific sub-watersheds associated with specific existing LANL facilities. The modeling covered an area of approximately 67 square miles.

58 GEOSCIENCES↗

Evaluating the Effectiveness of Soil Profile Rehabilitation for Pluvial Flood Mitigation Through Two-Dimensional Hydrodynamic Modeling

Pluvial flooding, driven by increasingly impervious surfaces and intense storm events, presents a growing challenge for urban areas worldwide. In Baltimore City, MD, USA, climate change, rapid urbanization, and aging stormwater infrastructure are exacerbating flooding impacts, resulting in significant socio-economic consequences. This study evaluated the effectiveness of a soil profile rehabilitation scenario using a 2D hydrodynamic modeling approach for the Tiffany Run watershed, Baltimore City. This study utilized different extreme storm events, a high-resolution (1 m) LiDAR Digital Terrain Model (DTM), building footprints, and hydrological soil data. These datasets were integrated into a fully coupled 2D hydrodynamic model, the City Catchment Analysis Tool (CityCAT), to simulate urban flood dynamics. The pre-soil rehabilitation simulation revealed a maximum water depth of 3.00 m in most areas, with hydrologic soil groups C and D, especially downstream of the study area. The post-soil rehabilitation simulation was targeted at vacant lots and public parcels, accounting for 33.20% of the total area of the watershed. This resulted in a reduced water depth of 2.50 m. Additionally, the baseline runoff coefficient of 0.49 decreased to 0.47 following the rehabilitation, and the model consistently recorded a peak runoff reduction rate of 4.10 across varying rainfall intensities. The validation using a contingency matrix demonstrated true-positive rates of 0.75, 0.50, 0.64, and 0 for the selected events, confirming the model’s capability at capturing real-world flood occurrences.

Baltimore City↗

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗

Climate-extreme modeling framework for sustainable flood management in the Arabian Peninsula

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.

CMIP6↗

Hydrological Modeling of the Jezero Crater Outlet-Forming Flood

Jezero crater is a site of prime scientific interest because it was a lake early in Mars history. Preserved clay- and carbonate-bearing sedimentary fans on Jezero's western and northwestern margin (Fig. 2) are accessible to future exploration. Geologic context [1] and stratigraphic analysis of the western fan strongly support the interpretation that these fans were deposited as deltas into the lake. This has helped establish Jezero as one of the final candidate landing sites for Mars 2020. The high level of certainty that Jezero was a lake results from the existence of its outlet valley, which required filling of the crater to form [e.g., 1,4]. Here, we specifically focus on how this outlet valley was carved by the dam breach flood that eroded the eastern crater rim. We have completed preliminary modeling in both 1D and 2D of the outlet's formation.

Geology↗

Calibration of Two-dimensional Floodplain Modeling in the Atchafalaya River Basin Using SAR Interferometry

Two-dimensional (2D) satellite imagery has been increasingly employed to improve prediction of floodplain inundation models. However, most focus has been on validation of inundation extent, with little attention on the 2D spatial variations of water elevation and slope. The availability of high resolution Interferometric Synthetic Aperture Radar (InSAR) imagery offers unprecedented opportunity for quantitative validation of surface water heights and slopes derived from 2D hydrodynamic models. In this study, the LISFLOOD-ACC hydrodynamic model is applied to the central Atchafalaya River Basin, Louisiana, during high flows typical of spring floods in the Mississippi Delta region, for the purpose of demonstrating the utility of InSAR in coupled 1D/2D model calibration. Two calibration schemes focusing on Manning s roughness are compared. First, the model is calibrated in terms of water elevations at a single in situ gage during a 62 day simulation period from 1 April 2008 to 1 June 2008. 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. The best-fit models show that the mean absolute errors are 3.8 cm for a single in situ gage calibration and 5.7 cm/46 days for InSAR water level calibration. The optimum values of Manning's roughness coefficients are 0.024/0.10 for the channel/floodplain, respectively, using a single in situ gage, and 0.028/0.10 for channel/floodplain the using SAR. Based on the calibrated water elevation changes, daily storage changes within the size of approx 230 sq km of the model area are also calculated to be of the order of 107 cubic m/day during high water of the modeled period. This study demonstrates the feasibility of SAR interferometry to support 2D hydrodynamic model calibration and as a tool for improved understanding of complex floodplain hydrodynamics

Hahn Chul Jung↗

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗

Intercomparison of flood inundation models across land use types and hydrological flood stages

Flood Inundation Mapping (FIM) model selection is a key operational decision because accurate, rapid mapping underpins early warning and resource allocation. FIM performance is context-dependent and can vary with hydrograph phase, land-use/land-cover (LULC), and the evaluation benchmark. Intercomparison studies typically assess a single near-peak snapshot against one reference dataset. Here, we provide a context-stratified intercomparison across (i) multiple hydrograph phases, (ii) LULC classes, and (iii) benchmark types, for five FIM approaches spanning a wide range of physical complexity and operational cost (TRITON, LISFLOOD-FP, HEC-RAS 2D, ARC-Curve2Flood, and OWP HAND-FIM). We use the Hurricane Matthew flood (2016) in the Neuse River Basin, North Carolina, USA, as a case study. Using high-resolution remote sensing-derived flood inundation maps, hand-labeled points, and building footprints, we assess model skill across two rising and two falling hydrograph limbs and across major LULC types. Results show that model rankings shift systematically across contexts: LISFLOOD-FP ranks highest in three of four flood phases, while TRITON leads during one rising limb phase; LISFLOOD-FP performs best in vegetated areas, whereas HEC-RAS improves relative performance in agricultural and urban areas; and benchmark choice influences conclusions, with LISFLOOD-FP performing best for flooded-building detection in the late falling limb, while TRITON ranks highest against hand-labeled points. We also report representative wall-clock runtimes for each workflow to provide use-case context for operational feasibility. Together, these results offer transferable guidance for model selection and for designing large-scale, benchmark-aware FIM intercomparison studies.

Nikrou, Parvaneh [University of Alabama]↗

A Preliminary Analysis of Precipitation Properties and Processes during NASA GPM IFloodS

The Iowa Flood Studies (IFloodS) is a NASA Global Precipitation Measurement (GPM) ground measurement campaign, which took place in eastern Iowa from May 1 to June 15, 2013. The goals of the field campaign were to collect detailed measurements of surface precipitation using ground instruments and advanced weather radars while simultaneously collecting data from satellites passing overhead. Data collected by the radars and other ground instruments, such as disdrometers and rain gauges, will be used to characterize precipitation properties throughout the vertical column, including the precipitation type (e.g., rain, graupel, hail, aggregates, ice crystals), precipitation amounts (e.g., rain rate), and the size and shape of raindrops. The impact of physical processes, such as aggregation, melting, breakup and coalescence on the measured liquid and ice precipitation properties will be investigated. These ground observations will ultimately be used to improve rainfall estimates from satellites and in particular the algorithms that interpret raw data for the upcoming GPM mission's Core Observatory satellite, which launches in 2014. The various precipitation data collected will eventually be used as input to flood forecasting models in an effort to improve capabilities and test the utility and limitations of satellite precipitation data for flood forecasting. In this preliminary study, the focus will be on analysis of NASA NPOL (S‐band, polarimetric) radar (e.g., radar reflectivity, differential reflectivity, differential phase, correlation coefficient) and NASA 2D Video Disdrometers (2DVDs) measurements. Quality control and processing of the radar and disdrometer data sets will be outlined. In analyzing preliminary cases, particular emphasis will be placed on 1) documenting the evolution of the rain drop size distribution (DSD) as a function of column melting processes and 2) assessing the impact of range on ground‐based polarimetric radar estimates of DSD properties.

Carey, Lawrence↗

Reactive transport modeling for supporting climate resilience at groundwater contamination sites

Abstract. Climate resilience is an emerging issue at contaminated sites and hazardous waste sites, since projected climate shifts (e.g., increased/decreased precipitation) and extreme events (e.g., flooding, drought) could affect ongoing remediation or closure strategies. In this study, we develop a reactive transport model (Amanzi) for radionuclides (uranium, tritium, and others) and evaluate how different scenarios under climate change will influence the contaminant plume conditions and groundwater well concentrations. We demonstrate our approach using a two-dimensional (2D) reactive transport model for the Savannah River Site F-Area, including mineral reaction and sorption processes. Different recharge scenarios are considered by perturbing the infiltration rate from the base case as well as considering cap-failure and climate projection scenarios. We also evaluate the uranium and nitrate concentration ratios between scenarios and the base case to isolate the sorption effects with changing recharge rates. The modeling results indicate that the competing effects of dilution and remobilization significantly influence pH, thus changing the sorption of uranium. At the maximum concentration on the breakthrough curve, higher aqueous uranium concentration implies that sorption is reduced with lower pH due to remobilization. To better evaluate the climate change impacts in the future, we develop the workflow to include the downscaled CMIP5 (Coupled Model Intercomparison Project) climate projection data in the reactive transport model and evaluate how residual contamination evolves through 2100 under four climate Representative Concentration Pathway (RCP) scenarios. The integration of climate modeling data and hydrogeochemistry models enables us to quantify the climate change impacts, assess which impacts need to be planned for, and therefore assist climate resiliency efforts and help guide site management.

54 ENVIRONMENTAL SCIENCES↗