Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Elements”

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

Rare earth element metasomatism in hydrothermal systems: The Willsboro-Lewis wollastonite ores, New York, USA

Wollastonite ores and garnet-pyroxene skarns in the Willsboro-Lewis district, New York, USA were formed in a complex hydrothermal system associated with the emplacement of a large anorthosite pluton. Contact-metamorphic marbles were replaced by wollastonite, garnet, and clinopyroxene during infiltration metasomatism involving large volumes of water of chiefly meteoric origin. Rare earth elements (REE) in these rocks show large departures from the protolith REE distribution, indicative of substantial REE mobility. Three types of chondrite-normalized REE distribution patterns are present. The most common, found in ores and skarns containing andradite-rich garnet, is convex-up in the light REE (LREE) with a maximum at Pr and a positive Eu anomaly. Europium anomalies and Pr/Yb ratios are correlated with X(Ad) in garnet. This pattern (type C) results from uptake of REE from hydrothermal fluids by growing crystals of calcsilicate minerals, principally andradite, with amounts of LREE controlled by the difference in ionic radius between Ca++ and REE3+ in garnet X sites. The Eu anomaly results either from prior interaction of the fluids with plagioclase-rich, Eu-positive anorthositic rocks in and near the ore zone, or by enrichment of divalent Eu on growth surfaces of garnet followed by entrapment, or both. Relative enrichment in heavy REE (type H) occurs in ores and skarn where calcsilicates, including grossularitic garnet, in contact-metamorphic marble have been concentrated by dissolution of calcite. In most cases a negative Eu anomaly is inherited from the marble protolith. Skarns containing titanite and apatite exhibit high total REE, relative light REE enrichment, and negative Eu anomalies (type L). These appear to be intrusive igneous rocks (ferrodiorites or anorthositic gabbros) that have been converted to skarn by Ca metasomatism. REE, sequestered in titanite, apatite, and garnet, preserve the approximate REE distribution pattern of the igneous protolith. Post-ore granulite facies metamorphism homogenized zoned mineral grains without causing complete intergranular reequilibration and does not appear to have significantly affected the whole-rock REE distributions. These results demonstrate that extensive REE metasomatism can occur in hydrothermal systems at shallow to intermediate depths and that REE geochemistry may be useful in discerning the origin of skarns and skarn-related ore deposits.

Geochimica et Cosmochimica Acta↗

Copper localization, elemental content, and thallus colour in the copper hyperaccumulator lichen Lecanora sierrae from California

An unusual dark blue-green lichen, Lecanora sierrae , was discovered over 30 years ago by Czehura near copper mines in the Lights Creek District, Plumas County, Northern California. Using atomic absorption spectroscopy, Czehura found that dark green lichen samples from Warren Canyon contained 4% Cu in ash and suggested that its colour was due to copper accumulation in the cortex. The present study addressed the hypothesis that the green colour in similar material we sampled from Warren Canyon in 2008, is caused by copper localization in the thallus. Optical microscopy and electron microprobe analysis of specimens of L. sierrae confirmed that copper localization took place in the cortex. Elemental analyses of L. sierrae and three other species from the same localities showed high enrichments of copper and selenium, suggesting that copper selenates or selenites might occur in these lichens and be responsible for the unusual colour.

Lichenologist↗

Concentration method for the spectrochemical determination of seventeen minor elements in natural water

A method for the quantitative spectrochemical determination of microgram amounts of 17 minor elements in water is given. The chelating reagents 8-quinolinol, tannic acid, and thionalide are utilized to concentrate traces (1 to 500 μ g.) of aluminum, cobalt, chromium, copper, iron, gallium, germanium, manganese, nickel, titanium, vanadium, bismuth, lead, molybdenum, cadmium, zinc, and beryllium. Indium is added as a buffer, and palladium is used as an internal standard. The ashed oxides of these 17 metals are subsequently subjected to direct current arcing conditions during spectrum analysis. The method can be used to analyze waters with dissolved solids ranging from less than 100 to more than 100,000 p.p.m. There is no limiting concentration range for the determination of the heavy metals since any volume of sample can be used that will contain a heavy metal concentration within the analytical range of the method. Both the chemical and spectrographic procedures are described, and precision and accuracy data are given.

