Geology Reports⌕ Search

SEARCH · Geology Reports

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 1,387 records · Page 77Linked to original sources

A comprehensive survey of faults, breccias, and fractures in and flanking the eastern Española Basin, Rio Grande rift, New Mexico

A comprehensive survey of geologic structures formed in the Earth’s brittle regime in the eastern Española Basin and flank of the Rio Grande rift, New Mexico, reveals a complex and protracted record of multiple tectonic events. Data and analyses from this representative rift flank-basin pair include measurements from 53 individual fault zones and 22 other brittle structures, such as breccia zones, joints, and veins, investigated at a total of just over 100 sites. Structures were examined and compared in poorly lithified Tertiary sediments, as well as in Paleozoic sedimentary and Proterozoic crystalline rocks. Data and analyses include geologic maps; field observations and measurements; orientation, kinematic, and paleostress analyses; statistical examination of fault trace lengths derived from aeromagnetic data; mineralogy and chemistry of host and fault rocks; and investigation of fault versus bolide-impact hypotheses for the origin of enigmatic breccias found in the Proterozoic basement rocks. Fault kinematic and paleostress analyses suggest a record of transitional, and perhaps partitioned, strains from the Laramide orogeny through Rio Grande rifting. Normal faults within Tertiary basin-fill sediments are consistent with more typical WNW-ESE Rio Grande rift extension, perhaps decoupled from bedrock structures due to strength contrasts favoring the formation of new faults in the relatively weak sediments. Analyses of the fault-length data indicate power-law length distributions similar to those reported from many geologic settings globally. Mineralogy and chemistry in Proterozoic fault-related rocks reveal geochemical changes tied to hydrothermal alteration and nearly isochemical transformation of feldspars to clay minerals. In sediments, faulted minerals are characterized by mechanical entrainment with minor secondary chemical changes. Enigmatic breccias in rift-flanking Proterozoic rocks are autoclastic and isochemical with respect to their protoliths and exist near shatter cones believed to be related to a previously reported pre-Pennsylvanian impact event. A weak iridium anomaly is associated with the breccias as well as adjacent protoliths, thus an impact shock wave cannot be ruled out for their origin. Major fault zones along the eastern rift-flank mountain front are discontinuous and unlikely to impede regional groundwater flow into Española Basin aquifers. The breccia bodies are not large enough to constitute aquifers, and no fault- or breccia-related geochemical anomalies were identified as potential contamination sources for ground or surface waters. The results of this work provide a broad picture of structural diversity and tectonic evolution along the eastern flank of the central Rio Grande rift and the adjacent Española Basin representative of the rift as a whole and many rifts worldwide.

New Mexico↗

Pore-pressure sensitivities to dynamic strains: observations in active tectonic regions

Triggered seismicity arising from dynamic stresses is often explained by the Mohr-Coulomb failure criterion, where elevated pore pressures reduce the effective strength of faults in fluid-saturated rock. The seismic response of a fluid-rock system naturally depends on its hydro-mechanical properties, but accurately assessing how pore-fluid pressure responds to applied stress over large scales in situ remains a challenging task; hence, spatial variations in response are not well understood, especially around active faults. Here I analyze previously unutilized records of dynamic strain and pore-pressure from regional and teleseismic earthquakes at Plate Boundary Observatory (PBO) stations from 2006 through 2012 to investigate variations in response along the Pacific/North American tectonic plate boundary. I find robust scaling-response coefficients between excess pore pressure and dynamic strain at each station that are spatially correlated: around the San Andreas and San Jacinto fault systems, the response is lowest in regions of the crust undergoing the highest rates of secular shear strain. PBO stations in the Parkfield instrument cluster are at comparable distances to the San Andreas fault (SAF), and spatial variations there follow patterns in dextral creep rates along the fault, with the highest response in the actively creeping section, which is consistent with a narrowing zone of strain accumulation seen in geodetic velocity profiles. At stations in the San Juan Bautista (SJB) and Anza instrument clusters, the response depends non-linearly on the inverse fault-perpendicular distance, with the response decreasing towards the fault; the SJB cluster is at the northern transition from creeping-to-locked behavior along the SAF, where creep rates are at moderate to low levels, and the Anza cluster is around the San Jacinto fault, where to date there have been no statistically significant creep rates observed at the surface. These results suggest that the strength of the pore pressure response in fluid-saturated rock near active faults is controlled by shear strain accumulation associated with tectonic loading, which implies a strong feedback between fault strength and permeability: dynamic triggering susceptibilities may vary in space and also in time.

