Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “depth migration”

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 199 records · Page 11

Risk assessment of wellbore leakage during underground hydrogen storage

The expansion of renewable energy sources would require large-scale energy storage options to overcome the intermittent nature of these sources. Underground hydrogen storage (UHS) in depleted hydrocarbon reservoirs offers a scalable and practical energy storage solution. These reservoirs are chosen for their availability and large capacity, but the unique properties of hydrogen raise concerns about potential leakage pathways, particularly through wellbores. In this study, we develop and apply, for the first time, reduced-order models (ROMs) specifically designed for efficient leakage risk prediction in UHS systems operating in depleted hydrocarbon reservoirs. Using 3,000 high-fidelity simulation scenarios, we examine the influence of 11 key parameters, including reservoir and aquifer depths, wellbore permeability and porosity, initial saturations of water, oil and gas fractions (hydrogen, light, intermediate, and heavy hydrocarbons), reservoir pressure multiplier, and the aquifer-to-reservoir volume ratio, to simulate leakage behavior over a 1,000-year timescale. We train ROMs using a two-step classification-regression approach, achieving R 2 values exceeding 99 % across all targets. These ROMs effectively capture the leakage evolution and identify critical controls of leakage, guiding the design of mitigation strategies. Results indicate that gas leakage occurs in about 27 % of scenarios as early as five years post-operation, reaching volumes of up to 106 ft3. Oil leakage is less frequent (~17 %) and typically begins decades later. Our findings also show that hydrogen often migrates first, owing to its smaller molecular size and higher buoyancy, followed by heavier hydrocarbons. Over time, these heavier components contribute significantly to the total leaked volume, reinforcing the need for targeted monitoring and remediation strategies. Our analysis highlights that deeper storage reservoirs, shallower aquifers, and low-permeability wellbores significantly reduce leakage risks. In conclusion, this work offers a robust framework for risk-informed UHS deployment, supporting energy security through reliable large-scale hydrogen storage while safeguarding environmental integrity.

08 HYDROGEN↗

Prototype Hanford Barrier Performance Monitoring Report: Fiscal Year 2024

This report provides an annual update on performance monitoring of the Prototype Hanford Barrier (PHB) at the Hanford Site. The PHB is part of a long-term study that serves as the scientific basis for many of the evapotranspiration-capillary barrier designs currently used globally and for future engineered barrier designs planned for remedial actions over waste and demolition sites at the Hanford Site. The PHB allows us to identify potential future issues and develop better monitoring techniques before these engineered barriers are constructed. Surface barriers like the PHB are essential for preventing water infiltration and curbing the spread of contaminants to groundwater. Effective monitoring of soil moisture levels above and below these barriers is crucial given that performance metrics could span up to 1,000 years. From July 2023 to June 2024, the water flux (as measured by tipping buckets) through the 2-m-thick silt loam layer of the PHB – a fine silty material that stores water under high capillary tension – remained well below the 0.5-mm-per-year performance threshold, demonstrating its effectiveness in preventing water penetration. In 2024, degraded tipping bucket gauges were replaced to ensure the accuracy of future data. Additionally, neutron probes revealed that the wetting front from the rainy season only penetrated to a maximum depth of 1.2 m into the silt loam. This further demonstrated the barrier's efficiency, as the wetting front did not fully penetrate the 2-m-thick capillary barrier. The western gravel slope of the PHB, composed of Hanford Site pit gravel of varying sizes, demonstrated low flux rates. In contrast, the eastern riprap slope exhibited higher flux rates. Zhang (2017) found that the riprap side slope of the PHB had the highest drainage rates in January and the lowest in late summer or early fall, indicating significant summer evaporation. Drainage rates increased significantly under enhanced precipitation conditions, far exceeding the design criterion, which could lead to water infiltration into the waste zone. The study introduced the “edge effect,” where elevated drainage rates from the riprap may migrate laterally beneath the barrier, compromising its ability to isolate waste, and recommended expanding the barrier and conducting further research to mitigate this issue. This report provides a review on surface barrier edge effects and their potential impact to subsurface contaminant migration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Time Lapse Charged Wellbore Casing Controlled Source Electromagnetic Surveys for Monitoring Injected Carbon Dioxide

