USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Earth Observer”

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 847 records · Page 47Linked to original sources
USGS2015

How to build and teach with QuakeCaster: an earthquake demonstration and exploration tool

QuakeCaster is an interactive, hands-on teaching model that simulates earthquakes and their interactions along a plate-boundary fault. QuakeCaster contains the minimum number of physical processes needed to demonstrate most observable earthquake features. A winch to steadily reel in a line simulates the steady plate tectonic motions far from the plate boundaries. A granite slider in frictional contact with a nonskid rock-like surface simulates a fault at a plate boundary. A rubber band connecting the line to the slider simulates the elastic character of the Earth’s crust. By stacking and unstacking sliders and cranking in the winch, one can see the results of changing the shear stress and the clamping stress on a fault. By placing sliders in series with rubber bands between them, one can simulate the interaction of earthquakes along a fault, such as cascading or toggling shocks. By inserting a load scale into the line, one can measure the stress acting on the fault throughout the earthquake cycle. As observed for real earthquakes, QuakeCaster events are not periodic, time-predictable, or slip-predictable. QuakeCaster produces rare but unreliable “foreshocks.” When fault gouge builds up, the friction goes to zero and fault creep is seen without large quakes. QuakeCaster events produce very small amounts of fault gouge that strongly alter its behavior, resulting in smaller, more frequent shocks as the gouge accumulates. QuakeCaster is designed so that students or audience members can operate it and record its output. With a stopwatch and ruler one can measure and plot the timing, slip distance, and force results of simulated earthquakes. People of all ages can use the QuakeCaster model to explore hypotheses about earthquake occurrence. QuakeCaster takes several days and about $500.00 in materials to build.

Open-File Report↗
USGS1998

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Development of partial rock veneers by root throw in a subalpine setting

Rock veneers stabilize hillslope surfaces, occur especially in areas of immature soil, and form through a variety of process sets that includes root throw. Near Westcliffe, Colorado, USA, data were collected from a 20 ?? 500 m transect on the east slope of the Sangre de Cristo Mountains. Ages of pit/mound complexes with rock fragments exposed at the surface by root throw ranged from recent (freshly toppled tree) to unknown (complete tree decay). Calculations based on dimensions of the pit/mound complexes, estimated time of free topppling, sizes of exposed rock fragments, and percentage rock covers at pit/mound complexes, as well as within the transect area, indicate that recent rates of root throw have resulted in only partial rock veneering since late Pleistocene deglaciation. Weathering of rock fragments prevent development of an extensive rock veneer and causes a balance, achieved within an estimated 700 years, between the rates of rock-fragment exposure by root throw and clast disintegration by chemical reduction. The estimated rate of rock-fragment reduction accounts for part of the fluvial sediment yields observed for forested subalpine areas of western North America. Copyright ?? 2005 John Wiley & Sons, Ltd.

Earth Surface Processes and Landforms↗
USGS1992

The Cretaceous-Tertiary extinction: A lethal mechanism involving anhydrite target rocks

The Chicxulub Crater, Yucatan, Mexico, is a leading contender as the site for the impact event that caused the Cretaceous-Tertiary (K-T) extinctions. A considerable thickness of anhydrite (CaSO 4 ) forms part of the target rock. High temperatures resulting from impact would drive SO 2 off from the anhydrite. Hundreds of billions of tonnes of sulfuric acid aerosol would thus enter the stratosphere and cause considerable cooling of the Earth's surface, decrease photosynthesis by orders of magnitude, deplete the ozone layer, and permit increased UV radiation to reach the Earth's surface. Finally, the aerosol would fall back to Earth as acid rain and devastate land and some lacustrine biota and near-surface marine creatures. The presence of anhydrite in the Chicxulub target rock may thus help explain the many extinctions observed at the K-T boundary.

Geochimica et Cosmochimica Acta↗
USGS2015

The effect of pressurized magma chamber growth on melt migration and pre-caldera vent locations through time at Mount Mazama, Crater Lake, Oregon

