Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “International Journal of Disaster Risk Reduction”

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,661 records · Page 11Linked to original sources
USGS2025

Manipulation of the Symbiodiniaceae microbiome confers multigenerational impacts on symbioses and reproductive ecology of its Exaiptasia diaphana host

Symbiodiniaceae-associated microbiota strongly influence cnidarian symbioses. We systematically reduced the bacterial and fungal communities associated with Symbiodiniaceae to study potential effects on the cnidarian holobiont Exaiptasia diaphana (Aiptasia). Clonal anemones were inoculated with xenic Breviolum minutum (SSB01) and microbiome-manipulated cultures after antibacterial or antifungal treatment. The asexual reproduction of pedal laceration allowed for three generations of clonal aposymbiotic Aiptasia to be utilised in this study, from the initial adult generation (G0), to the first (G1), and second (G2) generation. We inoculated small and large G1 Aiptasia with SSB01 algae and monitored onset of symbiosis, rate of algal proliferation, and holobiont characteristics. Sequencing the 16S and 18S rRNA gene regions identified significant differences in the bacterial and fungal communities of the G0 and G1 generations, alongside differences between the size classes of small and large G1 anemones. The microbiome of larger G1 individuals was distinct to the smaller G1 anemones, suggesting a microbiome maturation process. Control Breviolum minutum cultures exhibited a significantly greater proliferation rate in large G1 anemones when compared to antibacterial or antifungal treated cultures, whereas the opposite trend was documented in the small G1 anemones. Although no differences were observed between algal photochemical parameters, or the growth and polyp activity of G1 juveniles, we observed a significant influence in the production of G2 clones between treatments. Overall, we provide strong ecological implications of manipulating Symbiodiniaceae microbiome, not for the algae themselves, but for the maturation of the host Aiptasia, as well as for the cnidarian holobiont over multiple generations.

The ISME Journal↗
USGS2025

Using hydrogeologic context and water budgets to evaluate the potential for groundwater contributions to contaminants in Lake Superior

This study presents a synthesis of the hydrogeology in the U.S. Lake Superior watershed and the contribution of groundwater to the water budget of the U.S. Lake Superior basin. Much of the shoreline of Lake Superior in Minnesota and Michigan is composed of hydrogeologic units contributing very little direct groundwater discharge to the lake. Groundwater in watersheds adjacent to the lake typically flows in short, local flow systems characterized by thin glacial sediments with active groundwater flow in fractured bedrock within the top 60–90 m below land surface. The complex groundwater system in Wisconsin’s Bayfield Peninsula has the largest groundwater reservoir near the lake, characterized by thick sand and gravel glacial deposits and underlying sandstone aquifer. Although these thick sandy glacial deposits are not in direct contact with the lake at the shoreline, groundwater discharge may still be significant via subsurface exposures beyond the shoreline or flow through the underlying sandstone aquifer. Overall, most groundwater in the watershed is contributed as indirect base flow in streams around the lakeshore and comprises about 60 percent of the upland water budget. Direct groundwater flow to the shoreline contributes 2 to 9 percent of the inflow. Identifying possible contamination sources through direct sampling of groundwater would be an inefficient way to detect problems if sources are unknown, particularly for some chemicals of concern such as PFAS, pesticides, PCBs, chloride, and nutrients. Evaluating the chemical characteristics of contamination is also important to consider in evaluating how groundwater may contribute to pollution in Lake Superior.

Lake Superior↗
USGS2025

Ecohydrological response of a forested headwater catchment to a flash drought in the Southeastern U.S.

Flash droughts differ from traditionally defined droughts in their rapidity of intensification and often associated high vapor-pressure deficit. These droughts can lead to declines in streamflow and water table depth and induce water stress to vegetation at a greater rate than droughts that manifest over longer periods. However, little is known regarding the response of forested environments to flash drought because most studies of impacts have been conducted in agricultural settings. In this study we investigated water-use patterns of riparian trees using sap flow methods and examined the role of groundwater as a source of moisture over three periods that were delimited by antecedent soil moisture conditions. For a longer-term perspective we also examine monthly streamflow over the 35-year record. We observed that trees at only one monitoring plot showed a decrease in water use relative to evaporative demand during a flash drought. Total reverse sap flow (flow toward the roots rather than the canopy) greatly increased during the flash drought period, suggesting the likely occurrence of hydraulic redistribution to the excessively dry soils. Over the drought period groundwater became a more dominant source of moisture for sustaining forest water use. Monthly mean streamflow during the flash drought approached levels observed in past multiyear droughts. This is the first study, to our knowledge, to specifically investigate the response of multiple water budget components to flash drought in a humid forest. As more studies are conducted, a better understanding of the range of expected responses are likely to emerge.