Carbon dioxide (CO2) injection for enhanced oil recovery (EOR) and long-term carbon storage extends the life of an oil field while keeping CO2 out of the atmosphere. Monitoring CO2 migration during EOR operations is critical to maximizing production and detecting unexpected CO2 migration quickly. In this project, we monitor CO2 migration using time-lapse charged well casing controlled source electromagnetics (CSEM) at an active EOR site undergoing CO2 injection. Since the electrical conductivity of the reservoir changes as CO2 displaces other fluids, this novel survey technique aims to monitor CO2 injection effectively and unobtrusively from the surface via electromagnetic field measurements. Of four time-lapse field campaigns spaced over two years, we have completed two. Each campaign consists of a CSEM survey supplemented by magnetotelluric and time-domain electromagnetic surveys. In the CSEM survey, a transmitter delivers alternating current to the subsurface at frequencies from 0.125 Hz to 16 Hz using steel well casings as electrodes. Charged well casing surveys detect deeper conductivity contrasts than surveys that use only surface electrodes, allowing CO2 movement at the depth of the reservoir to be imaged more clearly. We transmit on two well pairs. One pair is oriented in the east-west direction, and the other is oriented in the northeast-southwest direction. Receiver stations measure orthogonal components of the electric and magnetic fields. The measurements are repeated at the same locations in each field campaign. Data differences between campaigns provide information about changes in the electrical conductivity of the subsurface. We present the latest results from this ongoing project, including modeling, survey design, data processing techniques, and preliminary results from the field campaigns.

Magnetic and electrical methods, Geomagnetic induc↗

Wilson Corners, Solid Waste Management Unit 001(SWMU 01) 2023 Annual Long-Term Monitoring Report

This report presents a summary of the long-term monitoring (LTM) activities that occurred in 2023 at Wilson Corners, Solid Waste Management Unit 001, at Kennedy Space Center (KSC), Florida. Annual LTM of groundwater is being conducted at the site. Based on results from groundwater sampling activities performed during the 2019 through 2020 LTM reporting period and the 2020 and 2021 DPT groundwater sampling, it was determined that the LTM sampling plan was no longer meeting the goal of LTM because delineation was not verified and the installation of an air sparge (AS) system to treat the area of the High Concentration Plume was recommended. The AS System was installed in late 2022 and early 2023. System start-up activities were initiated in April 2023. Following system startup, several site wells required retrofitting to equip wellheads for withstanding the air pressure released from air sparge wells during system operation. Some site wells also required repair or abandonment, and replacement. Survey of location and top-of-casing of newly installed monitoring wells was combined with scheduled AS system survey activities and was completed in January 2024. The activities presented in this report include the February and April 2023 LTM monitoring well installations; March and April 2023 LTM and performance monitoring well water level gauging and sampling; November 2023 LTM well retrofits and repairs; a summary of December 2023 LTM well abandonments and installations (complete site well abandonment activities will be presented under a separate cover); and January 2024 LTM well survey. During the March and April 2023 sampling events, the low-flow sampling method was used, and samples were analyzed for a select list of volatile organic compounds. In March 2023, groundwater flow for the site was generally to the west was generally consistent with historical observations at the site. The Low Concentration Plume (LCP) continues to extend both horizontally and vertically beyond the terminal depth of the current monitoring well network. Data, inclusive of the 2023 LTM and baseline performance monitoring sampling events, indicate that the LCP encompasses an estimated 19.5 acres, compared to the 2021 LCP footprint, inclusive of the 2020 and 2021 DPT sampling events of 20.7 acres. The vertical extent of VOCs was historically delineated by monitoring wells screened greater than 48 feet below land surface (bls). The results from the three vertical extent monitoring wells screened below 48 feet bls that were sampled during the 2023 LTM indicate that groundwater vinyl chloride concentrations in these three wells are greater than the GCTL. As presented in the 2021 Long-Term Monitoring Report (NASA 2022), the KSCRT agreed to delay deeper investigations in this area to prevent the creation of additional pathways for vertical migration. Based on groundwater sampling activities performed in 2023, recommendations are to perform the next annual LTM sampling event, scheduled for April 2024 and to conduct quarterly performance monitoring of the AS System. The current selection of monitoring wells in the recommended 2024 LTM plan will provide an adequate data set for monitoring groundwater plume behavior; however, the LTM monitoring well network will be evaluated and refined based on 2024 LTM and year one performance monitoring data.

King Linnea↗

Efficient U(VI) Reduction and Sequestration by Ti 2 CT x MXene

Although reduction of highly mobile U(VI) to less soluble U(IV) has been long considered an effective approach to in situ environmental remediation of uranium, candidate reducing agents are largely limited to Fe-based materials and microbials. The importance of titanium-containing compounds in natural uranium ore deposits suggests a role for titanium in uranium migration. In this report, for the first time, a two-dimensional transition metal carbide, Ti 2 CT x , is shown to efficiently remove uranium via a sorption-reduction strategy. Batch experiments demonstrate that Ti 2 CT x exhibits excellent U(VI) removal over a wide pH range, with an uptake capacity of 470 mg g –1 at pH 3.0. The mechanism for U(VI) to U(IV) reduction by Ti 2 CT x was deciphered by X-ray absorption spectroscopy and diffraction and photoelectron spectroscopy. The reduced U(IV) species at low pH is identified as mononuclear with bidendate binding to the MXene substrate. At near-neutral pH, nanoparticles of the UO 2+x phase adsorb to the substrate with some Ti 2 CT x transformed to amorphous TiO 2 . A subsequent in-depth study suggests Ti 2 CT x materials may be potential candidates for permeable reactive barriers in the treatment of wastewaters from uranium mining. Finally, this work highlights reduction-induced immobilization of U(VI) by Ti 2 CT x MXene including a pH-dependent reduction mechanism that might promote applications of titanium-based materials in the elimination of other oxidized contaminants.

