USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Applied Mechanics”

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 37 records · Page 2Linked to original sources
USGS2010

Modeling downstream fining in sand-bed rivers. I: Formulation

In this paper a numerical modeling formulation is presented for simulation of the development of the longitudinal profile and bed sediment distribution in sand-bed rivers. The objective of the model application, which is presented in the companion paper (Wright and Parker, 2005), is to study the development of two characteristics of large, low-slope, sand-bed rivers: (1) a downstream decrease in bed slope (i.e. concave upward longitudinal profile) and (2) a downstream decrease in characteristic bed sediment diameter (e.g. the median bed surface size D 50 ). Three mechanisms that lead to an upward concave profile and downstream fining are included in the modeling formulation: (1) a delta prograding into standing water at the downstream boundary, (2) sea-level rise, and (3) tectonic subsidence. In the companion paper (Wright and Parker, 2005) the model is applied to simulate the development of the longitudinal profile and downstream fining in sand-bed rivers flowing into the ocean during the past 5000 years of relatively slow sea-level rise.

Journal of Hydraulic Research↗
USGS1991

A model for hydrostatic consolidation of Pierre shale

This paper presents closed-form solutions for consolidation of transversely isotropic porous media under hydrostatic stress. The solutions are applied to model the time variation of pore pressure, volume strain and strains parallel and normal to bedding, and to obtain coefficients of consolidation and permeability, as well as other properties, and the bulk modulus resulting from hydrostatic consolidation of Pierre shale. It is found that the coefficients consolidation and permeability decrease and the bulk moduli increase with increasing confining pressure, reflecting the closure of voids in the rock. ?? 1991.

International Journal of Rock Mechanics and Mining↗
USGS2001

Impact of competitor species composition on predicting diameter growth and survival rates of Douglas-fir trees in southwestern Oregon

Mixed conifer and hardwood stands in southwestern Oregon were studied to explore the hypothesis that competition effects on individual-tree growth and survival will differ according to the species comprising the competition measure. Likewise, it was hypothesized that competition measures should extrapolate best if crown-based surrogates are given preference over diameter-based (basal area based) surrogates. Diameter growth and probability of survival were modeled for individual Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) trees growing in pure stands. Alternative models expressing one-sided and two-sided competition as a function of either basal area or crown structure were then applied to other plots in which Douglas-fir was mixed with other conifers and (or) hardwood species. Crown-based variables outperformed basal area based variables as surrogates for one-sided competition in both diameter growth and survival probability, regardless of species composition. In contrast, two-sided competition was best represented by total basal area of competing trees. Surrogates reflecting differences in crown morphology among species relate more closely to the mechanics of competition for light and, hence, facilitate extrapolation to species combinations for which no observations are available.

Canadian Journal of Forest Research↗
USGS2005

Are all prey created equal? A review and synthesis of differential predation on prey in substandard condition

Our understanding of predator-prey interactions in fishes has been influenced largely by research assuming that the condition of the participants is normal. However, fish populations today often reside in anthropogenically altered environments and are subjected to many kinds of stressors, which may reduce their ecological performance by adversely affecting their morphology, physiology, or behaviour. One consequence is that either the predator or prey, or both, may be in a substandard condition at the time of an interaction. We reviewed the literature on predator-prey interactions in fishes where substandard prey were used as experimental groups. Although most of this research indicates that such prey are significantly more vulnerable to predation, prey condition has rarely been considered in ecological theory regarding predator-prey interactions. The causal mechanisms for increased vulnerability of substandard prey to predation include a failure to detect predators, lapses in decision-making, poor fast-start performance, inability to shoal effectively, and increased prey conspicuousness. Despite some problems associated with empirical predator-prey studies using substandard prey, their results can have theoretical and applied uses, such as in ecological modelling or justification of corrective measures to be implemented in the wild. There is a need for more corroborative field experimentation, a better understanding of the causal mechanisms behind differential predation, and increased incorporation of prey condition into the research of predator-prey modellers and theoreticians. If the concept of prey condition is considered in predator-prey interactions, our understanding of how such interactions influence the structure and dynamics of fish communities is likely to change, which should prove beneficial to aquatic ecosystems.

