SEARCH · Geology Reports

Results for “The Energy Journal”

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 595 records · Page 33Linked to original sources

USGS2003

A shallow fault-zone structure illuminated by trapped waves in the Karadere-Duzce branch of the North Anatolian Fault, western Turkey

We discuss the subsurface structure of the Karadere-Duzce branch of the North Anatolian Fault based on analysis of a large seismic data set recorded by a local PASSCAL network in the 6 months following the Mw = 7.4 1999 Izmit earthquake. Seismograms observed at stations located in the immediate vicinity of the rupture zone show motion amplification and long-period oscillations in both P- and S-wave trains that do not exist in nearby off-fault stations. Examination of thousands of waveforms reveals that these characteristics are commonly generated by events that are well outside the fault zone. The anomalous features in fault-zone seismograms produced by events not necessarily in the fault may be referred to generally as fault-zone-related site effects. The oscillatory shear wave trains after the direct S arrival in these seismograms are analysed as trapped waves propagating in a low-velocity fault-zone layer. The time difference between the S arrival and trapped waves group does not grow systematically with increasing source-receiver separation along the fault. These observations imply that the trapping of seismic energy in the Karadere-Duzce rupture zone is generated by a shallow fault-zone layer. Traveltime analysis and synthetic waveform modelling indicate that the depth of the trapping structure is approximately 3-4 km. The synthetic waveform modelling indicates further that the shallow trapping structure has effective waveguide properties consisting of thickness of the order of 100 m, a velocity decrease relative to the surrounding rock of approximately 50 per cent and an S-wave quality factor of 10-15. The results are supported by large 2-D and 3-D parameter space studies and are compatible with recent analyses of trapped waves in a number of other faults and rupture zones. The inferred shallow trapping structure is likely to be a common structural element of fault zones and may correspond to the top part of a flower-type structure. The motion amplification associated with fault-zone-related site effects increases the seismic shaking hazard near fault-zone structures. The effect may be significant since the volume of sources capable of generating motion amplification in shallow trapping structures is large.

Geophysical Journal International↗

USGS2012

Acceleration spectra for subduction zone earthquakes

We estimate the source spectra of shallow earthquakes from digital recordings of teleseismic P wave groups, that is, P + p P + s P , by making frequency dependent corrections for the attenuation and for the interference of the free surface. The correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths. We apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from M S = 6.2 to 8.1, obtaining corrected P wave acceleration spectra on the frequency band from 0.01 to 2.0 Hz. Seismic moment estimates from surface waves and normal modes are used to extend these P wave spectra to the frequency band from 0.001 to 0.01 Hz. The acceleration spectra of moderate subduction zone earthquakes, that is, earthquakes whose seismic moments are less than 10 27 dyn cm, exhibit ω-square or Brune-type spectra, while the acceleration spectra of large subduction zone earthquakes, that is, earthquakes whose seismic moments are greater than 10 27 dyn cm, exhibit intermediate slopes where ü(ω) ∝ ω 5/4 for frequencies from 0.005 to 0.05 Hz. For this set of earthquakes, spectral shape appears to be a discontinuous function of seismic moment. Using reasonable assumptions for the phase characteristics, we transform the spectral shape observed for large earthquakes into the time domain to fit Ekström's (1987) moment rate functions for the M S = 8.1 Michoacan earthquake of September 19, 1985, and the M S = 7.6 Michoacan aftershock of September 21, 1985.

Journal of Geophysical Research Solid Earth↗

USGS2013

Scenarios of bioenergy development impacts on regional groundwater withdrawals