The pattern of eruptions at long-lived volcanic centers provides a window into the co-evolution of crustal magma transport, tectonic stresses, and unsteady magma generation at depth. Mount Mazama in the Oregon Cascades has seen variable activity over the last 400 ky, including the 50 km 3 climactic eruption at ca. 7.7 ka that produced Crater Lake caldera. The physical mechanisms responsible for the assembly of silicic magma reservoirs that are the precursors to caldera-forming eruptions are poorly understood. Here we argue that the spatial and temporal distribution of geographically clustered volcanic vents near Mazama reflects the development of a centralized magma chamber that fed the climactic eruption. Time-averaged eruption rates at Mount Mazama imply an order of magnitude increase in deep magma influx prior to the caldera-forming event, suggesting that unsteady mantle melting triggered a chamber growth episode that culminated in caldera formation. We model magma chamber–dike interactions over ∼50 ky preceding the climactic eruption to fit the observed distribution of surface eruptive vents in space and time, as well as petrologically estimated deep influx rates. Best fitting models predict an expanding zone of dike capture caused by a growing, oblate spheroidal magma chamber with 10–30 MPa of overpressure. This growing zone of chamber influence causes closest approaching regional mafic vent locations as well as more compositionally evolved Mazama eruptions to migrate away from the climactic eruptive center, returning as observed to the center after the chamber drains during the caldera-forming eruption.

Oregon↗
USGS2024

A great tsunami earthquake component of the 1957 Aleutian Islands earthquake

The great 1957 Aleutian Islands earthquake ruptured ∼1200 km of the plate boundary along the Aleutian subduction zone and produced a destructive tsunami across Hawaiʻi. Early seismic and tsunami analyses indicated that large megathrust fault slip was concentrated in the western Aleutian Islands, but tsunami waves generated by slip in the west cannot explain the large observed runup in Hawaiʻi far to the southeast. Recently mapped 1957 geologic deposits on eastern Aleutian Islands suggest occurrence of very large nearby slip. Jointly modeling tsunami runup along the eastern Aleutian and Hawaiian Islands together with tide gauge recordings across the Pacific resolves 12-26 m shallow slip along 600 km of the eastern Aleutian Islands in addition to modest, deeper western slip inferred from seismic records. The eastern near-trench slip results in an M W 8.3-8.6 tsunami earthquake component of the M W 8.6-8.8 rupture, comparable in size to the adjacent 1946 Aleutian tsunami earthquake to the east. The reexamination of the 1957 rupture confirms the tsunami hazards posed by the eastern Aleutian subduction zone to Hawaiʻi and lays the groundwork for investigation of large prehistoric earthquakes through modeling tsunami runup inferred from stratigraphic observations to constrain their rupture processes.

Alaska↗
USGS2004

Glacier ice mass fluctuations and fault instability in tectonically active Southern Alaska

Across the plate boundary zone in south central Alaska, tectonic strain rates are high in a region that includes large glaciers undergoing wastage (glacier retreat and thinning) and surges. For the coastal region between the Bering and Malaspina Glaciers, the average ice mass thickness changes between 1995 and 2000 range from 1 to 5 m/year. These ice changes caused solid Earth displacements in our study region with predicted values of -10 to 50 mm in the vertical and predicted horizontal displacements of 0-10 mm at variable orientations. Relative to stable North America, observed horizontal rates of tectonic deformation range from 10 to 40 mm/year to the north-northwest and the predicted tectonic uplift rates range from approximately 0 mm/year near the Gulf of Alaska coast to 12 mm/year further inland. The ice mass changes between 1995 and 2000 resulted in discernible changes in the Global Positioning System (GPS) measured station positions of one site (ISLE) located adjacent to the Bagley Ice Valley and at one site, DON, located south of the Bering Glacier terminus. In addition to modifying the surface displacements rates, we evaluated the influence ice changes during the Bering glacier surge cycle had on the background seismic rate. We found an increase in the number of earthquakes (ML???2.5) and seismic rate associated with ice thinning and a decrease in the number of earthquakes and seismic rate associated with ice thickening. These results support the hypothesis that ice mass changes can modulate the background seismic rate. During the last century, wastage of the coastal glaciers in the Icy Bay and Malaspina region indicates thinning of hundreds of meters and in areas of major retreat, maximum losses of ice thickness approaching 1 km. Between the 1899 Yakataga and Yakutat earthquakes (Mw=8.1, 8.1) and prior to the 1979 St. Elias earthquake (M s=7.2), the plate interface below Icy Bay was locked and tectonic strain accumulated. We used estimated ice mass change during the 1899-1979 time period to calculate the change in the fault stability margin (FSM) prior to the 1979 St. Elias earthquake. Our results suggest that a cumulative decrease in the fault stability margin at seismogenic depths, due to ice wastage over 80 years, was large, up to ???2 MPa. Ice wastage would promote thrust faulting in events such as the 1979 earthquake and subsequent aftershocks.