Journal of Fish Biology↗
USGS2006

An effective medium inversion algorithm for gas hydrate quantification and its application to laboratory and borehole measurements of gas hydrate-bearing sediments

The presence of gas hydrate in marine sediments alters their physical properties. In some circumstances, gas hydrate may cement sediment grains together and dramatically increase the seismic P- and S-wave velocities of the composite medium. Hydrate may also form a load-bearing structure within the sediment microstructure, but with different seismic wave attenuation characteristics, changing the attenuation behaviour of the composite. Here we introduce an inversion algorithm based on effective medium modelling to infer hydrate saturations from velocity and attenuation measurements on hydrate-bearing sediments. The velocity increase is modelled as extra binding developed by gas hydrate that strengthens the sediment microstructure. The attenuation increase is modelled through a difference in fluid flow properties caused by different permeabilities in the sediment and hydrate microstructures. We relate velocity and attenuation increases in hydrate-bearing sediments to their hydrate content, using an effective medium inversion algorithm based on the self-consistent approximation (SCA), differential effective medium (DEM) theory, and Biot and squirt flow mechanisms of fluid flow. The inversion algorithm is able to convert observations in compressional and shear wave velocities and attenuations to hydrate saturation in the sediment pore space. We applied our algorithm to a data set from the Mallik 2L–38 well, Mackenzie delta, Canada, and to data from laboratory measurements on gas-rich and water-saturated sand samples. Predictions using our algorithm match the borehole data and water-saturated laboratory data if the proportion of hydrate contributing to the load-bearing structure increases with hydrate saturation. The predictions match the gas-rich laboratory data if that proportion decreases with hydrate saturation. We attribute this difference to differences in hydrate formation mechanisms between the two environments.

Geophysical Journal International↗
USGS2001

An analytical solution for transient flow of Bingham viscoplastic materials in rock fractures

We present below an analytical solution to model the one-dimensional transient flow of a Bingham viscoplastic material in a fracture with parallel walls (smooth or rough) that is subjected to an applied pressure gradient. The solution models the acceleration and the deceleration of the material as the pressure gradient changes with time. Two cases are considered: A pressure gradient applied over a finite time interval and an applied pressure gradient that is constant over time. The solution is expressed in dimensionless form and can therefore be used for a wide range of Bingham viscoplastic materials. The solution is also capable of capturing the transition that takes place in a fracture between viscoplastic flow and rigid plug flow. Also, it shows the development of a rigid central layer in fractures, the extent of which depends on the fluid properties (viscosity and yield stress), the magnitude of the pressure gradient, and the fracture aperture and surface roughness. Finally, it is shown that when a pressure gradient is applied and kept constant, the solution for the fracture flow rate converges over time to a steady-state solution that can be defined as a modified cubic law. In this case, the fracture transmissivity is found to be a non-linear function of the head gradient. This solution provides a tool for a better understanding of the flow of Bingham materials in rock fractures, interfaces, and cracks. ?? 2001 Elsevier Science Ltd. All rights reserved.

International Journal of Rock Mechanics and Mining↗
USGS2008

Linking landscapes and habitat suitability scores for diadromous fish restoration in the susquehanna river basin

Dams within the Susquehanna River drainage, Pennsylvania, are potential barriers to migration of diadromous fishes, and many are under consideration for removal to facilitate fish passage. To provide useful input for prioritizing dam removal, we examined relations between landscape-scale factors and habitat suitability indices (HSIs) for native diadromous species of the Susquehanna River. We used two different methods (U.S. Fish and Wildlife Service method: Stier and Crance [1985], Ross et al. [1993a, 1993b, 1997], and Pardue [1983]; Pennsylvania State University method: Carline et al. [1994]) to calculate HSIs for several life stages of American shad Alosa sapidissima, alewives Alosa pseudoharengus, and blueback herring Alosa aestivalis and a single HSI for American eels Anguilla rostrata based on habitat variables measured at transects spaced every 5 km on six major Susquehanna River tributaries. Using geographical information systems, we calculated land use and geologic variables upstream from each transect and associated those data with HSIs calculated at each transect. We then performed canonical correlation analysis to determine how HSIs were linked to geologic and land use factors. Canonical correlation analysis identified the proportion of watershed underlain by carbonate rock as a positive correlate of HSIs for all species and life stages except American eels and juvenile blueback herring. We hypothesize that potential mechanisms linking carbonate rock to habitat suitability include increased productivity and buffering capacity. No other consistent patterns of positive or negative correlation between landscape-scale factors and HSIs were evident. This analysis will be useful for prioritizing removal of dams in the Susquehanna River drainage, because it provides a broad perspective on relationships between habitat suitability for diadromous fishes and easily measured landscape factors. This approach can be applied elsewhere to elucidate relationships between fine- and coarse-scale variables and suitability of habitat for fishes. ?? Copyright by the American Fisheries Society 2008.