Irrigation increases agricultural productivity, but it also stresses water resources (Huffaker and Hamilton 2007). Drought and the potential for drier conditions resulting from climate change could strain water supplies in landscapes where human populations rely on finite groundwater resources for drinking, agriculture, energy, and industry (IPCC 2007). For instance, in the North American Great Plains, rowcrops are utilized for livestock feed, food, and bioenergy production (Cassman and Liska 2007), and a large portion is irrigated with groundwater from the High Plains aquifer system (McGuire 2011). Under projected future climatic conditions, greater crop water use requirements and diminished groundwater recharge rates could make rowcrop irrigation less feasible in some areas (Rosenberg et al. 1999; Sophocleous 2005). The Rainwater Basin region of south central Nebraska, United States, is an intensively farmed and irrigated Great Plains landscape dominated by corn (Zea mays L.) and soybean (Glycine max L.) production (Bishop and Vrtiska 2008). Ten starch-based ethanol plants currently service the region, producing ethanol from corn grain (figure 1). In this study, we explore the potential of switchgrass (Panicum virgatum L.), a drought-tolerant alternative bioenergy feedstock, to impact regional annual groundwater withdrawals for irrigation under warmer and drier future conditions. Although our research context is specific to the Rainwater Basin and surrounding North American Great Plains, we believe the broader research question is internationally pertinent and hope that this study simulates similar research in other areas.

Nebraska↗

USGS2000

Characteristic length scale of input data in distributed models: implications for modeling grain size

The appropriate spatial scale for a distributed energy balance model was investigated by: (a) determining the scale of variability associated with the remotely sensed and GIS-generated model input data; and (b) examining the effects of input data spatial aggregation on model response. The semi-variogram and the characteristic length calculated from the spatial autocorrelation were used to determine the scale of variability of the remotely sensed and GIS-generated model input data. The data were collected from two hillsides at Upper Sheep Creek, a sub-basin of the Reynolds Creek Experimental Watershed, in southwest Idaho. The data were analyzed in terms of the semivariance and the integral of the autocorrelation. The minimum characteristic length associated with the variability of the data used in the analysis was 15 m. Simulated and observed radiometric surface temperature fields at different spatial resolutions were compared. The correlation between agreement simulated and observed fields sharply declined after a 10×10 m 2 modeling grid size. A modeling grid size of about 10×10 m 2 was deemed to be the best compromise to achieve: (a) reduction of computation time and the size of the support data; and (b) a reproduction of the observed radiometric surface temperature.

Idaho↗

USGS2018

Induced earthquake families reveal distinctive evolutionary patterns near disposal wells

The timing of events in seismic sequences can provide insights into the physical processes controlling fault slip. In southern Kansas, the rate of earthquakes rose rapidly starting in 2013 following expansion of energy production into the area, demanding the disposal of large volumes of wastewater into deep wells. Seismicity catalogs that are complete to low magnitudes can provide insights into the physical processes that induce seismicity near wastewater disposal. We develop a catalog of over 130,000 earthquakes recorded in southern Kansas from mid‐March 2014 through December 2017 by applying a matched filter algorithm to an original catalog of 5,831 template earthquakes. Detections have nearly identical waveforms to their associated template event and represent slip on nearly co‐located sections of a fault. We select template events with at least 100 associated detections and examine the characteristics of these prolific families of earthquakes. We find that families located close (<10 km) to areas with significant volumes of injected fluids have near‐Poissonian interevent times and the families remain active over longer durations. Families farther from high‐volume injection wells show strong clustering of interevent times and shorter sequence durations. We conclude that increasing pore fluid pressures from nearby disposal of large volumes of wastewater is the primary driver of these long duration episodes, with earthquake‐earthquake interactions driving sequences at greater distance from the wells.

Kansas↗

USGS2017

A comparison of three federal datasets for thermoelectric water withdrawals in the United States for 2010

Historically, thermoelectric water withdrawal has been estimated by the Energy Information Administration (EIA) and the U.S. Geological Survey's (USGS) water-use compilations. Recently, the USGS developed models for estimating withdrawal at thermoelectric plants to provide estimates independent from plant operator-reported withdrawal data. This article compares three federal datasets of thermoelectric withdrawals for the United States in 2010: one based on the USGS water-use compilation, another based on EIA data, and the third based on USGS model-estimated data. The withdrawal data varied widely. Many plants had three different withdrawal values, and for approximately 54% of the plants the largest withdrawal value was twice the smallest, or larger. The causes of discrepancies among withdrawal estimates included definitional differences, definitional noise, and various nondefinitional causes. The uncertainty in national totals can be characterized by the range among the three datasets, from 5,640 m 3 /s (129 billion gallons per day [bgd]) to 6,954 m 3 /s (158 bgd), or by the aggregate difference between the smallest and largest values at each plant, from 4,014 m 3 /s (92 bgd) to 8,590 m 3 /s (196 bgd). When used to assess the accuracy of reported values, the USGS model estimates identify plants that need to be reviewed.

