USGS ScienceSearch

SEARCH · USGS Science

Results for “Geochronology, Thermochronology and Time Scale Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

1,663 records · Page 88Linked to original sources

Incorporating location uncertainty improves inference with stop-level North American Breeding Bird Survey data

Ecological models should account for uncertainty to be most effective and useful. Yet, uncertainty from model covariates—unlike that from other sources, such as sampling error or process variability—is seldom explicitly incorporated. This can cause underestimates of uncertainty to cascade through model parameter estimates, predictions, and downstream uses. Burner et al. proposed a method for quantifying uncertainty in covariates and incorporating it into models using informative Bayesian priors. This method was applied to stop-level Breeding Bird Survey (BBS) analyses, where land cover uncertainty at each stop arises from substantial stop location uncertainty. A limited validation of model-estimated land cover, using stops with known locations, indicated the method’s potential effectiveness, but it was not rigorously evaluated. We conduct a robust simulation-based test, generating stop locations, extracting land cover, and simulating bird communities across 210 BBS routes in the upper Midwest. We compare 3 models: a “known” model with true land cover, a “naive” model assuming consistent 800-m stop spacing, and a “full” model using informative priors to estimate land cover. Species parameter estimates and predicted prevalence patterns across gradients in land cover from the full model approached those of the known model and were substantially closer to the true values used in simulations relative to those from the naive model. Naive model parameters were more biased relative to the other models, and credible intervals of predicted species prevalence rarely included the true simulated values. The full model also produced land cover covariate estimates closer to true simulation values relative to the mean informative priors. Our results show that, for the BBS, informative priors enable more accurate stop-level analyses despite location uncertainty. In contrast, naive models that ignore this uncertainty yield poor inferences. More broadly, we demonstrate empirically the utility of informative priors to account for covariate uncertainty in ecological models.

Michigan, Minnesota, Wisconson

Stratigraphic notes—Volume 1, 2022

