SEARCH · Geology Reports
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Reversal in estuarine sand supply driven by Holocene sea level rise: A model for sand transport in large structural estuaries, San Francisco Bay, California, USA
Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.
The middle Cretaceous greenhouse period experienced profound environmental change including episodes of enhanced global burial of organic carbon marked by carbon isotopic excursions (CIEs). However, the role and response of polar regions like the newly formed, partially enclosed Arctic Ocean Basin during middle Cretaceous carbon burial remains enigmatic. We present the first Arctic deepwater CIE record that characterizes conditions offshore of the Alaska margin north of 75°N paleolatitude. Organic carbon isotopes (δ 13 C org ) and 103–82 Ma ash zircon U-Pb dates from the distal Hue Shale record multiple Albian–Campanian CIEs during slow ∼3–15 m/Myr sediment accumulation rates. Average total organic carbon (TOC) increased substantially during large 2–3 ‰ CIEs of the ∼101 Ma Albian-Cenomanian boundary event (from 7 to 18 % TOC) and ∼94 Ma Cenomanian-Turonian boundary event (5 to 10 % TOC). Turonian TOC remained elevated (8–13 %) during high global sea levels and temperatures of the Cretaceous Thermal Maximum, followed by an increase from 7 to 11 % TOC during the ∼90 Ma late Turonian event 1.5 ‰ CIE. Average TOC subsequently decreased in the Coniacian–Campanian, but relative maxima occurred during subtle 0.5–1 ‰ CIEs interpreted as the ∼87 Ma late Coniacian event (increase from 4 to 7 % TOC), ∼85 Ma Horseshoe Bay event (3.5 to 4.5 % TOC), and ∼84 Ma Santonian-Campanian boundary event (3.5 to 5 % TOC). Increases in hydrogen index and productivity proxies (P, Ba, Nd) that accompanied each CIE episode with enhanced TOC suggest a strong link between marine productivity and organic carbon burial at short-term CIE timescales. However, long-term (>5–8 Myr) changes in trace metal redox (Mo, Fe, V) and salinity (B/Ga) proxies suggest shifts in prevailing environmental conditions at timescales longer than the CIEs. Late Albian–middle Turonian marine salinity occurred during euxinic (103–98 Ma) and suboxic (98–90 Ma) conditions with deposition interpreted to have occurred within and beneath an oxygen minimum zone, respectively. In contrast, late Turonian–early Campanian (90–82 Ma) freshening and restricted euxinic basin conditions may signal the start of widespread restriction known to characterize the Paleogene Arctic. Overall, these results highlight that middle Cretaceous Arctic deepwater remained a productive marine carbon sink coupled to the global carbon cycle despite evolving Arctic greenhouse conditions.
Two contrasting age models for initial mountain building in the northeastern (NE) Tibetan Plateau (Paleocene-early Eocene versus late Oligocene-early Miocene) have led to the debate on how the deformed continental lithosphere absorbs plate convergence in general. The initial compressional deformation in the West Qinling (WQL) of the NE Tibetan Plateau figures prominently in this ongoing debate. Here, apatite (U-Th)/He (AHe) thermochronology combined with geomorphological analysis are used to refine the onset of compressional deformation in the WQL. New AHe ages from two vertical transects and an updated reconstruction of an obliquely-tilted erosion surface document the accelerated exhumation in the northern WQL at 23-22 Ma, interpreted as the onset of north-vergent thrusting. The AHe results, together with sedimentary records in the intermontane and foreland basins, suggest that the entire WQL began experiencing compressional deformation in the late Oligocene-early Miocene. When integrated with previous studies, our findings show that the northern plateau boundary has not remained stationary since the collision, but has instead experienced ∼750 km of outward expansion during the late Oligocene to middle Miocene. This phase of rapid plateau growth is coeval with the ∼30–50 % reduction of the India-Eurasia convergence rate, which suggests that the increased gravitational potential energy of orogenic belts played a key role in plate motion changes.
We assess the growth of anomalously high relief on Denali, located in the Alaska Range, Alaska, and the tallest mountain in North America (6190 masl). Denali is 3000 m taller than most surrounding peaks. It lies inside a 19° restraining bend in the active Denali fault system that is moving at about 7 mm/yr, and thus there is a tectonic and structural driver for ongoing rock uplift. High relief around Denali is also due, in part, to its granitic rock type and low fracture density relative to adjacent metasediments. Here we show that unique climatic conditions at high elevations also contribute to the growth of relief. We examine 10Be concentrations in 1) three new gravel samples between 3500 and 5200 m elevation from sites unaffected by recent glacial erosion, 2) previously published samples from a sidewall of the Kahiltna Glacier from 2400 to 2800 masl, 3) previously published samples from medial moraines along the length of the Kahiltna Glacier from 1400 to ~500 masl, and 4) previously published alluvial samples along the Kahiltna River at an elevation of ~200 masl. These samples constitute a transect extending more than 5000 vertical meters, and the data establish that erosion rates decrease with elevation and contribute to the growth of relief. Erosion rates for the three new high-elevation samples are calculated to 4.6±0.6 mm/ka at 5.2 km, 28.6±3.7 mm/ka at 4.0 km, and 38±5 mm/ka at 3.5 km. Erosion rates at intermediate elevations, on the sidewall of the Kahiltna Glacier, range between 160-327 mm/ka. Along the medial moraines inferred erosion rates range between 140-537 mm/ka, and basin-wide erosion rates calculated from sediments in the river below the glacier range between 450-896 mm/ka. These differences in erosion rates can create relief of 3 km within 1-10 Ma, well within the estimated period of increase in rock uplift and exhumation on Denali over the last ~6 Ma. Meteorological data from 2130 to 5550 masl at 5 sites show it rarely gets above freezing above 4 km elevation, indicating that frost weathering currently plays a diminished role in erosion at high elevations. The immediate implication of this temperature and erosional correlation is an increase in relief. This is the first study to directly measure a significant decrease in erosion rates at high elevations in the relative absence of frost weathering. The results highlight the combined influence of rock type, glacial erosion, and permanent sub-zero temperatures on erosion rates. In combination with active faulting, the data explain the resultant increase in relief along the southern side of the Alaska Range over the past 100 ka.
The Coso geothermal field is a major geothermal power production site in the western United States. It has been observed that low-magnitude seismic events ( M < 3.71) are unevenly distributed in three distinct zones, namely, nearfield (<3 km), midfield (3–6 km), and farfield (> 6 km) from the Coso geothermal plant. These zones exhibit distinct changes in earthquake location before and during geothermal production episodes that began in 1986. After 1986, the midfield region of the main flank experiences a significantly lower seismicity rate than the surrounding areas before production episodes. During 2014–2019, the farfield earthquakes cluster in the eastern and western parts of the greater Coso area, which is discernably different from how those pre-production earthquake events were distributed along the conjugate NW-SE and SW-NW trending structures across the main flank. Here, we analyze the stage of stress with finite-element-based poroelastic simulations to illustrate how the spatiotemporal evolution of the seismicity is associated with the pattern of stress perturbations caused by fluid migration amid the operations of geothermal power plants. Generally, ∼70% of co-production seismicity is found in zones of increased Coulomb stress between 2014 and 2019 at >99% confidence. Meanwhile, the midfield zone of seismic paucity overlaps with the zone of decreasing pore-fluid pressure. Overall, the results provide a physical explanation of how decadal geothermal operations at Coso have perturbed stress-field changes and contributed to the evolving characteristic seismic pattern, shedding insights into assessing the seismic hazard in other geothermal settings.
Degassing of water from magmatic systems is key to transporting metals from magmas to form ore deposits, but elements like chlorine, through the formation of anion complexes, can be important in solubilizing and mobilizing these metals into water-rich fluids. Reconstructing the Cl systematics of evolving magmas is thus an important step towards understanding the origins of ore deposits, but the magmatic record is not well preserved because Cl can be lost during degassing. Here, we reconstruct the pre-degassing history of Cl in subduction zone (arc) magmas through amphiboles, which incorporate Cl directly into their crystal structures, preserving pre-eruptive magmatic signatures. Amphibole-reconstructed Cl contents indicate that magmatic differentiation can lead to a 4-fold increase in concentration due to Cl's incompatible behavior. The amphibole-reconstructed Cl contents of arc magmas are also significantly higher than values reported from melt inclusions, suggesting that many melt inclusions may have been trapped after magmas had already lost some Cl. We show that such Cl loss is likely associated with preferential partitioning of Cl into hydrous fluids degassed from the magma during crustal storage or ascent. The extent of Cl depletion can thus be used to estimate how much water was lost during early degassing. If Cl is important to certain ore deposits, magmatic water content may play an indirect role. Magmas too rich in water will lose water and hence Cl at greater depths, rendering such magmas less able to transport metals to the upper crust. By contrast, drier magmas may not produce enough Cl-rich fluids to mobilize metals. Thus, magmas with intermediate water contents may produce enough Cl-rich fluids at the right depths for certain types of ore deposits.
On Earth, silica-rich phases from opal to quartz are important indicators and tracers of geological processes. Hydrated silica, such as opal, is a particularly good matrix for the preservation of molecular and macroscopic biosignatures. Cherts, a type of silica-dominated rocks, provide a unique archive of ancient terrestrial life while quartz is the emblematic mineral of the Earth's continental crust. On Mars, hydrated silica has been detected in several locations based on remote sensing and rover-based studies. In the present article we report on the detection of cobbles made of hydrated silica (opal or chalcedony), as well as well-crystallized quartz. These detections were made with the SuperCam instrument onboard Perseverance (Mars 2020 mission), using a combination of LIBS, infrared and Raman spectroscopy. Quartz-dominated stones are detected unambiguously for the first time on the Martian surface, and based on grain size and crystallinity are proposed to be of hydrothermal origin. Although these rocks were all found as float, we propose that these detections are part of a common hydrothermal system, and represent different depths / temperatures of precipitation. This attests that hydrothermal processes were active in and around Jezero crater, possibly triggered by the Jezero crater-forming impact. These silica-rich rocks, in particular opaline silica, are very promising targets for sampling and return to Earth given their high biosignature preservation potential.
Many open-vent arc volcanoes display two modes in their continuous gas emissions, one with a characteristic CO 2 / S T ratio typical of periods of quiescent degassing and another punctuated by high CO 2 / S T gas emitted in the weeks before eruption, a recently recognized eruption precursor. In this study we explore the origin of the two modes of degassing revealed by time-series gas data at Turrialba volcano (Costa Rica) in the context of new melt inclusion (MI) data. To reconstruct the c[CO 2 ] of undegassed magma, we developed a rapid-quench piston-cylinder assembly to rehomogenize the vapor bubble commonly contained in MIs. We focus on olivine-hosted MIs from a mafic scoria sample erupted from Turrialba in 1864–1866. The reconstructed CO 2 contents in MIs decrease from ∼4000 to <1000 ppmw as S contents decrease from 3500 to <1000 ppmw. The highest reconstructed S and CO 2 in the MIs resulted in an initial magmatic CO 2 / S T ratio (molar) of 0.83. Informed by the MI data, we modeled the decompression degassing of Turrialba magma and vapor composition using the Sulfur_X and EVo models. Instead of being controlled by initial magmatic CO 2 /S T ratio as suggested by previous studies, we find that the quiescent gas emitted from Turrialba during 2014–2018 (CO 2 / S T = 2.3 ± 0.8, molar) appears to reflectequilibrium with magmas stored at 4–8 km (Sulfur_X) or 2 km (EVo) depth, when H 2 O is degassing extensively from the magma. A magma storage region at 4–8 km is also supported by seismic tomography. The second gas mode is noted by spikes in CO 2 / S T ∼ 7.9 ± 2 in the weeks prior to eruption. This gas reflects equilibrium with a magma at 12–18 km (Sulfur_X) or 4–8 km (EVo), where the ascending magma is saturated with a CO 2 -rich vapor. Thus, there are two important trans crustal depths beneath the volcano: one where the rate of H 2 O loss from the magma and thus magma viscosity increases, and one at greater depths where high CO 2 /S T vapor forms and may facilitate dike propagation. We interpret the shallower, H 2 O-loss region as the main site of magma stalling and storage, where quiescent gas is generated continuously. We interpret the greater depth (12–18 km) as the source of the precursory gas that precedes eruption, and where the mafic melt lastly equilibrated with a mush zone before ascending and triggering eruption weeks later. This hypothesis is ripe for testing at other volcanoes that exhibit two modes in gas geochemistry.
The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.
Sulfate is often touted as containing atmospheric oxygen whose isotopic signature can constrain redox, environmental conditions, and biological activity. Yet, the amount and isotopic fractionation associated with air-O 2 incorporation during sulfate formation is still debated, making its verification difficult. In this study, we identify a distinct, microbially dominated environment with the potential to preserve maximum signals of air-O 2 in sulfate. We report triple-oxygen isotope data for sulfate produced from pyrite oxidation in microbial and abiotic experiments, and from natural dissolved sulfate from the Rio Tinto, Spain, an acid mine drainage site. The oxygen isotope systematics of sulfate in these environments define a unique kinetic isotope effect associated with initial stage pyrite oxidation by Acidithiobacillus ferrooxidans that preserves >80 % oxygen from air-O 2 in sulfate. Unlike experiments, which evolve toward water-oxygen dominated sulfate on short time scales, Rio Tinto, Spain hosts a microbe rich environment with distinct geochemistry that maintains high O 2 -oxygen in sulfate. Therefore, in addition to containing isotopic records from water and air, sulfates can also contain a biosignature that is promising for understanding conditions on Mars and early Earth.
Based on observation and data from meteorites and in situ scientific missions, experiments as well as models, the Martian mantle is assumed to share some compositional and mineralogical affinity with the terrestrial mantle. However, there might be subtle differences like the Martian mantle being more ferroan. Yet, we do not have any direct analysis of a Martian mantle rock to confirm this assumption. NASA’s Perseverance rover found olivine-rich boulder-sized float rocks on the upper Jezero fan (Mars). These boulders have an ultramafic composition and their mineralogy is dominantly composed of Fo 73±3 olivine with high-Mg orthopyroxene, Cr-rich Ti-Fe oxides and minor plagioclase and high-Ca pyroxene. Microtextural and petrological analysis reveals that these minerals crystallized at equilibrium. In addition, these boulders are different from all the bedrocks analyzed by Perseverance along its traverse which are crustal igneous rocks and sediments. Comparing our data to Martian meteorites and available Mars bulk silicate models (BSM), we discuss that these boulders could represent primitive melts and/or lower crustal material, and we specifically hypothesize that they could be mantle peridotites. We propose that these putative mantle rocks could have been excavated by the succession of impacts from the shallow mantle or lower crust in the Isidis region where Jezero crater is located. These olivine-rich boulders could thereby constitute the first direct analysis of a Martian mantle rock.
The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10 5 –10 6 -year, rather than remaining in a long-term stationary configuration on ∼10 7 -year timescales.
Characterizing tectonic subsidence rates within depositional sequences provides direct insight into the driving mechanism(s) of accommodation in a basin. However, the temporal resolution of this record is often stymied by a lack of high-precision and high-resolution ages, which enable a more complete description of basin subsidence drivers. We explore the effect of high-precision and high-resolution ages in modeling accommodation for the Miocene La Jencia Basin of the central Rio Grande rift (RGR) and interpret driving mechanism behavior from these models (e.g., lithospheric thinning). We present a new geochronologic dataset of both laser ablation-inductively coupled plasma mass spectrometry (LA-ICPMS) and chemical abrasion-isotope dilution-thermal ionization mass spectrometry (CA-ID-TIMS) data and use these ages in Bayesian accommodation modeling. Models constrained by high-precision and high-accuracy TIMS ages yield peak tectonic subsidence rates exceeding 220 m/Myr, and an average Miocene subsidence of ∼120 m/Myr. While timing and magnitude vary, all models suggest two pulses of rapid Miocene tectonic subsidence, which we interpret to reflect basin-bounding fault movement. Prior to peak subsidence, there was an initial period of fault linkage and organization that occurred over <1–3 Myr that produced the basin-bounding La Jencia-Cerro Colorado fault zone. A comparison of published tectonic subsidence rates to those modeled here shows that while tectonic subsidence during continental rifting is highly variable, the La Jencia Basin rates appear relatively high. However, the significant difference between peak and average La Jencia Basin rates modeled here highlights the potential for underestimation of many records of tectonic subsidence due to a lack of high-precision and high-resolution age constraints. Furthermore, age data and modeling results presented here document fault movement and consequent rates of tectonic subsidence that lower-resolution data would not, providing a high-fidelity case study of continental rift basin development.
Removal of mafic-ultramafic lower crust (e.g., via delamination) is fundamental to making andesitic continental crust, yet direct evidence of this process remains elusive. A unique suite of garnet clinopyroxenite and hornblendite (arclogite) xenoliths from the northern volcanic zone of the Andes, erupted to the surface in the mid-Pleistocene, have bulk-rock geochemistry and osmium (Os) isotopic compositions that demonstrate their crustal affinity, yet have equilibration pressures and temperatures below the arc Mohorovičić discontinuity (~53 km) and some as deep as within the sub-arc mantle wedge (~105 km). Garnet websterites from the same xenolith suite, sourced from depths approaching the Wadati-Benioff zone (~140 km), have elemental and isotopic compositions indicative of a mantle origin, likely formed as products of peridotite-liquid reactions above the subducting slab. Variability in bulk-rock 143 Nd/ 144 Nd values and garnet oxygen isotope ratios ( 𝛿 18 O) for these samples is attributed to minor assimilation of subducted components and/or older crustal material, but assimilation alone cannot account for the highly radiogenic Os isotopic compositions. These results provide direct petrologic evidence for modern arclogite formation and foundering in the Andean orogen, the archetypal active continental subduction system. Elevated 187 Os/ 188 Os values of Mercaderes arclogites imply that recycling of gravitationally unstable arc cumulates during continental crust formation introduces highly radiogenic Os into the convective mantle, which has implications for mass fluxes across the crust-mantle boundary as well as the sources and evolution of mantle heterogeneity as seen in orogenic peridotite massifs and in the sources of oceanic basalts
Rheological data are crucial for the understanding of storage, transport, eruption, and emplacement of magma and lava. While viscosity estimates for lavas from Mauna Loa, the largest active volcano on Earth, exist, no direct rheological measurements have been done to date. This is especially surprising given the role Mauna Loa lavas and resulting landforms play in the interpretation of volcanic processes across our solar system and beyond. We present the first complete rheological characterization of Mauna Loa lava, based on the 2022 eruption. We constrain the melt viscosity and map the lava's rheology during crystallization at isothermal subliquidus conditions and at cooling at rates between 0.25 and 3.00 °C/minute. These experiments were performed at oxygen fugacities relevant to storage and eruption conditions at Mauna Loa (log fO 2 = -8.7). We integrate textural and thermal characteristics obtained from the experiments with those from water quenched natural samples (i.e. sampled from actively flowing lava and quenched within seconds to preserve the textural state of the lava while flowing) to reconstruct the evolution of lava rheology during emplacement. This integration suggests that the Mauna Loa lavas were emplaced at viscosities between 10 1.5 to 10 4.5 Pa s over the temperature range of 1150 – 1090 °C before entering a region of rapid viscosity increase and effective solidification caused by groundmass crystallization. We provide a detailed comparison of our data to rheological reconstructions of the 1984 eruption. The resulting rheological map can help guide physical property-based magma ascent and lava emplacement models, as well as the interpretation of flow morphologies and the conditions of formation of basaltic landforms on other planets.
n southwestern Sardinia, Italy, the Rio Naracauli drains a catchment that includes several abandoned mines. The drainage from the mines and associated waste rocks has led to extreme concentrations of dissolved Zn, but because of the near-neutral pH, concentrations of other metals remain low. In the reach from approximately 2300 to 3000 m downstream from the headwaters area, an amorphous Zn-silicate precipitates from the water. In this reach, concentrations of both Zn and silica remain nearly constant, but the loads (measured in mass/time) of both increase, suggesting that new Zn and silica are supplied to the stream, likely from emerging groundwater. Zinc isotope signatures of the solid are heavier than the dissolved Zn by about 0.5 permil in 66/64 Zn, suggesting that an extracellular biologically mediated adsorption process may be involved in the formation of the Zn-silicate.
Oil and natural gas have been the main sources of primary energy in the USA, providing 63% of the total energy consumption in 2011. Petroleum production, drilling operations, and improperly sealed abandoned wells have caused significant local groundwater contamination in many states, including at the USGS OSPER sites in Oklahoma. The potential for groundwater contamination is higher when producing natural gas and oil from unconventional sources of energy, including shale and tight sandstones. These reservoirs require horizontally-completed wells and massive hydraulic fracturing that injects large volumes (up to 50,000 m3/well) of high-pressured water with added proppant, and toxic organic and inorganic chemicals. Recent results show that flow back and produced waters from Haynesville (Texas) and Marcellus (Pennsylvania) Shale have high salinities (≥200,000 mg/L TDS) and high NORMs (up to 10,000 picocuries/L) concentrations. A major research effort is needed worldwide to minimize all potential environmental impacts, especially groundwater contamination and induced seismicity, when producing these extremely important new sources of energy.