Journal of the American Water Resources Associatio↗

USGS2012

Monitoring glacier surface seismicity in time and space using Rayleigh waves

Sliding glaciers and brittle ice failure generate seismic body and surface wave energy characteristic to the source mechanism. Here we analyze continuous seismic recordings from an array of nine short-period passive seismometers located on Bench Glacier, Alaska (USA) (61.033°N, 145.687°W). We focus on the arrival-time and amplitude information of the dominant Rayleigh wave phase. Over a 46-hour period we detect thousands of events using a cross-correlation based event identification method. Travel-time inversion of a subset of events (7% of the total) defines an active crevasse, propagating more than 200 meters in three hours. From the Rayleigh wave amplitudes, we estimate the amount of volumetric opening along the crevasse as well as an average bulk attenuation ( = 42) for the ice in this part of the glacier. With the remaining icequake signals we establish a diurnal periodicity in seismicity, indicating that surface run-off and subglacial water pressure changes likely control the triggering of these surface events. Furthermore, we find that these events are too weak (i.e., too noisy) to locate individually. However, stacking individual events increases the signal-to-noise ratio of the waveforms, implying that these periodic sources are effectively stationary during the recording period.

Journal of Geophysical Research F: Earth Surface↗

USGS2023

Evaluation of hydrodynamic mixing in an afterbay reservoir

This study focused on the mixing of a solute, assumed to be conservative, introduced to one arm of an afterbay reservoir, between Keswick and Shasta Dams on the Sacramento River near Redding, California. Rhodamine water tracer (WT) dye served as the solute in a field experiment, and was introduced over 4.5 h and monitored for 4 days by sondes moored in the reservoir. The scenario was modeled numerically using the Delft3D flexible mesh (FM) hydrodynamic and mixing model, with measured inflows, outflows, water level, water temperatures, and bathymetry as input. Manning’s n "> n and horizontal eddy viscosity served as the (constant) model calibration parameters, and each was adjusted an order of magnitude below the default values to force observed and modeled dye hydrographs to match in arrival time and duration. The low friction factor was concluded to be due to a combination of low flow speeds coupled with energy dissipation inherent to the model. The model and surface drifters equipped with dual-frequency Global Navigation Satellite System equipment revealed velocities in the 1 – 5    cm / s "> 1 – 5 cm / s range in much of the domain during the experiment. Simple analytical expressions were shown to be useful for estimating distance to full cross-sectional mixing, steady-state concentrations, and time to reach them, but the numerical model is required for investigation of the approach to steady state, and at locations where flows intersect. Time to steady-state concentrations was 1.5–13 days for 10 simulations that spanned a wide range of inflow conditions. Model sensitivity tests suggest that wind and heat fluxes were not important during the field study, but simulations of a summer scenario with small inflows of cold water upstream in warm weather should consider water temperature. Both field observations and numerical model results showed inflow to one arm of the reservoir reaching full cross-sectional mixing before plunging below the water surface near the intersection of this arm with the reservoir’s main stem. Model results are being used to guide management decisions related to inflows to the reservoir from a relic mining site that is also a USEPA Superfund site.

California↗

USGS2014

Age-Ratios and Condition of En Route Migrant Blackpoll Warblers in the British Virgin Islands