Analytical Chemistry↗

Effects of an extreme flood on trace elements in river water—From urban stream to major river basin

Major floods adversely affect water quality through surface runoff, groundwater discharge, and damage to municipal water infrastructure. Despite their importance, it can be difficult to assess the effects of floods on streamwater chemistry because of challenges collecting samples and the absence of baseline data. This study documents water quality during the September 2013 extreme flood in the South Platte River, Colorado, USA. Weekly time-series water samples were collected from 3 urban source waters (municipal tap water, streamwater, and wastewater treatment facility effluent) under normal-flow and flood conditions. In addition, water samples were collected during the flood at 5 locations along the South Platte River and from 7 tributaries along the Colorado Front Range. Samples were analyzed for 54 major and trace elements. Specific chemical tracers, representing different natural and anthropogenic sources and geochemical behaviors, were used to compare streamwater composition before and during the flood. The results differentiate hydrological processes that affected water quality: (1) in the upper watershed, runoff diluted most dissolved constituents, (2) in the urban corridor and lower watershed, runoff mobilized soluble constituents accumulated on the landscape and contributed to stream loading, and (3) flood-induced groundwater discharge mobilized soluble constituents stored in the vadose zone.

Colorado↗

Evaluating the behavior of gadolinium and other rare earth elements through large metropolitan sewage treatment plants

A primary pathway for emerging contaminants (pharmaceuticals, personal care products, steroids, and hormones) to enter aquatic ecosystems is effluent from sewage treatment plants (STP), and identifying technologies to minimize the amount of these contaminants released is important. Quantifying the flux of these contaminants through STPs is difficult. This study evaluates the behavior of gadolinium, a rare earth element (REE) utilized as a contrasting agent in magnetic resonance imaging (MRI), through four full-scale metropolitan STPs that utilize several biosolids thickening, conditioning, stabilization, and dewatering processing technologies. The organically complexed Gd from MRIs has been shown to be stable in aquatic systems and has the potential to be utilized as a conservative tracer in STP operations to compare to an emerging contaminant of interest. Influent and effluent waters display large enrichments in Gd compared to other REEs. In contrast, most sludge samples from the STPs do not display Gd enrichments, including primary sludges and end-product sludges. The excess Gd appears to remain in the liquid phase throughout the STP operations, but detailed quantification of the input Gd load and residence times of various STP operations is needed to utilize Gd as a conservative tracer.

Environmental Science & Technology↗

Metal and trace element burdens in two shorebird species at two sympatric wintering sites in southern California

In January and November, 1994, ten willets (Catoptrophorus semipalmatus) and ten black-bellied plovers (Pluvialis squatarola) were collected at a south-facing sandy beach at the North Island Naval Air Station (NASNI) and ten willets and three plovers were collected at the Tijuana Slough National Wildlife Refuge (TSNWR), 16 km to the southeast. Concentrations of 19 inorganic analytes were determined in the livers. Differences between sexes in diet or abilities to retain some analytes are indicated because male willets tend to have greater levels of Cd than females and male plovers tend to have greater levels of Fe, Hg, and Mn than females. Concentrations of Hg and Mg in willets decreased from January to November. If this decrease is a part of an annual cycle, then the birds may be acquiring those elements while at NASNI. The differences between sites shows higher levels of As and Se but lower levels of Mn in willets at NASNI while plovers have lower levels of Cu, Mg and Zn at NASNI. There were no correlations between Hg and Se concentrations in the livers of any sample group of birds.

California↗

Relative effect of temperature and pH on diel cycling of dissolved trace elements in Prickly Pear Creek, Montana

