SEARCH · Search USGS

Results for “Journal of Applied Geophysics”

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 235 records · Page 13Linked to original sources

USGS1997

Permeability and porosity of the Illinois UPH 3 drillhole granite and a comparison with other deep drillhole rocks

Permeability, porosity, and volumetric strain measurements were conducted on granite cores obtained at depths of 0.7 to 1.6 km from the Illinois UPH 3 drillhole at effective confining pressures from 5 to 100 MPa. Initial permeabilities were in the range of 10-17 to 10-19 m2 and dropped rapidly with applied pressure to values between 10-20 and 10-24 m2 at 100 MPa, typical of other deep granite core samples. These values are several decades lower than equivalent weathered surface granites at comparable effective confining pressures, where weathering products in cracks and pores inhibit crack closure with applied pressure. Permeabilities of the Illinois cores were inversely related to sample depth, suggesting that stress relief and thermal microfractures induced during core retrieval dominated the fluid flow. Thus these samples provide an upper bound on in situ matrix permeability values. A comparison of core permeability from UPH 3 and other deep drillholes shows that stress relief damage can often dominate laboratory permeability measurements. We conclude that it may be difficult to make meaningful estimates of in situ permeability based on either borehole samples (possible damage during retrieval) or surface-derived analogs (altered by weathering). Volumetric strain determined from porosity measurements was compared with differential strain analysis (DSA) data reported by other investigators on samples from the same depths in the drillhole. Our strain measurements (0.002 to 0.005 at 100 MPa) were nearly twice as large as the DSA values, probably because of the crack-enhancing effects of fluids present in our samples that are absent in the dry DSA cores, as well as other time-dependent deformation effects. This difference in observed strain magnitudes between the two measurement methods may be an important consideration if strain and/or porosity data from deep core samples are used in models of stress, fluid circulation, and excess fluid pressure generation in the midcrust.

Journal of Geophysical Research B: Solid Earth↗

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

The effect of stress changes on time-dependent earthquake probabilities for the central Wasatch Fault Zone, Utah, USA

Static and quasi-static Coulomb stress changes produced by large earthquakes can modify the probability of occurrence of subsequent events on neighboring faults. This approach is based on physical (Coulomb stress changes) and statistical (probability calculations) models, which are influenced by the quality and quantity of data available in the study region. Here, we focus on the Wasatch Fault Zone (WFZ), a well-studied active normal fault system having abundant geologic and paleoseismological data. Paleoseismological trench investigations of the WFZ indicate that at least 24 large, surface-faulting earthquakes have ruptured the fault’s five central, 35–59-km long segments since ~7 ka. Our goal is to determine if the stress changes due to the youngest paleoevents have significantly modified the present-day probability of occurrence of large earthquakes on each of the segments. For each segment, we modeled the cumulative (coseismic + postseismic) Coulomb stress changes (∆CFScum) due to earthquakes younger than the most recent event on the segment in question and applied the resulting values to the time-dependent probability calculations. Results from the Coulomb stress modeling suggest that the Brigham City, Salt Lake City, and Provo segments have accumulated ∆CFScum larger than 10 bars, whereas the Weber segment has experienced a stress decrease of 5 bars, in the scenario of recent rupture of the Great Salt Lake fault to the west. Probability calculations predict high probability of occurrence for the Brigham City and Salt Lake City segments, due to their long elapsed times (>1-2 ka) when compared to the Weber, Provo, and Nephi segments (< 1 ka). The range of calculated coefficients of variation (CV) has a large influence on the final probabilities, mostly in the case of the Brigham City segment. Finally, when the Coulomb stress and the probability models are combined, our results indicate that the ∆CFScum resulting from earthquakes postdating the youngest events on each of the five segments significantly affects the probability calculations for three of the segments: Brigham City, Salt Lake City, and Provo. The probability of occurrence of a large earthquake in the next 50 years on these three segments may therefore be underestimated if a time-independent approach, or a time-dependent approach that does not consider ∆CFS, is adopted.

Utah↗

USGS2016

Three-dimensional surface deformation derived from airborne interferometric UAVSAR: Application to the Slumgullion Landslide

