Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Basin”

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 901 records · Page 50

Large-Scale and Global Hydrology

Powered by the sun, water moves continuously between and through Earths oceanic, atmospheric, and terrestrial reservoirs. It enables life, shapes Earths surface, and responds to and influences climate change. Scientists measure various features of the water cycle using a combination of ground, airborne, and space-based observations, and seek to characterize it at multiple scales with the aid of numerical models. Over time our understanding of the water cycle and ability to quantify it have improved, owing to advances in observational capabilities, the extension of the data record, and increases in computing power and storage. Here we present some of the most recent estimates of global and continental ocean basin scale water cycle stocks and fluxes and provide examples of modern numerical modeling systems and reanalyses.Further, we discuss prospects for predicting water cycle variability at seasonal and longer scales, which is complicated by a changing climate and direct human impacts related to water management and agriculture. Changes to the water cycle will be among the most obvious and important facets of climate change, thus it is crucial that we continue to invest in our ability to monitor it.

land hydrology↗

Spatial and temporal agreement in climate model simulations of the Interdecadal Pacific Oscillation

Accelerated warming and hiatus periods in the long-term rise of Global Mean Surface Temperature (GMST) have, in recent decades, been associated with the Interdecadal Pacific Oscillation (IPO). Critically, decadal climate prediction relies on the skill of state-of-the-art climate models to reliably represent these low-frequency climate variations. We undertake a systematic evaluation of the simulation of the IPO in the suite of Coupled Model Intercomparison Project 5 (CMIP5) models. We track the IPO in pre-industrial (control) and all-forcings (historical) experiments using the IPO tripole index (TPI). The TPI is explicitly aligned with the observed spatial pattern of the IPO, and circumvents assumptions about the nature of global warming. We find that many models underestimate the ratio of decadal-to-total variance in sea surface temperatures (SSTs). However, the basin-wide spatial pattern of positive and negative phases of the IPO are simulated reasonably well, with spatial pattern correlation coefficients between observations and models spanning the range 0.4–0.8. Deficiencies are mainly in the extratropical Pacific. Models that better capture the spatial pattern of the IPO also tend to more realistically simulate the ratio of decadal to total variance. Of the 13% of model centuries that have a fractional bias in the decadal-to-total TPI variance of 0.2 or less, 84% also have a spatial pattern correlation coefficient with the observed pattern exceeding 0.5. This result is highly consistent across both IPO positive and negative phases. This is evidence that the IPO is related to one or more inherent dynamical mechanisms of the climate system.

54 ENVIRONMENTAL SCIENCES↗

The history of mare volcanism in the Orientale Basin: Mare deposit ages, compositions and morphologies

The eruptive history of mare basalts in the Orientale Basin has been studied, using Lunar Orbiter 4 high-resolution photographs, Zond 8 photographs, and recently acquired Galileo EM-1 multispectral images. This work represents a refined set of compositional data incorporating the use of a linear mixing model for mare compositions, crater count data, and a comprehensive morphologic analysis of Orientale Basin mare deposits. Evidence for multiple eruptive episodes has been found, with compositions ranging from medium- to high-Ti basalt (less than 4 to greater than 6 wt. percent TiO2). Eruptive styles included flood, rille-forming, and shield-forming eruptions. Impact crater densities of mare units in the Orientale Basin enable determination of the ages of these deposits, using the method of Neukum et al. Earliest eruptions of mare basalt in the basin occurred at greater than or equal to 3.80 Ga and the latest eruptions occurred at about 2.3-2.5 Ga. Hence, mare volcanism occurred over a period of nearly 1.5 Ga.

Kadel, S. D.↗

Computational Tools and Workflows for Quantitative Risk Assessment and Decision Support for Geologic Carbon Storage Sites: Progress and Insights from the U.S. DOE’s National Risk Assessment Partnership

