Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LAVA”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Performance Characteristics of a Boundary Layer Ingesting Inlet - Distortion Tolerant Fan Through CFD Code Coupling

A method of constructing a high-fidelity simulation using two NASA CFD codes, Launch Ascent and Vehicle Aerodynamics (LAVA), which is the general-purpose solver, and a turbomachinery specific code NASA-TURBO (TURBO), will be presented. Validation of the code coupling is performed through comparison to experimental data obtained in the 2017 wind tunnel test of NASA’s Boundary Layer Ingesting Inlet - Distortion Tolerant Fan (BLI2DTF). For the simulation, the external flow field is computed by LAVA utilizing a steady-state Reynolds-averaged Navier-Stokes (RANS) methodology. Data from LAVA is exchanged with TURBO, which simulated the propulsor utilizing an unsteady-RANS methodology. A sensitivity analysis is performed using averaged and instantaneous data transfer between the two codes, as well as varying the amount of simulation time TURBO computes between updates. The results of the coupled simulations appear to replicate the BLI2DTF experiment correctly. Static pressure taps leading into and around the nacelle, total pressure ratio and efficiency of the propulsor, and profiles of total temperature, pressure, and Mach number at the rear measurement plane, aft of the propulsor, are in good agreement with experimental data. These results provide confidence that the coupled simulation is behaving in a physically realizable manner.

CFD↗
NASA NTRS2024

Applications of ArcticDEM for measuring volcanic dynamics, landslides, retrogressive thaw slumps, snowdrifts, and vegetation heights

Topographical changes are of fundamental interest to a wide range of Arctic science disciplines faced with the need to anticipate, monitor, and respond to the effects of climate change, including geohazard management, glaciology, hydrology, permafrost, and ecology. This study demonstrates several geomorphological, cryospheric, and biophysical applications of ArcticDEM – a large collection of publicly available, time-dependent digital elevation models (DEMs) of the Arctic. Our study illustrates ArcticDEM's applicability across different disciplines and five orders of magnitude of elevation derivatives, including measuring volcanic lava flows, ice cauldrons, post-failure landslides, retrogressive thaw slumps, snowdrifts, and tundra vegetation heights. We quantified surface elevation changes in different geological settings and conditions using the time series of ArcticDEM. Following the 2014–2015 Bárðarbunga eruption in Iceland, ArcticDEM analysis mapped the lava flow field, and revealed the post-eruptive ice flows and ice cauldron dynamics. The total dense-rock equivalent (DRE) volume of lava flows is estimated to be (1431 ± 2) million m 3 . Then, we present the aftermath of a landslide in Kinnikinnick, Alaska, yielding a total landslide volume of (400 ± 8) × 103 m 3 and a total area of 0.025 km 2 . ArcticDEM is further proven useful for studying retrogressive thaw slumps (RTS). The ArcticDEM-mapped RTS profile is validated by ICESat-2 and drone photogrammetry resulting in a standard deviation of 0.5 m. Volume estimates for lake-side and hillslope RTSs range between 40,000 ± 9000 m 3 and 1,160,000 ± 85,000 m 3 , highlighting applicability across a range of RTS magnitudes. A case study for mapping tundra snow demonstrates ArcticDEM's potential for identifying high-accumulation, late-lying snow areas. The approach proves effective in quantifying relative snow accumulation rather than absolute values (standard deviation of 0.25 m, bias of −0.41 m, and a correlation coefficient of 0.69 with snow depth estimated by unmanned aerial systems photogrammetry). Furthermore, ArcticDEM data show its feasibility for estimating tundra vegetation heights with a standard deviation of 0.3 m (no bias) and a correlation up to 0.8 compared to the light detection and ranging (LiDAR). The demonstrated capabilities of ArcticDEM will pave the way for the broad and pan-Arctic use of this new data source for many disciplines, especially when combined with other imagery products. The wide range of signals embedded in ArcticDEM underscores the potential challenges in deciphering signals in regions affected by various geological processes and environmental influences.

Chunli Dai↗
DOE OSTI2023

Redox and mineral controls on Fe and Ti isotopic fractionations during calc-alkaline magmatic differentiation