Conference Paper↗
USGS2002

Draped aeromagnetic survey in Transantarctic Mountains over the area of the Butcher Ridge igneous complex showing extent of underlying mafic intrusion

A draped aeromagnetic survey over the area surrounding the Butcher Ridge igneous complex (BRIC), Transantarctic Mountains, was acquired in 1997–1998 as part of a larger Transantarctic Mountains Aerogeophysical Research Activity survey. The BRIC is a sill-like hypoabyssal intrusion ranging in composition from tholeiitic basalt to rhyolite. An 40 Ar/ 39 Ar age of 174 Ma and the chemical character of the basaltic rocks show the BRIC to be part of the widespread Jurassic Ferrar suite of continental tholeiitic rocks, that extends for 3500 km across Antarctica. The aeromagnetic survey shows a horseshoe-shaped pattern of anomalies reaching amplitudes as great as 1900 nT generally associated with the bedrock topography where it is exposed. It is apparent that the high-amplitude anomaly pattern is more extensive than the 10-km-long exposed outcrop, first crossed by a single 1990 aeromagnetic profile. The highest-amplitude anomalies appear south of the profile acquired in 1990 and extend out of the survey area. The new aeromagnetic data allow determination of the extent of the interpreted Butcher mafic(?) intrusion beneath exposures of Beacon sedimentary rock and ice in the area covered, as well as beneath the small BRIC exposure. The magnetic anomalies show a minimum area of 3000 km 2 , a much greater extent than previously inferred. Magnetic models indicate a minimum thickness of ∼1–2 km for a horizontal intrusion. However, nonunique models with magnetic layers decreasing in apparent susceptibility with depth are consistent with of a 4- to 8-km-thick layered intrusion. These magnetic models indicate progressively deeper erosion of the interpreted mafic-layered body from the south to north. The erosion has removed more magnetic upper layers that mask the magnetic effects of the lower less magnetic layers. The probable minimum volume of the intrusion in the area of the survey is ∼6000 km 3 . An alternate, but less likely, interpretation of a series of dikes can also fit the observed magnetic anomalies.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

The velocity field along the San Andreas Fault in central and southern California

The velocity field within a 100‐km‐broad zone centered on the San Andreas fault between the Mexican border and San Francisco Bay has been inferred from repeated surveys of trilateration networks in the 1973–1989 interval. The velocity field has the appearance of a shear flow that remains parallel to the local strike of the fault even through such major deflections as the big bend of the San Andreas fault in the Transverse Ranges of southern California. Across‐strike profiles of the fault‐parallel component of velocity exhibit the expected sigmoidal shape, whereas across‐strike profiles of the fault‐normal component of velocity are flat and featureless. No significant convergence upon the fault is observed even along the big bend sector of the fault. Simple dislocation models can explain most of the features of the observed velocity field, but those explanations are not unique. About 35 mm/yr of relative plate motion is accounted for within the span of the trilateration networks. Geologic studies indicate that the secular slip rate on the San Andreas fault is about 35 mm/yr. The agreement between these two estimates implies that most of the strain accumulation is elastic and will be recovered in subsequent earthquakes. The relative motion observed across the San Andreas fault (35 mm/yr) plus that observed across the Eastern California shear zone (8 mm/yr) accounts for most (43 mm/yr) of the observed North America‐Pacific relative plate motion (47 mm/yr).

California↗
USGS2024

Signatures of wave erosion in Titan’s coasts

The shorelines of Titan’s hydrocarbon seas trace flooded erosional landforms such as river valleys; however, it is unclear whether coastal erosion has subsequently altered these shorelines. Spacecraft observations and theoretical models suggest that wind may cause waves to form on Titan’s seas, potentially driving coastal erosion, but the observational evidence of waves is indirect, and the processes affecting shoreline evolution on Titan remain unknown. No widely accepted framework exists for using shoreline morphology to quantitatively discern coastal erosion mechanisms, even on Earth, where the dominant mechanisms are known. We combine landscape evolution models with measurements of shoreline shape on Earth to characterize how different coastal erosion mechanisms affect shoreline morphology. Applying this framework to Titan, we find that the shorelines of Titan’s seas are most consistent with flooded landscapes that subsequently have been eroded by waves, rather than a uniform erosional process or no coastal erosion, particularly if wave growth saturates at fetch lengths of tens of kilometers.