Diel (24 hr) cycles in dissolved metal and As concentrations have been documented in many northern Rocky Mountain streams in the U.S.A. The cause(s) of the cycles are unknown, although temperature- and pH-dependent sorption reactions have been cited as likely causes. A light/dark experiment was conducted to isolate temperature and pH as variables affecting diel metal cycles in Prickly Pear Creek, Montana. Light and dark chambers containing sediment and a strand of macrophyte were placed in the stream to simulate instream temperature oscillations. Photosynthesis-induced pH changes were allowed to proceed in the light chambers while photosynthesis was prevented in the dark chambers. Water samples were collected periodically for 22 hr in late July 2001 from all chambers and the stream. In the stream, dissolved Zn concentrations increased by 300% from late afternoon to early morning, while dissolved As concentrations exhibited the opposite pattern, increasing 33% between early morning and late afternoon. Zn and As concentrations in the light chambers showed similar, though less pronounced, diel variations. Conversely, Zn and As concentrations in the dark chambers had no obvious diel variation, indicating that light, or light-induced reactions, caused the variation. Temperature oscillations were nearly identical between light and dark chambers, strongly suggesting that temperature was not controlling the diel variations. As expected, pH was negatively correlated ( P < 0.01) with dissolved Zn concentrations and positively correlated with dissolved As concentrations in both the light and dark chambers. From these experiments, photosynthesis-induced pH changes were determined to be the major cause of the diel dissolved Zn and As cycles in Prickly Pear Creek. Further research is necessary in other streams to verify that this finding is consistent among streams having large differences in trace-element concentrations and mineralogy of channel substrate.

Montana↗

Finite element model predictions of static deformation from dislocation sources in a subduction zone: Sensitivities to homogeneous, isotropic, Poisson-solid, and half-space assumptions

Dislocation models can simulate static deformation caused by slip along a fault. These models usually take the form of a dislocation embedded in a homogeneous, isotropic, Poisson-solid half-space (HIPSHS). However, the widely accepted HIPSHS assumptions poorly approximate subduction zone systems of converging oceanic and continental crust. This study uses three-dimensional finite element models (FEMs) that allow for any combination (including none) of the HIPSHS assumptions to compute synthetic Green's functions for displacement. Using the 1995 Mw = 8.0 Jalisco-Colima, Mexico, subduction zone earthquake and associated measurements from a nearby GPS array as an example, FEM-generated synthetic Green's functions are combined with standard linear inverse methods to estimate dislocation distributions along the subduction interface. Loading a forward HIPSHS model with dislocation distributions, estimated from FEMs that sequentially relax the HIPSHS assumptions, yields the sensitivity of predicted displacements to each of the HIPSHS assumptions. For the subduction zone models tested and the specific field situation considered, sensitivities to the individual Poisson-solid, isotropy, and homogeneity assumptions can be substantially greater than GPS. measurement uncertainties. Forward modeling quantifies stress coupling between the Mw = 8.0 earthquake and a nearby Mw = 6.3 earthquake that occurred 63 days later. Coulomb stress changes predicted from static HIPSHS models cannot account for the 63-day lag time between events. Alternatively, an FEM that includes a poroelastic oceanic crust, which allows for postseismic pore fluid pressure recovery, can account for the lag time. The pore fluid pressure recovery rate puts an upper limit of 10-17 m2 on the bulk permeability of the oceanic crust. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Mafic alkaline magmatism and rare earth element mineralization in the Mojave Desert, California: The Bobcat Hills connection to Mountain Pass

Occurrences of alkaline and carbonatite rocks with high concentrations of rare earth elements (REE) are a defining feature of Precambrian geology in the Mojave Desert of southeastern California. The most economically important occurrence is the carbonatite stock at Mountain Pass, which constitutes the largest REE deposit in the United States. A central scientific goal is to understand the genesis of the carbonatite ore body in the context of widespread REE-rich igneous activity. A swarm of mafic alkaline (shonkinite) dikes has been mapped and sampled at Bobcat Hills, 65 km southeast of the Mountain Pass mine. Whole-rock geochemistry and zircon geochronology demonstrate a clear affinity to the ca. 1.4 Ga Mountain Pass intrusive system. Bobcat Hills dikes have comparably high REE concentrations (La ∼1,000× chondritic) and an error-weighted mean 207 Pb/ 206 Pb zircon crystallization age of 1,426 ± 2 Ma (2 σ ). Unlike the alkaline intrusions at Mountain Pass, which have abundant inherited zircon from Paleoproterozoic basement rocks and crustally influenced oxygen isotope compositions (δ 18 O zircon = 6.5–7.5‰), the Bobcat Hills dikes lack any evidence of crustal assimilation and have oxygen isotope values that overlap a mantle range (Bobcat Hills average δ 18 O zircon = 5.6 ± 0.3‰). The dikes were a high-temperature, early center of mafic alkaline magmatism in the Mojave Desert that serve as a snapshot of melt generation from a spatially extensive, metasomatized mantle source. We propose that modification of the crust over many tens of Myr at Mountain Pass created an environment that favored crustal assimilation and enabled ascent of late-stage, REE-rich carbonatite magmas.