North American Journal of Fisheries Management↗
USGS2018

A framework for identifying and characterising coral reef “oases” against a backdrop of degradation

Human activities have led to widespread ecological decline; however, the severity of degradation is spatially heterogeneous due to some locations resisting, escaping, or rebounding from disturbances. We developed a framework for identifying oases within coral reef regions using long‐term monitoring data. We calculated standardised estimates of coral cover ( z ‐scores) to distinguish sites that deviated positively from regional means. We also used the coefficient of variation (CV) of coral cover to quantify how oases varied temporally, and to distinguish among types of oases. We estimated “coral calcification capacity” (CCC), a measure of the coral community's ability to produce calcium carbonate structures and tested for an association between this metric and z ‐scores of coral cover. We illustrated our z ‐score approach within a modelling framework by extracting z ‐scores and CVs from simulated data based on four generalized trajectories of coral cover. We then applied the approach to time‐series data from long‐term reef monitoring programmes in four focal regions in the Pacific (the main Hawaiian Islands and Mo'orea, French Polynesia) and western Atlantic (the Florida Keys and St. John, US Virgin Islands). Among the 123 sites analysed, 38 had positive z ‐scores for median coral cover and were categorised as oases. Synthesis and applications . Our framework provides ecosystem managers with a valuable tool for conservation by identifying “oases” within degraded areas. By evaluating trajectories of change in state (e.g., coral cover) among oases, our approach may help in identifying the mechanisms responsible for spatial variability in ecosystem condition. Increased mechanistic understanding can guide whether management of a particular location should emphasise protection, mitigation or restoration. Analysis of the empirical data suggest that the majority of our coral reef oases originated by either escaping or resisting disturbances, although some sites showed a high capacity for recovery, while others were candidates for restoration. Finally, our measure of reef condition (i.e., median z ‐scores of coral cover) correlated positively with coral calcification capacity suggesting that our approach identified oases that are also exceptional for one critical component of ecological function.

Florida, Hawaii↗
USGS2016

A large refined catalog of earthquake relocations and focal mechanisms for the Island of Hawai'i and its seismotectonic implications

We present high-quality focal mechanisms based on a refined earthquake location catalog for the Island of Hawai'i, focusing on Mauna Loa and Kīlauea volcanoes. The relocation catalog is based on first-arrival times and waveform data of both compressional and shear waves for about 180,000 events on and near the Island of Hawai'i between 1986 and 2009 recorded by the seismic stations at the Hawaiian Volcano Observatory. We relocate all the earthquakes by applying ray tracing through an existing three-dimensional velocity model, similar event cluster analysis, and a differential-time relocation method. The resulting location catalog represents an expansion of previous relocation studies, covering a longer time period and consisting of more events with well-constrained absolute locations. The focal mechanisms are obtained based on the compressional-wave first-motion polarities and compressional-to-shear wave amplitude ratios by applying the HASH program to the waveform cross correlation relocated earthquakes. Overall, the good-quality (defined by the HASH parameters) focal solutions are dominated by normal faulting in our study area, especially in the active Ka'ōiki and Hīlea seismic zones. Kīlauea caldera is characterized by a mixture of approximately equal numbers of normal, strike-slip, and reverse faults, whereas its south flank has slightly fewer strike-slip events. Our relocation and focal mechanism results will be useful for mapping the seismic stress and strain fields and for understanding the seismic-volcanic-tectonic relationships within the magmatic systems.