In order to provide surface geodetic measurements with “landslide-wide” spatial coverage, we develop and validate a method for the characterization of 3-D surface deformation using the unique capabilities of the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) airborne repeat-pass radar interferometry system. We apply our method at the well-studied Slumgullion Landslide, which is 3.9 km long and moves persistently at rates up to ∼2 cm/day. A comparison with concurrent GPS measurements validates this method and shows that it provides reliable and accurate 3-D surface deformation measurements. The UAVSAR-derived vector velocity field measurements accurately capture the sharp boundaries defining previously identified kinematic units and geomorphic domains within the landslide. We acquired data across the landslide during spring and summer and identify that the landslide moves more slowly during summer except at its head, presumably in response to spatiotemporal variations in snowmelt infiltration. In order to constrain the mechanics controlling landslide motion from surface velocity measurements, we present an inversion framework for the extraction of slide thickness and basal geometry from dense 3-D surface velocity fields. We find that the average depth of the Slumgullion Landslide is 7.5 m, several meters less than previous depth estimates. We show that by considering a viscoplastic rheology, we can derive tighter theoretical bounds on the rheological parameter relating mean horizontal flow rate to surface velocity. Using inclinometer data for slow-moving, clay-rich landslides across the globe, we find a consistent value for the rheological parameter of 0.85 ± 0.08.

Colorado↗

USGS2026

Subsurface structure across the Tacoma Basin, Washington State, using trans-dimensional Bayesian inversion of fundamental mode spatial autocorrelation data

Spatial autocorrelation (SPAC), the azimuthal average of the normalized cross-correlation between equidistant station pairs deployed in a 2-D array, is widely used to image the subsurface structure. However, the rigorous estimate of subsurface structure and its uncertainties as a function of depth using SPAC data is challenging due to the nonlinear relation between the SPAC data and Earth structure as well as the trade-off between depth and velocity. Additionally, data noise is strongly correlated due to data processing (e.g. filtering, stacking from multiple time segments and azimuthal averaging). Most studies do not account for the correlated noise and fix the ratio of compressional-wave velocity ( V P ) to shear-wave velocity ( V s ) (i.e. V P / V s ratio) and the number of layers, both of which are typically unknown. To address these challenges, we develop a hierarchical trans-dimensional Bayesian inversion of fundamental mode of SPAC data that properly accounts for the correlated data noise, samples the V P / V s ratio and relaxes the number of layers (i.e. model parametrization) to be unknown in the inversion. We further examine the limitation of using only fundamental modes in the inversion. Our synthetic experiments show that the inversion recovers an incorrect model unless we sample the correlated noise and V P / V s ratio in the inversion. The inversion is then applied to SPAC data acquired at 19 sites across the Tacoma basin in Washington State to characterize the V s and the time-averaged V s over 30-m depth ( V s 30 ). Our results show that the V s 30 varies from ∼200 to 800 m s −1 . The V s 30 within the basin is higher in the middle and lower on the east and west sides. We find that these V s 30 values vary with geologic unit. The uncertainties for V s 30 are within 20 m s −1 in average except for the most eastern site TB28. Additionally, the uncertainties are greater for deeper depths beneath most of the sites as the sensitivity decreases as a function of depth. The Vs structure as a function of depth is also complex beneath some sites, possibly because the SPAC curves are affected by higher order Rayleigh modes that are not considered in the inversion. To better constrain the deeper V s structure, V s 30 and/or other average measures of V s over depth, additional constraints from complementary data, such as ellipticity or geologic data are needed. Moreover, our synthetic experiments show that higher order modes can have significant effect in the inversion results, particularly when there is a low-velocity layer.

Washington↗

USGS2003

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Compositional variations on the Moon: Recalibration of Galileo solid‐state imaging data for the Orientale region and farside

Updated radiometric calibration and systematic processing procedures for Galileo solid‐state imaging (SSI) data from the first (1990) Earth‐Moon encounter are presented. These procedures were applied to a whole‐disk imaging sequence of the Moon centered near Mare Orientale, called Lunmap 14 (L14). Processing of L14 data included radiometric calibration, subpixel coregistration, scattered light removal, geometric control and reprojection, photometric normalization, and calibration to Earth‐based spectra. Coregistration and scattered‐light removal procedures are improvements over the initial calibration of the SSI mosaics. The effects of scattered light correction are best seen using a whole‐Moon view such as L14; resolution of the debated amounts of light scattering from within or outside the camera field of view is not necessary. Scattered light removal particularly affects the 1‐μm spectral region and has implications for interpretation of mafic mineral signatures in mare deposits of the lunar limb and farside. Recalibrated spectra indicate that mare ponds of the limb show moderately deep 1‐μm absorptions, and thus mafic mineral contents, comparable to those of other nearside basalts. Mafic mineral contents of Schiller‐Shickard cryptomaria are higher than previously thought and are similar to some low‐Ti nearside basalts. Many of the recalibrated spectra from South Pole/Aitken are similar to those of Schiller‐Schickard cryptomaria, suggesting that many of these soils represent a mixed mare/highland lithology. The hypothesis that there is an olivine enrichment in southern South Pole/Aitken basin is not supported strongly by spectra shown here.