California↗

Multi-year tracing of spatial and temporal dynamics of post-fire aeolian sediment transport using rare earth elements provide insights into grassland management

Aeolian sediment transport occurs as a function of, and with feedback to ecosystem changes and disturbances. Many desert grasslands are undergoing rapid changes in vegetation, including the encroachment of woody plants, which alters fire regimes and in turn can change the spatial and temporal patterns of aeolian sediment transport. We investigated aeolian sediment transport and spatial distribution of sediment in the surface soil for 7 years following a prescribed fire using a multiple rare earth element (REE) tracer-based approach in a shrub-encroached desert grassland in the northern Chihuahuan desert. Results indicate that even though the aeolian horizontal sediment mass flux increased approximately three-fold in the first windy season in the burned areas compared to control areas, there were no significant differences after three windy seasons. The soil surface of bare microsites was the major contributor of aeolian sediments in unburned areas (87%), while the shrub microsites contributed the least (<2%) during the observation period. However, after the prescribed fire, the contribution of aeolian sediments from shrub microsites increased considerably (∼40%), indicating post-fire microsite-scale sediment redistribution. The findings of this study, which is the first to use multiple REE tracers for multi-year analysis of the spatial and temporal dynamics of aeolian sediment transport, illustrate how disturbance by prescribed fire can influence aeolian processes and alters dryland soil geomorphology in which distinct soils develop over time at very fine spatial scales of individual plants.

New Mexico↗

The McKinley Sequence of granitic rocks: A key element in the accretionary history of southern Alaska

The McKinley sequence of granitic rocks consists of several discrete plutons in the central Alaska Range. Most of these plutons crop out south of the Denali fault system (DFS) in the Talkeetna quadrangle. Plutons of the McKinley sequence largely intrude deformed upper Meszoic flysch between the DFS and the northern edges of Wrangellia and the Peninsular terrane, which jointly make up the Talkeetna superterrane. The average K-Ar age of biotite from nine granites of the McKinley sequence is 57.3 Ma; Rb-Sr data for whole rock samples indicate that the McKinley sequence cannot be older than 60 Ma. A selected suite of 20 samples of granite and granodiorite range in SiO 2 from 65.9 to 77.6%. All 20 samples are corundum normative, and 18 are moderately peraluminous. Initial 87 Sr/ 86 Sr ratios range from 0.7054 to 0.7085. The σ 18 O values range from +11.2 to +14.6‰. These high and variable Sr isotopic ratios, peraluminous nature, rare earth element patterns, and high σ 18 O values suggest that granitic rocks of the McKinley sequence crystallized from hybrid magmas produced by assimilation of sedimentary rocks by a mantle-derived melt. Mesozoic flysch is the likely source of the crustal component of the hybrid magmas. Geologic evidence suggests that the Talkeetna superterrane collided with stable Alaska after Early Cretaceous time. The flysch basin, lying south of stable Alaska, was closed by northward movement of the Talkeetna superterrane; maximum age for basin closure and terrane accretion is middle Cretaceous (Cenomanian). Paleomagnetic evidence indicates that all terranes north of the DFS have been part of stable Alaska since the Paleocene and that northward movement of Wrangellia was completed by 50 Ma. Granitic rocks of the McKinley sequence may be products of terrane accretion; the granitic rocks crystallized from hybrid magmas produced during terrane collision and deformation of the flysch basin. Isotopic ages of the McKinley sequence establish the time of final accretion of the Talkeetna superterrane as Paleocene.

Journal of Geophysical Research Solid Earth↗

Discrete element investigation of the influence of shallow soil density on the manifestations of strike-slip surface fault rupture

