Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Earth in Space”

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.

At least 649 records · Page 36Linked to original sources

Summary appraisals of the nation's ground-water resources – Pacific Northwest region

The Pacific Northwest Region's ground-water reservoirs are capable of providing large additional fresh-water supplies; these reservoirs become more important as undeveloped surface-storage sites and unapportioned surface-water supplies dwindle. Withdrawals of fresh water from all surface and underground sources are increasing; they may rise from the rate of 30 billion gallons per day in 1970 to about 60 billion gallons per day in 2020. By 1975 the withdrawal of ground water had increased 70 percent over the 1970 rate and accounted for 22 percent of total fresh-water withdrawal. Substantial increases in ground-water withdrawal must continue if projected water demands are to be met in the future. Large variations exist in the availability of and needs for ground water in the region, largely because of the great variety of landforms, climate, and earth materials. More than one-half the region is underlain by rock materials capable of yielding ground water to wells at moderate to large rates; six extensive areas are identified as major ground-water reservoirs. The most significant current problems pertaining to one or more of these ground-water reservoirs are: (1) Progressively declining water levels; (2) water-quality problems; (3) waterlogging; (4) inadequate information; and (5) competition for available supplies. In the future these same problems are expected to persist and generally worsen (especially water-quality deterioration) in most of the major reservoir areas. Management opportunities in the region include: (1) Development of new supplies and additional uses of ground water; (2) protection and enhancement of water quality; (3) reduction of waterlogging; (4) energy development from some ground-water reservoirs; (5) improving access to the ground water; (6) increased use of underground space for storage and disposal; and (7) greater use of advanced management and conservation techniques. Conjunctive use of surface and ground water to provide greater available supplies probably is the most promising water-management opportunity. However, if the full potential of the ground-water resources is to be realized, important constraints, including present water-right structures and serious deficiencies in information, must be overcome.

Idaho, Montana, Oregon, Washington↗

Methane hydrate formation in turbidite sediments of northern Cascadia, IODP Expedition 311

Expedition 311 of the Integrated Ocean Drilling Program (IODP) to northern Cascadia recovered gas-hydrate bearing sediments along a SW-NE transect from the first ridge of the accretionary margin to the eastward limit of gas-hydrate stability. In this study we contrast the gas gas-hydrate distribution from two sites drilled ~ 8??km apart in different tectonic settings. At Site U1325, drilled on a depositional basin with nearly horizontal sedimentary sequences, the gas-hydrate distribution shows a trend of increasing saturation toward the base of gas-hydrate stability, consistent with several model simulations in the literature. Site U1326 was drilled on an uplifted ridge characterized by faulting, which has likely experienced some mass wasting events. Here the gas hydrate does not show a clear depth-distribution trend, the highest gas-hydrate saturation occurs well within the gas-hydrate stability zone at the shallow depth of ~ 49??mbsf. Sediments at both sites are characterized by abundant coarse-grained (sand) layers up to 23??cm in thickness, and are interspaced within fine-grained (clay and silty clay) detrital sediments. The gas-hydrate distribution is punctuated by localized depth intervals of high gas-hydrate saturation, which preferentially occur in the coarse-grained horizons and occupy up to 60% of the pore space at Site U1325 and > 80% at Site U1326. Detailed analyses of contiguous samples of different lithologies show that when enough methane is present, about 90% of the variance in gas-hydrate saturation can be explained by the sand (> 63????m) content of the sediments. The variability in gas-hydrate occupancy of sandy horizons at Site U1326 reflects an insufficient methane supply to the sediment section between 190 and 245??mbsf. ?? 2008 Elsevier B.V.

Earth and Planetary Science Letters↗

Grooved terrain on Ganymede

