Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “volume distribution”

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 973 records · Page 54

Cometary Activity Begins at Kuiper Belt Distances: Evidence from C/2017 K2

We study the development of activity in the incoming long-period comet C/2017 K2 over the heliocentric distance range 9 ≲ r {sub H} ≲ 16 au. The comet continues to be characterized by a coma of submillimeter-sized and larger particles ejected at low velocity. In a fixed co-moving volume around the nucleus we find that the scattering cross section of the coma, C, is related to the heliocentric distance by a power law, C∝r{sub H}{sup -s}, with heliocentric index s = 1.14 ± 0.05. This dependence is significantly weaker than the r {sub H} {sup -2} variation of the insolation as a result of two effects. These are, first, the heliocentric dependence of the dust velocity and, second, a lag effect due to very slow-moving particles ejected long before the observations were taken. A Monte Carlo model of the photometry shows that dust production beginning at r {sub H} ~ 35 au is needed to match the measured heliocentric index, with only a slight dependence on the particle size distribution. Mass-loss rates in dust at 10 au are of order 10{sup 3} kg s{sup -1}, while loss rates in gas may be much smaller, depending on the unknown dust to gas ratio. Consequently, the ratio of the nongravitational acceleration to the local solar gravity, α', may, depending on the nucleus size, attain values of ~10{sup -7} ≲ α' ≲ 10{sup -5}, comparable to values found in short-period comets at much smaller distances. Nongravitational acceleration in C/2017 K2 and similarly distant comets, while presently unmeasured, may limit the accuracy with which we can infer the properties of the Oort cloud from the orbits of long-period comets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Crystallization Mechanism in Spark Plasma Sintered Bulk Metallic Glass Analyzed using Small Angle Neutron Scattering

Understanding the thermal stability of metallic glasses is critical to determining their safe temperatures of service. In this paper, the crystallization mechanism in spark plasma sintered Fe 48 Cr 15 Mo 14 Y 2 C15B6 metallic glass is established by analyzing the crystal size distribution using x-ray diffraction, transmission electron microscopy and in-situ small angle neutron scattering. Isothermal annealing at 700 °C and 725 °C for 100 min resulted in the formation of (Fe,Cr) 23 C 6 crystals, measured from transmission electron micrographs, to be from 10 to 30 nm. The small angle neutron scattering intensity measured in-situ, over a Q-range of 0.02 to 0.3 Å -1 , during isothermal annealing of the sintered samples, confirmed the presence of (Fe,Cr) 23 C 6 crystals. The measured scattering intensity, fitted by the maximum entropy model, over the Q-range of 0.02 to 0.06 Å -1 , revealed that the crystals had radii ranging from 3 to 18 nm. The total volume fraction of crystals were estimated to be 0.13 and 0.22 upon isothermal annealing at 700 °C and 725 °C for 100 min respectively. The mechanism of crystallization in this spark plasma sintered iron based metallic glass was established to be from pre-existing nuclei as confirmed by Avrami exponents of 0.25 ± 0.01 and 0.39 ± 0.01 at the aforesaid temperatures.

42 ENGINEERING↗

Sustainable Aviation Fuel from Hydrothermal Liquefaction of Wet Wastes

Hydrothermal liquefaction (HTL) uses heat and pressure to liquefy the organic matter in biomass/waste feedstocks to produce biocrude. When hydrotreated the biocrude is converted into transportation fuels including sustainable aviation fuel (SAF). Further, by liquifying the organic matter in wet wastes such as sewage sludge, manure, and food waste, HTL can prevent landfilling or other disposal methods such as anerobic digestion, or incineration. A significant roadblock to the development of a new route for SAF is the strict approval process, and the large volumes required (>400 L) for testing. Tier α and β testing can predict some of the properties required for ASTM testing with <400 mL samples. The current study is the first to investigate the potential for utilizing wet-waste HTL biocrude (WWHTLB) as an SAF feedstock. Herein, several WWHTLB samples were produced from food waste, sewage sludge, and fats, oils, and grease, and subsequently hydrotreated and distilled to produce SAF samples. The fuels (both undistilled and distilled samples) were analyzed via elemental and 2D-GC-MS. Herein, we report the Tier α and β analysis of an SAF sample derived originally from a WWHTLB. The results of this work indicate that the upgraded WWHTLB material exhibits key fuel properties, including carbon number distribution, distillation profile, surface tension, density, viscosity, heat of combustion, and flash point, which all fall within the required range for aviation fuel. WWHTLB has therefore been shown to be a promising candidate feedstock for the production of SAF.

