SEARCH · Search USGS

Results for “Journal of Geophysical Research - Atmospheres”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3Linked to original sources

USGS2001

Volcano monitoring using the Global Positioning System: Filtering strategies

Permanent Global Positioning System (GPS) networks are routinely used for producing improved orbits and monitoring secular tectonic deformation. For these applications, data are transferred to an analysis center each day and routinely processed in 24-hour segments. To use GPS for monitoring volcanic events, which may last only a few hours, real-time or near real-time data processing and subdaily position estimates are valuable. Strategies have been researched for obtaining station coordinates every 15 min using a Kalman filter; these strategies have been tested on data collected by a GPS network on Kilauea Volcano. Data from this network are tracked continuously, recorded every 30 s, and telemetered hourly to the Hawaiian Volcano Observatory. A white noise model is heavily impacted by data outages and poor satellite geometry, but a properly constrained random walk model fits the data well. Using a borehole tiltmeter at Kilauea's summit as ground-truth, solutions using different random walk constraints were compared. This study indicates that signals on the order of 5 mm/h are resolvable using a random walk standard deviation of 0.45 cm/√h. Values lower than this suppress small signals, and values greater than this have significantly higher noise at periods of 1–6 hours.

Hawai'i↗

USGS2018

Archie’s saturation exponent for natural gas hydrate in coarse-grained reservoirs

Accurately quantifying the amount of naturally occurring gas hydrate in marine and permafrost environments is important for assessing its resource potential and understanding the role of gas hydrate in the global carbon cycle. Electrical resistivity well logs are often used to calculate gas hydrate saturations, S h , using Archie's equation. Archie's equation, in turn, relies on an empirical saturation parameter, n . Though n = 1.9 has been measured for ice‐bearing sands and is widely used within the hydrate community, it is highly questionable if this n value is appropriate for hydrate‐bearing sands. In this work, we calibrate n for hydrate‐bearing sands from the Canadian permafrost gas hydrate research well, Mallik 5L‐38, by establishing an independent downhole S h profile based on compressional‐wave velocity log data. Using the independently determined S h profile and colocated electrical resistivity and bulk density logs, Archie's saturation equation is solved for n, and uncertainty is tracked throughout the iterative process. In addition to the Mallik 5L‐38 well, we also apply this method to two marine, coarse‐grained reservoirs from the northern Gulf of Mexico Gas Hydrate Joint Industry Project: Walker Ridge 313‐H and Green Canyon 955‐H. All locations yield similar results, each suggesting n ≈ 2.5 ± 0.5. Thus, for the coarse‐grained hydrate bearing ( S h > 0.4) of greatest interest as potential energy resources, we suggest that n = 2.5 ± 0.5 should be applied in Archie's equation for either marine or permafrost gas hydrate settings if independent estimates of n are not available.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Creep of water ices at planetary conditions: A compilation

Many constitutive laws for the flow of ice have been published since the advent of the Voyager explorations of the outer solar system. Conflicting data have occasionally come from different laboratories, and refinement of experimental techniques has led to the publication of laws that supersede earlier ones. In addition, there are unpublished data from ongoing research that also amend the constitutive laws. Here we compile the most current laboratory-derived flow laws for water ice phases I, II, III, V, and VI, and ice I mixtures with hard particulates. The rheology of interest is mainly that of steady state, and the conditions reviewed are the pressures and temperatures applicable to the surfaces and interiors of icy moons of the outer solar system. Advances in grain-size-dependent creep in ices I and II as well as in phase transformations and metastability under differential stress are also included in this compilation. At laboratory strain rates the several ice polymorphs are rheologically distinct in terms of their stress, temperature, and pressure dependencies but, with the exception of ice III, have fairly similar strengths. Hard particulates strengthen ice I significantly only at high particulate volume fractions. Ice III has the potential for significantly affecting mantle dynamics because it is much weaker than the other polymorphs and its region of stability, which may extend metastably well into what is nominally the ice II field, is located near likely geotherms of large icy moons.

Journal of Geophysical Research E: Planets↗

USGS2021

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗

USGS2006

Shrinking ponds in subarctic Alaska based on 1950-2002 remotely sensed images