The icy crust of Ganymede comprises bright and dark areas. Investigation of Voyager 1 and 2 images has shown that bright terrain is grooved and separates dark polygons of cratered terrain. The grooved terrain contains alternating ridges and grooves in straight and curvilinear sets, which are locally interrupted by smooth patches and swaths. Cratered terrain, where 'it occurs in small wedges and slivers, has a pervasive grain of narrowly spaced furrows, and thus is transitional to grooved terrain. An analysis of the morphology of terrain types, and of superposition and cross-cutting relations, suggests that grooved terrain grew at the expense of cratered terrain, that tracts of cratered terrain were converted into grooved terrain in situ , and that vertical tectonism and shear movements dominated in the restructuring of Ganymede's surface. It is postulated that during a period in the planet's history when the lithosphere was thin, upwelling convection currents caused incipient rifting accompanied by intensive normal faulting; where rifting went to completion, crustal segments separated, locally spread apart, and sheared past one another. In places subduction and compression may have occurred, but the evidence is inconclusive. Thus, the grooved terrain on Ganymede may record an early phase of ice-plate tectonics that caused rifting and drifting of the icy lithosphere, but, unlike silicate plate tectonics on Earth, may have resulted in only minor vertical turnover.

Icarus↗

Passive seismic monitoring of natural and induced earthquakes: Case studies, future directions and socio-economic relevance

An important discovery in crustal mechanics has been that the Earth’s crust is commonly stressed close to failure, even in tectonically quiet areas. As a result, small natural or man-made perturbations to the local stress field may trigger earthquakes. To understand these processes, Passive Seismic Monitoring (PSM) with seismometer arrays is a widely used technique that has been successfully applied to study seismicity at different magnitude levels ranging from acoustic emissions generated in the laboratory under controlled conditions, to seismicity induced by hydraulic stimulations in geological reservoirs, and up to great earthquakes occurring along plate boundaries. In all these environments the appropriate deployment of seismic sensors, i.e., directly on the rock sample, at the earth’s surface or in boreholes close to the seismic sources allows for the detection and location of brittle failure processes at sufficiently low magnitude-detection threshold and with adequate spatial resolution for further analysis. One principal aim is to develop an improved understanding of the physical processes occurring at the seismic source and their relationship to the host geologic environment. In this paper we review selected case studies and future directions of PSM efforts across a wide range of scales and environments. These include induced failure within small rock samples, hydrocarbon reservoirs, and natural seismicity at convergent and transform plate boundaries. Each example represents a milestone with regard to bridging the gap between laboratory-scale experiments under controlled boundary conditions and large-scale field studies. The common motivation for all studies is to refine the understanding of how earthquakes nucleate, how they proceed and how they interact in space and time. This is of special relevance at the larger end of the magnitude scale, i.e., for large devastating earthquakes due to their severe socio-economic impact.

Book chapter↗

Orogenic gold and geologic time: A global synthesis

Orogenic gold deposits have formed over more than 3 billion years of Earth's history, episodically during the Middle Archean to younger Precambrian, and continuously throughout the Phanerozoic. This class of gold deposit is characteristically associated with deformed and metamorphosed mid-crustal blocks, particularly in spatial association with major crustal structures. A consistent spatial and temporal association with granitoids of a variety of compositions indicates that melts and fluids were both inherent products of thermal events during orogenesis. Including placer accumulations, which are commonly intimately associated with this mineral deposit type, recognized production and resources from economic Phanerozoic orogenic-gold deposits are estimated at just over one billion ounces gold. Exclusive of the still-controversial Witwatersrand ores, known Precambrian gold concentrations are about half this amount. The recent increased applicability of global paleo-reconstructions, coupled with improved geochronology from most of the world's major gold camps, allows for an improved understanding of the distribution pattern of orogenic gold in space and time.

Ore Geology Reviews↗

Geomagnetism Program research plan, 2020–2024