Titanium and Fe isotopic compositions of lavas from a calc-alkaline differentiation suite and corresponding mineral separates from the Rindjani Volcano, Indonesia show that Fe and Ti isotopic fractionations between minerals and melts are lower than those recorded in other suites at all stages of differentiation. Here, the limited isotopic fractionation for Ti is likely due to low-Ti magnetite and clinopyroxene being the dominant carriers of Ti in Rindjani lavas, as these minerals are thought to have limited equilibrium Ti isotopic fractionation relative to silicate magmas. Other magmatic differentiation suites controlled by removal of Ti-rich magnetite and characterized by a lesser role of clinopyroxene have larger Ti isotopic fractionations. This effect is an indirect consequence of the elevated Fe 3+ /Fe 2+ ratio of calc-alkaline magmas such as Rindjani, which promotes Fe 3+ incorporation into magnetite at the expense of Fe 2+ -Ti 4+ pairs, such that increased oxygen fugacity will subdue Ti isotopic fractionation in global magmatic series. Similarly, we find negligible Fe isotopic fractionation in Rindjani bulk rocks and mineral separates, unlike previous studies. This is also likely due to the oxidized nature of the Rindjani differentiation suite, which leads to similar Fe 3+ /Fe 2+ ratios in melt and minerals and decreases overall mineral-melt Fe fractionation factors. Paired Ti and Fe isotopic analyses may therefore represent a powerful tool to assess oxygen fugacity during differentiation, independent from Fe 3+ determinations of erupted samples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Volatiles and Redox Along the East African Rift

Abstract The upper mantle under the Afar Depression in the East African Rift displays some of the slowest seismic wave speeds observed globally. Despite the extreme nature of the geophysical anomaly, lavas that erupted along the East African Rift record modest thermal anomalies. We present measurements of major elements, H 2 O, S, and CO 2 , and Fe 3+ /ΣFe and S 6+ /ΣS in submarine glasses from the Gulf of Aden seafloor spreading center and olivine‐, plagioclase‐, and pyroxene‐hosted melt inclusions from Erta Ale volcano in the Afar Depression. We combine these measurements with literature data to place constraints on the temperature, H 2 O, andfO 2 of the mantle sources of these lavas as well as the initial and final pressures of melting. The Afar mantle plume is C/FOZO/PHEM in isotopic composition, and we suggest that this mantle component is damp, with 852 ± 167 ppm H 2 O, not elevated infO 2 compared to the depleted MORB mantle, and has temperatures of ∼1401–1458°C. This is similar infO 2 and H 2 O to the estimates of C/FOZO/PHEM in other locations. Using the moderate H 2 O contents of the mantle together with the moderate thermal anomaly, we find that melting begins at around 93 km depth and ceases at around 63 km depth under the Afar Depression and at around 37 km depth under the Gulf of Aden, and that ∼1%–29% partial melts of the mantle can be generated under these conditions. We speculate that the presence of melt, and not elevated temperatures or high H 2 O contents, are the cause for the prominent geophysical anomaly observed in this region.

Geochemistry & Geophysics↗
DOE OSTI2022

Trans-crustal structural control of CO 2 -rich extensional magmatic systems revealed at Mount Erebus Antarctica

Erebus volcano, Antarctica, with its persistent phonolite lava lake, is a classic example of an evolved, CO 2 -rich rift volcano. Seismic studies provide limited images of the magmatic system. Here we show using magnetotelluric data that a steep, melt-related conduit of low electrical resistivity originating in the upper mantle undergoes pronounced lateral re-orientation in the deep crust before reaching shallower magmatic storage and the summit lava lake. The lateral turn represents a structural fault-valve controlling episodic flow of magma and CO 2 vapour, which replenish and heat the high level phonolite differentiation zone. This magmatic valve lies within an inferred, east-west structural trend forming part of an accommodation zone across the southern termination of the Terror Rift, providing a dilatant magma pathway. Unlike H 2 O-rich subduction arc volcanoes, CO 2 -dominated Erebus geophysically shows continuous magmatic structure to shallow crustal depths of < 1 km, as the melt does not experience decompression-related volatile supersaturation and viscous stalling.

58 GEOSCIENCES↗
DOE OSTI2021

Safely probing the chemistry of Chernobyl nuclear fuel using micro-focus X-ray analysis