Journal of Geophysical Research B: Solid Earth↗

A submarine landslide source for the devastating 1964 Chenega tsunami, southern Alaska

During the 1964 Great Alaska earthquake (M w 9.2), several fjords, straits, and bays throughout southern Alaska experienced significant tsunami runup of localized, but unexplained origin. Dangerous Passage is a glacimarine fjord in western Prince William Sound, which experienced a tsunami that devastated the village of Chenega where 23 of 75 inhabitants were lost – the highest relative loss of any community during the earthquake. Previous studies suggested the source of the devastating tsunami was either from a local submarine landslide of unknown origin or from coseismic tectonic displacement. Here we present new observations from high-resolution multibeam bathymetry and seismic reflection surveys conducted in the waters adjacent to the village of Chenega. The seabed morphology and substrate architecture reveal a large submarine landslide complex in water depths of 120–360 m. Analysis of bathymetric change between 1957 and 2014 indicates the upper 20–50 m ( ∼ 0.7 km 3 ) of glacimarine sediment was destabilized and evacuated from the steep face of a submerged moraine and an adjacent ∼ 21 km 2 perched sedimentary basin. Once mobilized, landslide debris poured over the steep, 130 m-high face of a deeper moraine and then blanketed the terminal basin ( ∼ 465 m water depth) in 11 ± 5 m of sediment. These results, combined with inverse tsunami travel-time modeling, suggest that earthquake- triggered submarine landslides generated the tsunami that struck the village of Chenega roughly 4 min after shaking began. Unlike other tsunamigenic landslides observed in and around Prince William Sound in 1964, the failures in Dangerous Passage are not linked to an active submarine delta. The requisite environmental conditions needed to generate large submarine landslides in glacimarine fjords around the world may be more common than previously thought.

Alaska↗

Precision of geodolite surveys: A reply to Jackson and Cheng

Jackson and Cheng have suggested that the changes in areal dilatation measured in Geodolite surveys by the U.S. Geological Survey may be simply an artifact of the measuring system. Although systematic error could conceivably account for the observed excursions in dilatation, we maintain that the specific criticisms by Jackson and Cheng are incorrect: the excursions in dilatation cannot be attributed to the offset correction nor to proportional error associated with temperature. The absence of both errors is demonstrated by using data that are particularly sensitive to the two effects.

Journal of Geophysical Research B: Solid Earth↗

A new dry hypothesis for the formation of Martian linear gullies

Long, narrow grooves found on the slopes of martian sand dunes have been cited as evidence of liquid water via the hypothesis that melt-water initiated debris flows eroded channels and deposited lateral levées. However, this theory has several short-comings for explaining the observed morphology and activity of these linear gullies. We present an alternative hypothesis that is consistent with the observed morphology, location, and current activity: that blocks of CO 2 ice break from over-steepened cornices as sublimation processes destabilize the surface in the spring, and these blocks move downslope, carving out levéed grooves of relatively uniform width and forming terminal pits. To test this hypothesis, we describe experiments involving water and CO 2 blocks on terrestrial dunes and then compare results with the martian features. Furthermore, we present a theoretical model of the initiation of block motion due to sublimation and use this to quantitatively compare the expected behavior of blocks on the Earth and Mars. The model demonstrates that CO 2 blocks can be expected to move via our proposed mechanism on the Earth and Mars, and the experiments show that the motion of these blocks will naturally create the main morphological features of linear gullies seen on Mars.

Icarus↗

Frictional strength heterogeneity and surface heat flow: Implications for the strength of the creeping San Andreas fault