This study investigates the influence of soil relative density on strike-slip surface fault rupture manifestations using three-dimensional numerical simulations performed with the discrete element method (DEM). The simulations capture the formation of distinctive fault strands within complex flower structures using tens of millions of grains. The tendency for dense soils to localize shear manifests as multiple localized shear bands within wide zones of deformation in strike-slip fault rupture, whereas diffuse shear deformation develops within narrow zones in loose soils. The spatial extents of soil deformation are consistent between simulations having similar relative densities but with different quantities of grains. However, individual shears are more distinguishable in assemblages of finer grains than in assemblages of coarser grains. The simulations show the progressive development of new shears within the bounds of previously developed shears. Shear activity transitions inward as fault activity diminishes along the outermost shears and continues along newly developed shears until a vertically dipping throughgoing shear structure develops that accommodates most of the fault displacement thereon. The throughgoing fault develops at smaller fault displacements in looser soils because the first shear rupture propagates closer to the vertical direction and does not undergo as much inward translation of shear activity, as is observed in denser soils. In all simulations, ground surface uplift develops between nonintersecting active shears, and ground surface subsidence tends to develop where new shears intersect previous shears. The surface traces in these simulations are shown to be consistent with analog models and case histories of surface fault rupture occurring in different shallow subsurface materials. Although computationally costly, these modeling results are valuable for providing a strong numerical supplement to traditional analog models used to represent the mechanics of strike-slip zones in soil, and they provide quantifiable stresses and large-strain deformations throughout the model domain.

Journal of Geotechnical and Geoenvironmental Engin↗

Identifying common decision problem elements for the management of emerging fungal diseases of wildlife

Emerging infectious diseases (EIDs) of wildlife have characteristics that make them difficult to manage, leading to reactive and often ineffective management strategies. Currently, two fungal pathogens, Pseudogymnoascus destructans (Pd) and Batrachochytrium salamandrivorans (Bsal), are causing declines in novel host species. To improve the application of management strategies addressing the risk of these pathogens to North American wildlife, we queried wildlife managers about their concerns regarding managing populations of bats and amphibians potentially impacted by Pd and Bsal. Using these responses, we identified aspects of each decision problem that were shared across pathogens, regions and agencies – and found similarities in decision-problem elements for disease management. Reframing management problems as decisions can enable managers to identify similarities across EIDs, i.e. uncertainties within management actions, and improve reactive responses if proactive management is not possible. Such an approach recognizes context-specific constraints and identifies relevant uncertainties that must be reduced in developing a response.

Society and Natural Resources↗

Coseismic and post-seismic gravity disturbance induced by seismic sources using a 2.5-D spectral element method

I present a prescription for computing free-air coseismic and post-seismic gravity changes induced by seismic sources in a viscoelastic earth model. I assume a spherical earth geometry and a 2.5-D calculation, that is, 3-D motions that satisfy the equations of quasi-static equilibrium on a 2-D viscoelastic structure. The prescription permits application to regional gravity computations where a 2-D structure adequately represents the structural heterogeneity. I use a hybrid approach where deformation is computed on a discretized domain and the resulting density perturbations are expanded with spherical harmonics to produce the free-air gravity field. Starting with a solution to the equations of quasi-static displacements in the Laplace transform domain for a given dislocation source, I solve Poisson’s equation using Lagrangian interpolation on spectral element nodes to compute the required deformation quantities that contribute to free-air gravity. A numerical inverse Laplace transform then yields time domain results. This methodology is tested with analytic solutions on a spherically stratified viscoelastic structure, then applied to evaluate the effect of a descending slab of relatively high viscosity on post-seismic gravity in a megathrust faulting setting.

Geophysical Journal International↗

Do life history traits influence patterns of maternal immune elements in New World blackbirds (Icteridae)?

Avian immunology developed originally by investigating domesticated poultry species (Galliformes), but in recent decades eco-immunological studies of wild bird species have revealed that avian immune systems are more diverse than initially assumed. This study compares six immunological elements in eggs of six species within the same family, the New World blackbirds (Icteridae),whose members differ most notably in two life history parameters, brood parasitism and body size. We measured the maternal immune investment of passive immune components in both yolk and albumen: lysozyme, ovotransferrin, and immunoglobulins (Igs), and LPS-specific Igs. We predicted that brood parasites would have higher levels of immune activity for both innate and adaptive immunity compared with non-brood parasites, and that increased body size could increase microbial exposure of larger animals, resulting in an increase in some adaptive immune responses, such as LPS-specific Igs. We found that brood parasites had significantly higher levels of Igs and lysozyme levels in albumen, but significantly lower levels of Igs in yolk compared with non-brood parasites. Igs in yolk scaled according to body size, with the smallest organisms (the brood parasites) having the lowest levels, and the largest organism (common grackle) having the highest. Our results confirm the findings of other studies of comparative immunity among species in a single taxon that (1) similarities in immune investment cannot be assumed among closely related species and (2) single measures of immune defense cannot be assumed to be indicators of a species’ overall immune strategy, as life history traits can differentially affect immune responses.