Detailed chemical analysis of the solidified molten fuel still residing in the stricken Chernobyl reactor unit 4 are inferred using multi-modal micro-focus X-ray analysis of a low-radioactivity proxy. A fascinating mixture of molten UO 2 , nuclear fuel cladding, concrete, stainless steel and other nuclear reactor components, these materials behaved like lava, solidifying to form a complex, highly radioactive glass-ceramic. Using element-specific chemical probes (micro-X-ray fluorescence and X-ray absorption spectroscopy), coupled with micro-diffraction analysis, the crystalline phase assemblage of simulants of these heterogeneous materials was established, which included “chernobylite” and a range of compositions in the (U 1-x Zr x )O 2 solid solution. Novel insight to nuclear accident fuel chemistry was obtained by establishing the oxidation state and local coordination of uranium not only in these crystalline phases, but uniquely in the amorphous fraction of the material, which varied depending on the history of the nuclear lava as it flowed through the reactor. This study demonstrates that micro-focus X-ray analysis of very small fractions of material can yield rich chemical information, which can be applied to nuclear-melt down materials to aid decommissioning and nuclear fuel management at nuclear accident sites.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Zinc isotope constraints on the cycling of carbon in the Bermuda mantle source

Volatile recycling and storage in the mantle transition zone (MTZ) is important for the refertilization of the upper mantle and is associated with the generation of high-µ (HIMU, where µ is 238 U/ 204 Pb) mantle. One way to probe the MTZ and the processes associated with mantle convection is to sample lavas that originate from the shallow mantle and were contaminated by upwelling from the MTZ, such as at the previously proposed shallow plume of Bermuda. Here we present the first δ 66 Zn isotopic compositions of Bermuda silica-undersaturated and silica-saturated lavas to explore the origin of the carbon-rich lithologies and the genesis of the large seamount found in the western North Atlantic Ocean. Contrasting with global δ 66 Zn data sets, our results (δ 66 Zn between 0.24 ± 0.04 and 0.41 ± 0.04) do not support direct sampling of recycled marine carbonates in the Bermuda HIMU mantle. Instead, we show that δ 66 Zn fractionation toward higher values is associated with magmatic processes and incorporation of carbon sourced from deep fluids associated with the formation of carbonatites. These carbon-rich fluids are likely sourced from the metasomatic reactions between the subducted cold slab of the Iapetus oceanic lithosphere ca. 500 Ma and the thickened continental lithospheric mantle of Pangea. Melting of this metasomatized mantle was triggered by the arrival of the Farallon slab to the eastern North American margin in the late Cenozoic via shallow convection.