Georgia↗
USGS2026

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↗
USGS2026

Transcriptional changes in wild Yukon River Chinook Salmon associated with Ichthyophonus infections

Objective We compared differentially expressed genes in Chinook Salmon Oncorhynchus tshawytscha with three divergent Ichthyophonus statuses (undetected, subclinical infections, or clinical disease; n = 100) to investigate associated transcriptomic responses. Disease associated with the fish parasite Ichthyophonus sp. was first diagnosed in adult Chinook Salmon from the Yukon River in the late 1980s and has subsequently been implicated in premature host mortality. Methods Chinook Salmon tissue sample collections and Ichthyophonus infection data were leveraged from a multi-agency collaboration during summer 2022 at three locations along the main-stem Yukon River that spanned nearly 2,000 km of freshwater migration (lower, middle, and upper river). We sequenced the transcriptome and compared this to infection status based on routine diagnostic procedures. Results Among the 17,569 genes for which messenger RNA was detected, we identified a transcription signature in the skeletal muscle that was associated with Ichthyophonus infections and included 53 differentially expressed genes. The differentially expressed genes and their pathways included those known for involvement in immune functions, energy synthesis, cellular breakdown, and reproduction—all processes that are known to be influenced by senescence during spawning migrations. Conclusions Results demonstrate a clear transcriptional difference between diseased fish (clinical disease group) and those in which Ichthyophonus was undetected, including identifying candidate markers for infection in this population. These results provide a foundation for development of nonlethal biomarkers to evaluate potential Ichthyophonus infections in Chinook Salmon based on gene transcription, protein products, or gene variants (e.g., polymorphisms).

Journal of Aquatic Animal Health↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS1976

Trace elements in some rocks and slag from the Yemen Arab Republic and their bearing on the iron ore near Sa'dah

Thirty-seven samples of rocks, saprolite, laterite, gossan, and slag were collected in the Yemen Arab Republic during June-July 1975 when field work was conducted to refine an interpretation of the regional geology made from LANDSAT-1 images by M. J. Grolier. These samples we reanalyzed in the laboratories of the U. S. Geological Survey, and the results were interpreted for use in programs of economic development as a guide for the evaluation of mineral resources in the Yemen Arab Republic. The work was undertaken as part of the program of the U. S. Agency for International Development with the cooperation of the Mineral and Petroleum Authority, Ministry of Economy, Yemen Arab Republic.

Open-File Report↗
USGS2024

GRAPES: Earthquake early warning by passing seismic vectors through the grapevine

Estimating an earthquake's magnitude and location may not be necessary to predict shaking in real time; instead, wavefield-based approaches predict shaking with few assumptions about the seismic source. Here, we introduce GRAph Prediction of Earthquake Shaking (GRAPES), a deep learning model trained to characterize and propagate earthquake shaking across a seismic network. We show that GRAPES’ internal activations, which we call “seismic vectors”, correspond to the arrival of distinct seismic phases. GRAPES builds upon recent deep learning models applied to earthquake early warning by allowing for continuous ground motion prediction with seismic networks of all sizes. While trained on earthquakes recorded in Japan, we show that GRAPES, without modification, outperforms the ShakeAlert earthquake early warning system on the 2019 M7.1 Ridgecrest, CA earthquake.

Shimane/HiroshimaPrefectures↗
USGS2024

Subduction zone geometry modulates the megathrust earthquake cycle: Magnitude, recurrence, and variability

Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably-dipping planar and variably-curved nonplanar megathrusts using the volumetric, high-order accurate code tandem to account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly-dipping megathrusts to host larger earthquakes than steeply-dipping ones. Under elevated pore pressure and less strongly velocity-weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly-dipping and sharply-curved megathrusts host multi-period supercycles of slow-to-fast, small-to-large slip events under higher effective stresses and more strongly velocity-weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second-order effects from structurally-derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co- and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat faults.

Journal of Geophysical Research: Solid Earth↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2025

Applying U-Th disequilbrium for dating siliceous sinters

Continental hydrothermal systems are critical avenues for the crustal transport of heat and mass captured for geothermal energy and mineral exploration. Thus, understanding their temporal evolution and longevity is important for resource characterization. Deposits of microlaminated siliceous sinter, common surface expressions of high temperature reservoirs (> 170 °C), have the potential to trace hydrothermal histories. Geothermal reservoirs are often located within uranium-bearing silicic volcanic rock where subsurface fluid-rock interactions extract U into hydrothermal fluids. U incorporated in the surface sinter deposit has the potential to provide a U— Th disequilibrium dating option. We focus on samples from El Tatio geyser field in the Altiplano of northern Chile, the largest geothermal system in the Andes. Our resulting 230 Th/U ages, along with the water and deposit elemental compositions, suggest concentrations of U and Th vary predictably along the sinter apron. While distal facies containing the highest U concentrations (> 50 μg/g) are least affected by detrital Th corrections, they can display suspected open-system behavior. In contrast, more medial facies, where bacterial mats and other porous textures are co mmonly concentrated, have only trace amounts of U (< 0.1 μg/g), which leads to unreliable or geologically improbable dates. Proximal facies tend to date most consistently. By comparing existing 14 C ages with 230 Th/U results, 230 Th/U ages tend to be younger than the 14 C ages, supporting the presence of a 14 C-dead carbon influence. New data confirm that the onset of geothermal activity at El Tatio goes back to the late Pleistocene.