Heat flow measurements along much of the San Andreas fault (SAF) constrain the apparent coefficient of friction (??app) of the fault to <0.2, much lower than laboratory-derived friction values for most geologic materials. However, heat flow data are sparse near the creeping section of the SAF, a frictional "asperity" where the fault slips almost exclusively by aseismic creep. We test the hypothesis that the creeping section has a substantially higher or lower ?? app than adjacent sections of the SAF. We use numerical models to explore the effects of faults with spatially and temporally heterogeneous frictional strength on the spatial distribution of surface heat flow. Heat flow from finite length asperities is uniformly lower than predicted by assuming an infinitely long fault. Over geologic time, lateral offset from strike-slip faulting produces heat flow patterns that are asymmetric across the fault and along strike. We explore a range of asperity sizes, slip rates, and displacement histories for comparing predicted spatial patterns of heat flow with existing measurements. Models with ??app ??? 0.1 fit the data best. For most scenarios, heat flow anomalies from a frictional asperity with ??app > 0.2 should be detectable even with the sparse existing observations, implying that ??app for the creeping section is as low as the surrounding SAF. Because the creeping section does not slip in large earthquakes, the mechanism controlling its weakness is not related to dynamic processes resulting from high slip rate earthquake ruptures. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Seismic reflection/refraction mapping of faulting and regional dips in the eastern Alaska Range

We present the results of a Trans‐Alaska Crustal Transect (TACT) investigation of the upper 2–5 km of the eastern Alaska Range in the vicinity of the Denali fault based on seismic reflection/refraction data, laboratory measurements of rock velocities, and structural mapping. The Denali fault is a major dextral slip structure mappable for more than 2000 km separating the Wrangellia and adjacent terranes to the south from the Yukon‐Tanana and adjacent terranes to the north. Geologic mapping suggests over 400 km of dextral slip has occurred on the fault, yet within the upper 1.5 km of the crust along the TACT corridor, basement rocks juxtaposed along the Denali fault reveal no significant seismic velocity differences, although the fault zone itself is associated with a minor lowering of velocity. The lack of seismic velocity contrast adjacent to the fault is in agreement with laboratory measurements of elastic wave velocities of samples from terranes bordering the fault. Laboratory measurements of elastic wave velocities of the metasedimentary mica‐quartz schists comprising the Yukon‐Tanana basement are highly anisotropic because of preferred orientation of mica and predict significant variations in velocity accompanying variations in foliation dip. Although other interpretations are possible, the northward shallowing of foliation dip of basement rocks in the Yukon‐Tanana terrane combined with the strong anisotropy associated with these highly foliated rocks can explain an observed northward increase in seismic velocity within this terrane. Seismic reflections from basement rocks within the Yukon‐Tanana terrane may originate from variations in anisotropy with depth and/or changes in composition reflecting different proportions of sandstone and shale in the protolith.

Alaska↗

Non-double-couple earthquakes. 1. Theory

Historically, most quantitative seismological analyses have been based on the assumption that earthquakes are caused by shear faulting, for which the equivalent force system in an isotropic medium is a pair of force couples with no net torque (a “double couple,” or DC). Observations of increasing quality and coverage, however, now resolve departures from the DC model for many earthquakes and find some earthquakes, especially in volcanic and geothermal areas, that have strongly non-DC mechanisms. Understanding non-DC earthquakes is important both for studying the process of faulting in detail and for identifying nonshear-faulting processes that apparently occur in some earthquakes. This paper summarizes the theory of “moment tensor” expansions of equivalent-force systems and analyzes many possible physical non-DC earthquake processes. Contrary to long-standing assumption, sources within the Earth can sometimes have net force and torque components, described by first-rank and asymmetric second-rank moment tensors, which must be included in analyses of landslides and some volcanic phenomena. Non-DC processes that lead to conventional (symmetric second-rank) moment tensors include geometrically complex shear faulting, tensile faulting, shear faulting in an anisotropic medium, shear faulting in a heterogeneous region (e.g., near an interface), and polymorphic phase transformations. Undoubtedly, many non-DC earthquake processes remain to be discovered. Progress will be facilitated by experimental studies that use wave amplitudes, amplitude ratios, and complete waveforms in addition to wave polarities and thus avoid arbitrary assumptions such as the absence of volume changes or the temporal similarity of different moment tensor components.