Mazza, Sarah E. [Smith College, Northampton, MA (U↗
DOE OSTI2022

Effects of Intrusive Confining Units on Groundwater Flow at Pahute Mesa, Nevada National Security Site

Geological, geophysical, hydrogeological, and hydrological evidence indicates that feeder dikes to lava flows play a role as confining units in the groundwater flow system underlying the Pahute Mesa area of underground nuclear testing at the Nevada National Security Site (NNSS). Of geological significance, the Pahute Mesa testing area is largely situated within the Silent Canyon caldera, which includes both tuffs and lavas that were fed by (and are connected to) intrusive rock units such as feeder dikes prior to and during volcanic eruptions (Orkild et al., 1968; Byers et al., 1976). Groundwater modeling results presented in this report support extension of “pool-dam” and “breach scenario” concepts presented in Jackson and Fenelon (2018) and Jackson et al. (2021) to low-permeability intrusive confining units (ICUs) such as feeder dikes that are expected to dissect and ring the caldera. As laterally extensive, steeply-dipping, and sub-parallel planar structures, ICUs can be expected to further compartmentalize groundwater flow in the volcanic aquifers systems at Pahute Mesa beyond the stratigraphically-defined tuff confining units (TCUs). ICUs provide an explanation for changes in water levels before and after the Cheshire test, detonated 2/14/76 (USDOE, 2015), and the unusually high horizontal gradient southwest of Cheshire. These ICU concepts are supported by groundwater modeling results presented in this report. The evidence and groundwater modeling results indicate that the downgradient extent of radionuclide contamination in some areas of Pahute Mesa could be overestimated by ignoring the confining effects of ICUs. On the other hand, contaminant extent could be underestimated by not considering transient flow effects and/or the re-equilibrated flow field caused by breaching of an ICU by an underground nuclear test.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Post-Closure Inspection Report for the Tonopah Test Range and Nevada Test and Training Range, Nevada for Calendar Year 2020 (Rev. 0)

This report provides the results of the annual post-closure use restriction (UR) inspections at corrective action sites (CASs) located on the Tonopah Test Range (TTR) and Nevada Test and Training Range (NTTR). This report covers post-closure UR inspections for calendar year 2020, and includes visual inspections and repair activities completed at the following corrective action units (CAUs): • CAU 407, Roller Coaster RadSafe Area (TTR) • CAU 424, Area 3 Landfill Complexes (TTR) • CAU 453, Area 9 UXO Landfill (TTR) • CAU 487, Thunderwell Site (TTR). Note: CAU 400, Bomblet Pit and Five Points Landfill (TTR) (CASs TA-19-001-05PT, Ordnance Disposal Pit; and TA-55-00-TAB2, Ordnance Disposal Pit) no longer requires inspection and reporting, as described in the Post-Closure Inspection Report for the Tonopah Test Range and Nevada Test and Training Range, Nevada for Calendar Year 2019. Site figures are included in Appendix A. The Post-Closure Inspection Plans and Record of Technical Change modifying the requirements for each UR are included in Appendix B. The Post-Closure Inspection Checklists are included in Appendix C. Field Notes are included in Appendix D. Appendix E, Inspection Photographs, does not include photographs because none were taken during the 2020 post-closure inspections. Appendix F, Post-Closure Vegetation Monitoring Report, does not include any data because vegetation monitoring was not conducted or required at any CAU in 2020. Visual inspections were conducted according to the post-closure requirements outlined in specific CAU and CAS URs, and approved by NDEP. The annual post-closure inspections for CASs associated with CAUs 407, 424, 453, and 487 were conducted on May 20, 2020. No maintenance or repair issues were noted at CAUs 487, or at CASs associated with CAU 424 (except CASs 03-08-002-A304 and 03-08-002-A308). Maintenance items and subsequent repairs include the following: • CAU 407, CAS TA-23-001-TARC (Roller Coaster RadSafe Area): A UR sign was loose and was rehung during the inspection. • CAU 424, CAS 03-08-002-A304 (Landfill Cell A3-4): The northeast monument needed a UR sign. The UR sign was attached to the monument on August 11, 2020. • CAU 424, CAS 03-08-002-A308 (Landfill Cell A3-8): The southeast monument needed additional lava rock. The lava rock was added during the inspection. • CAU 453, CAS 09-55-001-0952 (Area 9 Landfill): Two animal burrows were noted in the western portion of the northern and southern trenches (one burrow per trench). The animal burrows were repaired on August 11, 2020.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗
NASA NTRS1972

Lunar rilles and Hawaiian volcanic features - Possible analogues.

Study of the origin and development of lava conduits in terrestrial volcanic regions and the possible association of two important classes of lunar morphological features - i.e., channels and the sinuous rilles on the lunar maria - with stages in the evolution of lava conduits. The study was motivated by observations of some Hawaiian active volcanic regions. The origin of apparent flow channels on steep slopes on the moon is also considered.

Cruikshank, D. P.↗
NASA NTRS1973

Application of remote sensing in the study of vegetation and soils in Idaho

The author has identified the following significant results. Comparison of ERTS-1 imagery and USGS 1:250,000 scale maps of study areas with known ground points revealed significant map errors. These errors were sufficient to render impractical the projection of ERTS-1 imagery directly onto maps of the area. Marked differences were found in the delineation of ground features by different MSS bands. Generally, Band 4 was least useful, while Band 5 proved valuable for indicating patterns of native vegetation, cultivated areas - both dry and irrigated, lava fields, drainage basins, and deep bodies of water. Band 6 was better for landforms and drainages and for shallow bodies of water than Band 5 but inferior for indicating patterns in native vegetation and most types of cultivated land. Band 7 was best of all for indicating lava flows, water bodies, and landform features. Use of a additive color viewer-projector aided greatly in separation of images. A combination of Bands 5 and 7 with appropriate color filters proved best for separating most types of native vegetation and cultivated crops. Landform features and water bodies also showed well with this combination. The addition of Band 4 imagery to these further enhanced the identification of semi-dormant vegetation.

Tisdale, E. W.↗
NASA NTRS1972

Significant results from Apollo 14 lunar orbital photography.

Apollo 14 obtained 950 photographs from lunar orbit using the Hasselblad and Hycon cameras. The photographs reveal a number of new geologic features as well as previously unrecognized details of the morphology, structure, and stratigraphy of lunar surface units. The primary result is the verification of the extensive role of volcanism in the formation and modification of the lunar highlands, especially on the far side. Terra volcanism appears to be manifest in the formation of (1) constructional units of hilly and furrowed materials of regional extent as in the Kant Plateau in the central near-side highlands and northwest of the crater Pasteur near the eastern limb of the moon; (2) somewhat viscous lava flows and pools associated with fracture systems and/or what appear to be volcanic craters; (3) craters, crater chains, and irregular depressions, particularly on the lunar far side. The first photographs of a flow channel, a leveed sinuous rille that apparently originated by lava flowage on the surface, were obtained by Apollo 14. Another first is a high-resolution photograph of the interior of what appears to be the youngest lunar crater yet photographed in the 20- 40-km size range.

El-Baz, F.↗
NASA NTRS1973

A survey of the selenochemistry of major, minor and trace elements.

Average data for igneous and/or metaigneous rocks and soils from seven lunar sites are presented. There are compositional similarities between Apollo 11 and Luna 16 eastern maria, Ap 12 and 15 western maria, and between Ap 16 and L 20 highlands. Subtle differences do exist between the paired mare sites and the two highland sites and striking differences between the eastern and western maria. Chondritic normalized REE (rare earth element) patterns for igneous rocks and soils from all sites range from 7-350 generally with negative Eu anomalies. Anorthositic gabbroes to anorthosites, presumably highland material, exhibit a positive Eu anomaly. The REE patterns or Sr isotopic ratios suggest two lava flows each for the L 16 and Ap 14 sites, at least four lava flows for the Ap 11 and 12 site and about six for the Ap 15 site.

Schmitt, R. A.↗
NASA NTRS1973

Volcanism subprogram: Volcanological interpretation of the northern part of the Occidental Cordillera of Bolivia, utilizing ERTS imagery

The author has identified the following significant results. In the present study, 6 ERTS-1 images have been interpreted on a 1:1 million scale (black and white) with the respective field reconnaissance. The area studied is located in the region bordering with Chile and includes the western part of the Bolivian Altiplano, the volcano Cordillera (western cordillera) and the northern part of Chile to the Pacific Coast. The greater part of this region is formed by Pliocene/Pleistocene volcani rock, which is discordant with the Tertiary sediments with intercalations of calcareous tuff. The ERTS-1 imagery permits the tracing of regional boundaries of the great volcanic formations and the alinements of the volcanic bodies along the fault zones. They also permit a clear examination of the volcanic apparatus, including their secondary forms, such as lava flows, parasitic cones, and lava domes. Because of the great scale, it is not possible to identify either the small structures or those of low relief. On the basis of the interpretation of the images, it is possible to give an idea of the relative age of the volcanoes.

Brockmann, C. E.↗
NASA NTRS1974

Geologic applications of ERTS images on the Colorado Plateau, Arizona

Three areas in central and northern Arizona centered on the (1) Verde Valley, (2) Coconino Plateau, and (3) Shivwits Plateau were studied using ERTS photography. Useful applications results include: (1) upgrading of the existing state geologic map of the Verde Valley region; (2) detection of long NW trending lineaments in the basalt cap SE of Flagstaff which may be favorable locations for drilling for new water supplies; (3) tracing of the Bright Angel and Butte faults to twice their previously known length and correlating the extensions with modern seismic events, showing these faults to be present-day earthquake hazards; (4) discovering and successfully drilling perched sandstone aquifers in the Kaibab Limestone on the Coconino Plateau; and (5) determining the relationship between the Shivwits lavas and the formation of the lower Grand Canyon and showing that the lavas should be an excellent aquifer, as yet untapped.

Goetz, A. F. H.↗