09 BIOMASS FUELS↗

Lyman-$α$ forest holography: 3D predictions from 1D measurements

Cosmological analyses of Lyman-$α$ forest clustering rely on either one-dimensional correlations along individual sightlines or three-dimensional correlations between different sightlines. Because these observables probe the matter distribution on very different scales, they have traditionally been analyzed independently. In this work, we bridge this gap using ForestFlow, an emulator trained on a suite of cosmological hydrodynamical simulations that provides a unified description of Lyman-$α$ forest clustering from linear to nonlinear scales. This framework enables us to determine the range of three-dimensional clustering models compatible with the DESI one-dimensional flux power spectrum ($P_{\rm 1D}$). The resulting predictions successfully reproduce the large-scale clustering measured by the DESI BAO analysis and provide physically motivated priors on nonlinear clustering that are used in a companion paper presenting the full-shape analysis of the DESI DR2 Lyman-$α$ forest. We validate our methodology using the large-volume, high-resolution hydrodynamical simulation ACCEL-2, demonstrating excellent agreement across the full range of scales considered. Finally, we combine constraints from the $P_{\rm 1D}$ and BAO analyses on the parameter combinations $b_δσ_8$ and $b_ηf σ_8$, finding that the two probes provide comparable constraining power while exhibiting complementary parameter degeneracies. Our results establish a direct connection between one- and three-dimensional Lyman-$α$ forest measurements through ForestFlow, an approach we term Lyman-$α$ holography by analogy with the reconstruction of higher-dimensional structure from lower-dimensional information.

Chaves-Montero, J. [Barcelona, IFAE] (ORCID:000000↗

OPFLearnData: Dataset for Learning AC Optimal Power Flow

The datasets are resulting from OPFLearn.jl, a Julia package for creating AC OPF datasets. The package was developed to provide researchers with a standardized way to efficiently create AC OPF datasets that are representative of more of the AC OPF feasible load space compared to typical dataset creation methods. The OPFLearn dataset creation method uses a relaxed AC OPF formulation to reduce the volume of the unclassified input space throughout the dataset creation process. The dataset contains load profiles and their respective optimal primal and dual solutions. Load samples are processed using AC OPF formulations from PowerModels.jl. More information on the dataset creation method can be found in our publication, "OPF-Learn: An Open-Source Framework for Creating Representative AC Optimal Power Flow Datasets" and in the package website: https://github.com/NREL/OPFLearn.jl.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Bringing function to structure: Root–soil interactions shaping phosphatase activity throughout a soil profile in Puerto Rico

Abstract Large areas of highly productive tropical forests occur on weathered soils with low concentrations of available phosphorus (P). In such forests, root and microbial production of acid phosphatase enzymes capable of mineralizing organic phosphorus is considered vital to increasing available P for plant uptake. We measured both root and soil phosphatase throughout depth and alongside a variety of root and soil factors to better understand the potential of roots and soil biota to increase P availability and to constrain estimates of the biochemical mineralization within ecosystem models. We measured soil phosphatase down to 1 m, root phosphatase to 30 cm, and collected data on fine‐root mass density, specific root length, soil P, bulk density, and soil texture using soil cores in four tropical forests within the Luquillo Experimental Forest in Puerto Rico. We found that soil phosphatase decreased with soil depth, but not root phosphatase. Furthermore, when both soil and root phosphatase were expressed per soil volume, soil phosphatase was 100‐fold higher that root phosphatase. Both root and soil factors influenced soil and root phosphatase. Soil phosphatase increased with fine‐root mass density and organic P, which together explained over 50% of the variation in soil phosphatase. Over 80% of the variation in root phosphatase per unit root mass was attributed to specific root length (positive correlation) and available (resin) P (negative correlation). Synthesis : Fine‐root traits and soil P data are necessary to understand and represent soil and root phosphatase activity throughout the soil column and across sites with different soil conditions and tree species. These findings can be used to parameterize or benchmark estimates of biochemical mineralization in ecosystem models that contain fine‐root biomass and soil P distributions throughout depth.

59 BASIC BIOLOGICAL SCIENCES↗

Sample Preparation Methods for Targeted Single-Cell Proteomics

We compared three cell isolation and two proteomic sample preparation methods for single-cell and near-single-cell analysis. Whole blood was used to quantify hemoglobin (Hb) and glycated-Hb (gly-Hb) in erythrocytes using targeted mass spectrometry and stable isotope-labeled standard peptides. Each method differed in cell isolation and sample preparation as follows: 1) FACS and automated preparation in one-pot for trace samples (autoPOTS); 2) limited dilution via microscopy and a novel rapid one-pot sample preparation method that circumvented the need for the solid-phase extraction, low-volume liquid handling instrumentation and humidified incubation chamber; and 3) CellenONE-based cell isolation and the same one-pot sample preparation method used for limited dilution. Only the CellenONE device routinely isolated single-cells from which Hb was measured to be 540–660 amol per red blood cell (RBC), which was comparable to the calculated SI reference range for mean corpuscular hemoglobin (390–540 amol/RBC). FACSAria sorter and limited dilution could routinely isolate single-digit cell numbers, to reliably quantify CMV-Hb heterogeneity. Finally, we observed that repeated measures, using 5–25 RBCs obtained from N = 10 blood donors, could be used as an alternative and more efficient strategy than single RBC analysis to measure protein heterogeneity, which revealed multimodal distribution, unique for each individual.