Hawai'i↗
USGS1947

The sensitized luminescence of manganese-activated calcite

Synthetic manganese-activated calcites are shown to be practically inert to ultraviolet excitation in the range 2000-3500A, while they are luminescent under cathode-ray excitation. The incorporation of small amounts of an auxiliary impurity along with the manganese produces the strong response to ultraviolet radiation hitherto ascribed to CaCO3:Mn itself. Three such impurities have been studied: lead, thallium, and cerium. The first two induce excitation in the neighborhood of the mercury resonance line, while the cerium introduces a response principally to longer wave ultraviolet. The strong response to 2537A excitation shown by some natural calcites is likewise found to be due to the presence of lead along with the manganese, rather than to the manganese alone. The data do not warrant ascribing the longer wave-length ultraviolet-excited luminescence of all natural calcites to the action of an auxiliary impurity. The essential identity of the cathode-ray excited luminescence spectra of CaCO 3:Mn, CaCO3: (Pb+Mn), CaCO3:(Tl+Mn), and CaCO3:(Ce+Mn) with the 2537A-excited spectra of the latter three is evidence that the luminescent center in all cases is the manganese ion or the MnO6 group. It is shown that a "cascade" mechanism for the action of the auxiliary impurities, lead, thallium, and cerium, is incorrect; and that the phenomenon must be considered as a case of sensitized luminescence. Owing to the nature of cathode-ray excitation, the manganese activator can be excited by this agent even in the absence of a second impurity. For optical excitation, however, an absorption band for the ultraviolet must be established by building into the CaCO3:Mn a second impurity or "sensitizer.".

Journal of Applied Physics↗
USGS2012

Tectonic stresses in the lithosphere: constraints provided by the experimental deformation of rocks.

The strengths of rocks clearly place an upper limit on the stress that can be sustained by the upper half of the lithosphere. Laboratory data on rock rheology are generally lacking at intermediate temperatures and pressures on the important rock types expected in the lithosphere, so a definitive accounting of the strength distribution with depth in the upper lithosphere is presently unattainable. Analogies are drawn between the fragmentary strength data on silicates at intermediate temperature and the more extensive experimental data on marble and limestone, and several tentative conclusions are drawn: First, brittle processes, such as faulting and cataclasis, are expected to control rock strength at low pressures and temperatures. The strengths associated with these brittle mechanisms increase rapidly with increasing effective pressure and are relatively insensitive to temperature and strain rate. Second, the transitions between brittle and ductile processes occur at critical values of the least principal stress σ 3 . I suggest that the concept of the deformation mechanism map of Ashby (1972) be extended to brittle-ductile transitions by normalizing the applied differential stress σ by σ 3 , i.e., the transitions occur at critical values of σ/σ 3 . Third, the high temperature flow law of olivine, widely applied to problems involving flow in the asthenosphere, is inappropriate to the conditions of temperature found in the lithosphere, because different dislocation flow mechanisms dominate at low to intermediate temperatures. The fragmentary rheological data suggest the following rheological structure of the lithosphere where it is inelastically deforming: A rapid pressure-driven increase in rock strength with depth culminates with a shear strength maximum of up to 8 kbar at depths that depend on the state of stress and on the temperature distribution. A review of the mechanisms of weakening associated with water suggests that water weakening effects are probably not important in the oceanic lithosphere but are likely to be controlling in the continental crust.

Journal of Geophysical Research Solid Earth↗
USGS2013

Foreshocks during the nucleation of stick-slip instability