Over the past 50 years, Alaska has experienced a warming climate with longer growing seasons, increased potential evapotranspiration, and permafrost warming. Research from the Seward Peninsula and Kenai Peninsula has demonstrated a substantial landscape-level trend in the reduction of surface water and number of closed-basin ponds. We investigated whether this drying trend occurred at nine other regions throughout Alaska. One study region was from the Arctic Coastal Plain where depp permafrost occurs continuously across the landscape. The other eight study regions were from the boreal forest regions where discontinuous permafrost occurs. Mean annual precipitation across the study regions ranged from 100 to over 700 min yr-1. We used remotely sensed imagery from the 1950s to 2002 to inventory over 10,000 closed-basin ponds from at least three periods from this time span. We found a reduction in the area and number of shallow, closed-basin ponds for all boreal regions. In contrast, the Arctic Coastal Plain region had negligible change in the area of closed-basin ponds. Since the 1950s, surface water area of closed-basin ponds included in this analysis decreased by 31 to 4 percent, and the total number of closed-basin ponds surveyed within each study region decreased from 54 to 5 percent. There was a significant increasing trend in annual mean temperature and potential evapotranspiration since the 1950s for all study regions. There was no significant trend in annual precipitation during the same period. The regional trend of shrinking ponds may be due to increased drainage as permafrost warms, or increased evapotranspiration during a warmer and extended growing season. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗

USGS2003

Topographically driven groundwater flow and the San Andreas heat flow paradox revisited

Evidence for a weak San Andreas Fault includes (1) borehole heat flow measurements that show no evidence for a frictionally generated heat flow anomaly and (2) the inferred orientation of ??1 nearly perpendicular to the fault trace. Interpretations of the stress orientation data remain controversial, at least in close proximity to the fault, leading some researchers to hypothesize that the San Andreas Fault is, in fact, strong and that its thermal signature may be removed or redistributed by topographically driven groundwater flow in areas of rugged topography, such as typify the San Andreas Fault system. To evaluate this scenario, we use a steady state, two-dimensional model of coupled heat and fluid flow within cross sections oriented perpendicular to the fault and to the primary regional topography. Our results show that existing heat flow data near Parkfield, California, do not readily discriminate between the expected thermal signature of a strong fault and that of a weak fault. In contrast, for a wide range of groundwater flow scenarios in the Mojave Desert, models that include frictional heat generation along a strong fault are inconsistent with existing heat flow data, suggesting that the San Andreas Fault at this location is indeed weak. In both areas, comparison of modeling results and heat flow data suggest that advective redistribution of heat is minimal. The robust results for the Mojave region demonstrate that topographically driven groundwater flow, at least in two dimensions, is inadequate to obscure the frictionally generated heat flow anomaly from a strong fault. However, our results do not preclude the possibility of transient advective heat transport associated with earthquakes.

Journal of Geophysical Research B: Solid Earth↗

USGS2020

Plastic faulting in ice

Plastic faulting is a brittle‐like failure phenomenon exhibited by water ice and several other rock types under confinement. It is suspected to be the mechanism of deep earthquakes and extreme cases of shear localization in shallow rocks. Unlike ordinary Coulombic failure, plastic faulting is characterized by a pressure‐independent failure strength and fault plane oriented 45° to maximum principal stress. To research the question of how the instability initiates, we conducted over 50 constant‐displacement‐rate experiments on polycrystalline ice (phases Ih and II) near the brittle‐to‐ductile (B‐D) transition, at confining pressures P = 0–300 MPa, applied strain rates = 5 × 10 −5 – 7 × 10 −3 s −1 , temperatures T = 105–233 K, and mean grain sizes d = 0.25–1.18 mm. We find that (1) the width of the B‐D transition in variable space is vanishingly narrow, to the point of appearing as a crossover, (2) a plastic fault plane, once formed, is not a zone of subsequent weakness, (3) distributed ice I→II phase transformation in small amounts (<1 vol%) shows no causal relationship to subsequent failure, and (4) plastic faulting also occurs in ice II. We hypothesize that the elusive nucleating “trigger” parallels that of metals and ceramics undergoing severe plastic deformation, wherein transient local structural rearrangement occurs, in turn causing material strength to drop to a level sufficiently low, in a volume sufficiently large, that adiabatic instability is nucleated. Our results do not require and often are inconsistent with phase transformation. Plastic faulting may therefore be available to all solids undergoing severe deformation, and its appearance in so few is simply the result of insufficiently extreme conditions.