North Dakota, Puerto Rico, South Dakota↗

Rare earth element evidence for the petrogenesis of the banded series of the Stillwater Complex, Montana, and its anorthosites

A rare earth element (REE) study was made by isotope-dilution mass spectrometry of plagioclase separates from a variety of cumulates stratigraphically spanning the Banded series of the Stillwater Complex, Montana. Evaluation of parent liquid REE patterns, calculated on the basis of published plagioclase-liquid partition coefficients, shows that the range of REE ratios is too large to be attributable to fractionation of a single magma type. At least two different parental melts were present throughout the Banded series. This finding supports hypotheses of previous workers that the Stillwater Complex formed from two different parent magma types, designated the anorthositic- or A-type liquid and the ultramafic- or U-type liquid. On the basis of our data, one melt has a REE pattern with a distinctive shallow slope and is represented by samples from the thick, massive Anorthosite zones I and II (AN I and AN II) of the Middle Banded series. Although samples from AN I and AN II are separated by as much as 1400 m stratigraphically, they have remarkably similar calculated parent liquid characteristics, with (Ce/Sm)n = 1.7–1.9, (Nd/Sm)n = 1.3–1.4 and (Ce/Yb)n = 2.9–4.6 (where n denotes chondrite-normalized). These calculated liquids are probably close to representing A-type magma. In addition, plagioclase-bronzite cumulates from Norite zones I and II (N I and N II), although thought to be U-type cumulates, contain plagioclase that has A-type REE characteristics, implying that A-type magmas were injected into the magma chamber during formation of those zones. In contrast, calculated parent liquids of cumulus augite-bearing rocks have REE patterns that display distinctly steeper slopes than the A-type REE pattern. The extreme is the calculated parent liquid of a plagioclase-bronzite-augite cumulated with (Ce/Sm)n = 2.9, (Nd/Sm)n = 1.7, and (Ce/Yb)n = 10.1.

Journal of Petrology↗

Spectral element modelling of fault-plane reflections arising from fluid pressure distributions

The presence of fault-plane reflections in seismic images, besides indicating the locations of faults, offers a possible source of information on the properties of these poorly understood zones. To better understand the physical mechanism giving rise to fault-plane reflections in compacting sedimentary basins, we numerically model the full elastic wavefield via the spectral element method (SEM) for several different fault models. Using well log data from the South Eugene Island field, offshore Louisiana, we derive empirical relationships between the elastic parameters (e.g. P-wave velocity and density) and the effective-stress along both normal compaction and unloading paths. These empirical relationships guide the numerical modelling and allow the investigation of how differences in fluid pressure modify the elastic wavefield. We choose to simulate the elastic wave equation via SEM since irregular model geometries can be accommodated and slip boundary conditions at an interface, such as a fault or fracture, are implemented naturally. The method we employ for including a slip interface retains the desirable qualities of SEM in that it is explicit in time and, therefore, does not require the inversion of a large matrix. We performa complete numerical study by forward modelling seismic shot gathers over a faulted earth model using SEM followed by seismic processing of the simulated data. With this procedure, we construct post-stack time-migrated images of the kind that are routinely interpreted in the seismic exploration industry. We dip filter the seismic images to highlight the fault-plane reflections prior to making amplitude maps along the fault plane. With these amplitude maps, we compare the reflectivity from the different fault models to diagnose which physical mechanism contributes most to observed fault reflectivity. To lend physical meaning to the properties of a locally weak fault zone characterized as a slip interface, we propose an equivalent-layer model under the assumption of weak scattering. This allows us to use the empirical relationships between density, velocity and effective stress from the South Eugene Island field to relate a slip interface to an amount of excess pore-pressure in a fault zone. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