Reviews of Geophysics↗

High spectral resolution reflectance spectroscopy of minerals

The reflectance spectra of minerals are studied as a function of spectral resolution in the range from 0.2 to 3.0 μm. Selected absorption bands were studied at resolving powers (λ/Δλ) as high as 2240. At resolving powers of approximately 1000, many OH‐bearing minerals show diagnostic sharp absorptions at the resolution limit. At low resolution, some minerals may not be distinguishable, but as the resolution is increased, most can be easily identified. As the resolution is increased, many minerals show fine structure, particularly in the OH‐stretching overtone region near 1.4 μm. The fine structure can enhance the ability to discriminate between minerals, and in some cases the fine structure can be used to determine elemental composition. For example, in amphiboles and talcs, four absorption bands are observed in the samples analyzed in this study that are due to hydroxyl linked to Mg 3 , Mg 2 Fe, MgFe 2 , and Fe 3 sites. The band intensities have been shown by other investigators to give the Fe:Fe+Mg ratio from transmission spectra. This study shows that the same equations can be used to obtain the ratio from reflectance spectra of unprepared samples. High‐resolution reflectance Spectroscopy of minerals may prove to be a very important tool in the laboratory, in the field using field‐portable spectrometers, from aircraft, and from satellites looking at Earth or other planetary surfaces.

Journal of Geophysical Research Solid Earth↗

Location of long-period events below Kilauea Volcano using seismic amplitudes and accurate relative relocation

We present methods for improving the location of long-period (LP) events, deep and shallow, recorded below Kilauea Volcano by the permanent seismic network. LP events might be of particular interest to understanding eruptive processes as their source mechanism is assumed to directly involve fluid transport. However, it is usually difficult or impossible to locate their source using traditional arrival time methods because of emergent wave arrivals. At Kilauea, similar LP waveform signatures suggest the existence of LP multiplets. The waveform similarity suggests spatially close sources, while catalog solutions using arrival time estimates are widely scattered beneath Kilauea's summit caldera. In order to improve estimates of absolute LP location, we use the distribution of seismic amplitudes corrected for station site effects. The decay of the amplitude as a function of hypocentral distance is used for inferring LP location. In a second stage, we use the similarity of the events to calculate their relative positions. The analysis of the entire LP seismicity recorded between January 1997 and December 1999 suggests that a very large part of the LP event population, both deep and shallow, is generated by a small number of compact sources. Deep events are systematically composed of a weak high-frequency onset followed by a low-frequency wave train. Aligning the low-frequency wave trains does not lead to aligning the onsets indicating the two parts of the signal are dissociated. This observation favors an interpretation in terms of triggering and resonance of a magmatic conduit. Instead of defining fault planes, the precise relocation of similar LP events, based on the alignment of the high-energy low-frequency wave trains, defines limited size volumes. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Paleoclimate Records: Providing context and understanding of current Arctic change