Journal of Geophysical Research- Solid Earth↗

USGS2009

Responses of stream nitrate and dissolved organic carbon loadings to hydrological forcing and climate change in an upland forest of the northeast USA

[1] In coming decades, higher annual temperatures, increased growing season length, and increased dormant season precipitation are expected across the northeastern United States in response to anthropogenic forcing of global climate. We synthesized long-term stream hydrochemical data from the Sleepers River Research Watershed in Vermont, United States, to explore the relationship of catchment wetness to stream nitrate and DOC loadings. We modeled changes in growing season length and precipitation patterns to simulate future climate scenarios and to assess how stream nutrient loadings respond to climate change. Model results for the 2070–2099 time period suggest that stream nutrient loadings during both the dormant and growing seasons will respond to climate change. During a warmer climate, growing season stream fluxes (runoff +20%, nitrate +57%, and DOC +58%) increase as more precipitation (+28%) and quick flow (+39%) occur during a longer growing season (+43 days). During the dormant season, stream water and nutrient loadings decrease. Net annual stream runoff (+8%) and DOC loading (+9%) increases are commensurate with the magnitude of the average increase of net annual precipitation (+7%). Net annual stream water and DOC loadings are primarily affected by increased dormant season precipitation. In contrast, decreased annual loading of stream nitrate (−2%) reflects a larger effect of growing season controls on stream nitrate and the effects of lengthened growing seasons in a warmer climate. Our findings suggest that leaching of nitrate and DOC from catchment soils will be affected by anthropogenic climate forcing, thereby affecting the timing and magnitude of annual stream loadings in the northeastern United States.

Vermont↗

USGS2021

Cataloging tectonic tremor energy radiation in the Cascadia subduction zone

For the past ∼12 years the Pacific Northwest Seismic Network has been automatically detecting and locating tectonic tremor across the Cascadia subduction zone, resulting in a catalog of more than 500,000 tremor epicenters to date, which has served as a valuable resource for tremor and slip research. This manuscript presents an updated methodology for routine tremor detection in Cascadia and a new catalog of over 180,000 tremor epicenters including amplitudes detected along the subduction zone margin from 2017 to 2021. The events are detected via cross-correlation of continuous vertical envelope data of 128 stations from northern California to northern Vancouver Island. The modified approach results in less scatter and a 55% increase in detected epicenters than previously observed, as well as a newly identified tremor source offset updip from the main tremor and slip region at the southern edge of the subduction zone. Radiated seismic energy in the 1.5–5 Hz band is used to assign epicenters an energy magnitude ( M eL ), which is calibrated to the M L of local earthquakes. Southern Cascadia is most active, but the highest tremor energy rates occur in northern Cascadia. Tremor in central Cascadia is systematically weaker and less frequent. Individual epicenter magnitudes range from ∼0.5–2 and spatiotemporally cluster into 1,060 swarms with cumulative M eL ranging from ∼0.8 to 3.7. The swarms reflect underlying slow slip events and occur with an earthquake-like energy distribution with a b value ∼1. Tremor epicenters, however, follow a tapered Gutenberg-Richter distribution with high b values, suggesting individual tremor bursts and their constituent low-frequency earthquakes are fault-dimension limited.

California, Oregon, Washington↗

USGS2007

Mars global digital dune database and initial science results

