Search USGS⌕ Search

SEARCH · Search USGS

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

Molecular analysis of fungal communities and laccase genes in decomposing litter reveals differences among forest types but no impact of nitrogen deposition

The fungal community of the forest floor was examined as the cause of previously reported increases in soil organic matter due to experimental N deposition in ecosystems producing predominantly high-lignin litter, and the opposite response in ecosystems producing low-lignin litter. The mechanism proposed to explain this phenomenon was that white-rot basidiomycetes are more important in the degradation of high-lignin litter than of low-lignin litter, and that their activity is suppressed by N deposition. We found that forest floor mass in the low-lignin sugar-maple dominated system decreased in October due to experimental N deposition, whereas forest floor mass of high-lignin oak-dominated ecosystems was unaffected by N deposition. Increased relative abundance of basidiomycetes in high-lignin forest floor was confirmed by denaturing gradient gel electrophoresis (DGGE) and sequencing. Abundance of basidiomycete laccase genes, encoding an enzyme used by white-rot basidiomycetes in the degradation of lignin, was 5-10 times greater in high-lignin forest floor than in low-lignin forest floor. While the differences between the fungal communities in different ecosystems were consistent with the proposed mechanism, no significant effects of N deposition were detected on DGGE profiles, laccase gene abundance, laccase length heterogeneity profiles, or phenol oxidase activity. Our observations indicate that the previously detected accumulation of soil organic matter in the high-lignin system may be driven by effects of N deposition on organisms in the mineral soil, rather than on organisms residing in the forest floor. However, studies of in situ gene expression and temporal and spatial variability within forest floor communities will be necessary to further relate the ecosystem dynamics of organic carbon to microbial communities and atmospheric N deposition. ?? 2007 The Authors; Journal compilation ?? 2007 Society for Applied Microbiology and Blackwell Publishing Ltd.

Environmental Microbiology↗

The failure of earthquake failure models

In this study I show that simple heuristic models and numerical calculations suggest that an entire class of commonly invoked models of earthquake failure processes cannot explain triggering of seismicity by transient or "dynamic" stress changes, such as stress changes associated with passing seismic waves. The models of this class have the common feature that the physical property characterizing failure increases at an accelerating rate when a fault is loaded (stressed) at a constant rate. Examples include models that invoke rate state friction or subcritical crack growth, in which the properties characterizing failure are slip or crack length, respectively. Failure occurs when the rate at which these grow accelerates to values exceeding some critical threshold. These accelerating failure models do not predict the finite durations of dynamically triggered earthquake sequences (e.g., at aftershock or remote distances). Some of the failure models belonging to this class have been used to explain static stress triggering of aftershocks. This may imply that the physical processes underlying dynamic triggering differs or that currently applied models of static triggering require modification. If the former is the case, we might appeal to physical mechanisms relying on oscillatory deformations such as compaction of saturated fault gouge leading to pore pressure increase, or cyclic fatigue. However, if dynamic and static triggering mechanisms differ, one still needs to ask why static triggering models that neglect these dynamic mechanisms appear to explain many observations. If the static and dynamic triggering mechanisms are the same, perhaps assumptions about accelerating failure and/or that triggering advances the failure times of a population of inevitable earthquakes are incorrect.

Journal of Geophysical Research B: Solid Earth↗

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↗

Monitoring and modeling ice-rock avalanches from ice-capped volcanoes: A case study of frequent large avalanches on Iliamna Volcano, Alaska

Iliamna is an andesitic stratovolcano of the Aleutian arc with regular gas and steam emissions and mantled by several large glaciers. Iliamna Volcano exhibits an unusual combination of frequent and large ice-rock avalanches in the order of 1 ?? 106??m3 to 3 ?? 107??m3 with recent return periods of 2-4??years. We have reconstructed an avalanche event record for the past 45??years that indicates Iliamna avalanches occur at higher frequency at a given magnitude than other mass failures in volcanic and alpine environments. Iliamna Volcano is thus an ideal site to study such mass failures and its relation to volcanic activity. In this study, we present different methods that fit into a concept of (1) long-term monitoring, (2) early warning, and (3) event documentation and analysis of ice-rock avalanches on ice-capped active volcanoes. Long-term monitoring methods include seismic signal analysis, and space-and airborne observations. Landsat and ASTER satellite data was used to study the extent of hydrothermally altered rocks and surface thermal anomalies at the summit region of Iliamna. Subpixel heat source calculation for the summit regions where avalanches initiate yielded temperatures of 307 to 613??K assuming heat source areas of 1000 to 25??m2, respectively, indicating strong convective heat flux processes. Such heat flow causes ice melting conditions and is thus likely to reduce the strength at the base of the glacier. We furthermore demonstrate typical seismic records of Iliamna avalanches with rarely observed precursory signals up to two hours prior to failure, and show how such signals could be used for a multi-stage avalanche warning system in the future. For event analysis and documentation, space- and airborne observations and seismic records in combination with SRTM and ASTER derived terrain data allowed us to reconstruct avalanche dynamics and to identify remarkably similar failure and propagation mechanisms of Iliamna avalanches for the past 45??years. Simple avalanche flow modeling was able to reasonably replicate Iliamna avalanches and can thus be applied for hazard assessments. Hazards at Iliamna Volcano are low due to its remote location; however, we emphasize the transfer potential of the methods presented here to other ice-capped volcanoes with much higher hazards such as those in the Cascades or the Andes. ?? 2007 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