Journal of Geophysical Research: Planets↗

USGS2015

Point spread functions for earthquake source imaging: An interpretation based on seismic interferometry

Recently, various methods have been proposed and applied for earthquake source imaging, and theoretical relationships among the methods have been studied. In this study, we make a follow-up theoretical study to better understand the meanings of earthquake source imaging. For imaging problems, the point spread function (PSF) is used to describe the degree of blurring and degradation in an obtained image of a target object as a response of an imaging system. In this study, we formulate PSFs for earthquake source imaging. By calculating the PSFs, we find that waveform source inversion methods remove the effect of the PSF and are free from artifacts. However, the other source imaging methods are affected by the PSF and suffer from the effect of blurring and degradation due to the restricted distribution of receivers. Consequently, careful treatment of the effect is necessary when using the source imaging methods other than waveform inversions. Moreover, the PSF for source imaging is found to have a link with seismic interferometry with the help of the source-receiver reciprocity of Green’s functions. In particular, the PSF can be related to Green’s function for cases in which receivers are distributed so as to completely surround the sources. Furthermore, the PSF acts as a low-pass filter. Given these considerations, the PSF is quite useful for understanding the physical meaning of earthquake source imaging.

Geophysical Journal International↗

USGS2025

Identification of representative earthquakes for probabilistic tsunami hazard analysis (PTHA) using earthquake rupture forecasts and machine learning

As probabilistic tsunami hazard analysis (PTHA) focuses more on assessments for localized, populous regions, techniques are needed to identify a subsample of representative earthquake ruptures to make the computational requirements for producing high-resolution hazard maps tractable. Moreover, the greatest epistemic uncertainty in seismic PTHA is related to source characterization, which is often poorly defined and subjective. We address these two salient issues by applying streamlined earthquake rupture forecasts (ERFs), based on combinatorial optimization methods, to an unsupervised machine learning workflow for identifying representative ruptures. ERFs determine the optimal distribution of a millennia-scale sample of earthquakes by inverting the observed slip rate on major faults. We use two previously developed combinatorial optimization ERFs, integer programming and greedy sequential, to produce the optimal location of ruptures with seismic moments sampled from a regional Gutenberg–Richter magnitude–frequency distribution. These ruptures in turn are used to calculate peak nearshore tsunami amplitude, using computationally efficient tsunami Green's functions. An unsupervised machine learning workflow is then used to identify a small subsample of the earthquakes input to ERFs for onshore PTHA analysis. We eliminate epistemic uncertainty related to source distribution under traditional PTHA analysis; in its place, a quantifiable, less subjective and generally smaller uncertainty related to the input to ERFs is included. The Nankai subduction zone is used as a test case, where previous ERFs have been conducted. Results indicate that the locations of representative earthquakes are sensitive to choice of magnitude–area relation and to whether a minimum cumulative stress objective is imposed on the fault. In general, incorporating ERFs into PTHA provide a physically self-consistent method to incorporate fault slip information in determining representative earthquakes for onshore PTHA, eliminating a major source of epistemic uncertainty.

Nankai subduction zone↗

USGS2002

Topographic stress perturbations in southern Davis Mountains, west Texas 1. Polarity reversal of principal stresses

We have applied a previously developed analytical stress model to interpret subsurface stress conditions inferred from acoustic televiewer logs obtained in two municipal water wells located in a valley in the southern Davis Mountains near Alpine, Texas. The appearance of stress-induced breakouts with orientations that shift by 90° at two different depths in one of the wells is explained by results from exact solutions for the effects of valleys on gravity and tectonically induced subsurface stresses. The theoretical results demonstrate that above a reference depth termed the hinge point, a location that is dependent on Poisson's ratio, valley shape, and magnitude of the maximum horizontal tectonic stress normal to the long axis of the valley, horizontal stresses parallel to the valley axis are greater than those normal to it. At depths below this hinge point the situation reverses and horizontal stresses normal to the valley axis are greater than those parallel to it. Application of the theoretical model at Alpine is accommodated by the fact that nearby earthquake focal mechanisms establish an extensional stress regime with the regional maximum horizontal principal stress aligned perpendicular to the valley axis. We conclude that the localized stress field associated with a valley setting can be highly variable and that breakouts need to be examined in this context when estimating the orientations and magnitudes of regional principal stresses.