The 2005 Intergovernmental Panel on Climate Change (IPCC) Special Report on CCS raised the profile of CO2 capture and storage (CCS) as an important technology for reducing greenhouse gas (GHG) emissions. CCS is now recognized as a key component of most climate change mitigation scenarios. Since publication of that report the international research, development, and deployment (RD&D) community has advanced key technical aspects, clarified regulatory requirements, explored value chain and infrastructure solutions, and developed incentive paradigms to enable and promote large-scale deployment of CCS. These efforts have included research to better characterize geologic storage resources, to improve injection performance and storage efficiency, to assess and manage subsurface environmental risks, and to advance monitoring technologies to assure system conformance. These efforts have helped to build confidence in the viability of geologic carbon storage (GCS), but stakeholder concerns about long-term risks and liability associated with GCS remain a hurdle to broad acceptance and large-scale deployment of CCS. Since 2010, the U.S. DOE’s National Risk Assessment Partnership (NRAP) – a research collaboration between five contributing national laboratories – has worked to establish and demonstrate methods and tools to quantify and manage the subsurface environmental risks associated with GCS, amidst uncertainty. This work supports the Office of Fossil Energy and Carbon Management Carbon Transport and Storage Program’s goal of advancing safe and secure commercial-scale GCS deployment. To address the technical challenge of simulating the physical response of the GCS site to large-scale CO2 injection, NRAP has adopted an approach that relies on coupling computationally efficient reduced-order and/or data-driven proxy models of important system components (i.e., storage reservoir, sealing caprock, leakage pathways, intermediate formations, overlying groundwater aquifers, and the atmosphere) in integrated assessment framework. That integrated model of the physical system is complemented with fit-for purpose functionality to support site characterization and risk-related decisions. The recently released NRAP Phase II toolset includes the Open-Source Integrated Assessment Model (NRAP-Open-IAM) for evaluation of trends in leakage risk and potential impact, tools to support monitoring design optimization (Designs for Risk Evaluation and Management – DREAM v3.0 and Passive Seismic Monitoring Tool - PSMT), and tools for state of stress evaluation (State-of-Stress Analysis Tool - SOSAT) and forecasting induced seismicity risk. The NRAP team has also released a pair of reports describing conceptual workflows to incorporate physics-based, quantitative risk assessment into many of the design, planning, operation, and closure decisions for GCS projects. An online catalogue highlights published studies where these tools and methods are demonstrated. In this presentation, the utility of these products to assess risks and address key stakeholder questions will be highlighted through examples, and related insights about the safety and security of geologic carbon storage in qualified storage sites will be discussed. The prospect of rapid, large-scale deployment of GCS technology to aggressively reduce anthropogenic CO2 emissions requires careful consideration of interference between multiple commercial-scale storage projects within a basin. Going forward, NRAP is expanding and adapting site-scale risk quantification tools and methods to enable assessment of risks and inform management decisions for basin-scale deployment. Increasingly, this work will leverage next-generation approaches for surrogate modelling, fast prediction, and advanced visualization enabled by machine learning and artificial intelligence to promote virtual learning, scenario evaluation, and augment risk-based decision making.

quantitative risk assessment, geologic carbon stor↗

Amplified Mesoscale and Submesoscale Variability and Increased Concentration of Precipitation under Global Warming over Western North America

Abstract Cold-season precipitation statistics in simulations from the storm-resolving WRF Model at 6-km and 1-h resolution over western North America are analyzed. Pseudo–global warming future simulations for the 2041–80 period, constrained by GCMs under the RCP8.5 scenario, are compared to the 1981–2020 historical simulation. The analysis focuses on the dynamical properties of precipitation time series at subdaily scales and on the morphology of storms. The statistical distribution of precipitation intensities in each pixel of the simulation domain is characterized through nonparametric statistical indicators: frequency of wet hours, mean wet-hour precipitation intensity, and Gini coefficient as a measure of the temporal concentration of the precipitation volume. Additionally, the temporal and spatial Fourier power spectra of precipitation time series and precipitation fields are analyzed. The half-power period (HPP) and half-power wavelength (HPW) are defined as spectral measures of the characteristic scales of precipitation’s temporal and spatial patterns. The results show statistically significant increases in the mean wet-hour precipitation intensity and in the Gini coefficient in 99% of the pixels, indicating that the seasonal precipitation volume becomes more concentrated within a smaller number of hours with higher precipitation intensity. The statistics of change in the frequency of wet hours are more contrasted across the simulation domain. The changes are also reflected in the power spectra, which show the spatial and temporal variability increasing proportionally more with finer spatial and temporal scales and the HPW and HPP decreasing. These projected changes are expected to have consequences, not only in terms of hydrologic impacts but also in terms of the predictability of precipitation patterns. Significance Statement The precipitation characteristics of winter storms over the western United States and southwestern Canada are analyzed in future climate simulations for the 2041–80 period. As compared to present-day climate, the most intense parts of the storms are projected to produce a higher rainfall volume, with increased concentration over smaller areas and shorter time intervals. The propensity of rainfall intensity to vary rapidly over time will be enhanced in the future according to the simulations. These model predictions imply an increased risk of rapid flooding in small basins. They also suggest that predicting several hours ahead the time and location at which a storm will produce maximum rainfall may become more challenging in the future.