36 MATERIALS SCIENCE↗

Benchmark Exercise Report for Experimental Study of Bubble Scrubbing in Sodium Pool

Mechanistic source term (MST) analyses are likely to be an important part of advanced reactor licensing applications. For the purpose, an MST analysis code applicable to newly introduced advanced reactors, called SRT (Simplified Radionuclide Transport) code, has been developed by Argonne National Laboratory. SRT can track overall behaviors of radionuclides especially in metal fuel-based sodium fast reactors (SFRs) and microreactors. Throughout the simulation, migration inside fuel pins before failure, interaction with coolant (for SFR), removal/leakage in cover gas and containment (or confinement), and environmental dose impacts are considered alongside radioactive decay for short-lived nuclides. Among the postulated process, pool scrubbing phenomenon, especially under sodium pool condition, has been identified as high importance with limited supportive data. The phenomenon plays a crucial role in assessing the degree of radiological impacts as radioactive aerosols or vapors are efficiently and effectively removed during the process. To provide validation basis for SRT in assessing pool scrubbing performance inside sodium pools, the University of Wisconsin-Madison performed tests including extensive parametric effects. Separate effect tests were conducted to directly evaluate the SRT models and to estimate degree of contribution by each contributing factor. Specifically, bubble size, aerosol size, aerosol density, aerosol concentration, pool depth, system temperature, and bubble swarm effects were considered. According to the parametric effects, decontamination performance enhances with decreasing bubble size, large density, and deeper pool height. Aerosol concentration provides no effect for the whole range of interest, and pool temperature variation shows minor effects under the considered temperature condition. When multiple bubbles are injected generating a bubble swarm condition, DF performance further enhances by bubble interactions and turbulence characteristics. The measurement shows the exceptional importance of aerosol size range considered, with the lowest decontamination, where most radionuclides are expected to escape.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Enhancing trap and fault seal analyses by integrating observations from HR3D seismic data with well logs and conventional 3D seismic data, Texas inner shelf

An understanding of trap and fault seal quality is critical for assessing hydrocarbon prospectivity. To achieve this, modern analytical techniques leverage well data and conventional industry-standard 3D seismic data to evaluate the trap, and any faults displacing the reservoir and top seal intervals. Above all, geological interpretation provides the framework of trap and fault seal analyses, but can be hindered by the data resolution, quality and acquisition style of the conventional seismic data. Furthermore, limiting the analysis to only the petroleum system at depth may lead to erroneous perceptions because interpreting overburden features, such as shallow faults or gas chimneys, can provide valuable observations with respect to container performance, and can to help validate trap and fault seal predictions. Here, a supplement to conventional 3D data are high-resolution 3D seismic (HR3D) data, which provide detailed images of the overburden geology. This study utilizes an HR3D seismic volume in the San Luis Pass area of the Texas inner shelf, where shallow fault tips and a sizeable gas chimney are interpreted over an unsuccessful hydrocarbon prospect. Static post-drill fault seal and trap analyses suggest that the primary fault displacing the structural closure could have withheld columns of gas c. 100 m high, but disagree with our HR3D seismic interpretations and dry-well analyses. From our results, we hypothesize that tertiary gas migration through fault conduits reduced the hydrocarbon column in the prospective Early Miocene reservoir, and may have resulted from continued movement along the intersecting faults. Overall, this study reinforces the importance of understanding the overburden geology and geohistory of faulted prospects, and demonstrates the utility of pre-drill HR3D acquisition when conducting trap and fault seal analyses.

3D seismic↗

Electrical Resistivity Tomography Monitoring of In Situ Soil Flushing at the Hanford 100-K East Area: 100KE Soil Flushing Monitoring