At present, the Arctic Ocean is experiencing changes in ocean surface temperature and sea ice extent that are unprecedented in the era of satellite observations, which extend from the 1980s to the present (see sections 5c,d). To provide context for current changes, scientists turn to paleoclimate records to document and study anthropogenic influence and natural decadal and multidecadal climate variability in the Arctic system. Paleoceanographic records extend limited Arctic instrumental measurements back in time and are central to improving our understanding of climate dynamics and the predictive capability of climate models. By comparing paleoceanographic records with modern observations, scientists can place the rates and magnitudes of modern Arctic change in the context of those inferred from the geological record. Over geological time, paleoceanographic reconstructions using, for instance, marine sediment cores indicate that the Arctic has experienced huge sea ice fluctuations. These fluctuations range from nearly completely ice-free to totally ice-covered conditions. The appearance of ice-rafted debris and sea ice-dependent diatoms in Arctic marine sediments indicate that the first Arctic sea ice formed approxi-mately 47 million years ago (St. John 2008; Stickley et al. 2009; Fig. SB5.1), coincident with an interval of declining atmospheric carbon dioxide (CO2) concentration, global climate cooling, and expansion of Earths cryosphere during the middle Eocene. The development of year-round (i.e., perennial) sea ice in the central Arctic Ocean, similar to conditions that exist today, is evident in sediment records as early as 1418 million years ago (Darby 2008). These records suggest that transitions in sea ice cover occur over many millennia and often vary in concert with the waxing and waning of circum-Arctic land ice sheets, ice shelves, and long-term fluctuations in ocean and atmospheric temperature and atmospheric CO2 concentrations (Stein et al. 2012; Jakobsson et al. 2014). Over shorter time scales, shallow sediment records from Arctic Ocean continental shelves allow more detailed, higher-resolution (hundreds of years resolution) reconstructions of sea ice history extending through the Holocene (11 700 years ago to present), the most recent interglacial period. A notable feature of these records is an early Holocene sea ice minimum, corresponding to a thermal maximum (warm) period from 11 000 to 5000 years ago, when the Arctic may have been warmer and had less summertime sea ice than today (Kaufman et al. 2004). However, it is not clear that the Arctic was ice-free at any point during the Holocene (Polyak et al. 2010). High-resolution paleosea ice records from the western Arctic in the Chukchi and East Siberian Seas indicate that sea ice concentrations increased through the Holocene in concert with decreasing summer solar insolation (sunlight). Sea ice extent in this region also varied in response to the volume of Pacific water delivered via the Bering Strait into the Arctic Basin (Stein et al. 2017; Polyak et al. 2016). Records from the Fram Strait (Mller et al. 2012), Laptev Sea (Hrner et al. 2016), and Canadian Arctic Archipelago (Vare et al. 2009) also indicate a similar long-term expansion of sea ice and suggest sea ice extent in these regions is modulated by the varying influx of warm Atlantic water into the Arctic Basin (e.g., Werner et al. 2013). Taken together, available records support a circum-Arctic sea ice expansion during the late Holocene. A notably high-resolution summer sea ice history (<5-year resolution) has been established for the last 1450 years using a network of terrestrial records (tree ring , lake sediment, and ice core records) located around the margins of the Arctic Ocean (Kinnard et al. 2011). Results summarized in Fig. SB5.2 indicate a pronounced decline in summer sea ice extent beginning in the 20th century, with exceptionally low ice extent recorded since the mid-1990s, consistent with the satellite record (see section 5d). While several episodes of reduced and expanded sea ice extent occur in association with climate anomalies such as the Medieval Climate Warm Period (AD 8001300) and the Little Ice Age (AD 14501850), the magnitude and pace of the modern decline in sea ice is outside of the range of natural variability and unprecedented in the 1450-year reconstruction (Kinnard et al. 2011). A radiocarbon-dated driftwood record of the Ellesmere ice shelf in the Canadian High Arctic, the oldest landfast ice in the Northern Hemisphere, also demonstrates a substantial reduction in ice extents over the 20th century (England et al. 2017). A supporting sediment record indicates that inflowing Atlantic water in Fram Strait has warmed by 2C since 1900, driving break up and melt of sea ice (Spielhagen et al. 2011). Complementary mooring and satellite observations show the Atlantification of the eastern Arctic due to enhanced inflow of warm saline water through Fram Strait (Nilsen et al. 2016) and nutrient-rich Pacific water via the Bering has increased by more than 50% (Woodgate et al. 2012), further driving sea ice melt and warming seas. Similar high-resolution proxy records from Arctic regions also indicate that the modern rate of increasing annual surface air temperatures has not been observed over at least the last 2000 years (McKay and Kaufman 2014). Scientists conclude that broad-scale sea ice variations recorded in the paleo record were dominantly driven by changes in basin-scale changes in atmospheric circulation patterns, fluctuations in air temperature, strength of incoming solar radiation, and changes in the inflow of warm water via Pacific and Atlantic inflows (Polyak et al. 2010). There is general consensus that ice-free Arctic summers are likely before the end of the 21st century (e.g., Stroeve et al. 2007; Massonnet et al. 2012), while some climate model projections suggest ice-free Arctic summers as early as 2030 (Wang and Overland 2009). Paleoclimate studies and observational time series attribute the decline in sea ice extent and thickness over the last decade to both enhanced greenhouse warming and natural climate variability. While understanding the interplay of these factors is critical for future projections of Arctic sea ice and ecosystems, most observational time series records cover only a few decades. This highlights the need for additional paleoceanographic reconstructions across multiple spatial and temporal domains to better understand the drivers and implications of present and future Arctic Ocean change.