Climate change↗

Modeling of Peptides with Classical and Novel Machine Learning Force Fields: A Comparison

The replacement of classical force fields (FFs) with novel neural-network-based frameworks is an emergent topic in molecular dynamics (MD) simulations. In contrast to classical FFs, which have proven their capability to provide insights into complex soft matter systems at an atomistic resolution, the machine learning (ML) potentials have yet to demonstrate their applicability for soft materials. However, the underlying philosophy, which is learning the energy of an atom in its surrounding chemical environment, makes this approach a promising tool. In particular for the exploration of novel chemical compounds, which have not been considered in the original parametrization of classical FFs. In this article, we study the performance of the ANI-2x ML model and compare the results with those of two classical FFs, namely, CHARMM27 and the GROMOS96 43a1 FF. We explore the performance of these FFs for bulk water and two model peptides, trialanine and a 9-mer of the α-aminoisobutyric acid, in vacuum and water. The results for water describe a highly ordered water structure, with a structure similar to those using ab initio molecular dynamics simulations. The energy landscape of the peptides described by Ramachandran maps show secondary structure basins similar to those of the classical FFs but differ in the position and relative stability of the basins. Details of the sampled structures show a divergent performance of the different models, which can be related either to the short-ranged nature of the ML potentials or to shortcomings of the underlying data set used for training. These findings highlight the current state of the applicability of ANI-2x ML potential for MD simulations of soft matter systems. Simultaneously, they provide insights for future improvements of current ML potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Planetary Geophysics and Tectonics

The broad objective of this work is to improve understanding of the internal structures and thermal and stress histories of the solid planets by combining results from analytical and computational modeling, and geophysical data analysis of gravity, topography and tectonic surface structures. During the past year we performed two quite independent studies in the attempt to explain the Mariner 10 magnetic observations of Mercury. In the first we revisited the possibility of crustal remanence by studying the conditions under which one could break symmetry inherent in Runcorn's model of a uniformly magnetized shell to produce a remanent signal with a dipolar form. In the second we applied a thin shell dynamo model to evaluate the range of intensity/structure for which such a planetary configuration can produce a dipole field consistent with Mariner 10 results. In the next full proposal cycle we will: (1) develop numerical and analytical and models of thin shell dynamos to address the possible nature of Mercury s present-day magnetic field and the demise of Mars magnetic field; (2) study the effect of degree-1 mantle convection on a core dynamo as relevant to the early magnetic field of Mars; (3) develop models of how the deep mantles of terrestrial planets are perturbed by large impacts and address the consequences for mantle evolution; (4) study the structure, compensation, state of stress, and viscous relaxation of lunar basins, and address implications for the Moon s state of stress and thermal history by modeling and gravity/topography analysis; and (5) use a three-dimensional viscous relaxation model for a planet with generalized vertical viscosity distribution to study the degree-two components of the Moon's topography and gravity fields to constrain the primordial stress state and spatial heterogeneity of the crust and mantle.

Zuber, Maria↗

Integrating the Social, Hydrological and Ecological Dimensions of Freshwater Health: The Freshwater Health Index