Former operations in the 100-KE Area at the Hanford Site resulted in vadose zone hexavalent chromium that provides a source of groundwater contamination. Previous efforts to remediate vadose zone chromium involved excavation and offsite treatment of contaminated soils. Although much of the chromium was removed, contamination still exists in native soil between the water table and the bottom of the pit, which has since been backfilled. In-situ soil flushing was tested at the 100-KE Area to accelerate the removal of the remaining hexavalent chromium in the vadose zone. Soil flushing works by applying clean flush water at the surface, which mobilizes and transports chromium as it migrates downward to the water table. Once the chromium reaches the water table, it is hydraulically contained, extracted, and treated through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. Therefore, performance monitoring requires and understanding of flush water migration through the vadose zone. It is challenging to comprehensively assess soil flushing performance via borehole access due to the limited volume of investigation afforded by a given wellbore. As an alternative, 3D time-lapse electrical resistivity tomography (ERT) was tested as a method of monitoring the distribution of flush water over time within the vadose zone. ERT is a method of imaging the bulk electrical conductivity of the subsurface, which is highly dependent on water saturation levels in unconsolidated and unsaturated sediments. Thus, the timing and location of changes in bulk conductivity caused by soil flushing can be used to infer flush water migration pathways. In this report we present results of 3D time-lapse ERT imaging during two separate soil flushing campaigns conducted during the spring and summer months of 2022 and 2023. Results show generally that: 1) Pit backfill materials appear to nominally have larger permeability than native soils. Consequently, the boundary between pit backfill and native soil had a significant impact on flush water migration, causing some flush water to migrate along the pit boundary to the bottom of the pit. 2) Non-uniform flows, likely caused by variations in hydrogeologic properties, developed in the pit backfill materials, resulting in uneven flush water distribution on the southern margin of the soil flushing zone. 3) Redistribution of water at the interface between backfill materials and Hanford formation materials likely facilitated elevated and uniform distribution of flush water within native soils beneath the deeper parts of the pit boundary beneath the center of the flush zone, which presumably overlie soils with elevated chromium contamination. These areas appear to have been infiltrated by higher volumes of flush water than the northern and southern margins of the flush zone. 4) In comparison to 2022, high flush water application rates significantly improved flush water distribution throughout the target flushing zone. 5) Imaging resolution was limited to a depth of approximate 20 meters below the land surface, due primarily to limitations on the lateral extent of the surface ERT array. The ~10 m region of the vadose zone between approximately 130 meters elevation and the water table at approximately 120 meters elevation was unresolved.

100 Area↗

Lithologic controls on microfracturing from legacy underground nuclear explosions

Detection and verification of underground nuclear explosions (UNEs) can be improved with a better understanding of the nature and extent of explosion-induced damage in rock and the effect of this damage on radionuclide migration. Much of the previous work in this area has focused on centimeter- to meter-scale manifestations of damage, but to predict the effect of damage on permeability for radionuclide migration, observations at smaller scales are needed to determine deformation mechanisms. Based on studies of tectonic deformation in tuff, we expected that the heterogeneous tuff layers would manifest explosion-induced damage differently, with welded tuffs showing more fractures and nonwelded tuffs showing more deformation bands. In comparing post-UNE samples with lithologically matched pre-UNE equivalents, we observed damage in multiple lithologies of tuff through quantitative microfracture densities. We find that the texture (e.g., from deposition, welding, alteration, etc.) affects fracture densities, with stronger units fracturing more than weaker units. While we see no evidence of expected deformation bands in the nonwelded tuffs, we do observe, as expected, much larger microfracture densities at close range (<50 m) to the explosive source. We also observe a subtle increase in microfracture densities in post-UNE samples, relative to pre-UNE equivalents, in all lithologies and depths. The fractures that are interpreted to be UNE-induced are primarily transgranular and grain-boundary microfractures, with intragranular microfracture densities being largely similar to those of pre-UNE samples. This work has implications for models of explosion-induced damage and how that damage may affect flow pathways in the subsurface.

58 GEOSCIENCES↗

Wilson Corners SWMU 001 2015 Annual Long Term Monitoring Report Kennedy Space Center, Florida