Science Advances↗
USGS2020

The Modern Geological Survey; a model for research, innovation, synthjesis: A USGS perspective

Geological Surveys have long filled the role of providing Earth system science data and knowledge. These functions are increasingly complicated by accelerating environmental and societal change. Here we describe the USGS response to these evolving conditions. Underpinning the USGS approach is the recognition that many of the issues facing the U.S. and the world involve the interaction among geologic, hydrologic, and biologic processes, and how these interactions in turn affect society. Therefore, a goal of USGS planning is fostering interdisciplinary science. This focus is occurring in part through implementation of the recommendations of strategic planning teams. The USGS has also put in place groups building a broad information technology infrastructure as well as identifying and disseminating new Earth science research tools. In addition, the USGS has established an analysis and synthesis center that brings together groups of scientists who address interdisciplinary Earth system science issues. The goal is for these building blocks to evolve towards a comprehensive USGS data and knowledge platform; EarthMAP (Earth Monitoring, Assessment, and Projection). We also recognize that the modern geological survey must be a member of a community of geological surveys contributing data to a global database of 3-dimensional biogeophysical observations and interpretations.

Book chapter↗
USGS2006

Continuous borehole strain and pore pressure in the near field of the 28 September 2004 M 6.0 Parkfield, California, earthquake: Implications for nucleation, fault response, earthquake prediction and tremor

Near-field observations of high-precision borehole strain and pore pressure, show no indication of coherent accelerating strain or pore pressure during the weeks to seconds before the 28 September 2004 M 6.0 Parkfield earthquake. Minor changes in strain rate did occur at a few sites during the last 24 hr before the earthquake but these changes are neither significant nor have the form expected for strain during slip coalescence initiating fault failure. Seconds before the event, strain is stable at the 10 −11 level. Final prerupture nucleation slip in the hypocentral region is constrained to have a moment less than 2 × 10 12 N m ( M 2.2) and a source size less than 30 m. Ground displacement data indicate similar constraints. Localized rupture nucleation and runaway precludes useful prediction of damaging earthquakes. Coseismic dynamic strains of about 10 microstrain peak-to-peak were superimposed on volumetric strain offsets of about 0.5 microstrain to the northwest of the epicenter and about 0.2 microstrain to the southeast of the epicenter, consistent with right lateral slip. Observed strain and Global Positioning System ( gps ) offsets can be simply fit with 20 cm of slip between 4 and 10 km on a 20-km segment of the fault north of Gold Hill ( M 0 = 7 × 10 17 N m). Variable slip inversion models using gps data and seismic data indicate similar moments. Observed postseismic strain is 60% to 300% of the coseismic strain, indicating incomplete release of accumulated strain. No measurable change in fault zone compliance preceding or following the earthquake is indicated by stable earth tidal response. No indications of strain change accompany nonvolcanic tremor events reported prior to and following the earthquake.

California↗
USGS2014

1964 Great Alaska Earthquake: a photographic tour of Anchorage, Alaska

On March 27, 1964, at 5:36 p.m., a magnitude 9.2 earthquake, the largest recorded earthquake in U.S. history, struck southcentral Alaska (fig. 1). The Great Alaska Earthquake (also known as the Good Friday Earthquake) occurred at a pivotal time in the history of earth science, and helped lead to the acceptance of plate tectonic theory (Cox, 1973; Brocher and others, 2014). All large subduction zone earthquakes are understood through insights learned from the 1964 event, and observations and interpretations of the earthquake have influenced the design of infrastructure and seismic monitoring systems now in place. The earthquake caused extensive damage across the State, and triggered local tsunamis that devastated the Alaskan towns of Whittier, Valdez, and Seward. In Anchorage, the main cause of damage was ground shaking, which lasted approximately 4.5 minutes. Many buildings could not withstand this motion and were damaged or collapsed even though their foundations remained intact. More significantly, ground shaking triggered a number of landslides along coastal and drainage valley bluffs underlain by the Bootlegger Cove Formation, a composite of facies containing variably mixed gravel, sand, silt, and clay which were deposited over much of upper Cook Inlet during the Late Pleistocene (Ulery and others, 1983). Cyclic (or strain) softening of the more sensitive clay facies caused overlying blocks of soil to slide sideways along surfaces dipping by only a few degrees. This guide is the document version of an interactive web map that was created as part of the commemoration events for the 50th anniversary of the 1964 Great Alaska Earthquake. It is accessible at the U.S. Geological Survey (USGS) Alaska Science Center website: http://alaska.usgs.gov/announcements/news/1964Earthquake/ . The website features a map display with suggested tour stops in Anchorage, historical photographs taken shortly after the earthquake, repeat photography of selected sites, scanned documents, and small-scale maps, as well as links to slideshows of additional photographs and Google Street View™ scenes. Buildings in Anchorage that were severely damaged, sites of major landslides, and locations of post-earthquake engineering responses are highlighted. The web map can be used online as a virtual tour or in a physical self-guided tour using a web-enabled Global Positioning System (GPS) device. This publication serves the purpose of committing most of the content of the web map to a single distributable document. As such, some of the content differs from the online version.