59 BASIC BIOLOGICAL SCIENCES↗

High-fidelity micromechanical modeling of the effects of defects on damage and creep behavior in single tow ceramic matrix composite

Despite the superior properties of ceramic matrix composites (CMCs), their fabrication process generates inevitable defects with high density that significantly influence material integrity and residual useful life. Recent efforts have focused on CMC property prediction and investigation of their different inelastic mechanisms, but very limited work exists on understanding the influence of defects on inelastic responses of CMCs. Here, we introduce a three-dimensional micromechanics computational framework that includes experimentally informed material microstructure and architectural variabilities to investigate CMC response in the presence of manufacturing-induced defects. A developed microstructure generation algorithm is used to generate the representative volume elements in the micromechanics models from complex material morphology information obtained from extensive characterization studies of C/SiNC and SiC/SiNC CMCs. A fracture mechanics-informed matrix damage model is reformulated, in which the model considers the growth of porosity and microcracks in the as-received material. A progressive fiber damage model as well as Orthotropic viscoplasticity creep formulation are also utilized for the study. This methodology is implemented within the micromechanics framework to investigate the complex temperature- and time-dependent load transfer and damage mechanisms of CMCs under operation loading conditions. The developed framework explores the influence of as-received defects on the damage and creep behavior in service conditions. It also provides new insights into the effect of size, shape, distribution, and location of these defects on material response. The framework is then calibrated with unidirectional CMC minicomposite results available in the literature.

36 MATERIALS SCIENCE↗

Role of Parallel Solenoidal Electric Field on Energy Conversion in 2.5D Decaying Turbulence with a Guide Magnetic Field