We report on laboratory experiments which investigate interactions between aseismic slip, stress changes, and seismicity on a critically stressed fault during the nucleation of stick-slip instability. We monitor quasi-static and dynamic changes in local shear stress and fault slip with arrays of gages deployed along a simulated strike-slip fault (2 m long and 0.4 m deep) in a saw cut sample of Sierra White granite. With 14 piezoelectric sensors, we simultaneously monitor seismic signals produced during the nucleation phase and subsequent dynamic rupture. We observe localized aseismic fault slip in an approximately meter-sized zone in the center of the fault, while the ends of the fault remain locked. Clusters of high-frequency foreshocks (M w ~ −6.5 to −5.0) can occur in this slowly slipping zone 5–50 ms prior to the initiation of dynamic rupture; their occurrence appears to be dependent on the rate at which local shear stress is applied to the fault. The meter-sized nucleation zone is generally consistent with theoretical estimates, but source radii of the foreshocks (2 to 70 mm) are 1 to 2 orders of magnitude smaller than the theoretical minimum length scale over which earthquake nucleation can occur. We propose that frictional stability and the transition between seismic and aseismic slip are modulated by local stressing rate and that fault sections, which would typically slip aseismically, may radiate seismic waves if they are rapidly stressed. Fault behavior of this type may provide physical insight into the mechanics of foreshocks, tremor, repeating earthquake sequences, and a minimum earthquake source dimension.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

State of stress and modern deformation of the northern Basin and Range Province

Constraints on the current stress regime of the actively extending northern Basin and Range province are provided by deformation data (focal mechanisms and fault slip studies), hydraulic fracturing in situ stress measurements, borehole elongation (“breakouts”) analyses, and alignment of young volcanic vents. The integrated data indicate significant variations both in principal stress orientations and magnitudes. An approximately E-W least principal stress direction appears to characterize both the eastern and western margins of the Basin and Range province, whereas in the active interior parts of the province extension occurs in response to a least principal stress oriented NW to N60°W. The contrast in stress orientations between the province boundaries and in the interior suggests that along the margins the least principal stress direction may be locally controlled by the generally northerly trending profound lithospheric discontinuities associated with these margins. Active deformation along the southeastern and western province margins is characterized by a combination of strike-slip and normal faulting. Focal mechanisms along northeastern province margin (Wasatch front) and in central Nevada indicate a combination of normal and oblique-normal faulting. Temporal, regional, and depth-dependent variations in the relative magnitudes of the vertical and maximum horizontal stresses can explain much of the observed variations in deformation styles. However, some depth variation in faulting style inferred from focal mechanisms may be apparent and simply a function of the attitude of fault planes being reactivated. Evidence for significant temporal variation (or multiple cycles of variation) in relative stress magnitude comes from the Sierran front-Basin and Range boundary region where recent earthquakes are predominantly strike slip, whereas the profound relative vertical relief across the Sierra frontal fault zone in the last 9–10 m.y. implies a normal faulting stress regime. Using the best data on stress orientation, relative stress magnitudes are constrained from slip vectors of major earthquakes and young fault displacements. Analysis of well-constrained slip vectors in the Owens Valley, California, area indicate that large temporal variations in the magnitude of the approximately N-S oriented maximum horizontal stress are required to explain dominantly dip-slip and strike-slip offsets on subparallel faults. Similar faulting relations are observed throughout much of the boundary zone between the Basin and Range-Sierra Nevada (including the Walker Lane belt). Along the eastern province margin in the Wasatch front area in Utah, available data suggest that the maximum and minimum horizontal stresses may be approximately equal at depths of <4–5 km. Earthquake focal mechanisms in this area suggest more variability in relative magnitude of the two horizontal stresses with depth. Furthermore, superimposed sets of young fault striae along a segment of the Wasatch fault also indicate temporal variations of relative stress magnitudes. Sources of regional and temporal variations in the stress field may be linked to variable shear tractions applied to the base of the brittle crust related to intrusion, thermally induced flow, and the influence of the San Andreas plate boundary. Although difficult to date accurately, the fault slip data suggest that the temporal variations in relative magnitudes stress may occur on the time scale of both a single major earthquake cycle (1000–5000 years) and multiple earthquake cycles (10,000+ years).

Journal of Geophysical Research Solid Earth↗
USGS1981

A transient laboratory method for determining the hydraulic properties of 'tight' rocks-II. Application