The en route migration ecology of Blackpoll Warblers ( Setophaga striata ) is poorly understood, yet intriguing. Blackpoll Warblers undertake the longest open water migration of any wood warbler species, traveling from northeastern North America to South America, with the first potential landfall being the West Indies. This migration requires substantial energy reserves and subjects Blackpoll Warblers to unpredictable weather events, which may influence survival. Few studies have examined age ratios or condition of Blackpoll Warblers while the warblers are en route through the Caribbean region. I captured and banded Blackpoll Warblers in the British Virgin Islands over 10 consecutive autumn migrations. Ratios of hatch-year to adult Blackpoll Warblers were variable but averaged lower than the ratios reported at continental departure locations. Average mass of Blackpoll Warblers was less than that reported at continental departure locations, with 26% of adults and 40% of hatch-year birds below the estimated fat free mass; hatch-year birds were consistently in poorer condition than adults. Blackpoll Warblers captured in the British Virgin Islands were also in poorer condition than those reported from the Dominican Republic and Barbados; this may be because of the British Virgin Islands being the first landfall after the transatlantic crossing, whereas Blackpoll Warblers arriving at the other Caribbean study locations may have had opportunities for stopover prior to arrival or have departed from farther south on the continent. However, this suggests that the British Virgin Islands likely provide important stopover habitat as a first landfall location for Blackpoll Warblers arriving from the transatlantic migration route.

Guana Island↗

USGS2022

Pore systems and organic petrology of cretaceous Mowry and Niobrara source-rock reservoirs, Powder River Basin, Wyoming, USA

The Powder River Basin (PRB) is a world-class oil province, in large part thanks to contributions from premier source rocks, Cretaceous Mowry and Niobrara shales. Both formations are also unconventional reservoirs. A critical aspect of evaluating production potential and finding sweet spots is the nature of the pore systems in these fine-grained source-rock reservoirs. Variation by stratigraphic interval is important for selecting optimum target zones for horizontal wells. Understanding variation in pore type, size, and connectivity and relationships with mineralogy and fabric help in determining prospectivity in different parts of the basin. Deciphering controls on pore-system development helps predict intervals and locations of optimum reservoir quality. Imaging of Niobrara and Mowry samples from a range of thermal maturities provided observations and data on pore systems, organic matter (OM) types and associations with mineralogy and fabric, wettability, and microporosity associated with both diagenetic and detrital clays. Imaging techniques included scanning electron microscopy, organic petrography and correlative scanning electron microscopy, and mapping of mineralogy through energy dispersive spectroscopy. Mean solid bitumen (BR o ) and vitrinite reflectance (VR o ) values indicate all samples are in the oil window with values ranging from 0.52 to 1.15%. Organic fluorescence is prominent in amorphous OM, solid bitumen and some vitrinite in the early oil window. The fluorescence is extinguished at higher thermal maturity. Carbonate pellets (in Niobrara) mainly contain migrated solid bitumen and residual live oil and little or no terrigenous OM (vitrinite and inertinite). However, terrigenous OM is common in siliceous/argillaceous laminae in both formations, where it occurs with amorphous OM, some of which has converted in situ to a solid bitumen petroleum residue. One key finding is the widespread presence of migrated OM at very early oil window maturity. Distribution of such OM and associated wettability alteration is fabric-controlled, at all levels of thermal maturity studied. Clay morphology and abundance and supporting rigid mineral grain framework strongly influence pore development, preservation, and connectivity in both formations. Carbonate content is a good proxy for reservoir quality in Niobrara intervals due to association of porous solid bitumen with calcareous fecal pellets. High recrystallized microquartz content is associated with the best reservoir intervals in the Mowry.

Wyoming↗

USGS2018

Agassiz’s desert tortoise (Gopherus agassizii) activity areas are little changed after wind turbine-induced fires in California