Texas↗

USGS2012

Modeling thermal dynamics of active layer soils and near-surface permafrost using a fully coupled water and heat transport model

Thawing and freezing processes are key components in permafrost dynamics, and these processes play an important role in regulating the hydrological and carbon cycles in the northern high latitudes. In the present study, we apply a well-developed soil thermal model that fully couples heat and water transport, to simulate the thawing and freezing processes at daily time steps across multiple sites that vary with vegetation cover, disturbance history, and climate. The model performance was evaluated by comparing modeled and measured soil temperatures at different depths. We use the model to explore the influence of climate, fire disturbance, and topography (north- and south-facing slopes) on soil thermal dynamics. Modeled soil temperatures agree well with measured values for both boreal forest and tundra ecosystems at the site level. Combustion of organic-soil horizons during wildfire alters the surface energy balance and increases the downward heat flux through the soil profile, resulting in the warming and thawing of near-surface permafrost. A projection of 21st century permafrost dynamics indicates that as the climate warms, active layer thickness will likely increase to more than 3 meters in the boreal forest site and deeper than one meter in the tundra site. Results from this coupled heat-water modeling approach represent faster thaw rates than previously simulated in other studies. We conclude that the discussed soil thermal model is able to well simulate the permafrost dynamics and could be used as a tool to analyze the influence of climate change and wildfire disturbance on permafrost thawing.

Journal of Geophysical Research D: Atmospheres↗

USGS1996

Crustal structure of the Colorado Plateau, Arizona: Application of new long-offset seismic data analysis techniques

The Colorado Plateau is a large crustal block in the southwestern United States that has been raised intact nearly 2 km above sea level since Cretaceous marine sediments were deposited on its surface. Controversy exists concerning the thickness of the plateau crust and the source of its buoyancy. Interpretations of seismic data collected on the plateau vary as to whether the crust is closer to 40 or 50 km thick. A thick crust could support the observed topography of the Colorado Plateau isostatically, while a thinner crust would indicate the presence of an underlying low-density mantle. This paper reports results on long-offset seismic data collected during the 1989 segment of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that extended from the Transition Zone into the Colorado Plateau in northwest Arizona. We apply two new methods to analyze long-offset data that employ finite difference travel time calculations: (1) a first-arrival time inverter to find upper crustal velocity structure and (2) a forward-modeling technique that allows the direct use of the inverted upper crustal solution in modeling secondary reflected arrivals. We find that the crustal thickness increases from 30 km beneath the metamorphic core complexes in the southern Basin and Range province to about 42 km beneath the northern Transition Zone and southern Colorado Plateau margin. We observe some crustal thinning (to ∼37 km thick) and slightly higher lower crustal velocities farther inboard; beneath the Kaibab uplift on the north rim of the Grand Canyon the crust thickens to a maximum of 48 km. We observe a nonuniform crustal thickness beneath the Colorado Plateau that varies by ∼15% and corresponds approximately to variations in topography with the thickest crust underlying the highest elevations. Crustal compositions (as inferred from seismic velocities) appear to be the same beneath the Colorado Plateau as those in the Basin and Range province to the southwest, implying that the plateau crust represents an unextended version of the Basin and Range. Some of the variability in crustal structure appears to correspond to preserved lithospheric discontinuities that date back to the Proterozoic Era.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Wave-current interaction in Willapa Bay

This paper describes the importance of wave-current interaction in an inlet-estuary system. The three-dimensional, fully coupled, Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system was applied in Willapa Bay (Washington State) from 22 to 29 October 1998 that included a large storm event. To represent the interaction between waves and currents, the vortex-force method was used. Model results were compared with water elevations, currents, and wave measurements obtained by the U.S. Army Corp of Engineers. In general, a good agreement between field data and computed results was achieved, although some discrepancies were also observed in regard to wave peak directions in the most upstream station. Several numerical experiments that considered different forcing terms were run in order to identify the effects of each wind, tide, and wave-current interaction process. Comparison of the horizontal momentum balances results identified that wave-breaking-induced acceleration is one of the leading terms in the inlet area. The enhancement of the apparent bed roughness caused by waves also affected the values and distribution of the bottom shear stress. The pressure gradient showed significant changes with respect to the pure tidal case. During storm conditions the momentum balance in the inlet shares the characteristics of tidal-dominated and wave-dominated surf zone environments. The changes in the momentum balance caused by waves were manifested both in water level and current variations. The most relevant effect on hydrodynamics was a wave-induced setup in the inner part of the estuary.