Altiplano, El Tatio geyser field↗
USGS2024

Improving crop-specific groundwater use estimation in the Mississippi Alluvial Plain: Implications for integrated remote sensing and machine learning approaches in data-scarce regions

Study region The Mississippi Alluvial Plain (MAP) in the United States (US). Study focus Understanding local-scale groundwater use, a critical component of the water budget, is necessary for implementing sustainable water management practices. The MAP is one of the most productive agricultural regions in the US and extracts more than 11 km 3 /year for irrigation activities. Consequently, groundwater-level declines in the MAP region pose a substantial challenge to water sustainability, and hence, we need reliable groundwater pumping monitoring solutions to manage this resource appropriately. New hydrological insights for the region We incorporate remote sensing datasets and machine learning to improve an existing lookup table-based model of groundwater use previously developed by the U.S. Geological Survey (USGS). Here, we employ Distributed Random Forests, an ensemble machine learning algorithm to predict annual and monthly groundwater use (2014–2020) throughout this region at 1-km resolution, using pumping data from existing flowmeters in the Mississippi Delta. Our model compares favorably with the existing USGS model, with higher R 2 (0.51 compared to 0.42 in the previous model), and lower root mean square error (RMSE) and mean absolute error (MAE)— 0.14 m and 0.09 m, respectively in our model, compared to 0.15 m and 0.1 m in the previous model. Therefore, this work advances our ability to predict groundwater use in regions with scarce or limited in-situ groundwater withdrawal data availability.

Journal of Hydrology Regional Studies↗
USGS2026

Geophysical architecture and geochronology of the Neoarchean Mentor anorthosite intrusive complex, northwestern Minnesota: Largest anorthosite complex of the Superior Province?

The Neoarchean Mentor anorthosite intrusive complex (MAIC) lies within the Wawa subprovince in northwestern Minnesota, in a region where the Wawa, Quetico, and Wabigoon subprovinces are juxtaposed in close proximity. Archean rocks are entirely concealed, and interpretations are developed from aeromagnetic, gravity, and borehole samples. The MAIC includes both anorthosite (dense, weakly magnetized) and oxide-rich gabbro (strongly magnetized) lithologies. Anorthosite is coarse-grained to megacrystic, intensely altered, and locally brecciated. Pervasive epidote alteration enhanced the density of the anorthosite via introduction of higher-density mineral assemblages, explaining why the MAIC produces a significant gravity high. Oxide-rich gabbro forms a border phase of the MAIC and has potential for vanadium, chromium, and titanium mineralization, and produces a strong aeromagnetic high. The MAIC is interpreted to extend over an area of 640 km 2 , making it the largest known anorthosite complex of the Superior Province, as measured by preserved areal extent. Modeling indicates the MAIC extends more than 6 km into the subsurface. A new Pb–Pb zircon age of 2737.2 ± 4.5 Ma is interpretated as the crystallization age of anorthosite within the MAIC, showing that the MAIC formed well before the ca. 2690 Ma Shebandowanian orogeny, and raising new questions about correlations with other parts of the Wawa subprovince. A low-density batholith, here informally called the Fertile batholith, is interpreted to intrude the southern part of the MAIC. A new Pb–Pb zircon age of 2701.1 ± 6 Ma is interpreted as the magmatic age of the Fertile batholith.

Minnesota↗
USGS2025

New constraints on location and timing of the Great Lakes tectonic zone, central Upper Peninsula, Michigan, USA