Wind turbine-induced fires at a wind energy facility in California, USA, provided an opportunity to study the before and after effects of fire on a population of protected Agassiz’s desert tortoises ( Gopherus agassizii ) in the Sonoran Desert, a species and ecosystem poorly adapted to fire. We compared annual activity areas (AAs) of tortoises in 2011 and 2013, before and after two 2012 fires, with those of tortoises in adjacent areas unaffected by the same fires. Tortoises in both AAs affected by fire or unaffected by fire occupied the same general AAs in 2013, after the fires, as they did in 2011, before the fires. Some tortoises had both their 2011 and 2013 AAs completely or almost completely within the areas burned by the 2012 fires, despite the proximity of unburned habitat. None of the tortoises with 2011 AAs subsequently unaffected by the 2012 fires shifted their AAs into burned habitat in 2013. For the fire-affected group of tortoises, the mean percentages of 2011 and 2013 AAs burned by the 2012 fires were not significantly different, showing fidelity to the burned areas. Tortoises in both groups generally occupied consistent AAs, even post fire, placing them at potential risk of exposure to unfavourable burned habitat.

California↗

USGS2023

Crop water use dynamics over arid and semi-arid croplands in the lower Colorado River Basin

Numerous studies have evaluated the application of Remote Sensing (RS) techniques for mapping actual evapotranspiration (ETa) using Vegetation-Index-based (VI-based) and surface energy balance methods (SEB). SEB models computationally require a large effort for application. VI-based methods are fast and easy to apply and could therefore potentially be applied at high resolution; however, the accuracy of VI-based methods in comparison to SEB-based models remains unclear. We tested the ETa computed with the modified 2-band Enhanced Vegetation Index (METEVI2) implemented in the Google Earth Engine – for mapping croplands’ water use dynamics in the Lower Colorado River Basin. We compared METEVI2 with the well-established RS-based products of OpenET (Ensemble, eeMETRIC, SSEBop, SIMS, PT_JPL, DisALEXI and geeSEBAL). METEVI2 was then evaluated with measured ETa from four wheat fields (2017–2018). Results indicated that the monthly ETa variations for METEVI2 and OpenET models were comparable, though of varying magnitudes. On average, METEVI2 had the lowest difference rate from the average observed ETa with 17 mm underestimation, while SIMS had the highest difference rate (82 mm). Findings show that METEVI2 is a cost-effective ETa mapping tool in drylands to track crop water use. Future studies should test METEVI2’s applicability to croplands in more humid regions.

Lower Colorado River basin↗

USGS2016

Seismo-acoustic evidence for an avalanche driven phreatic eruption through a beheaded hydrothermal system: An example from the 2012 Tongariro eruption

The 6 August 2012 Te Maari eruption comprises a complex eruption sequence including multiple eruption pulses, a debris avalanche that propagated ~ 2 km from the vent, and the formation of a 500 m long, arcuate chasm, located ~ 300 m from the main eruption vent. The eruption included 6 distinct impulses that were coherent across a local infrasound network marking the eruption onset at 11:52:18 (all times UTC). An eruption energy release of ~ 3 × 1012 J was calculated using a body wave equation for radiated seismic energy. A similar calculation based on the infrasound record, shows that ~ 90% of the acoustic energy was released from three impulses at onset times 11:52:20 (~ 20% of total eruption energy), 11:52:27 (~ 50%), and 11:52:31 (~ 20%). These energy impulses may coincide with eyewitness accounts describing an initial eastward directed blast, followed by a westward directed blast, and a final vertical blast. Pre-eruption seismic activity includes numerous small unlocatable micro-earthquakes that began at 11:46:50. Two larger high frequency earthquakes were recorded at 11:49:06 and 11:49:21 followed directly by a third earthquake at 11:50:17. The first event was located within the scarp based on an arrival time location from good first P arrival times and probably represents the onset of the debris avalanche. The third event was a tornillo, characterised by a 0.8 Hz single frequency resonance, and has a resonator attenuation factor of Q ~ 40, consistent with a bubbly fluid filled resonator. This contrasts with a similar tornillo event occurring 2.5 weeks earlier having Q ~ 250–1000, consistent with a dusty gas charged resonator. We surmise from pre-eruption seismicity, and the observed attenuation change, that the debris avalanche resulted from the influx of fluids into the hydrothermal system, causing destabilisation and failure. The beheaded hydrothermal system may have then caused depressurisation frothing of the remaining gas charged system leading to the onset of explosive activity.