Bulletin American Meteorological Society↗

Disentangling stationary and dynamic estuarine fish habitat to inform conservation: Species-specific responses to physical habitat and water quality in San Francisco Estuary

Estuaries represent critical aquatic habitat that connects surface water distributed between Earth’s landmasses and oceans. They are dynamic transitional ecosystems, which provide important habitat for fishes and other aquatic organisms. Effective conservation of species inhabiting estuaries requires knowledge of the habitat features that drive their abundance and distribution. We sought to elucidate how stationary (i.e., wetlands, shoals, and channels) and dynamic (i.e., salinity, temperature, turbidity, and chlorophyll concentration) habitat features interact to drive distributions of individual fish species. The Pacific coast of the conterminous United States has over 400 estuaries of various types. The largest (historical surface area) is San Francisco Estuary, California. We conducted extensive field observations of fishes in the central San Francisco Estuary among stationary habitat types (i.e., wetland, shoal, and channel) over a 19-month period encompassing substantial variability in dynamic water quality conditions. Most of the species observed, especially native species of special management interest, were associated with tidal wetland habitat. Few species exhibited associations with water quality conditions driven by seasonal (temperature) or a combination of broad- and fine-scale ecosystem processes (salinity and turbidity). Our study provides (1) an empirical demonstration of how researchers can deal with the complex and dynamic expressions of habitat in estuarine systems to address urgent natural resource problems and (2) a clear demonstration of the urgent need for habitat restoration and its likely outcome in systems such as San Francisco Estuary. Restoration of suitable tidal wetland habitat on the West Coast of the United States is likely to be an effective conservation tool to support estuarine fishes given that over 90% of historical tidal wetland habitat in San Francisco Estuary and 85% of vegetated wetland habitat along the Pacific coast of the United States has been lost due to human modification.

California↗

Earlier plant growth helps compensate for reduced carbon fixation after 13 years of warming

1. Drylands play a dominant role in global carbon cycling and are particularly vulnerable to increasing temperatures, but our understanding of how dryland ecosystems will respond to climatic change remains notably poor. Considering that the area of drylands is projected to increase 11–23% by 2100, understanding the impacts of warming on the functions and services furnished by these arid and semiarid ecosystems has numerous implications. 2.In a unique 13‐year ecosystem warming experiment in a southwestern U.S. dryland, we investigated the consequences of rising temperature on Achnatherum hymenoides , a widespread, keystone grass species on the Colorado Plateau. We tracked individual‐ and population‐level responses to identify optimal strategies that may have been masked if considering only one level of plant response. 3.We found several factors combined to affect the timing and magnitude of plant responses during the 13 th year of warming. These included large warming‐induced biomass increases for individual plants, an 8.5‐day advancement in the growing season, and strong reductions in photosynthetic rates and population cover. 4.Importantly, we observed a lack of photosynthetic acclimation and, thus, a warming‐induced down regulation of photosynthetic rates. However, these physiological responses were concurrent with warmed‐plant increases in growing season length and investment in photosynthetic surfaces, demonstrating the species' ability to balance carbon fixation limitations with warming. 5.These results, which bring together ecophysiological, phenological, reproductive, and morphological assessments of plant responses to warming, suggest that the extent of change in A. hymenoides populations will be based upon numerous adaptive responses that vary in their direction and magnitude. Plant population responses to climatic warming remain poorly resolved, particularly for Earth's drylands, and our in situ experiment assessing multiple strategies offers a novel look into a warmer world.

Utah↗

Automated snow cover detection on mountain glaciers usingspaceborne imagery and machine learning