The Geomagnetism Program of the U.S. Geological Survey (USGS) monitors geomagnetic field variation through operation of a network of observatories across the United States and its territories, and it pursues scientific research needed to estimate and assess geomagnetic and geoelectric hazards. Over the next five years (2020–2024 inclusive) and in support of national and agency priorities, Geomagnetism Program research scientists plan to pursue an integrated set of research projects broadly encompassing empirical estimation and mapping of geomagnetic disturbance, modeling of solid-Earth conductivity structure and surface impedance, and mapping of magnetic-storm-induced geoelectric fields. Analyses are empirically based, relying on measured time series as well as statistical and numerical modeling of geomagnetic-monitoring data from ground-based observatories and surface-impedance tensors acquired during magnetotelluric surveys. The plan describes augmentation and development of the Geomagnetism Program's existing research portfolio, assuming present funding levels and staffing numbers. Because the projects are interdependent, they cannot be straightforwardly prioritized. They will all be pursued as resources and time permit; additional funding and staffing would enable the projects to be broadened and more rapidly completed. Where appropriate and subject to budgetary constraints and staffing numbers, research on specific projects might be accelerated or even judiciously expanded—some opportunities for expansion are discussed in this plan. Results will provide realistic illumination of the nature of the ground-level expression of space-weather disturbance, a subject of particular importance for projects focused on evaluating the vulnerability of electric-power-grid systems. This plan does not cover Geomagnetism Program operations, which are primarily concerned with the operation of magnetic observatories and, now, magnetotelluric surveys, although the context of such observatories and surveys is discussed. The research element of the program provides guidance for the expansion of program operations and research projects. In addition to the research projects summarized here, program scientists continue to provide leadership to the national and international geomagnetic, magnetotelluric, and space-weather communities.

Circular↗

A long-term geothermal observatory across subseafloor gas hydrates, IODP Hole U1364A, Cascadia accretionary prism

We report 4 years of temperature profiles collected from May 2014 to May 2018 in Integrated Ocean Drilling Program Hole U1364A in the frontal accretionary prism of the Cascadia subduction zone. The temperature data extend to depths of nearly 300 m below seafloor (mbsf), spanning the gas hydrate stability zone at the location and a clear bottom-simulating reflector (BSR) at ∼230 mbsf. When the hole was drilled in 2010, a pressure-monitoring Advanced CORK (ACORK) observatory was installed, sealed at the bottom by a bridge plug and cement below 302 mbsf. In May 2014, a temperature profile was collected by lowering a probe down the hole from the ROV ROPOS. From July 2016 through May 2018, temperature data were collected during a nearly two-year deployment of a 24-thermistor cable installed to 268 m below seafloor (mbsf). The cable and a seismic-tilt instrument package also deployed in 2016 were connected to the Ocean Networks Canada (ONC) NEPTUNE cabled observatory in June of 2017, after which the thermistor temperatures were logged by Ocean Networks Canada at one-minute intervals until failure of the main ethernet switch in the integrated seafloor control unit in May 2018. The thermistor array had been designed with concentrated vertical spacing around the bottom-simulating reflector and two pressure-monitoring screens at 203 and 244 mbsf, with wider thermistor spacing elsewhere to document the geothermal state up to seafloor. The 4 years of data show a generally linear temperature gradient of 0.055°C/m consistent with a heat flux of 61–64 mW/m 2 . The data show no indications of thermal transients. A slight departure from a linear gradient provides an approximate limit of ∼10 −10 m/s for any possible slow upward advection of pore fluids. In-situ temperatures are ∼15.8°C at the BSR position, consistent with methane hydrate stability at that depth and pressure.

Cascadia subduction zone↗

Thermal conductivity of hydrate-bearing sediments

A thorough understanding of the thermal conductivity of hydrate-bearing sediments is necessary for evaluating phase transformation processes that would accompany energy production from gas hydrate deposits and for estimating regional heat flow based on the observed depth to the base of the gas hydrate stability zone. The coexistence of multiple phases (gas hydrate, liquid and gas pore fill, and solid sediment grains) and their complex spatial arrangement hinder the a priori prediction of the thermal conductivity of hydrate-bearing sediments. Previous studies have been unable to capture the full parameter space covered by variations in grain size, specific surface, degree of saturation, nature of pore filling material, and effective stress for hydrate-bearing samples. Here we report on systematic measurements of the thermal conductivity of air dry, water- and tetrohydrofuran (THF)-saturated, and THF hydrate–saturated sand and clay samples at vertical effective stress of 0.05 to 1 MPa (corresponding to depths as great as 100 m below seafloor). Results reveal that the bulk thermal conductivity of the samples in every case reflects a complex interplay among particle size, effective stress, porosity, and fluid-versus-hydrate filled pore spaces. The thermal conductivity of THF hydrate–bearing soils increases upon hydrate formation although the thermal conductivities of THF solution and THF hydrate are almost the same. Several mechanisms can contribute to this effect including cryogenic suction during hydrate crystal growth and the ensuing porosity reduction in the surrounding sediment, increased mean effective stress due to hydrate formation under zero lateral strain conditions, and decreased interface thermal impedance as grain-liquid interfaces are transformed into grain-hydrate interfaces.