This document presents the findings of the 2015 Long Term Monitoring (LTM) that was completed at the Wilson Corners site, located at the National Aeronautics and Space Administration John F. Kennedy Space Center, Florida. The objectives of the 2015 LTM event were to evaluate the groundwater flow direction and gradient, to monitor the vertical and horizontal extent of the volatile organic compounds (VOCs; including the upgradient and sidegradient extents, which are monitored every five years), and to monitor select locations internal to the dissolved groundwater plume. The 2015 LTM event included several upgradient and sidegradient monitoring wells that are not sampled annually to verify the extent of VOCs in this portion of the site. The December 2015 LTM groundwater sampling event included, depth to groundwater measurements, 40 VOC samples collected using passive diffusion bags, and one VOC sample collected using low-flow techniques. Additionally, monitoring well MW0052DD was overdrilled and abandoned using rotasonic drilling techniques. The following conclusions can be made based on the 2015 LTM results: groundwater flow is generally to the west with northwest and southwest flow components from the water table to approximately 55 feet below land surface (ft BLS); peripheral monitoring wells generally delineate VOCs to groundwater cleanup target levels (GCTLs) except for monitoring wells MW0088, MW0090, MW0095, and NPSHMW0039, which had vinyl chloride (VC) concentrations near the GCTL and MW0062, which had trichloroethene (TCE), cis-1,2-dichloroethenen (cDCE), and VC concentrations above natural attenuation default concentrations (NADCs); VOCs in interior downgradient wells generally fluctuate within historic ranges except for monitoring wells in the north-northwest portion of the site, which have increasing VC concentrations indicating potential plume migration and expansion; Historically, the vertical extents of the VOCs were delineated by monitoring wells screened greater than 60 ft BLS (MW0083 through MW0086, and MW0078). The 2015 LTM results indicate that concentrations of daughter product cDCE is greater than the NADC in MW0078 and that cDCE and VC are greater than NADCs in MW0130. TCE was greater than the GCTL in monitoring well MW130 and not detected above method detection limits (9 micrograms per Liter) in monitoring well MW0078. No GCTL exceedances were identified in monitoring wells MW0083 or MW0086; the dissolved plume footprint appears generally stable, though not fully delineated by monitoring well data in the northeast portion of the site; and 2015 LTM results generally support the existing Conceptual Site Model. Geosyntec recommends modifying the LTM program, collecting a verification sample from monitoring well MW0062, and performing a direct push technology instigation in the northnortheastern portion of the site. A cluster of monitoring wells (MW0132 [2 to 12 ft BLS], MW0133 [15 to 25 ft BLS], and MW0134 [29 to 34 ft BLS]) is proposed to delineate impacts in the northnortheast portion of the site. Geosyntec recommends installation of a vertical extent well in the center of the site (Hot Spot 2 area) post-remediation implementation.

groundwater↗

Pyroscopegridding: Geo-Leo Aerosol Data Fusion (an Open-Source Package)

The retrieval of aerosol optical depths (AODs) from sun-synchronous polar orbiting satellites, such as MODISs, VIIRSs, OMI, TROPOMI, etc., has been widely adopted as a method for obtaining information regarding particulate matter (PM) and related atmospheric processes. However, the advent of recently launched geostationary satellites, such as GOES-16/17/18, Himawari-8/9, and Meteosat Third Generation (MTG), has led to an increased temporal resolution of AOD observations (order of 10 minutes), resulting in typically one or more images per hour during daylight hours, compared to the once-per-day observations obtained from LEO satellites. By integrating these observations, the diurnal cycle of global AOD can be characterized at local, regional, and global scales. The scientific community is still examining the novel data from geostationary satellite observations and evaluating methods for effectively merging these observations with differing spatial and temporal resolutions. This presents a significant ""Big Data"" challenge, encompassing not only data storage, but also data discoverability, accessibility, and migration within cloud computing environments. This study presents our attempts at fusing Level 2 aerosol data from six satellites, three of which are geostationary (GOES-16/17 and Himawari-8) and three of which are polar orbiting (TERRA/MODIS, AQUA/MODIS, and SNPP-VIIRS), using the Dark Target aerosol retrieval algorithm. The ability to fuse remote sensing products on demand into desired temporal and spatial domains empowers researchers and practitioners to more efficiently work with satellite and sensor data. It is our hope that through making our open-source package and accompanying functionality available, the scientific community will have improved access to aerosol data processing resources.

Jennifer Wei↗

Plate Margin Deformation and Active Tectonics Along the Northern Edge of the Yakutat Terrane in the Saint Elias Orogen, Alaska and Yukon, Canada

The northwest directed motion of the Pacific plate is accompanied by migration and collision of the Yakutat terrane into the cusp of southern Alaska. The nature and magnitude of accretion and translation on upper crustal faults and folds is poorly constrained, however, due to pervasive glaciation. In this study we used high-resolution topography, geodetic imaging, seismic, and geologic data to advance understanding of the transition from strike-slip motion on the Fairweather fault to plate margin deformation on the Bagley fault, which cuts through the upper plate of the collisional suture above the subduction megathrust. The Fairweather fault terminates by oblique-extensional splay faulting within a structural syntaxis, allowing rapid tectonic upwelling of rocks driven by thrust faulting and crustal contraction. Plate motion is partly transferred from the Fairweather to the Bagley fault, which extends 125 km farther west as a dextral shear zone that is partly reactivated by reverse faulting. The Bagley fault dips steeply through the upper plate to intersect the subduction megathrust at depth, forming a narrow fault-bounded crustal sliver in the obliquely convergent plate margin. Since . 20 Ma the Bagley fault has accommodated more than 50 km of dextral shearing and several kilometers of reverse motion along its southern flank during terrane accretion. The fault is considered capable of generating earthquakes because it is linked to faults that generated large historic earthquakes, suitably oriented for reactivation in the contemporary stress field, and locally marked by seismicity. The fault may generate earthquakes of Mw <= 7.5.