In Part I a general analytical solution for the transient pulse test was presented. Part II presents a graphical method for analyzing data from a test to obtain the hydraulic properties of the sample. The general solution depends on both hydraulic conductivity and specific storage and, in theory, analysis of the data can provide values for both of these hydraulic properties. However, in practice, one of two limiting cases may apply in which case it is possible to calculate only hydraulic conductivity or the product of hydraulic conductivity times specific storage. In this paper we examine the conditions when both hydraulic parameters can be calculated. The analyses of data from two tests are presented. In Appendix I the general solution presented in Part I is compared with an earlier analysis, in which compressive storage in the sample is assumed negligible, and the error in calculated hydraulic conductivity due to this simplifying assumption is examined.

International Journal of Rock Mechanics and Mining↗
USGS2013

Conditions favouring Bromus tectorum dominance of endangered sagebrush steppe ecosystems

1. Ecosystem invasibility is determined by combinations of environmental variables, invader attributes, disturbance regimes, competitive abilities of resident species and evolutionary history between residents and disturbance regimes. Understanding the relative importance of each factor is critical to limiting future invasions and restoring ecosystems. 2. We investigated factors potentially controlling Bromus tectorum invasions into Artemisia tridentata ssp. wyomingensis communities across 75 sites in the Great Basin. We measured soil texture, cattle grazing intensity, gaps among perennial plants and plant cover including B. tectorum, biological soil crusts (BSCs) and bare soil. Using a priori knowledge, we developed a multivariate hypothesis of the susceptibility of Artemisia ecosystems to B. tectorum invasion and used the model to assess the relative importance of the factors driving the magnitude of such invasions. 3. Model results imply that bunchgrass community structure, abundance and composition, along with BSC cover, play important roles in controlling B. tectorum dominance. Evidence suggests abundant bunchgrasses limit invasions by limiting the size and connectivity of gaps between vegetation, and BSCs appear to limit invasions within gaps. Results also suggest that cattle grazing reduces invasion resistance by decreasing bunchgrass abundance, shifting bunchgrass composition, and thereby increasing connectivity of gaps between perennial plants while trampling further reduces resistance by reducing BSC. 4. Synthesis and applications. Grazing exacerbates Bromus tectorum dominance in one of North America's most endangered ecosystems by adversely impacting key mechanisms mediating resistance to invasion. If the goal is to conserve and restore resistance of these systems, managers should consider maintaining or restoring: (i) high bunchgrass cover and structure characterized by spatially dispersed bunchgrasses and small gaps between them; (ii) a diverse assemblage of bunchgrass species to maximize competitive interactions with B. tectorum in time and space; and (iii) biological soil crusts to limit B. tectorum establishment. Passive restoration by reducing cumulative cattle grazing may be one of the most effective means of achieving these three goals.

Journal of Applied Ecology↗
USGS1998

Three-dimensional seismic structure and moment tensors of non-double-couple earthquakes at the Hengill-Grensdalur volcanic complex, Iceland

The volcanic and geothermal areas of Iceland are rich sources of non-double-couple (non-DC) earthquakes. A state-of-the-art digital seismometer network deployed at the Hengill–Grensdalur volcanic complex in 1991 recorded 4000 small earthquakes. We used the best recorded of these to determine 3-D V P and V P /V S structure tomographically and accurate earthquake moment tensors. The V P field is dominated by high seismic wave speed bodies interpreted as solidified intrusions. A widespread negative (−4 per cent) V P /V S anomaly in the upper 4km correlates with the geothermal field, but is too strong to be caused solely by the effect of temperature upon liquid water or the presence of vapour, and requires in addition mineralogical or lithological differences between the geothermal reservoir and its surroundings. These may be caused by geothermal alteration. Well-constrained moment tensors were obtained for 70 of the best-recorded events by applying linear programming methods to P - and S -wave polarities and amplitude ratios. About 25 per cent of the mechanisms are, within observational error, consistent with DC mechanisms consistent with shear faulting. The other 75 per cent have significantly non-DC mechanisms. Many have substantial explosive components, one has a substantial implosive component, and the deviatoric component of many is strongly non-DC. Many of the non-DC mechanisms are consistent, within observational error, with simultaneous tensile and shear faulting. However, the mechanisms occupy a continuum in source-type parameter space and probably at least one additional source process is occurring. This may be fluid flow into newly formed cracks, causing partial compensation of the volumetric component. Studying non-shear earthquakes such as these has great potential for improving our understanding of geothermal processes and earthquake source processes in general.