Journal of Geophysical Research B: Solid Earth↗

Modified GIC estimation using 3-D Earth conductivity

Geomagnetically induced currents (GICs) are quasi-direct current (DC) electric currents that flow in technological conductors during geomagnetic storms. Extreme GICs are hazardous to man-made infrastructure. GICs enter and exit the technological systems, such as the electric power grid, at grounding points, and their magnitudes depend on the currents that flow underground. They are, therefore, a function of the Earth's electrical conductivity, represented at ground level as Earth impedances, as well as the resistance parameters of the power network. Traditional GIC estimation practices are based on Earth impedances obtained from laterally homogeneous or piecewise layered-Earth models. We refer to these methods, collectively, as the 1-D approximation. However, GIC hazard mitigation can be improved with more accurate GIC modeling that takes the spatially heterogeneous Earth's conductivity into account. Here, we propose a modified approximation for GIC estimation that is very similar to the 1-D approximation but is instead derived from empirical 3-D Earth impedances. Our formulation sets up the computation of static, frequency-dependent power line telluric response functions, which, once computed, may be considered part of the power grid system model. These response functions may then be used for historical scenario analysis of GIC hazards and for simplified real-time, albeit approximate, GIC estimation in a power grid. This modest modification to the simpler local field formulation approach avoids real-time integration of geoelectric fields along power lines while taking the realistic 3-D Earth into account in a rigorous manner. Once implemented, the method provides a power grid operator with the benefits of convenience and computational speed for a first look real-time operational GIC hazard assessment. We estimate that the proposed modified 3-D GIC modeling approach produces GIC values that are well within 50% of those obtained with the full-scale power line integration of spatially variable geoelectric fields, for storms comparable in scale to the 2003 Halloween storm, all geological structures, and power lines located in the contiguous United States and other low- to middle-latitude regions.

Space Weather↗

A trait-based framework for predicting when and where microbial adaptation to climate change will affect ecosystem functioning

As the earth system changes in response to human activities, a critical objective is to predict how biogeochemical process rates (e.g. nitrification, decomposition) and ecosystem function (e.g. net ecosystem productivity) will change under future conditions. A particular challenge is that the microbial communities that drive many of these processes are capable of adapting to environmental change in ways that alter ecosystem functioning. Despite evidence that microbes can adapt to temperature, precipitation regimes, and redox fluctuations, microbial communities are typically not optimally adapted to their local environment. For example, temperature optima for growth and enzyme activity are often greater than in situ temperatures in their environment. Here we discuss fundamental constraints on microbial adaptation and suggest specific environments where microbial adaptation to climate change (or lack thereof) is most likely to alter ecosystem functioning. Our framework is based on two principal assumptions. First, there are fundamental ecological trade-offs in microbial community traits that occur across environmental gradients (in time and space). These trade-offs result in shifting of microbial function (e.g. ability to take up resources at low temperature) in response to adaptation of another trait (e.g. limiting maintenance respiration at high temperature). Second, the mechanism and level of microbial community adaptation to changing environmental parameters is a function of the potential rate of change in community composition relative to the rate of environmental change. Together, this framework provides a basis for developing testable predictions about how the rate and degree of microbial adaptation to climate change will alter biogeochemical processes in aquatic and terrestrial ecosystems across the planet.

Biogeochemistry↗

Multidecadal geomorphic evolution of a profoundly disturbed gravel-bed river system—a complex, nonlinear response and its impact on sediment delivery