Alaska↗
USGS2006

The sand seas of titan: Cassini RADAR observations of longitudinal dunes

The most recent Cassini RADAR images of Titan show widespread regions (up to 1500 kilometers by 200 kilometers) of near-parallel radar-dark linear features that appear to be seas of longitudinal dunes similar to those seen in the Namib desert on Earth. The Ku-band (2.17-centimeter wavelength) images show ∼100-meter ridges consistent with duneforms and reveal flow interactions with underlying hills. The distribution and orientation of the dunes support a model of fluctuating surface winds of ∼0.5 meter per second resulting from the combination of an eastward flow with a variable tidal wind. The existence of dunes also requires geological processes that create sand-sized (100- to 300-micrometer) particulates and a lack of persistent equatorial surface liquids to act as sand traps.

Science↗
USGS2001

hypoDD-A Program to Compute Double-Difference Hypocenter Locations

HypoDD is a Fortran computer program package for relocating earthquakes with the double-difference algorithm of Waldhauser and Ellsworth (2000). This document provides a brief introduction into how to run and use the programs ph2dt and hypoDD to compute double-difference (DD) hypocenter locations. It gives a short overview of the DD technique, discusses the data preprocessing using ph2dt, and leads through the earthquake relocation process using hypoDD. The appendices include the reference manuals for the two programs and a short description of auxiliary programs and example data. Some minor subroutines are presently in the c language, and future releases will be in c. Earthquake location algorithms are usually based on some form of Geiger’s method, the linearization of the travel time equation in a first order Taylor series that relates the difference between the observed and predicted travel time to unknown adjustments in the hypocentral coordinates through the partial derivatives of travel time with respect to the unknowns. Earthquakes can be located individually with this algorithm, or jointly when other unknowns link together the solutions to indivdual earthquakes, such as station corrections in the joint hypocenter determination (JHD) method, or the earth model in seismic tomography. The DD technique (described in detail in Waldhauser and Ellsworth, 2000) takes advantage of the fact that if the hypocentral separation between two earthquakes is small compared to the event-station distance and the scale length of velocity heterogeneity, then the ray paths between the source region and a common station are similar along almost the entire ray path (Fréchet, 1985; Got et al., 1994). In this case, the difference in travel times for two events observed at one station can be attributed to the spatial offset between the events with high accuracy. DD equations are built by differencing Geiger’s equation for earthquake location. In this way, the residual between observed and calculated travel-time difference (or double-difference) between two events at a common station are a related to adjustments in the relative position of the hypocenters and origin times through the partial derivatives of the travel times for each event with respect to the unknown. HypoDD calculates travel times in a layered velocity model (where velocity depends only on depth) for the current hypocenters at the station where the phase was recorded. The double-difference residuals for pairs of earthquakes at each station are minimized by weighted least squares using the method of singular value decomposition (SVD) or the conjugate gradients method (LSQR, Paige and Saunders, 1982). Solutions are found by iteratively adjusting the vector difference between nearby hypocentral pairs, with the locations and partial derivatives being updated after each iteration. Details about the algorithm can be found in Waldhauser and Ellsworth (2000). When the earthquake location problem is linearized using the double-difference equations, the common mode errors cancel, principally those related to the receiver-side structure. Thus we avoid the need for station corrections or high-accuracy of predicted travel times for the portion of the raypath that lies outside the focal volume. This approach is especially useful in regions with a dense distribution of seismicity, i.e. where distances between neighboring events are only a few hundred meters. The improvement of double-difference locations over ordinary JHD locations is shown in Figure 1 for about 10,000 earthquakes that occurred during the 1997 seismic crisis in the Long Valley caldera, California. While the JHD locations (left panel) show a diffuse picture of the seismicity, double-difference locations (right panel) bring structural details such as the location of active fault planes into sharp focus.