The Great Lakes tectonic zone (GLTZ) forms the boundary between the Wawa–Abitibi and Minnesota River Valley subprovinces within the Archean Superior Province. The GLTZ is concealed for all of its 1100 km length, except for a segment in the central Upper Peninsula of Michigan. There, it is exposed as a northwest-striking mylonite zone along a 11 km segment, extending to the onlap of Paleozoic rocks to the east. Farther east, its location has been unknown. Here, we use aeromagnetic and gravity data to develop interpretations of the expression of the GLTZ and to define its extent under cover. Aeromagnetic gradients over the mylonite zone are interpreted to be produced by structurally juxtaposed rocks with varying magnetizations. Gravity data show a regional gradient along the GLTZ, produced by the juxtaposition of a dense greenstone belt on the north against lower-density gneisses and granites on the south. The GLTZ is interpreted to extend ∼55 km under cover to the east. The GLTZ is terminated on the east by the buried eastern arm of the ca. 1100 Ma Midcontinent Rift. An undeformed granitic dike that cuts the mylonitic foliation produces a U–Pb apatite age of 2523 ± 33 Ma, implying no major post-Archean shearing occurred, and is at odds with previous interpretations of major Proterozoic reactivation. A granite intrusion in the Minnesota River Valley subprovince produces a Pb–Pb zircon age of 2606.9 ± 3.6/7.4 Ma. This suggests that magmatism related to the Sacred Heart orogeny, previously known in Minnesota, extended to Michigan.

Michigan↗
USGS2025

Impact of gas/liquid phase change of CO2 during injection for sequestration

CO 2 sequestration in deep saline formations is an effective and important process to control the rapid rise in CO 2 emissions. The process of injecting CO 2 requires reliable predictions of the stress in the formation and the fluid pressure distributions – particularly since monitoring of the CO 2 migration is difficult – to mitigate leakage, prevent induced seismicity, and analyze wellbore stability. A key aspect of CO 2 is the gas–liquid phase transition at the temperatures and pressures of relevance to leakage and sequestration, which has been recognized as being critical for accurate predictions but has been challenging to model without ad hoc empiricisms. This paper presents a robust multiphase thermodynamics-based poromechanics model to capture the complex phase transition behavior of CO 2 and predict the stress and pressure distribution under super- and sub- critical conditions during the injection process. A finite element implementation of the model is applied to analyze the behavior of a multiphase porous system with CO 2 as it displaces the fluid brine phase. We find that if CO 2 undergoes a phase transition in the geologic reservoir, the spatial variation of the density is significantly affected, and the migration mobility of CO 2 decreases in the reservoir. A key feature of our approach is that we do not a priori assume the location of the CO 2 gas/liquid interface – or even if it occurs at all – but rather, this is a prediction of the model, along with the spatial variation of the phase of CO 2 and the change of the saturation profile due to the phase change.

Journal of the Mechanics and Physics of Solids↗
USGS2025

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado↗
USGS2026

Creation of aerial camera system boresight installations in Bangor, Maine; Monroe, Louisiana; and Bassville Park, Florida

Aerial camera systems rely on boresights for defining internal and external metrics. The U.S. Fish and Wildlife Service’s Migratory Bird Surveys Branch required boresights to calibrate their camera systems for creating georeferenced imagery that can then be used in geographic information systems. The U.S. Geological Survey created high-accuracy boresight installations in Bangor, Maine, in July 2024; Monroe, Louisiana, in September 2024; and Bassville Park, Florida, in January 2025 to provide this metric-defining capacity to the U.S. Fish and Wildlife Service. This U.S. Geological Survey Data Report details the methods and results of the boresight installation process and provides camera operators with the data needed to use each of the installations.

Florida, Louisiana, Maine↗
USGS2024

Remote sensing evapotranspiration in ensemble-based framework to enhance cascade routing and re-infiltration concept in integrated hydrological model applied to support decision making

Integrated hydrological models (IHMs) help characterize the complexity of surface–groundwater interactions. The cascade routing and re-infiltration (CRR) concept, recently applied to a MODFLOW 6 IHM, improved conceptualization and simulation of overland flow processes. The CRR controls the transfer of rejected infiltration and groundwater exfiltration from upslope areas to adjacent downslope areas where that water can be evaporated, re-infiltrated back to subsurface, or discharged to streams as direct runoff. The partitioning between these three components is controlled by uncertain parameters that must be estimated. Thus, by quantifying and reducing those uncertainties, next to uncertainties of the other model parameters (e.g. hydraulic and storage parameters), the reliability of the CRR is improved and the IHM is better suited for decision support modelling, the two key objectives of this work. To this end, the remotely sensed MODIS-ET product was incorporated into the calibration process for complementing traditional hydraulic head and streamflow observations. A total of approximately 150,000 observations guided the calibration of a 13-year MODFLOW 6 IHM simulation of the Sardon catchment (Spain) with daily stress periods. The model input uncertainty was represented by grid-cell-scale parameterization, yielding approximately 500,000 unknown input parameters to be conditioned. The calibration was carried out through an iterative ensemble smoother. Incorporating the MODIS-ET data improved the CRR implementation, and reduced uncertainties associated with other model parameters. Additionally, it significantly reduced the uncertainty associated with net recharge, a critical flux for water management that cannot be directly measured and rather is commonly estimated by IHM simulations.

Sardon catchment↗