A 2.5-km3 debris avalanche during the 1980 eruption of Mount St. Helens reset the fluvial landscape of upper North Fork Toutle River valley. Since then, a new drainage network has formed and evolved. Cross-section surveys repeated over nearly 40 years at 16 locations along a 20-km reach of river valley document channel evolution, geomorphic processes, and their impacts on sediment delivery. We analyzed spatial and temporal changes in channel morphology using two new metrics: 1) a shape index that defines the degree of U-shaped or V-shaped valley geometry; and 2) an alluvial phase-space diagram that relates bed degradation or aggradation between consecutive surveys to increases or decreases in cross-section area. Metric relations reveal more diverse channel evolution than originally described by a simple, linear-response model of sequential channel initiation and incision; aggradation and widening; and subsequent episodic scour and fill with little change in bed elevation. Instead, vertical and lateral adjustments have been crucial processes intertwined throughout channel evolution. Channel evolution has followed a distinctly nonlinear and non-sequential trajectory, migrating through several phase spaces and involving varied combinations of (1) degradation and aggradation with widening and narrowing, (2) bed-level fluctuations with little change in cross-section area, and (3) changes in cross-section area with little change of bed elevation. Persistent channel widening and reworking of the channel bed presently drive elevated sediment delivery from this basin. Elevated sediment delivery is likely to persist until valley-floor widths greatly exceed that of the channel-migration corridor, and/or channel banks and valley walls stabilize.

Washington↗

On the use of volumetric strain meters to infer additional characteristics of short-period seismic radiation

Volumetric strain meters (Sacks-Evertson design) are installed at 15 sites along the San Andreas fault system, to monitor long-term strain changes for earthquake prediction. Deployment of portable broadband, high-resolution digital recorders (GEOS) at several of the sites extends the detection band for volumetric strain to periods shorter than 5 × 10 −2 sec and permits the simultaneous observation of seismic radiation fields using conventional short-period pendulum seismometers. Simultaneous observations establish that the strain detection bandwidth extends from periods greater than 10 7 seconds to periods near 5 × 10 −2 sec with a dynamic range exceeding 140 dB. Measurements of earth-strain noise for the period band, 10 7 to 10 −2 sec, show that ground noise, not instrument noise, currently limits the measurement of strain over a bandwidth of more than eight orders of magnitude in period. Comparison of the short-period portion of earth-strain, noise spectra (20 to 5 × 10 −2 sec) with average spectra determined from pendulum seismometers, suggest that observed noise is predominantly dilatational energy. Recordings of local and regional earthquakes indicate that dilatometers respond to P energy but not direct shear energy and that straingrams can be used to resolve superimposed reflected P and S waves for inference of wave characteristics not permitted by either sensor alone. Simultaneous measurements of incident P - and S -wave amplitudes are used to introduce a technique for single-station estimates of wave field inhomogeneity, free-surface reflection coefficients and local material P velocity. Estimates of these parameters derived for the North Palm Springs earthquake ( M w 5.9) respectively for an incident P wave of 29° are −85°, 1.71, 2.9 km/sec, and for an incident S wave of 17° are 79°, 0.85, 2.9 km/sec. The empirical estimates of reflection coefficients are consistent with model estimates derived using an anelastic half-space model with incident inhomogeneous wave fields.

California↗

Toward a near-lossless image compression strategy for the NASA/USGS Landsat Next mission