Washington↗

USGS2003

When and where the aftershock activity was depressed: Contrasting decay patterns of the proximate large earthquakes in southern California

Seismic quiescence has attracted attention as a possible precursor to a large earthquake. However, sensitive detection of quiescence requires accurate modeling of normal aftershock activity. We apply the epidemic-type aftershock sequence (ETAS) model that is a natural extension of the modified Omori formula for aftershock decay, allowing further clusters (secondary aftershocks) within an aftershock sequence. The Hector Mine aftershock activity has been normal, relative to the decay predicted by the ETAS model during the 14 months of available data. In contrast, although the aftershock sequence of the 1992 Landers earthquake (M = 7.3), including the 1992 Big Bear earthquake (M = 6.4) and its aftershocks, fits very well to the ETAS up until about 6 months after the main shock, the activity showed clear lowering relative to the modeled rate (relative quiescence) and lasted nearly 7 years, leading up to the Hector Mine earthquake (M = 7.1) in 1999. Specifically, the relative quiescence occurred only in the shallow aftershock activity, down to depths of 5-6 km. The sequence of deeper events showed clear, normal aftershock activity well fitted to the ETAS throughout the whole period. We argue several physical explanations for these results. Among them, we strongly suspect aseismic slips within the Hector Mine rupture source that could inhibit the crustal relaxation process within "shadow zones" of the Coulomb's failure stress change. Furthermore, the aftershock activity of the 1992 Joshua Tree earthquake (M = 6.1) sharply lowered in the same day of the main shock, which can be explained by a similar scenario.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

The strength and rheology of methane clathrate hydrate

Methane clathrate hydrate (structure I) is found to be very strong, based on laboratory triaxial deformation experiments we have carried out on samples of synthetic, high-purity, polycrystalline material. Samples were deformed in compressional creep tests (i.e., constant applied stress, ??), at conditions of confining pressure P = 50 and 100 MPa, strain rate 4.5 ?? 10-8 ??? ?? ??? 4.3 ?? 10-4 s-1, temperature 260 ??? T ??? 287 K, and internal methane pressure 10 ??? PCH4 ??? 15 MPa. At steady state, typically reached in a few percent strain, methane hydrate exhibited strength that was far higher than expected on the basis of published work. In terms of the standard high-temperature creep law, ?? = A??ne-(E*+PV*)/RT the rheology is described by the constants A = 108.55 MPa-n s-1, n = 2.2, E* = 90,000 J mol-1, and V* = 19 cm3 mol-1. For comparison at temperatures just below the ice point, methane hydrate at a given strain rate is over 20 times stronger than ice, and the contrast increases at lower temperatures. The possible occurrence of syntectonic dissociation of methane hydrate to methane plus free water in these experiments suggests that the high strength measured here may be only a lower bound. On Earth, high strength in hydrate-bearing formations implies higher energy release upon decomposition and subsequent failure. In the outer solar system, if Titan has a 100-km-thick near-surface layer of high-strength, low-thermal conductivity methane hydrate as has been suggested, its interior is likely to be considerably warmer than previously expected.

Journal of Geophysical Research B: Solid Earth↗

USGS2010

Locations and magnitudes of historical earthquakes in the Sierra of Ecuador (1587–1996)