Bruhn, Ronald L.↗

Alteration of Granite Surface with the Formation of Calcium-Silicate-Hydrate by Alkaline Groundwater - 19165

Under the hyper alkaline Ca-rich condition caused by cementitious materials used for the geological disposal system, host rocks containing various silicate minerals will alter with the dissolution of the minerals, and the expansion of flow paths may accelerate the migration of radionuclides. On the other hand, the formation of secondary minerals under such condition will contribute to the confinement for radionuclides. As a typical secondary mineral, we focused on calcium-silicate-hydrate (C-S-H). In particular, this study examined the C-S-H formation on the granite surface under the hyper-alkaline Ca-rich condition by the immersion tests and the estimation based on the film theory. In the immersion tests, granite chips immersed in a Ca(OH)2 or a NaOH solution adjusted pH to 12.3 over 60 days. From the results, the Si concentration in the Ca(OH)2 solution was much lower than that in the NaOH solution, and the Ca concentration decreased in the case of Ca(OH)2 solution. These indicates the C-S-H formation around the granite surface. Furthermore, assuming that the C-S-H layer was formed on the granite surface, the depth point of the C-S-H layer, where C-S-H was formed by reacting Ca ions with silicic acid leached from granite, was estimated by the film theory. The estimation result showed the C-S-H layer would grow with the C-S-H formation at the surface contacting granite and C-S-H. These results in this study suggest that micro flow-paths in the host rocks will be clogged by the continuous growth of the C-S-H layer on the surface of silicate minerals. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ma'adim Vallis Estuarine Delta in Elysium Basin and Its Relevance as a Landing Site for Exobiology Exploration on Mars