As orbiting Earth imaging platforms carry more complex and capable instruments, efficient methods are needed to reduce the time and cost associated with storing and downlinking greater volumes of image data. The upcoming NASA/USGS Landsat Next mission, with an increase in spatial and spectral resolution over previous Landsat missions, is no exception. Landsat Next will produce nearly six times the amount of image data per day over either of the current Landsat 8 or Landsat 9 observatories. Near-lossless compression, where the image after compression is not identical to the original image, allows for the efficient storage and transmission of all image data while meeting the mission’s global coverage, temporal revisit frequency, and science measurement and performance requirements. Although the Landsat user community is understandably cautious about lossy compression, it is possible to constrain the maximum loss, or error, introduced during compression, ensuring that any added error remains within the intrinsic noise level of the instrument. The Consultative Committee for Space Data Systems image compression standard, CCSDS 123.0-B-2, was chosen for the Landsat Next mission because it is an internationally supported standard suited for integration with space hardware, and it allows control over the magnitude and distribution of compression error. Using several proxy datasets as a surrogate for Landsat Next image data, an investigation was performed to determine a preliminary set of parameter values that would keep the added compression error within acceptable limits. The results of these studies demonstrate that near-lossless image compression can be utilized by the Landsat Next instruments to store and downlink all science data without compromising image quality or mission requirements.

Remote Sensing of Environment↗

Sorted bed forms as self-organized patterns: 2. complex forcing scenarios

We employ a numerical model to study the development of sorted bed forms under a variety of hydrodynamic and sedimentary conditions. Results indicate that increased variability in wave height decreases the growth rate of the features and can potentially give rise to complicated, a priori unpredictable, behavior. This happens because the system responds to a change in wave characteristics by attempting to self-organize into a patterned seabed of different geometry and spacing. The new wavelength might not have enough time to emerge before a new change in wave characteristics occurs, leading to less regular seabed configurations. The new seabed configuration is also highly dependent on the preexisting morphology, which further limits the possibility of predicting future behavior. For the same reasons, variability in the mean current magnitude and direction slows down the growth of features and causes patterns to develop that differ from classical sorted bed forms. Spatial variability in grain size distribution and different types of net sediment aggradation/degradation can also result in the development of sorted bed forms characterized by a less regular shape. Numerical simulations qualitatively agree with observed geometry (spacing and height) of sorted bed forms. Also in agreement with observations is that at shallower depths, sorted bed forms are more likely to be affected by changes in the forcing conditions, which might also explain why, in shallow waters, sorted bed forms are described as ephemeral features. Finally, simulations indicate that the different sorted bed form shapes and patterns observed in the field might not necessarily be related to diverse physical mechanisms. Instead, variations in sorted bed form characteristics may result from variations in local hydrodynamic and/or sedimentary conditions.

Journal of Geophysical Research F: Earth Surface↗

Aeromagnetic map of Mountain Pass and vicinity, California and Nevada

Magnetic investigations of Mountain Pass and vicinity were begun as part of an effort to study regional crustal structures as an aid to understanding the geologic framework and mineral resources of the eastern Mojave Desert. The study area, which straddles the state boundary between southeastern California and southern Nevada, encompasses Mountain Pass, which is host to one of the world’s largest rare earth element carbonatite deposits. The deposit is found along a north-northwest-trending, fault-bounded block that extends along the eastern parts of the Clark Mountain Range, Mescal Range, and Ivanpah Mountains. This Paleoproterozoic block is composed of a 1.7-Ga metamorphic complex of gneiss and schist that underwent widespread metamorphism and associated plutonism during the Ivanpah orogeny. The Paleoproterozoic rocks were intruded by a Mesoproterozoic (1.4 Ga) ultrapotassic alkaline intrusive suite and carbonatite body. The intrusive rocks include, from oldest to youngest, shonkinite, mesosyenite, syenite, quartz syenite, potassic granite, carbonatite, carbonatite dikes, and late shonkinite dikes. Generally speaking, magnetic anomalies reflect lateral changes in subsurface magnetization that can be used to infer subsurface geologic structure, revealing variations in lithology and delineating geologic features such as faults, plutons, volcanic rocks, calderas, and sedimentary basins. A regional aeromagnetic map was derived from statewide aeromagnetic maps of California and Nevada that were compiled from numerous surveys flown at various flightline altitudes and spacings. This compilation, although composed of surveys acquired using different specifications, allows seamless interpretation of magnetic anomalies across survey boundaries. In addition, a high-resolution aeromagnetic survey was flown by helicopter over parts of the Clark Mountain Range, Mescal Range, and Ivanpah Mountains. The resulting mapped magnetic anomalies show in much greater detail the complex subsurface structures in the Mountain Pass area.