The excitation and characteristic frequency of the long-period volcanic event: An approach based on an inhomogeneous autoregressive model of a linear dynamic system

We present a method to quantify the source excitation function and characteristic frequencies of long-period volcanic events. The method is based on an inhomogeneous autoregressive (AR) model of a linear dynamic system, in which the excitation is assumed to be a time-localized function applied at the beginning of the event. The tail of an exponentially decaying harmonic waveform is used to determine the characteristic complex frequencies of the event by the Sompi method. The excitation function is then derived by operating an AR filter constructed from the characteristic frequencies to the entire seismogram of the event, including the inhomogeneous part of the signal. We apply this method to three long-period events at Kusatsu-Shirane Volcano, central Japan, whose waveforms display simple decaying monochromatic oscillations except for the beginning of the events. We recover time-localized excitation functions lasting roughly 1 s at the start of each event and find that the estimated functions are very similar to each other at all the stations of the seismic network for each event. The phases of the characteristic oscillations referred to the estimated excitation function fall within a narrow range for almost all the stations. These results strongly suggest that the excitation and mode of oscillation are both dominated by volumetric change components. Each excitation function starts with a pronounced dilatation consistent with a sudden deflation of the volumetric source which may be interpreted in terms of a choked-flow transport mechanism. The frequency and Q of the characteristic oscillation both display a temporal evolution from event to event. Assuming a crack filled with bubbly water as seismic source for these events, we apply the Van Wijngaarden-Papanicolaou model to estimate the acoustic properties of the bubbly liquid and find that the observed changes in the frequencies and Q are consistently explained by a temporal change in the radii of the bubbles characterizing the bubbly water in the crack.

Journal of Geophysical Research B: Solid Earth↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Numerical modeling of injection, stress and permeability enhancement during shear stimulation at the Desert Peak Enhanced Geothermal System

Creation of an Enhanced Geothermal System relies on stimulation of fracture permeability through self-propping shear failure that creates a complex fracture network with high surface area for efficient heat transfer. In 2010, shear stimulation was carried out in well 27-15 at Desert Peak geothermal field, Nevada, by injecting cold water at pressure less than the minimum principal stress. An order-of-magnitude improvement in well injectivity was recorded. Here, we describe a numerical model that accounts for injection-induced stress changes and permeability enhancement during this stimulation. In a two-part study, we use the coupled thermo-hydrological-mechanical simulator FEHM to: (i) construct a wellbore model for non-steady bottom-hole temperature and pressure conditions during the injection, and (ii) apply these pressures and temperatures as a source term in a numerical model of the stimulation. In this model, a Mohr-Coulomb failure criterion and empirical fracture permeability is developed to describe permeability evolution of the fractured rock. The numerical model is calibrated using laboratory measurements of material properties on representative core samples and wellhead records of injection pressure and mass flow during the shear stimulation. The model captures both the absence of stimulation at low wellhead pressure (WHP ≤1.7 and ≤2.4 MPa) as well as the timing and magnitude of injectivity rise at medium WHP (3.1 MPa). Results indicate that thermoelastic effects near the wellbore and the associated non-local stresses further from the well combine to propagate a failure front away from the injection well. Elevated WHP promotes failure, increases the injection rate, and cools the wellbore; however, as the overpressure drops off with distance, thermal and non-local stresses play an ongoing role in promoting shear failure at increasing distance from the well.

Nevada↗

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↗

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↗

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↗

Foraging optimally for home ranges

Economic models predict behavior of animals based on the presumption that natural selection has shaped behaviors important to an animal's fitness to maximize benefits over costs. Economic analyses have shown that territories of animals are structured by trade-offs between benefits gained from resources and costs of defending them. Intuitively, home ranges should be similarly structured, but trade-offs are difficult to assess because there are no costs of defense, thus economic models of home-range behavior are rare. We present economic models that predict how home ranges can be efficient with respect to spatially distributed resources, discounted for travel costs, under 2 strategies of optimization, resource maximization and area minimization. We show how constraints such as competitors can influence structure of homes ranges through resource depression, ultimately structuring density of animals within a population and their distribution on a landscape. We present simulations based on these models to show how they can be generally predictive of home-range behavior and the mechanisms that structure the spatial distribution of animals. We also show how contiguous home ranges estimated statistically from location data can be misleading for animals that optimize home ranges on landscapes with patchily distributed resources. We conclude with a summary of how we applied our models to nonterritorial black bears (Ursus americanus) living in the mountains of North Carolina, where we found their home ranges were best predicted by an area-minimization strategy constrained by intraspecific competition within a social hierarchy. Economic models can provide strong inference about home-range behavior and the resources that structure home ranges by offering falsifiable, a priori hypotheses that can be tested with field observations.

Journal of Mammalogy↗