The debouche of Ma'adim Vallis in the Elysium Basin generated a transitional transported sediment structure, which planimetric shape is controlled by the enclosing topography of a deep reentrant gulf of the Basin into the highland. We defined it as an estuarine delta. The location and the importance of this estuarine delta is supported by the theoretical model of graded profile constructed for Ma'adim Vallis, and by two approaches: (i) the reconstruction of Ma'adim Vallis downstream course from Gusev to Elysium Basin, and (ii) the survey of the sediment deposit in the alleged estuary. The longitudinal graded profile of Ma'adim Vallis finds its base-level in the Elysium Basin, at a about 1000 m elevation, which is in agreement with the observed Basin shoreline. This model is supported by observational evidence of flow between the northern rim of Gusev crater, and the Elysium Basin shoreline. This downstream course of Ma'adim Vallis can be divided into three hydrogeologic regions. into three hydrogeologic regions. (a) The first region is a flooded plain (Zephiria Mensae), consisting in chaotic terrain formed by highland rocks, and disintegrated lava of the western flank of Apollinaris. Morphologic indicators of the flood process are: (1) the sediment deposit over the Gusev crater northern rim that reflects the overspilling of the crater-lake water through a 40-km wide gap provided by an ancient impact crater, (2) the tear-drop shaped feature on the northeastern flank of Apollinaris Patera, and (3) the chaotic terrain that suggest the emergence of ground water generated by the seepage of the crater lake through high-permeable broken rampart material. This underground water circulation sustained by the hydrostatic pressure of the crater-lake has likely generated a hydrothermal system in the volcanic environment of Apollinaris Patera. The stratigraphy of the flooded area is identified as Hesperian age, with occurrences of Noachian hilly individual features, and as Amazonian flooded plain and chaotic material, (b) The second region is located on the western flank of Apollinaris Patera. It is surrounded by relics of deep valleys that suggest a former downstream course of Ma'adim Vallis. The geologic setting of this region (Lucus Planum) is interpreted to be an Amazonian formation composed by the middle and lower members of the Medusae Formation., c) The third region corresponds to the convergence of the west and east branches of Ma'adirn Vallis into a deep re-entrant wide gulf that penetrates about 100 km into the highland. This topographic depression is delineated by the 1000 in elevation contour. This gulf has formed an estuarine configuration centered at 3S/190W within the Elysium Basin. This configuration has favored the formation of a estuarine sedimentary delta, because of topographically controlled lateral migration. This estuarine structure is strongly dominated by the incoming supply of Ma'adim Vallis fluvial sediment extracted from Zephiria Mensae and Lucus Planum. The obtuse-angle geometry of the estuary increases the sedimentation rate, which is higher than in the course of the channel. The sediment deposition process is governed by the estuarine water circulation. The inflowing loaded fluvial water enters the estuary as a bottom current, and mixes with the relatively less-loaded water of the receiving basin. When they mixed. the inflowing fluvial material, and the landward basin circulating water generate an accumulation of highly-diversified estuarine deposit stratification. This accumulation of material is mostly centered in the transitional zone of the delta. The sediment trapping efficiency of the estuary is function of the energy balance between the inflowing fluvial water, and the ingoing basin current. The submergence of the delta by the rising of the water-level increases the estuary water-depth, and consequently the sediment entrapment is favored. The locus of sediment accumulation moves landward in the zone of inflowing fluvial water. This results in the rising of the channel base-level, thus in the increase of the length of the longitudinal graded-profile. The sediment deposit facies of the zone A shows a generally smooth surface. The longitudinal deposit is bordered by alluvial terraces that reflect the variations of the channel level. The waning of the Elysium Basin caused the erosion of the Basin estuarine zone by small channels, this episode being characterized by dissected tear-drop shaped mesa-like morphologies in the delta. Our estuarine delta model predicts a lithostratigraphic depositional sequence associated with the water submergence and the transgression of Elysium Basin. The thickness of the estuarine sediment corresponds to the Elysium Basin levels changes relatively to the bed floor of the estuary, The depositional sequence of Ma'adim Vallis are described: (1) a pro-current filled region (A), where fluvial are longitudinally accumulated by the inflowing water, (2) inverse current from Elysium Basin (B), where fluvial and lacustrine sediments are accumulated, and (3) zone of current equilibrium (C), where the sediments are distributed as a shoreline at the boundary of the estuarine delta. The estuary sedimentology dynamics collects and keeps the record of the geologic unit material crossed by Ma'adim Vallis, and those of the lakebed deposit of Elysium Basin. The predicted mixed stratigraphic sequence from fluvial and lacustrine sediment makes this site an exceptional environment to concentrate potential multi-origin biologic records. We envision four possible strategies to explore this sedimentologic record: (1) longitudinal surface and subsurface traverses in region A to investigate outcrop levees, (2) exploration of the mesa walls in region B, (3) deep drilling hole lodging of the sequential deposits in the zones A and B, and (4) surface and subsurface exploration of the shoreline delta. The expected results for each of these strategies are: (1) in the deepest layers of region A are predicted frequent and abundant coarse material, sandy lenses lamination grading downward from sand to cobbles. Volcanic debris from the Noachian crustal Plateau unit material, hydrothermal altered rocks, carbonates, Hesperian and possibly Amazonian volcanic material, from Apollinaris Patera, altered rocks and carbonates from Zephiria Mensae are expected. As a favorable environment for inception of life, possible biological records are expected in transported rock, (2) At the surface, and subsurface (<=100 m), large deposits sandy to silted material from Elysium paleolake basin mixed with fine-grained sediments from Ma'adim Vallis are expected mostly in the upstream part of region B, (3) on the shoreline of the estuarine delta, abundant fine material from Elysium paleolake basin (evaporites, carbonates), mostly Amazonian in age are expected. The Ma'adini estuary is a favorable landing site for all the above mentioned science aspects, and .for its location. The site lies near the equator, which is favorable for the rover solar power supply, and at 1000m elevation, which is a favorable configuration for the descent system braking. Another advantage is the extent of the area of high scientific interest (33,000 sq km), which provides a good ellipse, and potential long study traverses.

Grin, E. A.↗

Model for radon diffusion through the lunar regolith.

Description of a model for radon diffusion through the lunar regolith in which the atom migrates by random walk. The regolith is represented by a system of randomly oriented baffles in which the mean distance which the atom travels between two collisions takes on the role of a mean free path. The effective mean time between two collisions depends on two entities: the actual mean time-of-flight and the mean sticking time on grain surfaces for one collision. The latter depends strongly on the temperature and the heat of adsorption of radon on regolith materials. Both the mean free path as well as the heat of adsorption are either poorly known or unknown for the lunar regolith; hence these quantities are treated as free parameters. Because of the greatly different mean lifetimes against radioactive decay of Rn219, Rn220, and Rn222, the regolith acts as a powerful 'filter' for these species. Rn222 escape is significant (32%) even for a mean free path of 1 micron, a heat of adsorption of 7.0 kcal/mole and a regolith depth of 4 m. Calculations of radon escape from a 4 m thick regolith, using mean free paths of 1, 10, and 80 microns and heats of adsorption of 4.0, 5.2, and 7.0 kcal/mole show that the Rn222/Rn220 escape ratio can be as small as 7.7 and as large as, or larger than 47. The small value of 7.7 is of particular interest, because it is nearly equal to the escape ratio inferred by Turkevich et al. (1970) from their Surveyor 5 results.

Friesen, L. J.↗

Design and Modeling of a Coded-Interleaving System in the Presence of Fading