Tracking the extent of seasonal snow on glaciers over time is critical for assessing glacier vulnerability and the response of glacierized watersheds to climate change. Existing snow cover products do not reliably distinguish seasonal snow from glacier ice and firn, preventing their use for glacier snow cover detection. Despite previous efforts to classify glacier surface facies using machine learning on local scales, currently there is no published comparison of machine learning models for classifying glacier snow cover across different satellite image products. We present an automated snow detection workflow for mountain glaciers using supervised machine-learning-based image classifiers and Landsat 8 and 9, Sentinel-2, and PlanetScope satellite imagery. We develop the image classifiers by testing numerous machine learning algorithms with training and validation data from the U.S. Geological Survey Benchmark Glacier Project glaciers. The workflow produces daily to twice monthly time series of several glacier mass balance and snowmelt indicators (snow-covered area, accumulation area ratio, and seasonal snow line) from 2013 to present. Workflow performance is assessed by comparing automatically classified images and snow lines to manual interpretations at each glacier site. The image classifiers exhibit overall accuracies of 92%–98%, K scores of 84%–96%, and F scores of 93%–98% for all image products. The median difference between automatically and manually delineated median snow line altitudes is 31m (IQR of 73to0m)across all image products. The Sentinel-2 classifier (support vector machine) produces the most accurate glacier mass balance and snowmelt indicators and distinguishes snow from ice and f irn the most reliably. Although they are less accurate, the Landsat- and PlanetScope-derived estimates greatly enhance the temporal coverage of observations. The transient accumulation area ratio produces the least noisy time series, making it the most reliable indicator for characterizing seasonal snow trends. The temporally detailed accumulation area ratio time series reveal that the timing of minimum snow cover conditions varies by up to a month between Arctic (63°N) and midlatitude (48°N) sites, underscoring the potential for bias when estimating glacier minimum snow cover conditions from a single late-summer image. Widespread application of our automated snow detection workflow has the potential to improve regional assessments of glacier mass balance, land ice representations within Earth system models, water resources, and the impacts of climate change on snow cover across broad spatial scales.

The Cryosphere↗

Reexamination of the magnitudes for the 1906 and 1922 Chilean earthquakes using Japanese tsunami amplitudes: Implications for source depth constraints

Far-field tsunami records from the Japanese tide gauge network allow the reexamination of the moment magnitudes ( M w ) for the 1906 and 1922 Chilean earthquakes, which to date rely on limited information mainly from seismological observations alone. Tide gauges along the Japanese coast provide extensive records of tsunamis triggered by six great ( M w >8) Chilean earthquakes with instrumentally determined moment magnitudes. These tsunami records are used to explore the dependence of tsunami amplitudes in Japan on the parent earthquake magnitude of Chilean origin. Using the resulting regression parameters together with tide gauge amplitudes measured in Japan we estimate apparent moment magnitudes of M w 8.0–8.2 and M w 8.5–8.6 for the 1906 central and 1922 north-central Chile earthquakes. The large discrepancy of the 1906 magnitude estimated from the tsunami observed in Japan as compared with those previously determined from seismic waves ( M s 8.4) suggests a deeper than average source with reduced tsunami excitation. A deep dislocation along the Chilean megathrust would favor uplift of the coast rather than beneath the sea, giving rise to a smaller tsunami and producing effects consistent with those observed in 1906. The 1922 magnitude inferred from far-field tsunami amplitudes appear to better explain the large extent of damage and the destructive tsunami that were locally observed following the earthquake than the lower seismic magnitudes ( M s 8.3) that were likely affected by the well-known saturation effects. Thus, a repeat of the large 1922 earthquake poses seismic and tsunami hazards in a region identified as a mature seismic gap.

Journal of Geophysical Research B: Solid Earth↗

Crustal migration of CO2-rich magmatic fluids recorded by tree-ring radiocarbon and seismicity at Mammoth Mountain, CA, USA