This is the first volume in the U.S. Geological Survey (USGS) series of reports on stratigraphy entitled “Stratigraphic Notes,” which consists of short papers that highlight stratigraphic studies, changes in stratigraphic nomenclature, and explanations of stratigraphic names and concepts used on published geologic maps. “Stratigraphic Notes” is a long-term (multiyear), multivolume publication containing articles that address updates or revisions to stratigraphic nomenclature (and whose content ultimately will be incorporated by National Geologic Map Database personnel into Geolex, https://ngmdb.usgs.gov/Geolex/ ). We welcome papers for the “Stratigraphic Notes” series from geoscientists of the USGS, of State Geological Surveys, and from academicians. Papers can be submitted for publication in “Stratigraphic Notes” by contacting the USGS Geologic Names Committee ( gnc@usgs.gov ). As new “Stratigraphic Notes” volumes are published, links to the volumes will be posted at https://doi.org/10.3133/pp1879 . This first volume ("Stratigraphic notes—Volume 1, 2022") includes articles that provide guidance for those who wish to submit papers to “Stratigraphic Notes,” as well as information on how to make your manuscripts compliant for geologic names reviews and how to organize your paper’s content to facilitate inclusion of new or revised names in Geolex. This volume also includes some specific guidance on conducting geologic names reviews of geologic and hydrogeologic reports.

Professional Paper

Potential for leakage among principal aquifers in the Memphis area, Tennessee

The principal aquifers in the Memphis area consist primarily of sand or sand and gravel, and the confining beds consist of clay, silt, sand, and lignite. The Jackson Formation and upper part of the Claiborne Group serve as the confining bed separating the water table aquifers from the Memphis Sand; the Flour Island Formation separates the Memphis Sand from the Fort Pillow Sand. Differences in total hydraulic head among the principal aquifers in the Memphis urban area result in vertical hydraulic gradients which create a potential for inter-aquifer exchange of water. Throughout this area, the gradient is downward from the water table aquifers to the Memphis Sand. In the central part of the Memphis urban area, the vertical hydraulic gradient is upward from the Fort Pillow Sand to the Memphis Sand, and in the eastern and western parts, it is downward from the Memphis Sand to the Fort Pillow Sand. The vertical distribution of carbon-14 data for water from the fluvial deposits, Memphis Sand, and Fort Pillow Sand shows an increase in the relative age of the water with depth. The areal distribution of carbon-14 data for water from the upper part of the Memphis Sand indicates that relatively recent water has been brought into the major cone of depression in the potentiometric surface of the Memphis Sand, either by horizontal movement or from downward vertical leakage. The normal, near-surface geothermal gradient in the Memphis area was determined to be 0.6 C/100 ft. Deviations from the normal geothermal gradient, in areas affected by intense pumping from the Memphis Sand, indicate that downward vertical leakage occurs from the water table aquifers through the Jackson-upper Claiborne confining bed to the Memphis Sand. The velocity of downward vertical leakage of water from the Memphis Sand through the Flour Island confining bed to the Fort Pillow Sand was determined to be 0.0066 ft/day by analysis of borehole temperature data from an observation well in the northeastern part of the Memphis area. From this velocity and the head difference between the Memphis Sand and the Fort Pillow Sand at this locality, the hydraulic conductivity of the Flour Island confining bed was determined to be 0.00114 ft/day. (Lantz-PTT)

Tennessee

Analysis of factors affecting plume remediation in a sole-source aquifer system, southeastern Nassau County, New York

Several plumes of dissolved, chlorinated solvents, including trichloroethylene, have been identified in a sole-source aquifer near the former Northrop Grumman Bethpage Facility and Naval Weapons Industrial Reserve Plant sites in southeastern Nassau County, New York. Past investigations have documented that the groundwater contamination originated from this industrial area and now extends to the south, in the direction of groundwater flow. The intermixed plumes are commonly referred to as the “Navy Grumman groundwater plume.” Detailed groundwater-flow modeling was needed for the New York State Department of Environmental Conservation (NYSDEC) to evaluate design options necessary for the construction, operation, optimization, maintenance, and monitoring of a groundwater extraction and treatment cleanup plan selected in a December 2019 Amended Record of Decision by the NYSDEC to comprehensively address these plumes. Consequently, the NYSDEC began a cooperative study with the U.S. Geological Survey in 2020 to better understand the local hydrogeologic framework using two independent approaches to characterize aquifer heterogeneity and update an existing regional groundwater-flow model to provide transient boundary conditions for new inset groundwater-flow models of the plume area. We developed these detailed inset models for the two independent aquifer characterizations using history-matching techniques coupled with a novel approach to risk-based management optimization of the remedial design. We also used the updated regional model to assess this optimized groundwater extraction and treatment design for potential saltwater intrusion. The ensembles of parameters resulting from history matching provided a platform with which to evaluate capture by water-supply and remedial wells using particle-tracking techniques. Using the ensemble to select a risk stance, we performed multiobjective optimization to identify various configurations of remedial pumping that are consistent with external constraints and that favor potentially competing objectives. Multiple solutions provide tradeoffs that NYSDEC can consider. In general, pumping redistribution may help to prevent further contamination migration downgradient. These and other study results are intended to support decisions for the remedial design focused on the local area encompassing the full extent of the Navy Grumman groundwater plume.

New York

Predictable seismic cycles result from structural rupture barriers on oceanic transform faults

Earthquakes of magnitude ( M ) >5.5 on oceanic transform faults (OTFs) repeatedly rupture the same locked patches, sometimes quasiperiodically. These patches are separated by “barriers” that halt earthquake propagation and slip mostly aseismically. However, the physical processes governing this systematic behavior remain unclear. We analyzed two barriers along the Gofar transform fault that have arrested ~15 M 6 earthquakes over the past three decades. Ocean bottom seismometer data indicate that the barriers hosted intense microseismicity before the mainshocks and comprise multistrand faults and transtensional stepovers with 100- to 400-m lateral offset. These characteristics contradict earthquake rupture termination models invoking velocity-strengthening friction or large geometric steps and instead point to damage-enhanced porosity and dilatancy-strengthening mechanisms. By isolating rupture segments, the barriers regulate the quasiperiodic recurrence of OTF earthquakes.

Science

Characterizing sedimentary organic carbon in a hydrothermal spreading center, the Escanaba Trough

Sediments in critical marine mineral environments are of wide importance due to their preservation of both marine minerals and organic carbon (OC) stocks. However, OC storage and cycling is often overlooked in mineral system studies. This work characterizes sedimentary OC within the Escanaba Trough, a hydrothermal sulfide system off the coast of northern California. By utilizing ROV-based push coring, we collected sediments near and far from high temperature, low temperature, and inactive vents. We applied a multipronged organic geochemical approach, measuring bulk sediment, OC fractions of varying labilities, and biomarkers to tease apart the storage, source, and cycling of OC within this complex system. In contrast to past work indicating a primarily terrestrial source to deeper, Pleistocene sediments, our results suggest a primarily marine source in surface sediments at least 50 m away from active venting. Near active venting, we see evidence of locally produced chemosynthetic OC in addition to the background marine source. This chemosynthetic OC appears to be rapidly remineralized and supports more complex deep-sea food webs compared to hemipelagic sites. Still, the greatest contribution of labile OC was observed at inactive vent sites, which we suggest is due to hydrothermal fluid alteration of more labile OC at actively venting sites. Depending on fluid temperature, hydrothermally altered OC is either preserved in the sediments as petroleum or migrated with fluid flow. Considering inactive sites are those most likely to be targeted by potential mining, we suggest additional studies could verify if these results apply at other sedimented seafloor massive sulfide systems.

California

Thermal detectability of subsurface water ice on Mars: A comparative analysis for the Subsurface Water Ice Mapping (SWIM) Project

We have developed a new global map of shallow ground ice distribution, SWIM23, based on Mars Global Surveyor Thermal Emission Spectrometer data and made systematic comparisons between this new map and two similar, previously developed data products. We have explored the origin of differences between the three ice maps by detailing technical and procedural differences in their development, by making global pixel-by-pixel comparisons, and by carrying out a series of one-dimensional thermal simulations to explore fundamental physical limitations of thermal ice-detection techniques. These efforts and the production of a composite thermal ice-consistency map supported integration of multiple geophysical data products relevant to ice detection in the upper meter of the Martian regolith into the larger Mars Subsurface Water Ice Mapping project. Our work also highlights fundamental physical limitations to thermal ice detection as a technique, particularly the rapid fall-off in ice detection sensitivity at depths >30 cm, which produces maximum uncertainty in the presence and depth of ice within regions preferred for potential human landing sites. A future Mars orbiter mission designed to detect ice and support crewed landing site selection in the midlatitude region should give payload priority to an instrument capable of probing the 1–5 m depth range (i.e., a high-frequency radar), over a next-generation thermal spectrometer, which is unlikely to offer clarity on ice table depths or lateral continuity of the ice table in the locations of highest interest.

Planetary Science Journal

Subsurface water ice mapping on Mars: A probabilistic approach

Subsurface water ice deposits on Mars are an important resource for potential future human exploration. They are also an indicator of the planet’s past climate. However, the distribution of subsurface water ice in Mars’s midlatitudes is uncertain because spacecraft imagery cannot directly observe subsurface ice in most cases. Various spacecraft remote sensing instruments are sensitive to subsurface water ice, including thermal imaging spectrometers, radar sounders, and neutron spectrometers. Geomorphic analyses of images can also implicate subsurface ice. Building upon the data products from the Mars Subsurface Water Ice Mapping project, we provide a probabilistic framework to jointly interpret existing data and estimate the likelihood of subsurface water ice in the Martian midlatitudes between 60 ∘ S and 60 ∘ N with uncertainty. Broadly, we find that near-surface ice is likely present poleward of ∼45 ∘ in both the northern and southern hemispheres. However, closer to the equator, existing remote sensing data cannot uniquely constrain the presence of subsurface water ice. Our probabilistic results provide a framework for quantifying the abundance of ice on Mars, and our uncertainty estimates allow future analysis and exploration to target regions of high uncertainty.

Planetary Science Journal

The distribution of bromide in water in the Floridan aquifer system, Duval County, northeastern Florida

Although Duval County, Florida, has ample ground-water resources for public supply, the potential exists for a problem with excessive disinfectant by-products. These disinfectant by-products result from the treatment of raw water containing low concentrations of bromide and naturally occurring organic compounds. Because of this potential problem, the relation of bromide concentrations to aquifer tapped, well location and depth, and chemical characteristics of water in the Floridan aquifer system underlying Duval County were studied to determine if these relations could be applied to delineate water with low-bromide concentrations for future supplies. In 1992, water samples from 106 wells that tap the Floridan aquifer system were analyzed for bromide and major dissolved constituents. A comparison of bromide concentrations from the 1992 sampling with data from earlier studies (1979-80) indicates that higher bromide concentrations were detected during the earlier studies. The difference between the old and new data is probably because of a change in analytical methodology in the analysis of samples. Bromide concentrations exceeded the detection limit (0.10 milligrams per liter) in water from 28 of the 106 wells (26 percent) sampled in 1992. The maximum concentration was 0.56 milligrams per liter. There were no relations between bromide and major dissolved constituents, well depth, or aquifer tapped that would be useful for determining bromide concentrations. Areal patterns of bromide occurrence are not clearly defined, but areas with relatively high bromide concentrations tend to be located in a triangular area near the community of Sunbeam, Florida, and along the St. Johns River throughout Duval County.

Florida

Prediction of regional broadband strong ground motions using a teleseismic source model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, earthquake

To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( ⁠M w ⁠ ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.

Papanoa

Updating the crustal fault model for the 2023 National Seismic Hazard Model for Alaska

We present the crustal fault model for Alaska, based on geologic observations, as a primary input for the 2023 revision of the U.S. Geological Survey National Seismic Hazard Model. We update the 2013 Alaska Quaternary fault and fold database to produce a simplified model of 105 fault sections and four fault zone polygons with basic geologic parameters including slip sense and rate. Significant updates include the following: (1) a slip rate of ∼53 mm/year on the Queen Charlotte Fault indicating it accommodates all of the plate boundary motion; (2) quantified slip rates on megathrust splay faults in the southern Prince William Sound region and near Kodiak Island; (3) improved details of structures in the Chugach-St. Elias orogen; (4) revision of the Castle Mountain Fault from right-lateral slip to a predominantly reverse fault; (5) improved Interior Alaska tectonic models that clarify relationships between the Denali, Totschunda, and thrust faults on both sides of the Alaska Range; (6) identified large earthquake sources in the eastern Brooks Range; and (7) omission of the Chatham Strait section of the Denali Fault. The fault model underscores that the collision of the Yakutat microplate is the dominant driver of active crustal faulting in most of Alaska.

Alaska, Yukon

Mapping bedrock outcrops in the Sierra Nevada Mountains (California, USA) using machine learning

Accurate, high-resolution maps of bedrock outcrops can be valuable for applications such as models of land–atmosphere interactions, mineral assessments, ecosystem mapping, and hazard mapping. The increasing availability of high-resolution imagery can be coupled with machine learning techniques to improve regional bedrock outcrop maps. In the United States, the existing 30 m U.S. Geological Survey (USGS) National Land Cover Database (NLCD) tends to misestimate extents of barren land, which includes bedrock outcrops. This impacts many calculations beyond bedrock mapping, including soil carbon storage, hydrologic modeling, and erosion susceptibility. Here, we tested if a machine learning (ML) model could more accurately map exposed bedrock than NLCD across the entire Sierra Nevada Mountains (California, USA). The ML model was trained to identify pixels that are likely bedrock from 0.6 m imagery from the National Agriculture Imagery Program (NAIP). First, we labeled exposed bedrock at twenty sites covering more than 83 km 2 (0.13%) of the Sierra Nevada region. These labels were then used to train and test the model, which gave 83% precision and 78% recall, with a 90% overall accuracy of correctly predicting bedrock. We used the trained model to map bedrock outcrops across the entire Sierra Nevada region and compared the ML map with the NLCD map. At the twenty labeled sites, we found the NLCD barren land class, even though it includes more than just bedrock outcrops, accounted for only 41% and 40% of mapped bedrock from our labels and ML predictions, respectively. This substantial difference illustrates that ML bedrock models can have a role in improving land-cover maps, like NLCD, for a range of science applications.

California

Landslide-channel feedbacks amplify channel widening during floods

Channel widening is a major hazard during floods, particularly in confined mountainous catchments. However, channel widening during floods is not well understood and not always explained by hydraulic variables alone. Floods in mountainous regions often coincide with landslides triggered by heavy rainfall, yet landslide-channel interactions during a flood event are not well known or documented. Here we demonstrate with an example from the Great Colorado Flood in 2013, a 1000 year precipitation event, how landslide-channel feedbacks can substantially amplify channel widening and flood risk. We use a combination of DEM differencing, field analysis, and multiphase flow modeling to document landslide-channel interaction during the flood event in which sediment delivered by landslides temporarily dammed the channel before failing and generating substantial channel widening. We propose that such landslide-flood interactions will become increasingly important to account for in flood hazard assessment as flooding and landsliding both increase with extreme rainfall under climate change.

Colorado

Cotton farming pesticides affect ileal microbiota activity expressions of virulome but not resistome or metabolic pathways in a sedentary wild passerine

The increased use of agrochemicals to enhance crop production has had detrimental environmental effects including implication in the sharp decline of North American farmland-breeding birds. Here, using a combination of deep shotgun metatranscriptomics and pesticide exposure data, we sought to assess whether exposure to cotton ( Gossypium spp.) production had a differential effect on ileum multi-kingdom microbial activity, metabolism, anti-microbial resistance, and virulence factors of sedentary northern mockingbirds ( Mimus polyglottos ) sampled from two cotton-producing areas (16 birds in total) and one uncultivated area (7 birds) in Texas, USA. Both Shannon Index values (Adj. r 2 = 0.174, F (1,21) = 5.633, p = 0.027) and a Mantel test (Spearman ρ = 0.184, p = 0.013) supported a relationship between metabolically active microbiota Bray–Curtis dissimilarities and differences in pesticide mixtures among study areas. Virulence factor richness (Adj. r 2 = 0.182, F (1,21) = 5.890, p = 0.024), Shannon Index (Adj. r 2 = 0.231, F (1,21) = 7.612, p = 0.012), and load (sum of virulence factor abundances; Adj. r 2 = 0.160, F (1,21) = 5.194, p = 0.033) were related to total pesticide load (total quantity of pesticides). We found no pesticide effects on expression of either antimicrobial resistance genes or metabolic pathways.

Texas

Accurate simulation of flow through dipping aquifers with MODFLOW 6 using enhanced cell connectivity

In simulations of groundwater flow through dipping aquifers, layers of model cells are often “deformed” to follow the top and bottom elevations of the aquifers. When this approach is used in MODFLOW, adjacent cells within the same model layer are vertically offset from one another, and the standard conductance-based (two-point) formulation for flow between cells does not rigorously account for these offsets. The XT3D multi-point flow formulation in MODFLOW 6 is designed to account for geometric irregularities in the grid, including vertical offsets, and to provide accurate results for both isotropic and anisotropic groundwater flow. A recent study evaluated the performance of the standard formulation and XT3D using a simple, synthetic benchmark model of a steeply dipping aquifer. Although XT3D generally improved the accuracy of flow simulations relative to the standard formulation as expected, neither formulation produced accurate flows in cases that involved large vertical offsets. In this paper, we explain that the inability of XT3D to produce accurate flows in the steeply dipping aquifer benchmark was not due to an inherent limitation of the flow formulation, but rather to the limited cell connectivity inherent in the most commonly used discretization packages in MODFLOW 6. Furthermore, we demonstrate that XT3D is able to produce the expected accuracy when adequate cell connectivity is introduced using MODFLOW's unstructured grid type and the aquifer is discretized vertically using at least two model layers.

Groundwater

From model to action: Identifying the gaps in coastal marsh models for decision making

Coastal salt marsh modeling projects are designed to answer questions about where salt marshes may migrate, how they can migrate, and how they will be impacted by changing water levels. While this information is vital for climate-resilient marsh conservation, translating these findings into actionable results remains a challenge. The authors examined the actionability of marsh modeling efforts in the United States by analyzing data collected from semi-structured interviews of marsh model users ( n = 24) across two ongoing projects along the U.S. East Coast and the Gulf of America (Gulf of Mexico). By qualitatively analyzing the interview data, the authors found that the tasks of practitioners who use coastal marsh model outputs fall into three major themes: (1) marsh restoration, (2) planning with uncertainty, and (3) conserving habitat for marsh-reliant species. For each of these themes, the authors identify unmet needs, including high spatial resolution information for local planning, accessible descriptions of uncertainty to increase user confidence, and the incorporation of human dimensions data (e.g., human alterations to the landscape such as impoundments and culverts) for a comprehensive understanding of the coastal salt marsh. Marsh modeling projects should strategize how to fulfil these unmet user needs so that marsh modeling outputs can better support decisions related to marsh conservation and restoration.

Gulf of America, United States East Coast

A crosswalk of the 2015 World Terrestrial Ecosystems to the International Union for the Conservation of Nature Global Ecosystem Typology Framework

To support ecosystem mapping and accounting applications, we aligned the 2015 U.S. Geological Survey/Esri/The Nature Conservancy-World Terrestrial Ecosystems (WTEs) with the International Union for Conservation of Nature Global Ecosystem Typology (GET) framework. This process, known as “crosswalking,” enabled the development of a global map of GET level 3 Ecosystem Functional Groups (EFGs) at a 250-meter spatial resolution. Crosswalking involved manually assigning 1,781 biogeographically stratified WTEs to their most probable EFG based on similarities in climate, terrain, vegetation, and geographic distribution. We compared attributes of the WTE dataset with summary characteristics of the EFGs. The resulting crosswalked global map of International Union for Conservation of Nature GET ecosystems is intended to be useful for standardizing ecosystem classification and reporting under frameworks such as the Kunming-Montreal Global Biodiversity Framework and the United Nations System of Environmental-Economic Accounting. We discuss key challenges in reconciling non-identical classifications, such as many-to-one relationships and variation in data quality.

Open-File Report

Representation of surface-water flows using Gradient-Related Discharge in an Everglades Network

The Everglades Depth Estimation Network interpolates water-level gage data to produce daily water-level elevations for the Everglades in south Florida. These elevations were used to estimate flow vectors (gradients and directions) and volumetric flow rates using the Gradient-Related Discharge in an Everglades Network (GARDEN) application developed by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers. Flow rates in both the east-west and north-south directions were computed on a 400-meter square grid using modified parameters in the Manning’s equation. The frictional resistance parameter in the Manning’s equation was calibrated to measured flow rates at coastal creeks fed by Everglades Depth Estimation Network boundary flows. Levees and other features that act as barriers to flow were defined as “no-flow” grid cells where vectors were set to zero. The flow volume magnitudes were calibrated with 2020 daily values of coastal river flows, and verification was performed using 2021 data. Within a given day, the measured coastal river flows fluctuate more than the GARDEN boundary flows because of tidal and wind forcings. Because the GARDEN boundary flows were the upstream water source for the coastal rivers, calibration focused on matching average daily flow volumes rather than daily fluctuations. The Pearson’s correlation coefficient is 0.766 for the 2020 calibration period and 0.566 for the 2021 verification period. Applying GARDEN to periods with hydraulic-control-structure releases allows the propagation of structure flows to be seen in the daily flow-vector maps along with the multiday response of flows farther downgradient. Flow vectors may be overestimated near control structures because of difficulties in resolving the water gradient downstream from the structure. Flow vectors farther from the structure are more accurate than those near the structure.

Florida