The whole territory of Ecuador is exposed to seismic hazard. Great earthquakes can occur in the subduction zone (e.g. Esmeraldas, 1906, M w 8.8), whereas lower magnitude but shallower and potentially more destructive earthquakes can occur in the highlands. This study focuses on the historical crustal earthquakes of the Andean Cordillera. Several large cities are located in the Interandean Valley, among them Quito, the capital (∼2.5 millions inhabitants). A total population of ∼6 millions inhabitants currently live in the highlands, raising the seismic risk. At present, precise instrumental data for the Ecuadorian territory is not available for periods earlier than 1990 (beginning date of the revised instrumental Ecuadorian seismic catalogue); therefore historical data are of utmost importance for assessing seismic hazard. In this study, the Bakun & Wentworth method is applied in order to determine magnitudes, locations, and associated uncertainties for historical earthquakes of the Sierra over the period 1587–1976. An intensity-magnitude equation is derived from the four most reliable instrumental earthquakes ( M w between 5.3 and 7.1). Intensity data available per historical earthquake vary between 10 (Quito, 1587, Intensity ≥VI) and 117 (Riobamba, 1797, Intensity ≥III). The bootstrap resampling technique is coupled to the B&W method for deriving geographical confidence contours for the intensity centre depending on the data set of each earthquake, as well as confidence intervals for the magnitude. The extension of the area delineating the intensity centre location at the 67 per cent confidence level (±1σ) depends on the amount of intensity data, on their internal coherence, on the number of intensity degrees available, and on their spatial distribution. Special attention is dedicated to the few earthquakes described by intensities reaching IX, X and XI degrees. Twenty-five events are studied, and nineteen new epicentral locations are obtained, yielding equivalent moment magnitudes between 5.0 and 7.6. Large earthquakes seem to be related to strike slip faults between the North Andean Block and stable South America to the east, while moderate earthquakes ( M w ≤ 6) seem to be associated with to thrust faults located on the western internal slopes of the Interandean Valley.

Geophysical Journal International↗

USGS2015

Limited role for thermal erosion by turbulent lava in proximal Athabasca Valles, Mars

The Athabasca Valles flood lava is among the most recent (<50 Ma) and best preserved effusive lava flows on Mars and was probably emplaced turbulently. The Williams et al. (2005) model of thermal erosion by lava has been applied to what we term “proximal Athabasca,” the 75 km long upstream portion of Athabasca Valles. For emplacement volumes of 5000 and 7500 km 3 and average flow thicknesses of 20 and 30 m, the duration of the eruption varies between ~11 and ~37 days. The erosion of the lava flow substrate is investigated for three eruption temperatures (1270°C, 1260°C, and 1250°C), and volatile contents equivalent to 0–65 vol % bubbles. The largest erosion depths of ~3.8–7.5 m are at the lava source, for 20 m thick and bubble-free flows that erupted at their liquidus temperature (1270°C). A substrate containing 25 vol % ice leads to maximum erosion. A lava temperature 20°C below liquidus reduces erosion depths by a factor of ~2.2. If flow viscosity increases with increasing bubble content in the lava, the presence of 30–50 vol % bubbles leads to erosion depths lower than those relative to bubble-free lava by a factor of ~2.4. The presence of 25 vol % ice in the substrate increases erosion depths by a factor of 1.3. Nevertheless, modeled erosion depths, consistent with the emplacement volume and flow duration constraints, are far less than the depth of the channel (~35–100 m). We conclude that thermal erosion does not appear to have had a major role in excavating Athabasca Valles.

Journal of Geophysical Research E: Planets↗

USGS2012

Landscape controls on total and methyl Hg in the Upper Hudson River basin, New York, USA

Approaches are needed to better predict spatial variation in riverine Hg concentrations across heterogeneous landscapes that include mountains, wetlands, and open waters. We applied multivariate linear regression to determine the landscape factors and chemical variables that best account for the spatial variation of total Hg (THg) and methyl Hg (MeHg) concentrations in 27 sub-basins across the 493 km 2 upper Hudson River basin in the Adirondack Mountains of New York. THg concentrations varied by sixfold, and those of MeHg by 40-fold in synoptic samples collected at low-to-moderate flow, during spring and summer of 2006 and 2008. Bivariate linear regression relations of THg and MeHg concentrations with either percent wetland area or DOC concentrations were significant but could account for only about 1/3 of the variation in these Hg forms in summer. In contrast, multivariate linear regression relations that included metrics of (1) hydrogeomorphology, (2) riparian/wetland area, and (3) open water, explained about 66% to >90% of spatial variation in each Hg form in spring and summer samples. These metrics reflect the influence of basin morphometry and riparian soils on Hg source and transport, and the role of open water as a Hg sink. Multivariate models based solely on these landscape metrics generally accounted for as much or more of the variation in Hg concentrations than models based on chemical and physical metrics, and show great promise for identifying waters with expected high Hg concentrations in the Adirondack region and similar glaciated riverine ecosystems.

New York↗