Unrest at Mammoth Mountain over the past several decades, manifest by seismicity, ground deformation, diffuse CO 2 emissions, and elevated 3 He/ 4 He ratios in fumarolic gases has been driven by the release of CO 2 -rich fluids from basaltic intrusions in the middle to lower crust. Recent unrest included the occurrence of three lower-crustal (32–19 km depth) seismic swarms beneath Mammoth Mountain in 2006, 2008 and 2009 that were consistently followed by peaks in the occurrence rate of shallow (≤10 km depth) earthquakes. We measured 14 C in the growth rings (1998–2012) of a tree growing in the largest (∼0.3 km 2 ) area of diffuse CO 2 emissions on Mammoth Mountain (the Horseshoe Lake tree kill; HLTK) and applied atmospheric CO 2 concentration source area modeling to confirm that the tree was a reliable integrator of magmatic CO 2 emissions over most of this area. The tree-ring 14 C record implied that magmatic CO 2 emissions from the HLTK were relatively stable from 1998 to 2009, nearly doubled from 2009 to 2011, and then declined by the 2012 growing season. The initial increase in CO 2 emissions was detected during the growing season that immediately followed the largest (February 2010) peak in the occurrence rate of shallow earthquakes. Migration of CO 2 -rich magmatic fluids may have driven observed patterns of elevated deep, then shallow seismicity, while the relationship between pore fluid pressures within a shallow (upper 3 km of crust) fluid reservoir and permeability structure of the reservoir cap rock may have controlled the temporal pattern of surface CO 2 emissions.

California↗

Competition between dissolved organic matter and freshwater plankton control methylmercury isotope fractionation during uptake and photochemical demethylation

Isotope fractionation related to photochemical reactions and planktonic uptake at the base of the food web is a major uncertainty in the biological application of mercury (Hg) stable isotopes. In freshwater systems, it is unclear how competitive interactions among methylmercury (MeHg), dissolved organic matter (DOM), and phytoplankton govern the magnitude of mass-dependent and mass-independent fractionation. This study investigated how DOM alters rates of planktonic MeHg uptake and photodegradation and corresponding Hg isotope fractionation in the presence of freshwater phytoplankton species, Raphidocelis subcapitata . Outdoor sunlight exposure experiments utilizing R. subcapitata were performed in the presence of different DOM samples using environmentally relevant ratios of MeHg-DOM thiol groups. The extent of Δ 199 Hg in phytoplankton incubations (2.99‰ St. Louis River HPOA, 1.88‰ Lake Erie HPOA) was lower compared to paired abiotic control experiments (4.29 and 2.86‰, respectively) after ∼30 h of irradiation, resulting from cell shading or other limiting factors reducing the extent of photodemethylation. Although the Δ 199 Hg/Δ 201 Hg ratio was uniform across experiments (∼1.4), Δ 199 Hg/δ 202 Hg slopes varied dramatically (from −0.96 to 15.4) across incubations with R. subcapitata and DOM. In addition, no evidence of Hg isotope fractionation was observed within R. subcapitata cells. This study provides a refined examination of Hg isotope fractionation markers for key processes occurring in the lower food web prior to bioaccumulation, critical for accurately accounting for the photochemical processing of Hg isotopes across a wide spectrum of freshwater systems.

ACS Earth and Space Chemistry↗

The accommodation of relative motion at depth on the San Andreas fault system in California

Plate motion below the seismogenic layer along the San Andreas fault system in California is generally assumed to occur by aseismic slip along a deeper extension of the fault. It is also possible that below the seismogenic layer, deformation is distributed laterally over a zone. Several observed features of the San Andreas fault in California have implications about the mode of accommodation of relative motion along the plate boundary beneath the seismogenic zone: the shallow depth of all earthquakes in California, the depth to which coseismic slip occurred during the 1906 San Francisco earthquake, the broad zone of strain accumulation, the broad heat flow anomaly, and the existence of widely separated parallel faults. The observations strongly imply that below the seismogenic zone, relative motion is distributed over a zone and occurs by inelastic flow rather than by aseismic slip on discrete fault planes. The existence of multiple faults further suggests that tractions at the base of the brittle layer are significant over time periods of years to hundreds of years.

Journal of Geophysical Research Solid Earth↗