A new Mars Global Digital Dune Database (MGD 3 ) constructed using Thermal Emission Imaging System (THEMIS) infrared (IR) images provides a comprehensive and quantitative view of the geographic distribution of moderate‐ to large‐size dune fields (area >1 km 2 ) that will help researchers to understand global climatic and sedimentary processes that have shaped the surface of Mars. MGD 3 extends from 65°N to 65°S latitude and includes ∼550 dune fields, covering ∼70,000 km 2 , with an estimated total volume of ∼3,600 km 3 . This area, when combined with polar dune estimates, suggests moderate‐ to large‐size dune field coverage on Mars may total ∼800,000 km 2 , ∼6 times less than the total areal estimate of ∼5,000,000 km 2 for terrestrial dunes. Where availability and quality of THEMIS visible (VIS) or Mars Orbiter Camera narrow‐angle (MOC NA) images allow, we classify dunes and include dune slipface measurements, which are derived from gross dune morphology and represent the prevailing wind direction at the last time of significant dune modification. For dunes located within craters, the azimuth from crater centroid to dune field centroid (referred to as dune centroid azimuth) is calculated and can provide an accurate method for tracking dune migration within smooth‐floored craters. These indicators of wind direction are compared to output from a general circulation model (GCM). Dune centroid azimuth values generally correlate to regional wind patterns. Slipface orientations are less well correlated, suggesting that local topographic effects may play a larger role in dune orientation than regional winds.

Journal of Geophysical Research E: Planets↗

USGS2020

Geological and thermal control of the hydrothermal system in northern Yellowstone Lake: Inferences from high resolution magnetic surveys

A multiscale magnetic survey of the northern basin of Yellowstone Lake was undertaken in 2016 as part of the Hydrothermal Dynamics of Yellowstone Lake Project (HD‐YLAKE)—a broad research effort to characterize the cause‐and‐effect relationships between geologic and environmental processes and hydrothermal activity on the lake floor. The magnetic survey includes lake surface, regional aeromagnetic, and near‐bottom autonomous underwater vehicle (AUV) data. The study reveals a strong contrast between the northeastern lake basin, characterized by a regional magnetic low punctuated by stronger local magnetic lows, many of which host hydrothermal vent activity, and the northwestern lake basin with higher‐amplitude magnetic anomalies and no obvious hydrothermal activity or punctuated magnetic lows. The boundary between these two regions is marked by a steep gradient in heat flow and magnetic values, likely reflecting a significant structure within the currently active ~20‐km‐long Eagle Bay‐Lake Hotel fault zone that may be related to the ~2.08‐Ma Huckleberry Ridge caldera rim. Modeling suggests that the broad northeastern magnetic low reflects both a shallower Curie isotherm and widespread hydrothermal activity that has demagnetized the rock. Along the western lake shoreline are sinuous‐shaped, high‐amplitude magnetic anomaly highs, interpreted as lava flow fronts of upper units of the West Thumb rhyolite. The AUV magnetic survey shows decreased magnetization at the periphery of the active Deep Hole hydrothermal vent. We postulate that lower magnetization in the outer zone results from enhanced hydrothermal alteration of rhyolite by hydrothermal condensates while the vapor‐dominated center of the vent is less altered.

Wyoming↗

USGS2012

Structure of the lower crust beneath the Carolina Trough, U.S. Atlantic continental margin

Data from three large-offset seismic profiles provide information on the crustal structure beneath the Carolina trough. The profiles, obtained by the U.S. Geological Survey, the Naval Oceanographic Research Development Agency, and the Scripps Institution of Oceanography in 1985, were oriented parallel to the trough and were located (1) seaward of the East Coast Magnetic Anomaly (ECMA), which is generally thought to represent the boundary between oceanic and continental crust; (2) along the axis of the trough between the ECMA and the hinge zone, which is thought to reflect the landward limit of highly stretched and altered transitional crust; and (3) along the Carolina platform landward of the basement hinge zone on crust thought to have been thinned only slightly during rifting. These data constrain the velocity structure of the lower crust and provide evidence for a thick lens of high-velocity (>7.1 km/s) lower crustal material that extends beneath the Carolina trough and the adjacent ocean basin. This lens reaches a maximum thickness of about 13 km beneath the deepest part of the trough, thins to about 5 km seaward of the ECMA, and is either very thin or absent landward of the hinge zone. It is interpreted to represent material that was underplated beneath and/or intruded into the crust during the late stage of continental rifting and that led to an anomalously thick plutonic layer during the early seafloor spreading phase. These data thus support the recent conclusions of White et al. (1987b) and Mutter et al. (1988) that the initiation of seafloor spreading is attended in many, if not most, cases by the generation of an anomalously large volume of melt.