California, Nevada↗

Mineral deposits of the Mesoproterozoic Midcontinent Rift System in the Lake Superior region – Metallogeny of the prolifically mineralized Keweenawan LIP

The Keweenawan large igneous province (LIP) of the Midcontinent Rift System (MRS) of North America is perhaps the most prolifically and diversely mineralized LIP known on Earth (Nicholson et al., 1992). The MRS is an approximately 2,200 km curvilinear continental rift that stretches from Kansas northeast to the Lake Superior region where it turns southeast and extends through lower Michigan (Fig. 1). Rocks of the MRS host a varied suite of magmatic and hydrothermal mineral deposits in the Lake Superior region of the United States and Canada where rift rocks are exposed at or near the surface. Historically, hydrothermal deposits, such as Michigan’s native Cu deposits and the White Pine sediment-hosted stratiform Cu deposit, were major MRS metal producers. On-going exploration for and potential development of Cu-Ni sulfide deposits hosted by the Duluth Complex of Minnesota and the opening of the Eagle Ni mine in Michigan indicate an expanding interest in MRS magmatic deposits. Many of the MRS hydrothermal and magmatic mineral deposits are significant past, present, and likely future providers of critical minerals. We have placed these deposits into a space and time metallogenic framework (Woodruff et al., 2020a) that is summarized here.

Iowa, Kansas, Minnesota, Michigan, Nebraska, Ontar↗

InSAR Evidence for an active shallow thrust fault beneath the city of Spokane Washington, USA

In 2001, a nearly five month long sequence of shallow, mostly small magnitude earthquakes occurred beneath the city of Spokane, a city with a population of about 200,000, in the state of Washington. During most of the sequence, the earthquakes were not well located because seismic instrumentation was sparse. Despite poor-quality locations, the earthquake hypocenters were likely very shallow, because residents near the city center both heard and felt many of the earthquakes. The combination of poor earthquake locations and a lack of known surface faults with recent movement make assessing the seismic hazards related to the earthquake swarm difficult. However, the potential for destruction from a shallow moderate-sized earthquake is high, for example Christchurch New Zealand in 2011, so assessing the hazard potential of a seismic structure involved in the Spokane earthquake sequence is important. Using interferometric synthetic aperture radar (InSAR) data from the European Space Agency ERS2 and ENVISAT satellites and the Canadian Space Agency RADARSAT-1, satellite we are able to show that slip on a shallow previously unknown thrust fault, which we name the Spokane Fault, is the source of the earthquake sequence. The part of the Spokane Fault that slipped during the 2001 earthquake sequence underlies the north part of the city, and slip on the fault was concentrated between ~0.3 and 2 km depth. Projecting the buried fault plane to the surface gives a possible surface trace for the Spokane Fault that strikes northeast from the city center into north Spokane.

Washington↗

Compact seismicity bursts have different characteristics from regional seismicity

Earthquakes tend to cluster, developing into sequences driven by stress perturbations and transient fault-zone processes. Depending on the driving process, earthquake sequences show differing behaviors. This variability challenges our ability to observe or distinguish these driving processes in high resolution. Here we systematically identify seismicity bursts throughout southern California using new statistical methods and examine their causes with other independent geophysical observations. Seismicity bursts are defined as statistically significant seismicity rate anomalies. Our approach has the advantage of being data driven, depending on few earthquake occurrence assumptions. We find abundant seismicity bursts across southern California, most frequently occurring along the San Jacinto Fault and in the Salton Sea and Coso geothermal fields. These seismicity bursts are highly compact in space and time, often encompassed by a 5 km radius and have durations less than 10 hr. Many of the seismicity bursts have their largest earthquake near the beginning of the sequence, but the largest earthquake is usually not an obvious mainshock. We utilize a variety of independent geophysical data sets to analyze the characteristics of the seismicity bursts, finding that the seismicity bursts have low b-values, low spectral stress drops, and varied stress ratios compared to regional seismicity. These differences suggest that seismicity bursts are driven by transient processes acting frequently across fault networks.

California↗