Open-File Report↗
USGS2025

Petrogenesis and mineralization potential of spinifex komatiitic basalts in the Bradley Peak greenstone terrane, Wyoming Province

Komatiitic volcanic rocks are important hosts of Ni sulfide mineralization and record early Earth evolution; however, those in the well-studied Archean Wyoming Province have received little attention. Here, we elucidate the timing and petrogenesis of the Bradley Peak komatiitic volcanic rocks using field and textural observations, geochronology, and geochemistry. Detrital and igneous zircon U-Pb ages for two samples from previously undated units support published age determinations, placing the eruption age at 2.72 Ga. Stratigraphy of the volcanic flows was mapped and 36 samples including cumulates, greenschists, and spinifex-textured rocks were collected. Whole-rock geochemistry was used to classify the spinifex-textured samples as Al-undepleted komatiitic basalts (11–17 wt% MgO). Platinum-group element concentrations (n = 25) are like those in global Al-undepleted komatiitic basalts, and PGE/Ti ratios do not indicate the volcanic flows likely host sulfide mineralization. Initial ε Nd values of −0.5 to +4.7 (n = 16), indicate that these lavas were derived from a depleted mantle source and have negligible evolved crust contamination. The primary magma to the komatiitic basalt flows is estimated to have had 19 wt% MgO and be derived from ∼15 to 25 % mantle partial melting at 3–4 GPa. Trace element chemistry and thermodynamic modeling suggest the primary melt assimilated local banded iron formation. Although the Bradley Peak komatiitic basalts do not contain positive evidence of magmatic sulfide deposits, depleted Au in the flows suggests they could be source rocks for nearby orogenic gold deposits.

Wyoming↗
USGS2019

Dome formation on Ceres by sold-state flow analogous to terrestrial salt tectonics

The dwarf planet Ceres’s outer crust is a complex, heterogeneous mixture of ice, clathrates, salts and silicates. Numerous large domes on Ceres’s surface indicate a degree of geological activity. These domes have been attributed to cryovolcanism, but that is difficult to reconcile with Ceres’s small size and lack of long-lived heat sources. Here we alternatively propose that Ceres’s domes form by solid-state flow within the compositionally heterogeneous crust, a mechanism directly analogous to salt tectonics on Earth. We use numerical simulations to illustrate that differential loading of a crust with compositional heterogeneity on a scale of tens of kilometres can produce dome-like features of scale similar to those observed. The mechanism requires the presence of low-viscosity and low-density, possibly ice-rich, material in the upper 1–10 km of the subsurface. Such substantial regional heterogeneity in Ceres’s crustal composition is consistent with observations from the National Aeronautics and Space Administration’s Dawn mission. We conclude that deformation analogous to that in terrestrial salt tectonics is a viable alternative explanation for the observed surface morphologies, and is consistent with Ceres being both cold and geologically active.

Nature Geoscience↗
USGS2013

Projected future changes in vegetation in western North America in the 21st century

Rapid and broad-scale forest mortality associated with recent droughts, rising temperature, and insect outbreaks has been observed over western North America (NA). Climate models project additional future warming and increasing drought and water stress for this region. To assess future potential changes in vegetation distributions in western NA, the Community Earth System Model (CESM) coupled with its Dynamic Global Vegetation Model (DGVM) was used under the future A2 emissions scenario. To better span uncertainties in future climate, eight sea surface temperature (SST) projections provided by phase 3 of the Coupled Model Intercomparison Project (CMIP3) were employed as boundary conditions. There is a broad consensus among the simulations, despite differences in the simulated climate trajectories across the ensemble, that about half of the needleleaf evergreen tree coverage (from 24% to 11%) will disappear, coincident with a 14% (from 11% to 25%) increase in shrubs and grasses by the end of the twenty-first century in western NA, with most of the change occurring over the latter half of the twenty-first century. The net impact is a ~6 GtC or about 50% decrease in projected ecosystem carbon storage in this region. The findings suggest a potential for a widespread shift from tree-dominated landscapes to shrub and grass-dominated landscapes in western NA because of future warming and consequent increases in water deficits. These results highlight the need for improved process-based understanding of vegetation dynamics, particularly including mortality and the subsequent incorporation of these mechanisms into earth system models to better quantify the vulnerability of western NA forests under climate change.

Journal of Climate↗