We perform 2.5D particle-in-cell simulations of decaying turbulence in the presence of a guide (out-of-plane) background magnetic field. The fluctuating magnetic field initially consists of Fourier modes at low wavenumbers (long wavelengths). With time, the electromagnetic energy is converted to plasma kinetic energy (bulk flow+thermal energy) at the rate per unit volume of J· E for current density J and electric field E. Such decaying turbulence is well known to evolve toward a state with strongly intermittent plasma current. Here we decompose the electric field into components that are irrotational, E ir , and solenoidal (divergence-free), E so . E ir is associated with charge separation, and J · E ir is a rate of energy transfer between ions and electrons with little net change in plasma kinetic energy. Therefore, the net rate of conversion of electromagnetic energy to plasma kinetic energy is strongly dominated by J · E so , and for a strong guide magnetic field, this mainly involves the component E so,∥ parallel to the total magnetic field B. We examine various indicators of the spatial distribution of the energy transfer rate J ∥ · E so,∥ , which relates to magnetic reconnection, the best of which are (1) the ratio of the out-of-plane electric field to the in-plane magnetic field, (2) the out-of-plane component of the nonideal electric field, and (3) the magnitude of the estimate of current helicity.

79 ASTRONOMY AND ASTROPHYSICS↗

First Measurement of Inclusive Muon Neutrino Charged Current Triple Differential Cross Section on Argon

The field of accelerator neutrino experiments is entering an era of precision oscillation measurements where the remaining unknown neutrino measurements will be determined. The upcoming DUNE and Hyper-K experiments aim to determine the neutrino mass hierarchy and degree of Charge-Parity (CP) violation in the neutrino sector, providing potential insight on the matter-antimatter imbalance observed in the universe. However, these experiments require highly accurate measurements, and neutrino cross section modeling uncertainties may limit their capabilities. Cross section measurements at current- generation experiments can aid the development of neutrino interaction models to reduce these uncertainties. This is especially true for measurements of neutrino energy, as it drives neutrino oscillations and is of key importance to oscillation experiments. The MicroBooNE experiment uses a Liquid Argon Time Projection Chamber (LArTPC) to produce neutrino-argon cross sections as one of its physics goals. The MicroBooNE detector’s fully active volume, precision reconstruction, and calorimetry information are leveraged in the Wire-Cell analysis to produce a muon neutrino selection that is 92% pure while maintaining 68% efficiency. A reconstruction chain featuring a fully 3D charge reconstruction and a graph-based particle trajectory fit are used to produce accurate measurements of lepton kinematics as well as visible hadronic energy produced in a neutrino interaction. This thesis presents the first neutrino-argon triple-differential cross section measurement, targeting inclusive charged-current final states. Wiener SVD unfolding is used to produce a measurement over neutrino energy, muon momentum, and muon scattering angle. A series of constrained goodness of fit tests are used to demonstrate the validity of MicroBooNE’s model in describing the distribution of reconstructed kinematics seen in data to ensure the accuracy of unfolding. The validated unfolding to neutrino energy represents a step forward in the field of neutrino cross sections, and demonstrates the capabilities of the LArTPC detector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multi-modal chemical characterization of highly viscous submicrometer organic particles

Distinguishing highly viscous organic particles within complex mixtures of atmospheric aerosol and accurate descriptions of their composition, size distributions, and mixing states are challenges at the forefront of aerosol measurement science and technology. Here, we present results obtained from complementary single-particle measurement techniques employed for the in-depth characterization of highly viscous particles. In this study, we demonstrate advantages and synergy of this multi-modal particle characterization approach based on the analysis of individual viscous particles formed in the air-discharged waste produced by a common sewer pipe rehabilitation technology. Using oil immersion flow microscopy, we investigate particle size distributions and morphology of colloidal components present in field-collected aqueous waste condensates. We compare these results with corresponding measurements of viscous particles formed in drying droplets of the aerosolized discharged waste. The colloidal components and viscous particles were found to be approximately 10 µm and 0.5 µm, respectively. The aerosolized viscous particles exhibited a spherical morphology, while the colloidal particles appeared noticeably fractal, resembling fragments of a cured composite material. Chemical imaging of the viscous particles collected on substrates was performed using scanning electron microscopy and soft X-ray spectro-microscopy techniques. Through these methods, comprehensive description of these particles emerged, confirming their high solid-like viscosity, wide-ranging sizes, diverse carbon speciation with high degrees of oxygenation, and high organic volume fractions. The aerosolized viscous particles were further characterized using high-throughput single particle mass spectrometry. This technique provides real-time measurements of composition, size, and morphological metrics for large numbers of individual particles, enabling the identification of their distinct mass spectrometric signatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimizing Fast Charging and Wetting in Lithium-Ion Batteries with Optimal Microstructure Patterns Identified by Genetic Algorithm