Degradation of freshwater ecosystems and the services they provide is a primary cause of increasing water insecurity, raising the need for integrated solutions to freshwater management. While methods for characterizing the multi-faceted challenges of managing freshwater ecosystems abound, they tend to emphasize either social or ecological dimensions and fall short of being truly integrative. This paper suggests that management for sustainability of freshwater systems needs to consider the linkages between human water uses, freshwater ecosystems and governance. We present a conceptualization of freshwater resources as part of an integrated social-ecological system and propose a set of corresponding indicators to monitor freshwater ecosystem health and to highlight priorities for management. We demonstrate an application of this new framework -the Freshwater Health Index (FHI) - in the Dongjiang River Basin in southern China, where stakeholders are addressing multiple and conflicting freshwater demands. By combining empirical and modeled datasets with surveys to gauge stakeholders' preferences and elicit expert information about governance mechanisms, the FHI helps stakeholders understand the status of freshwater ecosystems in their basin, how ecosystems are being manipulated to enhance or decrease water-related services, and how well the existing water re-source management regime is equipped to govern these dynamics over time. This framework helps to operationalize a truly integrated approach to water resource management by recognizing the interplay between governance, stakeholders, freshwater ecosystems and the services they provide

ecosystem services↗

Surface Quantitative Precipitation Estimates (SQUIRE) of Snow Water Equivalent from the Surface Atmospheric Integrated Field Laboratory

The upper Colorado River basin is the primary source of water for 40 million people. With declining snowpack in the basin, forecasting hydrological budgets in the Southwest United States is more important than ever. However, due in part, to a lack of reliable observations of precipitation in complex terrain, hydrological models struggle to assess and forecast snowpack snow water equivalent (SWE) in the upper Colorado River basin (UCRB). Therefore, the need for more reliable SWE forecasts in the UCRB motivated the U.S. Department of Energy Atmospheric Radiation Measurement Facility’s Surface Atmospheric Integrated Field Laboratory (SAIL) that occurred from June 2021 to June 2023. During SAIL, the X-band precipitation radar from Colorado State University conducted volume scans sampling the precipitation properties over the UCRB. The ARM facility developed a gridded Surface Quantitative Precipitation Estimates (SQUIRE) product from the radar observations. To do this, various daily SWE estimates from radar using the radar reflectivity factor Z e and specific differential phase K dp were compared against ground-based precipitation gauges. SWE in precipitation calculated from Wolfe and Snider’s S–Z e estimator was in best agreement with the rain gauges for the days when SWE < 12 mm. For days with SWE > 12 mm, the WSR-88D Intermountain West relationship had the best agreement with the precipitation gauges. Airborne snow depth observations show that SQUIRE captures regions of orographic enhancement in the mountains to the west and northwest of the SAIL study area, indicating that the scientific community should focus on understanding and ultimately simulating orographic atmospheric precipitation processes to improve UCRB snowpack SWE assessment and forecasting.

Hydrology↗

AK112: Full Waveform Inversion Tomography of Alaska Improves Waveform Fits While Imaging Crustal, Mantle, and Slab Structure

We report a full waveform inversion tomography model of Alaska and the surrounding regions, inferring radially anisotropic shear and isotropic compressional wavespeeds by fitting complete waveforms from 120 regional earthquakes. Our multiscale approach inverted time–frequency phase misfits (maximum period of 100 s), starting with a minimum period of 40 s and ending at 20 s in 7 stages and 112 total iterations. The model (AK112) was evaluated by computing the misfits for 36 independent validation events. We find that misfit reductions were large and equal (∼55%) for both the inversion and validation data sets, providing confidence in the model. AK112 also provides much better waveform fits compared to other reported models for the region, including an isotropic version of itself, highlighting the importance of anisotropy. The model resolves known crustal, upper mantle, and slab structure to depths of 100 km with new detail: sedimentary basins in the Alaskan Shelf, Cook Inlet, and Colville basins, among others; discontinuous lithospheric structure across major terrane boundaries; and subducting slab geometry and back‐arc volcanic sources. In addition to tectonic interpretations, the model enables full waveform simulations for long‐period earthquake ground motions and source characterization (e.g., moment tensor and finite‐fault inversion).