Geophysical Journal International↗
USGS1990

Theory and application of an approximate model of saltwater upconing in aquifers

Motion and mixing of salt water and fresh water are vitally important for water-resource development throughout the world. An approximate model of saltwater upconing in aquifers is developed, which results in three non-linear coupled equations for the freshwater zone, the saltwater zone, and the transition zone. The description of the transition zone uses the concept of a boundary layer. This model invokes some assumptions to give a reasonably tractable model, considerably better than the sharp interface approximation but considerably simpler than a fully three-dimensional model with variable density. We assume the validity of the Dupuit-Forchheimer approximation of horizontal flow in each layer. Vertical hydrodynamic dispersion into the base of the transition zone is assumed and concentration of the saltwater zone is assumed constant. Solute in the transition zone is assumed to be moved by advection only. Velocity and concentration are allowed to vary vertically in the transition zone by using shape functions. Several numerical techniques can be used to solve the model equations, and simple analytical solutions can be useful in validating the numerical solution procedures. We find that the model equations can be solved with adequate accuracy using the procedures presented. The approximate model is applied to the Smoky Hill River valley in central Kansas. This model can reproduce earlier sharp interface results as well as evaluate the importance of hydrodynamic dispersion for feeding salt water to the river. We use a wide range of dispersivity values and find that unstable upconing always occurs. Therefore, in this case, hydrodynamic dispersion is not the only mechanism feeding salt water to the river. Calculations imply that unstable upconing and hydrodynamic dispersion could be equally important in transporting salt water. For example, if groundwater flux to the Smoky Hill River were only about 40% of its expected value, stable upconing could exist where hydrodynamic dispersion into a transition zone is the primary mechanism for moving salt water to the river. The current model could be useful in situations involving dense saltwater layers.

Journal of Hydrology↗
USGS2013

Seismotectonic framework of the 2010 February 27 M w 8.8 Maule, Chile earthquake sequence

After the 2010 M w 8.8 Maule earthquake, an international collaboration involving teams and instruments from Chile, the US, the UK, France and Germany established the International Maule Aftershock Deployment temporary network over the source region of the event to facilitate detailed, open-access studies of the aftershock sequence. Using data from the first 9-months of this deployment, we have analyzed the detailed spatial distribution of over 2500 well-recorded aftershocks. All earthquakes have been relocated using a hypocentral decomposition algorithm to study the details of and uncertainties in both their relative and absolute locations. We have computed regional moment tensor solutions for the largest of these events to produce a catalogue of 465 mechanisms, and have used all of these data to study the spatial distribution of the aftershock sequence with respect to the Chilean megathrust. We refine models of co-seismic slip distribution of the Maule earthquake, and show how small changes in fault geometries assumed in teleseismic finite fault modelling significantly improve fits to regional GPS data, implying that the accuracy of rapid teleseismic fault models can be substantially improved by consideration of existing fault geometry model databases. We interpret all of these data in an integrated seismotectonic framework for the Maule earthquake rupture and its aftershock sequence, and discuss the relationships between co-seismic rupture and aftershock distributions. While the majority of aftershocks are interplate thrust events located away from regions of maximum co-seismic slip, interesting clusters of aftershocks are identified in the lower plate at both ends of the main shock rupture, implying internal deformation of the slab in response to large slip on the plate boundary interface. We also perform Coulomb stress transfer calculations to compare aftershock locations and mechanisms to static stress changes following the Maule rupture. Without the incorporation of uncertainties in earthquake locations, just 55 per cent of aftershock nodal planes align with faults promoted towards failure by co-seismic slip. When epicentral uncertainties are considered (on the order of just ±2–3 km), 90 per cent of aftershocks are consistent with occurring along faults demonstrating positive stress transfer. These results imply large sensitivities of Coulomb stress transfer calculations to uncertainties in both earthquake locations and models of slip distributions, particularly when applied to aftershocks close to a heterogeneous fault rupture; such uncertainties should therefore be considered in similar studies used to argue for or against models of static stress triggering.

Maule↗