To sustain the high-rate current required for fast charging electric vehicle batteries, electrodes must exhibit sufficiently high effective ionic diffusion. Additionally, to reduce battery manufacturing costs, wetting time must decrease. Both of these issues can be addressed by structuring the electrodes with mesoscale pore channels. However, their optimal spatial distribution, or patterns, is unknown. Herein, a genetic algorithm has been developed to identify these optimal patterns using a CPU-cheap proxy distance-based model to evaluate the impact of the added pore networks. Both coin-cell and pouch cell form factors have been considered for the wetting analysis, with their respective electrolyte infiltration mode. Regular hexagonal and mud-crack-like patterns, respectively, for fast charging and fast wetting were found to be optimal and have been compared with pre-determined, easier to manufacture, patterns. The model predicts that using cylindrical channels arranged in a regular hexagonal pattern is ∼6.25 times more efficient for fast charging as compared to grooved lines with both structuring strategies being restricted to a 5% electrode total volume loss. The model also shows that only a very limited electrode volume loss (1%–2%) is required to dramatically improve the wetting (5–20 times) compared to an unstructured electrode.

25 ENERGY STORAGE↗

Observation of Cosmic-Ray Anisotropy in the Southern Hemisphere with 12 yr of Data Collected by the IceCube Neutrino Observatory

We analyzed the 7.92 × 10 11 cosmic-ray-induced muon events collected by the IceCube Neutrino Observatory from 2011 May 13, when the fully constructed experiment started to take data, to 2023 May 12. This data set provides an up-to-date cosmic-ray arrival direction distribution in the Southern Hemisphere with unprecedented statistical accuracy covering more than a full period length of a solar cycle. Improvements in Monte Carlo event simulation and better handling of year-to-year differences in data processing significantly reduce systematic uncertainties below the level of statistical fluctuations compared to the previously published results. We confirm the observation of a change in the angular structure of the cosmic-ray anisotropy between 10 TeV and 1 PeV, more specifically in the 100–300 TeV energy range. For the first time, we analyzed the angular power spectrum at different energies. The observed variations of the power spectra with energy suggest relatively reduced large-scale features at high energy compared to those of medium and small scales. The large volume of data enhances the statistical significance at higher energies, up to the PeV scale, and smaller angular scales, down to approximately 6° compared to previous findings.

79 ASTRONOMY AND ASTROPHYSICS↗

$\bar{\nu}_\mu$ charged-current $\pi^0$ data release

Data release for the NOvA muon antineutrino charged-current (CC) pi^0 cross section presented in arXiv:2511.05807. The signal for this analysis is defined as muon antineutrino CC interactions in the fiducial volume of the NOvA near detector (a 2.7 m × 2.7 m × 9.0 m region) that produce at least one pi^0 in the final state emerging from the nucleus, within the phase space of muon momentum [0.5, 2.5) GeV/c and muon angle [0, 60) degree, as described in arXiv:2511.05807. The released zip file contains two files: NOvA_NumubarCCPi0_DataRelease.root README.txt The ROOT file includes the cross-section results as well as the statistical and systematic covariance matrices for each variable used in this analysis. The README provides a detailed description of the contents of the data release. Official Flux: The flux used in this analysis is available from the NOvA Public Docs: https://publicdocs.fnal.gov/cgi-bin/ShowDocument?docid=8. File structure --- The ROOT file contains the following TDirectories: pi0p - pi^0 momentum distributions pi0dir - pi^0 angular distributions muonp - muon momentum distributions muondir - muon angular distributions Q2 - reconstructed Q^2 distributions Wmass - reconstructed W_mass distributions Each directory contains three histograms: xsec (TH1D): Cross section result cov_stat (TH2D): Statistical covariance matrix cov_syst (TH2D): Systematic covariance matrix Usage notes: - xsec gives the measured differential cross section w.r.t. the corresponding variable. - cov_stat and cov_syst provide the full covariance matrices. - The bin definitions and kinematic phase spaces follow those used in arXiv:2511.05807. Citation --- If you use these data, please cite: NOvA Collaboration, arXiv:2511.05807.