Carolina Trough↗

USGS2022

Quantifying interdependencies in geyser eruptions at the Upper Geyser Basin, Yellowstone National Park

The Upper Geyser Basin at Yellowstone National Park (Wyoming, USA) harbors the greatest concentration of geysers worldwide. Research suggests that individual geysers are not isolated but rather are hydraulically connected in the subsurface with other geysers and thermal springs. To quantify such connections, we combined techniques from machine learning, causal inference, and dynamical systems to characterize the collective eruptive behavior of a set of 10 geysers over 18 months (April 2007 – September 2008) focusing on geyser-geyser interactions. Model predictions were up to 15 times more accurate when we sought to predict a geyser's eruption time series based on outflow channel temperatures from the network than based on its own time series alone, suggesting the existence of a complex interconnected subsurface groundwater system. On average, cone-type geysers had larger impacts on other geysers than did fountain-type geysers. Similarly, cone-type geysers were on average more insulated from other geysers. However, substantial unexplained variation remained after considering the cone versus fountain dichotomy. Distance between geysers also affected interactions: nearby geysers had stronger effects on focal geysers than did geysers located farther away. Collectively, results support the hypothesis of geyser interdependence at timescales of 5 min–10 days. Our analyses highlight the existence of quantifiable geyser-to-geyser interactions that can be resolved through pairwise and system-level analyses. These findings emphasize the subsurface interconnectedness of thermal features, provide information relevant to visitor experiences in Yellowstone National Park, and suggest strategies for exploring patterns of interdependence that may exist among other episodic geological phenomena.

Wyoming↗

USGS2014

Temporal variability of carbon and nutrient burial, sediment accretion, and mass accumulation over the past century in a carbonate platform mangrove forest of the Florida Everglades.

The objective of this research was to measure temporal variability in accretion and mass sedimentation rates (including organic carbon (OC), total nitrogen (TN), and total phosphorous (TP)) from the past century in a mangrove forest on the Shark River in Everglades National Park, USA. The 210 Pb Constant Rate of Supply model was applied to six soil cores to calculate annual rates over the most recent 10, 50, and 100 year time spans. Our results show that rates integrated over longer timeframes are lower than those for shorter, recent periods of observation. Additionally, the substantial spatial variability between cores over the 10 year period is diminished over the 100 year record, raising two important implications. First, a multiple-decade assessment of soil accretion and OC burial provides a more conservative estimate and is likely to be most relevant for forecasting these rates relative to long-term processes of sea level rise and climate change mitigation. Second, a small number of sampling locations are better able to account for spatial variability over the longer periods than for the shorter periods. The site average 100 year OC burial rate, 123 ± 19 (standard deviation) g m -2 yr -1 , is low compared with global mangrove values. High TN and TP burial rates in recent decades may lead to increased soil carbon remineralization, contributing to the low carbon burial rates. Finally, the strong correlation between OC burial and accretion across this site signals the substantial contribution of OC to soil building in addition to the ecosystem service of CO 2 sequestration.

Florida↗

USGS2019

Identifying the molecular signatures of agricultural expansion in Amazonian headwater streams

Agricultural impacts on aquatic ecosystems are well-studied, however, most research has focused on temperate regions, whereas the forefront of agricultural expansion is currently in the tropics. At the vanguard of this growth is the boundary between the Amazon and Cerrado biomes in Brazil, driven primarily by expansion of soybean and corn croplands. Here we examine the impacts of cropland expansion on receiving lowland Amazon Basin headwater streams in terms of dissolved organic carbon (DOC) concentration and dissolved organic matter (DOM) composition via ultra-high resolution mass spectrometry. Streams draining croplands had lower DOC concentrations and DOM molecular signatures enriched in N- and S-containing formula in comparison to forested streams. Cropland streams were also enriched in aliphatic, peptide-like and highly unsaturated and phenolic (low O/C) compound categories in comparison to forest streams (enriched in polyphenolics, condensed aromatics and highly unsaturated and phenolic (high O/C) compound categories) indicative of the shifting of sources from organic-rich surface soils and litter layers to autochthonous and more microbial biomass. Distinct molecular assemblages were strongly correlated with cropland and forest catchments, highlighting headwater streams as sentinels for detecting change. On investigation of unique molecular formulae present in only cropland sites, four cropland markers provided the ability to track agricultural impacts in the region. Overall, these patterns indicate reduced organic matter inputs in croplands and greater microbial degradation at these sites leading to declining DOC concentrations, and DOM of more microbial character in receiving streams that is more biolabile, with clear ramifications for downstream ecology and biogeochemical cycles.