Journal of Volcanology and Geothermal Research↗

USGS1987

A Pennsylvanian-age terrestrial storm deposit: using plant fossils to characterize the history and process of sediment accumulation

A thin black shale overlying the B-coal underclay (in the Middle Pennsylvanian post-Pottsville strata of the Bernice Basin) contains a compression flora composed of large, prostrate, unidirectionally oriented lycopod trunks and randomly oriented pteridosperm stems. Analyses of modern log accumulations indicate that unidirectional trunk orientations can be produced by riverflood currents, volcanic blasts, and most high-energy windstorms. Since there are neither fluvial sediments nor ash deposits associated with the Bernice assemblage, this deposit is believed to have been formed by high-energy winds. Furthermore, this deposit is interpreted to be in situ because storm winds (and volcanic blasts) rarely have sufficient energy for the physical transport of large, intact tree trunks. The sedimentary history of the B-coal underclay can be determined from the successional changes in the species and plant part compositions (leaves, seeds, branches, trunks, etc.) of the preserved plant material. The underclay is an accretionary floodplain soil which accumulated as discrete increments during episodic floods. The sediments deposited with each flood incorporated the litter layer of the lycopod-pteridosperm forest which occupied this site. Ordinarily, the flood water would recede, and renewed root growth would destroy the primary sedimentary structures and the newly incorporated organic material. Because the bedding and forest litter are preserved in the top 5 cm of the underclay, root growth and silt deposition must have been terminated by the last flooding event. The site eventually became permanently inundated, and an organic-rich mud began to accumulate in the resulting floodplain lake. The lycopod-pteridosperm forest drowned and, at some later time, was blown down into the lake. The trunks are preserved on a single bedding plane in a 2-cm-thick, organic-rich lacustrine black shale. Continued organic accumulation in the lake resulted in the accumulation of a hypautochthonous peat which eventually was colonized by a peat-forming flora.

Journal of Sedimentary Petrology↗

USGS2021

Subsurface swimming and stationary diving are metabolically cheap in adult Pacific walruses (Odobenus rosmarus divergens)

Walruses rely on sea-ice to efficiently forage and rest between diving bouts while maintaining proximity to prime foraging habitat. Recent declines in summer sea ice have resulted in walruses hauling out on land where they have to travel farther to access productive benthic habitat while potentially increasing energetic costs. Despite the need to better understand the impact of sea ice loss on energy expenditure, knowledge about metabolic demands of specific behaviours in walruses is scarce. In the present study, 3 adult female Pacific walruses ( Odobenus rosmarus divergens ) housed in professional care participated in flow-through respirometry trials to measure metabolic rates while floating inactive at the water surface during a minimum of 5 min, during a 180 s stationary dive, and while swimming ∼90 m horizontally underwater. Metabolic rates during stationary dives (3.82±0.56 l O 2 min −1 ) were lower than those measured at the water surface (4.64±1.04 l O 2 min −1 ), which did not differ from rates measured during subsurface swimming (4.91±0.77 l O 2 min −1 ). Thus, neither stationary diving nor subsurface swimming resulted in metabolic rates above those exhibited by walruses at the water surface. These results suggest that walruses minimize their energetic investment during underwater behaviours as reported for other marine mammals. Although environmental factors experienced by free-ranging walruses (e.g. winds or currents) likely affect metabolic rates, our results provide important information for understanding how behavioural changes affect energetic costs and can be used to improve bioenergetics models aimed at predicting the metabolic consequences of climate change on walruses.

Journal of Experimental Biology↗

USGS2005

Perceptions of earthquake and tsunami issues in U.S. Pacific Northwest port and harbor communities