Rodgers, Arthur [Lawrence Livermore National Labor↗

Materials and formation of the Imbrium basin

A study of Imbrium basin deposits has been conducted, including geologic mapping, analysis of orbital geochemical data, earth-based near-IR spectroscopy, and analytical modeling. The morphologic facies of the deposits display an unusual bilateral symmetry with respect to the main rim of the basin, which is thought to be 1160 km in diameter and composed of Apennine-Carpathian-Alpes mountain ranges. The Imbrium ring structure consists of six concentric rings ranging from 550-3200 km in diameter. The evolution of the target for the Imbrium basin impact is discussed. It is suggested that the Imbrium impact, which occurred 3.85 Gyrs ago, formed a basin predominantly by a proportional-growth crater-forming mechanism.

Spudis, Paul D.↗

Deep Learning Experiments for Tropical Cyclone Intensity Forecasts

Reducing tropical cyclone (TC) intensity forecast errors is a challenging task that has interested the operational forecasting and research community for decades. To address this, we developed a deep learning (DL)-based multilayer perceptron (MLP) TC intensity prediction model. The model was trained using the global Statistical Hurricane Intensity Prediction Scheme (SHIPS) predictors to forecast the change in TC maximum wind speed for the Atlantic basin. In the first experiment, a 24-h forecast period was considered. To overcome sample size limitations, we adopted a leave one year out (LOYO) testing scheme, where a model is trained using data from all years except one and then evaluated on the year that is left out. When tested on 2010–18 operational data using the LOYO scheme, the MLP outperformed other statistical–dynamical models by 9%–20%. Additional independent tests in 2019 and 2020 were conducted to simulate real-time operational forecasts, where the MLP model again outperformed the statistical–dynamical models by 5%–22% and achieved comparable results as HWFI. The MLP model also correctly predicted more rapid intensification events than all the four operational TC intensity models compared. In the second experiment, we developed a lightweight MLP for 6-h intensity predictions. When coupled with a synthetic TC track model, the lightweight MLP generated realistic TC intensity distribution in the Atlantic basin. Therefore, the MLP-based approach has the potential to improve operational TC intensity forecasts, and will also be a viable option for generating synthetic TCs for climate studies.

58 GEOSCIENCES↗

An Inversion of Gravity and Topography for Mantle and Crustal Structure on Mars

Analysis of the gravity and topography of Mars presently provides our primary quantitative constraints on the internal structure of Mars. We present an inversion of the long-wavelength (harmonic degree less than or equal to 10) gravity and topography of Mars for lateral variations of mantle temperature and crustal thickness. Our formulation incorporates both viscous mantle flow (which most prior studies have neglected) and isostatically compensated density anomalies in the crust and lithosphere. Our nominal model has a 150-km-thick high-viscosity surface layer over an isoviscous mantle, with a core radius of 1840 km. It predicts lateral temperature variations of up to a few hundred degrees Kelvin relative to the mean mantle temperature, with high temperature under Tharsis and to a lesser extent under Elysium and cool temperatures elsewhere. Surprisingly, the model predicts crustal thinning beneath Tharsis. If correct, this implies that thinning of the crust by mantle shear stresses dominates over thickening of the crust by volcanism. The major impact basins (Hellas, Argyre, Isidis, Chryse, and Utopia) are regions of crustal thinning, as expected. Utopia is also predicted to be a region of hot mantle, which is hard to reconcile with the surface geology. An alternative model for Utopia treats it as a mascon basin. The Utopia gravity anomaly is consistent with the presence of a 1.2 to 1.6 km thick layer of uncompensated basalt, in good agreement with geologic arguments about the amount of volcanic fill in this area. The mantle thermal structure is the dominant contributor to the observed geoid in our inversion. The mantle also dominates the topography at the longest wavelengths, but shorter wavelengths (harmonic degrees greater than or equal to 4) are dominated by the crustal structure. Because of the uncertainty about the appropriate numerical values for some of the model's input parameters, we have examined the sensitivity of the model results to the planetary structural model (core radius and core and mantle densities), the mantle's viscosity stratification, and the mean crustal thickness. The model results are insensitive to the specific thickness or viscosity contrast of the high-viscosity surface layer and to the mean crustal thickness in the range 25 to 100 km. Models with a large core radius or with an upper mantle low-viscosity zone require implausibly large lateral variations in mantle temperature.

Kiefer, Walter S.↗

A review of SWAT applications, performance and future needs for simulation of hydro-climatic extremes

Hydro-climatic extremes, such as droughts and floods, have most likely increased due to climatic change and could lead to severe impacts on socio-economic, structural and environmental sectors. With nearly 4000 publications, the Soil and Water Assessment Tool (SWAT) is clearly one of the most extensively used ecohydrological models worldwide. The model has been widely used for projecting the impacts of future hydro-climatic changes, but application for extreme streamflow conditions is still rarely reported. To date, SWAT application reviews have focused on compilations of SWAT studies for specific or relatively new applications such as eco-hydrological modelling, ecosystem services, sub-daily simulations, and pesticide fate and transport simulations. However, no existing SWAT review studies have focused on simulation of hydro-climatic extremes. Furthermore, this research aims to bridge this gap by compiling and reviewing the findings of studies reporting SWAT hydro-climatic extremes including highlighting the performance and future research needs. A total of 111 articles have been identified since 1999; most of these studies were conducted in the United States and China. These articles can be divided into extreme flow assessments, drought studies, flood studies, drought and flood studies, SWAT coupling with other models, and SWAT improvements. Most of the extreme performance assessment studies reported “satisfactory ”performance, with a particular emphasis on peak flow comparisons. Future research needs regarding this topic include: (1) a unified SWAT extreme performance assessment framework; (2) SWAT improvements that result in improved replication of peak and low flows; (3) reliability assessment of global and satellite products for SWAT extreme simulations; (4) bias correction of CMIP6 and regional climate projections; (5) comparison of SWAT + and SWAT for extreme flow simulations in different types of basins; (6) development of an extreme flow module within an overall SWAT modelling system; and (7) integration of artificial intelligence within SWAT modelling.

54 ENVIRONMENTAL SCIENCES↗

Does Ocean Color Data Assimilation Improve Estimates of Global Ocean Inorganic Carbon?

Ocean color data assimilation has been shown to dramatically improve chlorophyll abundances and distributions globally and regionally in the oceans. Chlorophyll is a proxy for phytoplankton biomass (which is explicitly defined in a model), and is related to the inorganic carbon cycle through the interactions of the organic carbon (particulate and dissolved) and through primary production where inorganic carbon is directly taken out of the system. Does ocean color data assimilation, whose effects on estimates of chlorophyll are demonstrable, trickle through the simulated ocean carbon system to produce improved estimates of inorganic carbon? Our emphasis here is dissolved inorganic carbon, pC02, and the air-sea flux. We use a sequential data assimilation method that assimilates chlorophyll directly and indirectly changes nutrient concentrations in a multi-variate approach. The results are decidedly mixed. Dissolved organic carbon estimates from the assimilation model are not meaningfully different from free-run, or unassimilated results, and comparisons with in situ data are similar. pC02 estimates are generally worse after data assimilation, with global estimates diverging 6.4% from in situ data, while free-run estimates are only 4.7% higher. Basin correlations are, however, slightly improved: r increase from 0.78 to 0.79, and slope closer to unity at 0.94 compared to 0.86. In contrast, air-sea flux of C02 is noticeably improved after data assimilation. Global differences decline from -0.635 mol/m2/y (stronger model sink from the atmosphere) to -0.202 mol/m2/y. Basin correlations are slightly improved from r=O.77 to r=0.78, with slope closer to unity (from 0.93 to 0.99). The Equatorial Atlantic appears as a slight sink in the free-run, but is correctly represented as a moderate source in the assimilation model. However, the assimilation model shows the Antarctic to be a source, rather than a modest sink and the North Indian basin is represented incorrectly as a sink rather than the source indicated by the free-run model and data estimates.

Gregg, Watson↗

mosartwmpy: A Python implementation of the MOSART-WM coupled hydrologic routing and water management model

mosartwmpy is a Python implementation of the Model for Scale Adaptive River Transport with Water Management (MOSART-WM). This new version retains the functionality of the legacy model (written in FORTRAN) while providing new features to enhance user experience and extensibility. MOSART is a large-scale river-routing model used to study riverine dynamics of water, energy, and biogeochemistry cycles across local, regional, and global scales (Li et al., 2013). The WM component introduced by Voisin et al. (2013) represents river regulation through reservoir storage and release operations, diversions from reservoir releases, and allocation to sectoral water demands. Each reservoir release is independently calibrated using long-term mean monthly inflow into the reservoir, long-term mean monthly demand associated with this reservoir, and reservoir goals (flood control, irrigation, recreation, etc.). Generic monthly pre-release rules and storage targets are set up for individual reservoirs; however, those releases are updated annually for inter-annual variability (dry or wet year) and daily for environmental constraints such as flow minimum release and minimum/maximum storage levels. The WM model allows an evaluation of the impact of water management over multiple river basins at once (global, continental scales) and with consistent representation of human operations over the full domain.

97 MATHEMATICS AND COMPUTING↗

Limits on Runoff Episode Duration for Early Mars: Integrating Lake Hydrology and Climate Models

Fluvio-lacustrine features on the martian surface attest to a climate that was radically different in the past. Since climate models have difficulty sustaining a liquid hydrosphere at the surface, multiple cycles of runoff episodes may have characterized the ancient Mars climate. A fundamental question thus remains: what was the duration of these runoff-producing episodes? Here we use morphometric measurements from newly identified coupled lake systems, containing both an open- and a closed-basin lake (n=7). We combined hydrological balances with precipitation outputs from climate models, and found that breaching runoff episodes likely lasted 102−105yr; other episodes may have been shorter but could not be longer. Runoff episode durations are model-dependent and spatiallyvariable, and no27climate model scenario can satisfy a unique duration for all coupled systems. In the near future, these quantitative constraints on early Mars lake persistence may be tested through in situ observations from Perseverance rover. Plain Language Summary Ancient river and lake beds on the martian surface point to a climate that was radically different in the past, unlike the dry, cold desert that Mars is today. It experienced significant liquid water flowing across the surface (i.e., ‘runoff’) and hosted long-lived bodies of water (lakes and seas). However, scientists that generate models to simulate the ancient Mars climate have difficulty sustaining this proposed liquid water cycle at the surface. One possible hypothesis is that multiple cycles of individual liquid water episodes may have characterized the ancient martian climate (i.e. runoff episodes). A fundamental question thus remains: what was the duration of these individual runoff-producing episodes? Here we use an approach that takes advantage of unique martian lake system, ‘coupled’ lake systems, and their shapes. We combined these measurements with precipitation results from previous climate models, and found that runoff episodes likely lasted 100-10,000 years, requiring a global climate mechanism that generated repeated cycles of this duration. Ultimately, our results can be tested through observations made by NASA’s Perseverance rover, shedding light on the uninterrupted duration of surface water on Mars and its potential to host life in the past.

Gaia Stucky de Quay↗

Are we modeling the properties of unconventional shales correctly?

Technological advancements have led to impressive growth in hydrocarbon (HC) production from unconventional shale reservoirs. However, major limiting factors in shale gas development are i) low recovery of gas in place (<20%), (ii) a rapid decline in well productivity, iii) release of contaminants, despite the use of advanced hydraulic fracturing fluids and multiple fracturing stages. This is mainly due to a lack of understanding of the nanoscale physicochemical properties of shales, especially that of kerogen, the macromolecule that is not only the source of the majority of the HC’s in shales but also holds most of the HC’s in adsorbed state. Over the last few years, a large number of studies have been published on simulating the physicochemical properties of shale utilizing molecular models of kerogen. However, the molecular models of kerogen input into these simulations are based on the kerogen “type” derived from a very limited number of shale samples. In this paper, we examine the variations that can exist in kerogen structure within a particular kerogen “type” across different shale basins and single shale basin at similar thermal maturity levels. We propose that using kerogen “type” based structural models for molecular simulations could lead to inaccurate estimation of HC reserves, HC recovery, HC production, frackability, and quantity and quality of produced water. Our study highlights the need for developing a better classification of kerogen based on its “molecular structure” instead of “type” for a more accurate prediction of physicochemical properties of shales.

04 OIL SHALES AND TAR SANDS↗