Amazon River basin, Tanguro Ranch, Xinga River↗

USGS2012

Characterizing mercury concentrations and fluxes in a Coastal Plain watershed: Insights from dynamic modeling and data

Mercury (Hg) is one of the leading water quality concerns in surface waters of the United States. Although watershed-scale Hg cycling research has increased in the past two decades, advances in modeling watershed Hg processes in diverse physiographic regions, spatial scales, and land cover types are needed. The goal of this study was to assess Hg cycling in a Coastal Plain system using concentrations and fluxes estimated by multiple watershed-scale models with distinct mathematical frameworks reflecting different system dynamics. We simulated total mercury (Hg T , the sum of filtered and particulate forms) concentrations and fluxes from a Coastal Plain watershed (McTier Creek) using three watershed Hg models and an empirical load model. Model output was compared with observed in-stream Hg T . We found that shallow subsurface flow is a potentially important transport mechanism of particulate Hg T during periods when connectivity between the uplands and surface waters is maximized. Other processes (e.g., stream bank erosion, sediment re-suspension) may increase particulate Hg T in the water column. Simulations and data suggest that variable source area (VSA) flow and lack of rainfall interactions with surface soil horizons result in increased dissolved Hg T concentrations unrelated to DOC mobilization following precipitation events. Although flushing of DOC-Hg T complexes from surface soils can also occur during this period, DOC-complexed Hg T becomes more important during base flow conditions. TOPLOAD simulations highlight saturated subsurface flow as a primary driver of daily Hg T loadings, but shallow subsurface flow is important for Hg T loads during high-flow events. Results suggest limited seasonal trends in Hg T dynamics.

South Carolina↗

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↗

USGS2013

Inferring fault rheology from low-frequency earthquakes on the San Andreas

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor (NVT) on the San Andreas fault in central California show strong sensitivity to shear stress induced by the daily tidal cycle. LFEs occur at all levels of the tidal shear stress and are in phase with the very small, ~400 Pa, stress amplitude. To quantitatively explain the correlation, we use a model from the existing literature that assumes the LFE sources are small, persistent regions that repeatedly fail during shear of a much larger scale, otherwise aseismically creeping fault zone. The LFE source patches see tectonic loading, creep of the surrounding fault which may be modulated by the tidal stress, and direct tidal loading. If the patches are small relative to the surrounding creeping fault then the stressing is dominated by fault creep, and if patch failure occurs at a threshold stress, then the resulting seismicity rate is proportional to the fault creep rate or fault zone strain rate. Using the seismicity rate as a proxy for strain rate and the tidal shear stress, we fit the data with possible fault rheologies that produce creep in laboratory experiments at temperatures of 400 to 600°C appropriate for the LFE source depth. The rheological properties of rock-forming minerals for dislocation creep and dislocation glide are not consistent with the observed fault creep because strong correlation between small stress perturbations and strain rate requires perturbation on the order of the ambient stress. The observed tidal modulation restricts ambient stress to be at most a few kilopascal, much lower than rock strength. A purely rate dependent friction is consistent with the observations only if the product of the friction rate dependence and effective normal stress is ~ 0.5 kPa. Extrapolating the friction rate strengthening dependence of phyllosilicates (talc) to depth would require the effective normal stress to be ~50 kPa, implying pore pressure is lithostatic. If the LFE source is on the order of tens of meters, as required by the model, rate-weakening friction rate dependence (e.g., olivine) at 400 to 600°C requires that the minimum effective pressure at the LFE source is ~ 2.5 MPa.

California↗