Although there is considerable energy focused on assessing natural hazards associated with earthquakes and tsunamis in the U.S. Pacific Northwest, little has been done to understand societal vulnerability to these hazards. Part of understanding societal vulnerability includes assessing the perceptions and priorities of public sector individuals with traditional emergency management responsibilities and of private citizens who could play key roles in community recovery. In response to this knowledge gap, we examine earthquake and tsunami perceptions of stakeholders and decision makers from coastal communities in the U.S. Pacific Northwest, focusing on perceptions of (1) regional hazards and societal vulnerability, (2) the current state of readiness, and (3) priorities for future hazard adjustment efforts. Results of a mailed survey suggest that survey participants believe that earthquakes and tsunamis are credible community threats. Most communities are focusing on regional mitigation and response planning, with less effort devoted to recovery plans or to making individual organizations more resilient. Significant differences in expressed perceptions and priorities were observed between Oregon and Washington respondents, mainly on tsunami issues. Significant perception differences were also observed between private and public sector respondents. Our results suggest the need for further research and for outreach and planning initiatives in the Pacific Northwest to address significant gaps in earthquake and tsunami hazard awareness and readiness.

International Journal of Mass Emergencies and Disa↗

USGS1996

The cyclic and fractal seismic series preceding an mb 4.8 earthquake on 1980 February 14 near the Virgin Islands

Seismic activity in the 10 months preceding the 1980 February 14, m b 4.8 earthquake in the Virgin Islands, reported on by Frankel in 1982, consisted of four principal cycles. Each cycle began with a relatively large event or series of closely spaced events, and the duration of the cycles progressively shortened by a factor of about 3/4. Had this regular shortening of the cycles been recognized prior to the earthquake, the time of the next episode of seismicity (the main shock) might have been closely estimated 41 days in advance. That this event could be much larger than the previous events is indicated from time-to-failure analysis of the accelerating rise in released seismic energy, using a non-linear time- and slip-predictable foreshock model. Examination of the timing of all events in the sequence shows an even higher degree of order. Rates of seismicity, measured by consecutive interevent times, when plotted on an iteration diagram of a rate versus the succeeding rate, form a triangular circulating trajectory. The trajectory becomes an ascending helix if extended in a third dimension, time. This construction reveals additional and precise relations among the time intervals between times of relatively high or relatively low rates of seismic activity, including period halving and doubling. The set of 666 time intervals between all possible pairs of the 37 recorded events appears to be a fractal; the set of time points that define the intervals has a finite, non-integer correlation dimension of 0.70. In contrast, the average correlation dimension of 50 random sequences of 37 events is significantly higher, close to 1.0. In a similar analysis, the set of distances between pairs of epicentres has a fractal correlation dimension of 1.52. Well-defined cycles, numerous precise ratios among time intervals, and a non-random temporal fractal dimension suggest that the seismic series is not a random process, but rather the product of a deterministic dynamic system.

Geophysical Journal International↗

USGS1998

Nonlethal evaluation of the physiological health of unionid mussels: Method for biopsy and glycogen analysis

In long-lived unionid mussels, many short-term measures of growth are of limited value. Changes in physiological condition may be an early indication of stress, because the increased energy demand associated with stress often results in a depletion of glycogen reserves, the principal storage form of carbohydrates in unionid mussels. Our goal was to nonlethally extract tissue from freshwater mussels and then to develop a rapid and dependable method for the analysis of glycogen in the tissue extracts. A biopsy technique was developed to remove between 5 and 10 mg of foot tissue in Amblema plicata plicata. The survival rate did not differ between biopsied and non-biopsied mussels during a 581-d observation period, demonstrating that the biopsy technique will allow nonlethal evaluation of the physiological condition of individual mussels through measurement of changes in contaminant, genetic, and biochemical indicators in tissue. We also modified the standard alkaline digestion and phenol-sulfuric acid analysis of glycogen for use on the small samples of biopsied tissue and to reduce analysis time and cost. We present quality control data, including method detection limits and estimates of precision and bias. The modified analytical method is rapid and accurate and has a method detection limit of 0.014 mg glycogen. Glycogen content in the biopsied samples was well above the method detection limit; it ranged from 0.09 to 0.36 mg, indicating that the method should be applicable to native mussels.

Journal of the North American Benthological Societ↗