Wu, Wanwei [Pittsburgh U.] (ORCID:0000000326327215↗

Fast HARDI Uncertainty Quantification and Visualization with Spherical Sampling

In this paper, we study uncertainty quantification and visualization of orientation distribution functions (ODF), which corresponds to the diffusion profile of high angular resolution diffusion imaging (HARDI) data. The shape inclusion probability (SIP) function is the state‐of‐the‐art method for capturing the uncertainty of ODF ensembles. The current method of computing the SIP function with a volumetric basis exhibits high computational and memory costs, which can be a bottleneck to integrating uncertainty into HARDI visualization techniques and tools. We propose a novel spherical sampling framework for faster computation of the SIP function with lower memory usage and increased accuracy. In particular, we propose direct extraction of SIP isosurfaces, which represent confidence intervals indicating spatial uncertainty of HARDI glyphs, by performing spherical sampling of ODFs. Our spherical sampling approach requires much less sampling than the state‐of‐the‐art volume sampling method, thus providing significantly enhanced performance, scalability, and the ability to perform implicit ray tracing. Our experiments demonstrate that the SIP isosurfaces extracted with our spherical sampling approach can achieve up to 8164× speedup, 37282× memory reduction, and 50.2% less SIP isosurface error compared to the classical volume sampling approach. We demonstrate the efficacy of our methods through experiments on synthetic and human‐brain HARDI datasets.

97 MATHEMATICS AND COMPUTING↗

Accelerated Data Analytics for Power-System Time-Series (ADAPT) [Final Project Report]

Over the past decade, electric utilities have made significant progress in deploying networks of phasor measurement units (PMUs). The high-speed measurements provided by PMUs are valuable for offline analysis, but the sheer volume of data presents a challenge for utilities. The Accelerated Data Analytics for Power-System Time-Series (ADAPT) project was created to address this challenge by developing SciSync, an open-source software tool that enables utilities to rapidly read, process, analyze, and review grid measurements. SciSync was designed to provide the capabilities of Archive Walker, a research tool developed by PNNL, along with the reliability and performance of commercial-grade software.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Modeling discontinuous dynamic recrystallization containing second phase particles in magnesium alloys utilizing phase field method

Pre-existing second phase particles can induce particle stimulated nucleation (PSN) and pinning effect on grain boundaries during DRX process of magnesium alloys. The interaction among pinning effect, PSN mechanism and strain-induced grain boundary migration (SIBM) nucleation mechanism imposes challenge in quantitative study of microstructure evolution. A phase field model that describes discontinuous dynamic recrystallization process considering second phase particles (PF-DDRX-SP) which mainly distribute along grain boundaries is developed. The second phase order parameter with diffusion interface is introduced to the free energy function. With considering the pinning effect of second phase particles on grain boundaries, the grain boundary energy couples the interaction between second phase and grains. The effect of second phase on dislocation density evolution is considered in the model, and the PSN and SIBM mechanism are both executed by limiting the nucleation position to the grain boundary and the second phase interface. For the validation, the initial topology consistent with initial microstructure is applied to the simulation of DRX kinetics and flow behavior of AZ80 magnesium alloy, and great reliability and accuracy of developed model is indicated. Furthermore, through the model, the microstructure evolution with different particle size and with different volume fraction of round particle during DRX process were simulated. Finally, the results show that particles with higher dispersion reduce the number of DRX nuclei, which indicate that second phase particles suppress the nucleation at grain boundaries because of the pinning effect. More importantly, the proposed model could also quantitatively predict the relationships between the parameters of second phase particles and DRX behaviors, and enable to optimize the initial second phase structure in a uniform grain structure during thermomechanical process.

36 MATERIALS SCIENCE↗

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