Future deep space links are migrating towards higher frequency bands such as Ka band and optical. These links are susceptible to non linear effects (e.g. cloud coverage, atmospheric turbulence, scintillation, antenna mis-pointing, jitter, relative motion, etc.) that cause various degrees of fading loss over various time scales. We consider a channel experiencing block fading where the channel coherence time is a lot longer than the length of a codeword. As long as the fading duration is not excessive, an interleaver can be used in a coded-communication system that randomizes the fades among symbols of multiple codewords. This achieves an “averaging” effect that spreads the signal energy more evenly among the codewords, and improves the likelihood that each of them can be decoded correctly. For simplicity, we consider a block interleaver where adjacent symbols of a codeword have a fixed delay (or separation) with each other, and investigate the error-rate performance of a coded-interleaving system (with fixed and finite interleaving depth) in the presence of block fading. Unlike prior research which assumes theoretical channel coherence models, we aim to build upon our analysis using empirical spacecraft signal measurements obtained with the Deep Space Network, and to use the fundamental theories of probability and statistics (e.g. variations of Law of Large Numbers) to quantify the error-rate performance of the coded-interleaving system. We propose to utilize these results in the design of a two-dimensional coded-interleaving system for fading channels.

Sanchez-Net, Marc↗

Deformation of pores in response to uniaxial and hydrostatic stress cycling in Marcellus Shale: Implications for gas recovery

One of the main challenges during gas production from shale reservoirs is low recovery rate. One contributing factor to this outcome is an insufficient understanding of pore systems, especially pore behaviour following changes in reservoir conditions or resulting from gas production practices. Because the pressure in the producing well can be controlled, understanding the effects of pressure variation on the pore size distribution and methane trapping is necessary to help design optimal conditions to improve the gas recovery rate. This work is the first systematic study of sub-millimeter pore deformation in shale caused by uniaxial and hydrostatic stress up to 100 MPa. Overmature samples from the Middle Devonian Marcellus Shale were analyzed using neutron scattering (SANS and USANS) techniques to interpret the response of nanopores to stress cycling of magnitude and duration compatible with the hydraulic fracturing treatments. Experiments reported here are performed at a series of uniaxial pressures up to 100 MPa and at hydrostatic pressures of deuterated methane 0 and 50 MPa. Since at the original depth of the shale samples’ burial of 2184 m the hydrostatic pressure is approximately 27 MPa and the lithostatic pressure is about 55 MPa, the experimental conditions reasonably well simulate the reservoir pressure regime. Our SANS and USANS results show that different pore sizes are affected by uniaxial stress in different ways. Specifically, in the pore size range from 1 nm to 800 nm, a decrease of pore density with pressure is observed, with the most depleted being mesopores of about 100 nm in diameter. The observed decrease is likely related to deformation of kerogen, followed by a loss of pore nano-volume, as well as methane trapped in the micropores. For pores larger than 5 μm, USANS data suggest that the negative trend is reversed at above 74 MPa, and the number density of large macropores may increase with increased stress even above the original value. The increased number of macropores at high pressure may create new interconnected conduits for gas migration, resulting in a better recovery rate. Another important finding of this study is an irreversible rearrangement of pore size distribution taking place after pressure cycling. Furthermore, this irreversible reorganization of pore size distribution should be taken into account during management of well production to maximize recovery rate.

03 NATURAL GAS↗

Thermal Tides During the 2001 Martian Global-Scale Dust Storm

The 2001 (Mars Year 25) global dust storm radically altered the dynamics of the Martian atmosphere. Using observations from the Thermal Emission Spectrometer onboard the Mars Global Surveyor spacecraft and Mars WRF general circulation model simulations, we examine the changes to thermal tides and planetary waves caused by the storm. We find that the extratropical diurnal migrating tide is dramatically enhanced during the storm, particularly in the southern hemisphere, reaching amplitudes of more than 20 K. The tropical diurnal migrating tide is weakened to almost undetectable levels. The diurnal Kelvin waves are also significantly weakened, particularly during the period of global expansion at Ls=200deg-210deg. In contrast, the westward propagating diurnal wavenumber 2 tide strengthens to 4-8 K at altitudes above 30km. The wavenumber 1 stationary wave reaches amplitudes of 10-12 K at 50deg-70degN, far larger than is typically seen during this time of year. The phase of this stationary wave and the enhancement of the diurnal wavenumber 2 tide appear to be responses to the high-altitude westward propagating equatorial wavenumber 1 structure in dust mixing ratio observed during the storm in previous works. This work provides a global picture of dust storm wave dynamics that reveals the coupling between the tropics and high-latitude wave responses. We conclude that the zonal distribution of thermotidal forcing from atmospheric aerosol concentration is as important to understanding the atmospheric wave response as the total global mean aerosol optical depth.

Thermal tides↗