|
Diciembre 2007
Percolation conditions in fractured hard rocks: A numerical
approach using the three-dimensional binary fractal fracture
network (3D-BFFN) model
Authors: Shinji Nakaya, Kiminori Nakamura et al
Link: Click Aquí
Abstract
We numerically investigate fracture connectivity and percolation
conditions in fractured hard rocks using a three-dimensional
binary fractal fracture network (3D-BFFN) model based on three
fractal geometric parameters: the fractal dimensions (D2) of
the spatial distribution of fractures, the exponent of the power-law
cumulative fracture length distribution (a), and the maximum
fracture length (l max) normalized by the domain length (L),
l max/L. Numerical results clarify that the percolation threshold
in 3D-BFFN models is strongly controlled by fractal geometric
parameters and is independent of any anisotropy in the orientations
T. In addition, when a < 1.8 and
|
l
max/L < 1.0, percolation seldom occurs independently of D2
and T. In the current study the analytical solution of percolation
probability (P) is presented as a function of the three fractal
parameters within the 3D-BFFN model. Application of the 3D-BFFN
model to seismogenic fractures determined from the earthquake
catalogue in an offshore volcanic region between Miyake-jima Island
(MI) and Kozu-shima Island (KI) off the Izu Peninsula, Japan,
suggests that P is mainly affected by the error involved in determining
a during actual surveys. Otherwise, P provides a useful index
for determining whether a three-dimensional domain percolates
in fracture networks in fractured hard rocks. The basis of this
approach is the observation from fracture network connections
that domains with P > 0.55 are percolated domains. The zone
of percolation within seismogenic fracture networks between MI
and KI reveals that the networks formed from seismic-swarm-related
seismogenic fractures over a 7-week period related to the intrusion
of a dyke, inferred previously from seismicity and deformation
data. |
Diciembre 2007
Asperity generating upper crustal sources revealed by b value
and isostatic residual anomaly grids in the area of Antofagasta,
Chile
Authors: M. Sobiesiak, U. Meyer et al
Link: Click Aquí
Abstract
In our study we show that the locations of largest coseismic
slip (asperities) on the fault plane of the M w = 8.0 1995 Antofagasta
earthquake in Northern Chile can be mapped by the spatial distribution
of the seismic b value obtained from the aftershock sequence
of the megathrust earthquake. These areas of high seismic moment
release and concurrent high-b values are congruent with anomalies
of the gravity isostatic residual (IR) field in the Antofagasta
region. They are superimposed on the seismogenic part of the
north Chilean subduction zone where the
|
strongest coupling of
the upper and lower plate is expected. The IR anomalies are interpreted
to be caused by large Jurassic-Early Cretaceous batholiths which
intruded into the upper crust. The observed positive correlations
between high seismic b values, IR anomalies, and geologic structures
enable us to propose a mechanical model for the generation of
the asperities in the Antofagasta region. We suggest that the
batholiths in conjunction with buoyant forces acting on the subducted
slab of the Nazca plate are responsible for locking the interface
where the asperities are located. This implies long-term conditions
for the existence of the asperity generating tectonic situation.
Concequently, the asperities around Antofagasta could be stationary
features, at least for several seismic cycles. Hence we propose
that the IR anomalies along the north Chilean convergent margin
can be used as an indicator for high moment release and slip in
future large earthquakes. |
Diciembre
2007
Characterizing earthquake recurrence parameters for offshore
faults in the low-strain, compressional Kapiti-Manawatu Fault
System, New Zealand
Authors: Scott D. Nodder, Geoffroy Lamarche et al
Link: Click Aquí
Abstract
Seafloor fault scarps and near-surface deformation of late Quaternary
seismic reflectors occur along the eastern margin of the Wanganui
Basin, 200 km behind the active Hikurangi subduction front,
southern North Island, New Zealand. The offshore scarps are
associated with the low-strain, compressional Kapiti-Manawatu
Fault System (KMFS), which comprises high-angle (>60°)
reactivated reverse and normal faults oriented NE-SW, highly
oblique to the coast. Seafloor scarps range from <10 to 50
km in length with vertical seafloor offsets of 2 to 30 m. The
longest structure is the Mascarin Fault, with
|
maximum
late Quaternary vertical slip rates of 3 mm a-1 (where a is years).
Other faults in the KMFS have typical rates of <1 mm a-1, comparable
to long-term estimates. Three zones of recent deformation are
identified: faults in the north and south of the KMFS are characterized
by high scarps (>10 m high) and moderate to long fault seafloor
rupture lengths, and those in central parts of the fault system
are characterized by low scarps (<5 m high), variable slip
rates, and short to moderate fault lengths. Empirical equations
indicate that KMFS faults may generate earthquakes with moment
magnitudes (Mw) of 5.7-7.5 (mean 6.9 ± 0.3, ±1 standard
deviation, for sources with Mw = 6.5). Estimated recurrence intervals
that are generally >10,000 a, suggest that the seismic hazard
of the Kapiti-Manawatu region is relatively low. Incorporation
of these new geological data, however, is likely to increase slightly
the expected seismic hazard in southern North Island. The method
of determining the earthquake recurrence parameters of offshore
faults has potentially wider applications elsewhere. |
Diciembre 2007
Fractal characterization of fracture networks: An improved
box-counting technique
Authors: Ankur Roy, Edmund Perfect et al
Link: Click Aquí
Abstract
Box counting is widely used for characterizing fracture networks
as fractals and estimating their fractal dimensions (D). If
this analysis yields a power law distribution given by N r -D
, where N is the number of boxes containing one or more fractures
and r is the box size, then the network is considered to be
fractal. However, researchers are divided in their opinion about
which is the best box-counting algorithm to use, or whether
fracture networks are indeed fractals. A synthetic fractal fracture
network with a known D value was used to develop a new algorithm
for the box-counting method that returns improved estimates
of D. The method is based on identifying the lower limit of
fractal behavior (r cutoff) using the condition ds/dr ? 0, where
s is the standard deviation from a linear regression equation
fitted to log(N) versus log(r) with data for r < r cutoff
sequentially excluded. A set of 7 nested fracture maps from
the Hornelen Basin, Norway was used to test the improved method
and demonstrate its accuracy for natural patterns. We also reanalyzed
a suite of 17 fracture trace maps that had previously been evaluated
for their fractal nature. The improved estimates of D for these
maps ranged from 1.56 ± 0.02 to 1.79 ± 0.02, and
were much greater than the original estimates. These higher
D values imply a greater degree of fracture connectivity and
thus increased propensity for fracture flow and the transport
of miscible or immiscible chemicals.
|
Diciembre
2007
Microearthquake streaks and seismicity triggered by slow earthquakes
on the mobile south flank of Kilauea Volcano, Hawai'i
Authors: Cecily J. Wolfe, Benjamin A. Brooks et al
Link: Click Aquí
Abstract
We perform waveform cross correlation and high precision relocation
of both background seismicity and seismicity triggered by periodic
slow earthquakes at Kilauea Volcano's mobile south flank. We
demonstrate that the triggered seismicity dominantly occurs
on several preexisting fault zones at the Hilina region. Regardless
of the velocity model employed, the relocated earthquake epicenters
and triggered seismicity localize onto distinct fault zones
that form streaks aligned with the slow earthquake surface displacements
determined from GPS. Due to the unknown effects of velocity
heterogeneity and nonideal station coverage, our relocation
analyses cannot distinguish whether some of these fault zones
occur within the volcanic crust at shallow depths or whether
all occur on the decollement between the volcano and preexisting
oceanic crust at depths of ~8 km. Nonetheless, these Hilina
fault zones consistently respond to stress perturbations from
nearby slow earthquakes.
|
Diciembre
2007
Repeating earthquake finite source models: Strong asperities
revealed on the San Andreas Fault
Authors: Douglas Dreger, Robert M. Nadeau et al
Link: Click Aquí
Abstract
We investigate the rupture process of a sequence of repeating
Mw 2.1 earthquakes on the San Andreas Fault in Parkfield spanning
the occurrence of the September 28, 2004 mainshock by inverting
seismic moment rate functions obtained from empirical Green's
function deconvolution. The results show that these events have
extremely concentrated slip patches with radii on the order
of 10-20 m, with peak slip between 8.4 and 11.4 cm. The rupture
speed and rise time are consistent with values of larger earthquakes.
The spatial distribution of stress drop for the events shows
low average values 2.5-5.6 MPa and very large peak values of
66.7-93.9 MPa. The results show that strong asperities can exist
at small scales on an otherwise weak fault, and helps reconcile
differences between traditional spectra-based and tectonic loading
methods for determining the stress drop of small repeating earthquakes.
|
Diciembre
2007
Regional and teleseismic double-difference earthquake relocation
using waveform cross-correlation and global bulletin data
Authors: Felix Waldhauser and David Schaff
Link: Click Aquí
Abstract
We have developed a double-difference algorithm to relocate
earthquakes recorded at global seismic networks, using differential
arrival times for first and later arriving regional and global
phases to invert for the vectors connecting the hypocenters.
Differential times are formed from global seismic bulletins
and are accurately measured on similar seismograms by time domain
waveform cross correlation. We evaluate the performance of this
spherical, multiphase double-difference algorithm using three-dimensional
regional-scale synthetic data and two sets of earthquake data
in different tectonic settings. The first includes 3783 intermediate
depth earthquakes that occurred between 1964 and 2000 in the
subducting Nazca plate beneath northern Chile, where the relocated
seismicity confirms a narrowly spaced double seismic zone previously
imaged with temporary local seismic data. Residual statistics
and comparison with accurately known locations indicate mean
relative location errors at the 90% confidence level of 2.4
km laterally and 1.8 km vertically. Later events typically constrained
by cross-correlation data have errors of 1.6 km laterally and
1.4 km vertically. The second data set includes 75 crustal earthquakes
in the 1999 Izmit and Düzce, Turkey, aftershock sequences,
where the double-difference solutions image orientation and
dip of individual fault segments that are consistent with focal
mechanisms and near-surface information. Fault complexity likely
causes a low level of waveform similarity in this aftershock
sequence and thus generates fewer correlated events compared
to the Chile earthquakes. Differences between the double-difference
locations and corresponding locations in global seismicity catalogs
(Earthquake Data Report, EDR; International Seismological Centre,
ISC; Engdahl-Hilst-Buland, EHB) are typically greater than 10
km. We evaluate the potential of cross-correlation and double-difference
methods to improve hypocenter locations on a global scale.
|
Diciembre 2007
Gas hydrate concentration and characteristics within Hydrate
Ridge inferred from multicomponent seismic reflection data
Authors: Dhananjay Kumar, Mrinal K. Sen et al
Link: Click Aquí
Abstract
A seismic experiment composed of streamer and ocean bottom seismometer
(OBS) surveys was conducted in the summer of 2002 at southern
Hydrate Ridge, offshore Oregon, to map the gas hydrate distribution
within the hydrate stability zone. Gas hydrate concentrations
within the reservoir can be estimated with P wave velocity (V
p ); however, we can further constrain gas hydrate concentrations
using S wave velocity (V s ), and use V s through its relationship
to V p (V p /V s ) to reveal additional details such as gas
hydrate form within the matrix (i.e., hydrate cements the grains,
becomes part of the matrix frame or floats in pore
|
space). Both V p and
V s can be derived simultaneously by inverting multicomponent
seismic data. In this study, we use OBS data to estimate seismic
velocities where both gas hydrate and free gas are present in
the shallow sediments. Once V p and V s are estimated, they are
simultaneously matched with modeled velocities to estimate the
gas hydrate concentration. We model V p using an equation based
on a modification of Wood's equation that incorporates an appropriate
rock physics model and V s using an empirical relation. The gas
hydrate concentration is estimated to be up to 7% of the rock
volume, or 12% of the pore space. However, V p and V s do not
always fit the model simultaneously. V p can vary substantially
more than V s . Thus we conclude that a model, in which higher
concentrations of hydrate do not affect shear stiffness, is more
appropriate. Results suggest gas hydrates form within the pore
space of the sediments and become part of the rock framework in
our survey area. |
Diciembre
2007
Spatial mapping of the b value at Mount Etna, Italy, using
earthquake data recorded from 1999 to 2005
Authors: M. Murru, R. Console et al
Link: Click Aquí
Abstract
The spatial pattern of the b value of the frequency-magnitude
relation has been analyzed using gridding techniques beneath
Mount Etna, Italy. A regional data set of 2900 events with M
d (duration magnitude) =1.5 up to 15 km depth occurring between
August 1999 and December 2005 has been used. Two regions with
an abnormally high b value have been found, one centered beneath
the southern part of the Valle del Bove, above the 6 km below
sea level (bsl) deep basement, and the other beneath the summit
region 2 km bsl
|
east
of the Central Craters. We can infer that these high b value anomalies
are regions of increased crack density, and/or high pore pressure,
related to the presence of nearby magma storage. This interpretation
is supported by all the available geophysical evidence, such as
tomographic studies and geodetic deformation measurements. The
data set has also been subdivided into five periods, corresponding
to different phases of volcanic activity: 2001 preeruption, 2001
eruptive, 2002-2003 preeruption, 2002-2003 eruptive, and 2002-2003
posteruption. The minimum magnitude of completeness, M c , and
the b value were computed for each period. A volume of anomalously
high b values can be observed in each of these periods (except
for the 2002-2003 preeruption interval). This approach has allowed
the detection of the transient presence of magmatic intrusions
during the various periods evaluated. |
Diciembre 2007
Repeating earthquakes and seismic potential along the northern
Longitudinal Valley fault of eastern Taiwan
Authors: Ruey-Juin Rau, Kate Huihsuan Chen et al
Link: Click Aquí
Abstract
The northern Longitudinal Valley fault in eastern Taiwan creeps
at the surface with a small rate of ~1 cm/yr but slips in large
earthquakes. To
|
improve seismic hazard assessment, it is important
to comprehend the slip deficit rate distribution at depth. We
discovered 25 M L 2.1-4.6 repeating earthquakes in this area
and inverted GPS measurements for producing an image of the
along-strike spatial distribution of deep fault slip rates.
The repeating events are located at the depths of 10-22 km with
24.9-77.5 mm/yr slip rates, which are comparable with the GPS-derived
slip rates of 47.5 ± 5.8 mm/yr at similar depth ranges.
Based on distribution of GPS-derived slip deficits, since 1951,
the northern Longitudinal Valley fault has become capable of
releasing stored strain in a future M w = 7.3 earthquake.
|
|
Diciembre 2007
Interacciones entre
terremotos
y actividad volcánica
Autores: Nicolas Lemarchand y Jean-Robert
Grasso
Link: Click Aquí
Abstract
Utilizando el catálogo mundial 1973-2005 para sismicidad
M >= 4.8 y erupciones volcánicas con VEI >= 0,
encontramos un incremento del inicio de la actividad volcánica
el día del terremoto. Este resultado emerge del estudio
de las series de tiempo de razones diarias de erupción
relativas al tiempo to del terremoto, respecto de la sismicidad
global incluída en el catálogo. La correlación
es fuerte para los pares terremoto-actividad volcánica
(TV), siempre que el volcán se ubique dentro de una distancia
al epicentro menor o igual a diez veces la longitud de ruptura.
Los resultados muestran que para M >= 4.8 y VEI >= 0,
los pares TV son tan importantes como los procesos de interacción
con M >= 7 and VEI >= 2-5 previamente reportados. La agrupación
temporal para los pares TV no está necesaria-mente referida
a to. La correlación se distribuye en torno de to
con una dispersión de 6 a 10 días (antes y después
del terremoto) y sigue una ley de potencias . Estos
patrones, los cuales son similares a las leyes directas e inversas
de Omori observadas para terremotos tectónicos, corresponden
a la primera evidencia que muestra que las erupciones volcánicas
son contempo-ráneas al rompimiento de la corteza terrestre.
Las agrupaciones espaciotemporales que observamos nos permiten
concluir lo siguiente:
i) Es incorrecto que el disparo del terremoto es el único
proceso que influye en el par TV
ii) Se aportan pruebas que nos indican un acoplamiento regional
tectónico en el daño de la corteza terrestre.
|
Diciembre 2007

Interactions between earthquakes and
volcano activity
Authors: Nicolas Lemarchand and Jean-Robert Grasso
Link: Click Aquí
Abstract
Using the 1973-2005 worldwide catalogues for M >= 4.8 seismicity
and VEI >= 0 volcano eruptions, we find a significant, when
tested against catalogue randomizations, increase of eruption
onsets on the earthquake day. This result emerges from stacking
time series of daily eruption rates relatively to earthquake
time to, over the whole seismicity catalogue. It is stronger
for earthquake-volcano pairs for which the volcano is within
ten rupture size from the epicenter. These results show that
M >= 4.8 - VEI >= 0 earthquake-volcano pairs are as important
for interaction processes as the M >= 7 and VEI >= 2-5
pairs previously reported to interact. The clustering in time
for earthquake-eruption pairs is not bounded to t 0. It remains
above the background noise 6-10 days before and after t 0, and
follows a power law distribution. These patterns, which are
similar to the direct and inverse Omori's laws observed for
tectonic earthquakes, are the first evidence for the volcano
eruptions to be contemporary of a stochastic brittle damage
in the earth crust. The clustering we observe in time and space
(i) rejects the earthquake triggering as the single process
that drive earthquake-volcano interactions; (ii) supports a
regional tectonic coupling in the earth crust damage.
|
Diciembre 2007
On the preparatory processes of the M6.6 earthquake of June
17th, 2000, in Iceland
Authors: M. Bonafede, C. Ferrari et al
Link: Click Aquí
Abstract
A model is proposed to explain the spatial distribution of foreshocks
of the June 17th 2000, M s 6.6 earthquake in the South Iceland
Seismic Zone (SISZ) and the high stress drop of the
|
mainshock. Fluids of
magmatic origin, ascending at near-lithostatic pressure through
a low permeability layer perturb the regional stress field, inhibiting
fluid flow laterally, where a high strength asperity is left.
The asperity is modeled as elastic, embedded within a medium with
low effective rigidity. Regional stresses due to tectonic motions
are perturbed by the presence of the asperity, enhancing the production
of hydrofractures and foreshocks in the NW and SE quadrants and
increasing considerably the shear stress within the asperity,
leading to the NS striking mainshock. |
Diciembre
2007
Search for past life on Mars: Physical and chemical characterization
of minerals of biotic and abiotic origin: 2. Aragonite
Authors: F. Stalport, P. Coll et al
Link: Click Aquí
Abstract
One of the major objectives of the future Martian surface probes
will be to reveal a past or present biological activity. We
propose that biominerals could have recorded such an activity
at Mars, and thus could be interesting targets for these missions.
Therefore, we try to find a method
|
capable
to discriminate biominerals from their geochemical counterparts.
With this aim, various terrestrial aragonites of biotic and abiotic
origins were studied as reference minerals, because they could
have also been produced at Mars. Their thermal properties were
studied with differential thermal analysis, and then compared.
The results show that biotic aragonites thermally decompose at
temperatures at least 20°C lower than the temperatures of
decomposition of abiotic aragonites. Therefore, the temperatures
of thermal degradation of such biominerals could be a relevant
parameter to find a past biological activity at Mars, and differential
thermal analysis could be useful for situ astrobiological exploration
of Mars. |
Diciembre 2007
Predictive model for permeability reduction by small wetting
phase saturations
Authors: Siyavash Motealleh and Steven L. Bryant
Link: Click Aquí
Abstract
Laboratory experiments show that the rate of decrease in gas
phase permeability as a function of water saturation can be
large at small saturations. Drying out the water saturation
during gas injection has been observed to reduce absolute permeability.
This paper proposes a quantitative grain-scale explanation for
these observations. We also show why the effect is
|
magnified in rocks having
lower porosity. To compute phase geometry and permeability, we
use a physically representative network model. The network is
extracted from a model rock, built from a dense random packing
of spheres modified geometrically to simulate rock-forming processes.
At small saturations the wetting phase exists largely in the form
of pendular rings held at grain contacts. The rings decrease the
void area available for flowing nonwetting phase. Because the
hydraulic conductance of the throat varies with the square of
the void area (other factors being equal), the effect on permeability
is disproportionate to the volume occupied by the rings. The same
approach quantifies the reduction in permeability by salt precipitation
during drying. |
Diciembre
2007
Magnetic characteristics of synthetic pseudo-single-domain
and multi-domain greigite (Fe3S4)
Authors: Liao Chang, Andrew P. Roberts et al
Link: Click Aquí
Abstract
We report the magnetic properties of pure synthetic pseudo-single-domain
(PSD) and multi-domain (MD) greigite and the grain size dependence
of the magnetic properties of greigite for the first time. The
dominantly PSD-like and MD-like behavior are demonstrated by
hysteresis, first-order reversal curve diagrams, low -
|
temperature
cycling (LTC) of room temperature saturation isothermal remanent
magnetization (SIRM) and low-temperature SIRM warming curves.
Variations in a range of magnetic properties clearly correlate
with grain size. Characteristic PSD/MD behavior is preserved at
low temperatures, which, coupled with the small decrease in remanence
during warming, rule out the presence of substantial superparamagnetic
behavior in the studied samples. LTC-SIRM measurements indicate
a continuous demagnetization of remanence during cooling. Knowledge
of this expanded range of magnetic properties of greigite should
be widely useful in environmental magnetic and paleomagnetic studies.
|
Noviembre
2007
Predictibilidad: Recientes investigaciones a partir de la Teoría
de la Información
Authors: Timothy DelSole and Michael K. Tippett
Link: Click Aquí
Abstract
Este paper resume un marco conceptual para la investigación
de la predictibilidad a partir de la Teoría de la Información.
Este marco conecta y unifica una gran variedad de métodos
estadísticos utilizados tradicionalmente en análisis
de predictibilidad, incluyendo regresión lineal, análisis
de correlación canónica, descomposición
de valores singulares, análisis de discriminante y asimilación
de datos. La base de este marco es un procedimiento llamado
Análisis de Componentes Predecibles (PrCA). El PrCA realiza
una descomposión óptima de las variables que se
pueden utilizar para realizar predicitibilidad, así como
el Análisis de Componentes Principales descompone óptimamente
las variables a partir de su varianza. Para distribuciones normales
se obtiene los mismos componentes predecibles, ya sea que se
optimice la información predecible, la dispersión
de la entropía relativa, la información mutua,
el error de Mahalanobis, la razón entre la señal
promedio y el ruido, el error cuadrático medio normalizado
o la anomalía de correlación. Para distribuciones
normales conectadas, el PrCA es equivalente al análisis
de correlación canónico entre pronóstico
y observaciones. El operador de regresión que mapea de
observaciones a pronósticos juega un importante rol en
este marco, dejando los vectores singulares del operador como
los componentes predecibles, al tiempo que los valores singulares
son las correlaciones canónicas. Esta última correspondencia
ocurre sólo si los vectores singulares son calculados
mediante la norma de Mahalanobis, un resultado que aclara el
rol de la norma en la predictibilidad. En modelos estocásticos
lineales el hecho que impone la minimización de la predicitibilidad
es lo único que define la dinámica "blanca"
del operador normal. Esta condición para predictibilidad
mínima es invariante frente a las transformaciones lineales
y es equivalente a un balance detallado. El marco también
inspira algunos nuevos acercamientos para contabilizar las deficiencias
de los modelos de pronóstico y estimar distribuciones
a partir de un muestreo finito.
|
Noviembre 2007
Predictability: Recent
insight from information theory
Authors: Timothy DelSole and Michael K.
Tippett
Link: Click Here
Abstract
This paper summarizes a framework for investigating
predictability based on information theory. This framework connects
and unifies a wide variety of statistical methods traditionally
used in predictability analysis, including linear regression,
canonical correlation analysis, singular value decomposition,
discriminant analysis, and data assimilation. Central to this
framework is a procedure called predictable component analysis
(PrCA). PrCA optimally decomposes variables by predictability,
just as principal component analysis optimally decomposes variables
by variance. For normal distributions the same predictable components
are obtained whether one optimizes predictive information, the
dispersion part of relative entropy, mutual information, Mahalanobis
error, average signal to noise ratio, normalized mean square
error, or anomaly correlation. For joint normal distributions,
PrCA is equivalent to canonical correlation analysis between
forecast and observations. The regression operator that maps
observations to forecasts plays an important role in this framework,
with the left singular vectors of this operator being the predictable
components and the singular values being the canonical correlations.
This correspondence between predictable components and singular
vectors occurs only if the singular vectors are computed
using Mahalanobis norms, a result that sheds light on the
role of norms in predictability. In linear stochastic models
the forcing that minimizes predictability is the one that renders
the "whitened" dynamical operator normal. This condition
for minimum predictability is invariant to linear transformation
and is equivalent to detailed balance. The framework also inspires
some new approaches to accounting for deficiencies of forecast
models and estimating distributions from finite samples.
|
Noviembre
2007
Transition warming and cooling remanences in magnetite
Author: David J. Dunlop
Link: Click Here
Abstract
Insight into the size and morphology of assemblages of magnetite
particles can be gained by comparing temperature variations
of remanence or susceptibility after zero-field cooling (ZFC)
and after field cooling (FC) through the Verwey transition around
T v = 120 K. At 10 K a sample is demagnetized following ZFC,
but in the FC initial state before warming the sample has a
transition cooling remanence (TrCRM) acquired in crossing T
v. The matching transition warming remanence (TrWRM) acquired
as a result of heating a demagnetized sample from low temperature
across T v is often called inverse thermoremanent magnetization
(ITRM). In TrCRM experiments, initially demagnetized samples
were cooled in a 2 mT field from 300 K to 10 K. Magnetization
M was measured at 1 K to 5 K
|
intervals,
the highest-resolution data being taken between 140 K and 90 K.
The field was zeroed at 10 K, and the TrCRM was monitored during
zero-field warming back to 300 K. The properties of TrCRMs are
generally similar to those of TrWRMs produced by heating a ZFC
sample in a 2 mT field from 10 K. In 10 of 12 samples (grain sizes
from 0.6 to 135 µm), M of monoclinic magnetite produced
by field cooling through T v exactly equals M of cubic magnetite
produced by field warming through T v, even though the ultimate
TrCRM and TrWRM values when H ? 0 are entirely different. Mirror-image
symmetry was observed between in-field warming curves tracking
the acquisition of TrWRM and zero-field warming curves of TrCRM
between 10 and 300 K. The symmetry, with increases in the field-on
M curves mirroring decreases in the field-off M r curves, was
almost perfect from 10 to 110 K. Approximate symmetry was also
observed between in-field cooling curves tracking TrCRM production
and zero-field cooling curves of TrWRM between 300 K and T v.
Detailed study of the properties and mechanism(s) of transition
remanences will help clarify why the ZFC/FC method is diagnostic
in some instances and not in others. |
Noviembre 2007
Magnetofossil spike during the Paleocene-Eocene thermal maximum:
Ferromagnetic resonance, rock magnetic, and electron microscopy
evidence from Ancora, New Jersey, United States
Authors: Robert E. Kopp, Timothy D. Raub et al
Link: Click Here
Abstract
Previous workers identified a magnetically anomalous clay layer
deposited on the northern United States Atlantic Coastal Plain
during the Paleocene-Eocene thermal maximum (PETM). The finding
inspired the highly controversial hypothesis that a cometary
impact triggered the PETM. Here we present ferromagnetic resonance
(FMR), isothermal and anhysteretic remanent magnetization, first-order
reversal curve, and
|
transmission electron
microscopy analyses of late Paleocene and early Eocene sediments
in drill core from Ancora, New Jersey. A novel paleogeographic
analysis applying a recent paleomagnetic pole from the Faeroe
Islands indicates that New Jersey during the initial Eocene had
a ~6°-9° lower paleolatitude (~27.3° for Ancora) and
a more zonal shoreline trace than in conventional reconstructions.
Our investigations of the PETM clay from Ancora reveal abundant
magnetite nanoparticles bearing signature traits of crystals produced
by magnetotactic bacteria. This result, the first identification
of ancient biogenic magnetite using FMR, argues that the anomalous
magnetic properties of the PETM sediments are not produced by
an impact. They instead reflect environmental changes along the
eastern margin of North America during the PETM that led to enhanced
production and/or preservation of magnetofossils. |
|
|
Noviembre
2007
A three-dimensional crustal seismic velocity model for southern
California from a composite event method
Authors: Guoqing Lin, Peter M. Shearer et al
Link: Click Here
Abstract
We present a new crustal seismic velocity model for southern
California derived from P and S arrival times from local earthquakes
and explosions. To reduce the volume of data and ensure a more
uniform source distribution, we compute "composite event"
picks for 2597 distributed master events that include pick information
for other events within spheres of 2 km radius. The approach
reduces random picking error and maximizes the number of S wave
picks. To constrain absolute event locations and shallow
|
velocity structure, we also use times from controlled sources,
including both refraction shots and quarries. We implement the
SIMULPS tomography algorithm to obtain three-dimensional (3-D)
V p and V p /V s structure and hypocenter locations of the composite
events. Our new velocity model in general agrees with previous
studies, resolving low-velocity features at shallow depths in
the basins and some high-velocity features in the midcrust. Using
our velocity model and 3-D ray tracing, we relocate about 450,000
earthquakes from 1981 to 2005. We observe a weak correlation between
seismic velocities and earthquake occurrence, with shallow earthquakes
mostly occurring in high P velocity regions and midcrustal earthquakes
occurring in low P velocity regions. In addition, most seismicity
occurs in regions with relatively low V p /V s ratios, although
aftershock sequences following large earthquakes are often an
exception to this pattern. |
Noviembre 2007
Energy-Constrained Recharge, Assimilation, and Fractional Crystallization
(EC-RA?FC): A Visual Basic computer code for calculating trace
element and isotope variations of open-system magmatic systems
Authors: Wendy A. Bohrson, Frank J. Spera et al
Link: Click Here
Abstract
Volcanic and plutonic rocks provide abundant evidence for complex
processes that occur in magma storage and transport systems.
The fingerprint of these processes, which include fractional
crystallization, assimilation, and magma recharge, is captured
in petrologic and geochemical characteristics of suites of cogenetic
rocks. Quantitatively evaluating the relative contributions
of each process requires integration of mass, species, and energy
constraints, applied in a self-consistent way. The energy-constrained
model Energy-Constrained Recharge, Assimilation, and Fractional
Crystallization (EC-Ra FC) tracks the trace element and isotopic
evolution of a magmatic system (melt + solids) undergoing simultaneous
fractional crystallization,
|
recharge, and assimilation. Mass, thermal, and
compositional (trace element and isotope) output is provided
for melt in the magma body, cumulates, enclaves, and anatectic
(i.e., country rock) melt. Theory of the EC computational method
has been presented by Spera and Bohrson (2001, 2002, 2004),
and applications to natural systems have been elucidated by
Bohrson and Spera (2001, 2003) and Fowler et al. (2004). The
purpose of this contribution is to make the final version of
the EC-RA?FC computer code available and to provide instructions
for code implementation, description of input and output parameters,
and estimates of typical values for some input parameters. A
brief discussion highlights measures by which the user may evaluate
the quality of the output and also provides some guidelines
for implementing nonlinear productivity functions. The EC-RA?FC
computer code is written in Visual Basic, the programming language
of Excel. The code therefore launches in Excel and is compatible
with both PC and MAC platforms. The code is available on the
authors' Web sites http://magma.geol.ucsb.edu/and http://www.geology.cwu.edu/ecrafc)
as well as in the auxiliary material.
|
Noviembre
2007
Isotopic fractionation of methane and ethane hydrates between
gas and hydrate phases
Authors: Akihiro Hachikubo, Tomoko Kosaka et al
Link: Click Here
Abstract
Isotopic fractionation of carbon and hydrogen in methane and
ethane during the formation of gas hydrates was investigated.
The gas hydrate samples were experimentally prepared in a
|
pressure
cell and isotopic compositions of both residual and hydrate-bound
gases were measured. dD of hydrate-bound molecules of methane
and ethane hydrates was several per mil lower than that of residual
gas molecules in the formation processes, while there was no difference
in the case of d 13C. These isotopic differences in dD are enough
small for discussing the source types of hydrate-bound gases using
the d 13C-dD diagram of Whiticar et al. [1986]. These results
may provide useful insight into the formation process of gas hydrates.
|
Noviembre 2007
Powering Mercury's dynamo
Author: J.-P. Williams, O. Aharonson et al
Link: Click Here
Abstract
The presence of the global magnetic field of Mercury has implications
for the interior structure of the planet and its thermal evolution.
We use a thermal evolution model to explore the conditions under
which excess entropy is available to drive a convective dynamo.
The current state of the core is
|
strongly affected by
its sulfur concentration and the viscosity of the overlying mantle.
A present-day dynamo is difficult to achieve. The minimum rate
of entropy production required to drive a dynamo is attained in
only the most optimistic models, and requires present-day mantle
convection. An additional entropy source such as the addition
of a radiogenic heat source in the core increases the probability
of a present-day dynamo. Given the uncertainty, more specific
characterization of the planet's interior and magnetic field is
required to alleviate ambiguities in the original Mariner 10 observations.
|
Noviembre
2007
Multistep ahead streamflow forecasting: Role of calibration
data in conceptual and neural network modeling
Authors: Elena Toth and Armando Brath
Link: Click Here
Abstract
When choosing the rainfall-runoff modeling approach to be integrated
in a river flow forecasting system, two crucial issues are the
minimum data requirement for calibration purposes and the reliability
of the predictions over different time horizons (lead-times).
The paper presents an investigation of the real-time
|
forecasting
ability of a conceptual and a neural network model, comparing
the performances obtainable for increasing lead-times and analyzing
the influence of the amount of the calibration data over two real-data
case studies. Neural networks proved to be an excellent tool for
the real-time rainfall-runoff simulation of continuous periods
(including low, average and peak flows), provided that an extensive
set of hydro-meteorological data is available for calibration
purposes. On the other hand, the comparison highlights that a
conceptual formulation may allow a significant forecasting improvement
in comparison with the data-driven approach when focusing on the
prediction of flood events and especially in case of a limited
availability of calibration data. |
Noviembre 2007
Stochastic simulation and spatial estimation with multiple
data types using artificial neural networks
Authors: Lance E. Besaw and Donna M. Rizzo
Link: Click Here
Abstract
A novel data-driven artificial neural network (ANN) that quantitatively
combines large numbers of multiple types of soft data is presented
for performing stochastic simulation and/or spatial estimation.
A counterpropagation ANN is extended with a radial basis function
to estimate parameter fields that reproduce the spatial structure
exhibited in autocorrelated parameters. Applications involve
using three geophysical properties measured on a slab of Berea
sandstone and the delineation of landfill leachate
|
at a site in the Netherlands
using electrical formation conductivity as our primary variable
and six types of secondary data (e.g., hydrochemistry, archaea,
and bacteria). The ANN estimation fields are statistically similar
to geostatistical methods (indicator simulation and cokriging)
and reference fields (when available). The method is a nonparametric
clustering/classification algorithm that can assimilate significant
amounts of disparate data types with both continuous and categorical
responses without the computational burden associated with the
construction of positive definite covariance and cross-covariance
matrices. The combination of simplicity and computational speed
makes the method ideally suited for environmental subsurface characterization
and other Earth science applications with spatially autocorrelated
variables. |
Noviembre
2007
Teleseismic P wave imaging of the 26 December 2004 Sumatra-Andaman
and 28 March 2005 Sumatra earthquake ruptures using the Hi-net
array
Authors: Miaki Ishii, Peter M. Shearer et al
Link: Click Here
Abstract
Seismograms from a dense, high-quality seismic network in Japan
are used to investigate the characteristics of the 26 December
2004 Sumatra-Andaman and the 28 March 2005 Sumatran earthquakes.
The onset of the P waveforms are aligned through cross correlation,
and a simple concept of back-projecting seismic energy to a
grid of potential source locations is applied. The waveform
alignment removes the effects due to lateral variations in wave
speed between the hypocenter and each station. To better approximate
the effects of three- dimensional heterogeneity for paths originating
from grid points away from the hypocenter, cross-
|
correlation
results of the P waveforms from aftershocks are introduced. This
additional information leads to improved resolution of smaller-scale
features near many of the aftershocks by reducing wavefront distortion.
The back-projection analysis provides a quick assessment of the
spatiotemporal extent and variability of relative high-frequency
energy release, which can be translated into an estimate of the
moment magnitude, as well as an unparalleled view of high-frequency
rupture propagation. The results are, in general, consistent with
those obtained from more involved source inversion methods. The
2004 Sumatra-Andaman earthquake released most energy in a region
northwest of the Sumatra island and the rupture extended to the
northern Andaman islands, about 1300 km from the epicenter. This
northern portion of the rupture radiated a considerable amount
of energy, but there is little evidence of slow slip. The 2005
event is imaged to have bilateral rupture with northwestern slip
occurring for about 50 s before it moved to the southeast of the
epicenter. |
Octubre 2007
Rupture characterization and aftershock relocations for the
1994 and 2006 tsunami earthquakes in the Java subduction zone
Authors: Susan L. Bilek and E. Robert Engdahl
Link: Click Here
Abstract
Two tsunami earthquakes occurred in the Java subduction zone
within the last 12 years, providing multiple well-recorded tsunami
events for analysis. The June 2, 1994 Mw = 7.8 and July 17,
2006 Mw = 7.7 reverse mechanism earthquakes produced tsunami
with locally large (>8 m)
|
runups, were deficient
in high frequency energy relative to long period radiation, had
primarily normal faulting aftershocks, and occurred near subducting
seamounts. Previous models of the 1994 event and aftershocks suggested
extension in the outer rise due to slip at depth at a subducted
seamount. Here we relocate aftershocks for the 1994 and 2006 events
and compare to slip distributions and bathymetry. Most 1994 aftershocks
occurred within the slab updip of the mainshock, consistent with
the extensional hypothesis, but the 2006 aftershocks have a more
complex pattern. Regional subducting features as asperities and
barriers may influence both events, but in different ways. |
Octubre
2007
Radon daily signals in the Elat Granite, southern Arava, Israel
Authors: G. Steinitz, O. Piatibratova et al
Link: Click Here
Abstract
High time resolution monitoring of radon (=222Rn) in three boreholes,
4, 10 and 53 m deep, along a 0.6 km transect is carried out
in massive granite in southern Israel. Three components of variation
occur in the measured signal (MS) - seasonal radon (SR - periodic),
multiday (MD), and daily radon (DR - periodic). Temporal variation
of the components suggests an association between the overall
level of the long-term variation and the amplitude of the daily
variation. The daily mean level of radon and the daily standard
deviation vary periodically throughout the year. Time offsets
occur among
|
time
series of the MS and were investigated also for the MD and DR
components, using consecutive 20-day intervals spanning +900 days.
The resulting time series show that systematic time offsets occur,
whereby the radon signal always occurs first at the easternmost
site. The MD shows a gradually varying lag of 0-12 h, and the
DR a stable 1-3 h lag. Spectral analysis shows that diurnal (24-h)
and semidiurnal (12-h) periodic components characterize the DR.
The amplitudes of these components exhibit regular temporal variation
having a seasonal pattern. The ratios of co-occurring amplitudes
of these components define a linear pattern indicating a fundamental
statistical property in the frequency domain of the radon time
series. The results indicate that unrecognized dynamic processes
are driving the radon signal in the subsurface regime of the pluton,
suggesting new prospects for radon behavior in the frame of interacting
geodynamic and Earth-Sun system related processes. |
Octubre 2007
Representing aquifer architecture in macrodispersivity models
with an analytical solution of the transition probability matrix
Authors: Zhenxue Dai, Andrew Wolfsberg et al
Link: Click Here
Abstract
The multi-dimensional transition probability model represents
hydrofacies architecture in modeling aquifer heterogeneity.
The structure of the aquifer architecture is mathematically
characterized by a canonical representation of the transition
probability matrix, solved by an eigenvalue decomposition method.
Whereas the eigenvalue decomposition has been numerically solved
previously, we show here that it can be analytically solved
under the assumptions that
|
cross-transition probabilities
are dictated by facies proportions and that the juxtapositional
tendencies of the facies are symmetric. Although limited by the
assumptions, analytical solutions provide more immediate insights
about the relationships between transition probability and facies
proportion and mean length. The analytical solution is first tested
by comparison with the numerical solutions and then used to represent
hydrofacies architecture within expressions for the spatial covariance
of conductivity and the macrodispersivity. The relationship between
the longitudinal macrodispersivity and integral scale, the indicator
correlation length, and the facies proportion is represented in
an equation for estimating the field-scale dispersivity. An example
is used to show how sedimentary structures, conductivity contrasts,
and facies mean lengths affect the scales of the macrodispersivity.
|
Octubre
2007
Resistivity-based monitoring of biogenic gases in peat soils
Authors: Lee Slater, Xavier Comas et al
Link: Click Here
Abstract
Biogenic free-phase gas (FPG) formation was induced in a peat
block (dimensions 0.28 × 0.21 × 0.21 m) extracted
from a peatland in Maine. Electrical resistivity (ER), surface
deformation, and methane (CH4) flux from the peat surface was
monitored over a 48-day period during which the temperature
remained constant at 21 ± 1°C. ER measurements were
made on 5 vertical electrode arrays, each containing 20 electrodes
spaced at 0.01-m intervals. Surface deformation was monitored
using 30 elevation rods equally spaced across the surface of
the block, and average CH4 flux from the peat surface estimated
by integration of measurements obtained with a portable gas
detector over a 20-min time period. Pore water conductivity
was recorded at three depths (0.06, 0.09, and 0.15 m) at a single
point in the block. ER measurements were inverted for the ratio
resistivity change relative to a background data set and corrected
for changes in pore fluid conductivity, permitting an estimate
of equivalent change in gas content assuming (1) insignificant
surface conduction, (2) porosity changes estimated from peat
surface expansion, and (3)
|
an
Archie saturation exponent n of 1.3 based on results from a parallel
block experiment. The resistivity ratios reveal a pattern of FPG
evolution consistent with surface deformation and CH4 flux data.
During the first part of the experiment (approximately the first
24 days), a gradual buildup in FPG within a layer 4-6 cm below
the peat surface (water table) occurs concurrent with modest surface
deformation and low CH4 fluxes. In contrast, during the latter
half of the experiment (approximately 25-48 days), a complex pattern
of more pronounced gas buildup and release at multiple depths
occurs concurrent with large rates of surface deformation and
higher CH4 fluxes. The experiment demonstrates that ER monitoring
is a viable geophysical technology for imaging and monitoring
biogenic gas fluxes in peat soils. Here the resistivity clearly
shows that FPG is preferentially generated in layers about 0.04-0.06
m below the peat surface and that the buildup of gas is spatially
nonuniform even in a relatively small peat block. Furthermore,
the experimental results suggest that factors other than temperature
and atmospheric pressure must control biogenic gas accumulation
and release. As the method is readily deployable at the field
scale, possibly in an autonomous monitoring mode, resistivity
measurements may permit significant improvements in understanding
of carbon gas generation and release from northern peatlands.
|
Octubre 2007
A program for magnetic susceptibility-equivalent pore conversion
Authors: J. Jezek and F. Hrouda
Link: Click Here
Abstract
Pore magnetic anisotropy can be used to estimate the average
geometry of void spaces in rocks in the form of the equivalent
pore (EP) ellipsoid. Direct computation of EP from
|
measured magnetic susceptibility
is impossible. We present a method and a Matlab program for automatic
magnetic susceptibility-equivalent pore conversion. Input data
are the magnetic parameters (P and T, or L and F) representing
the bulk magnetic anisotropy, and the intrinsic susceptibility
of the fluid used in the measurement. EP is estimated iteratively
by a repeated look-up table procedure using P and T values computed
in a coarse grid of EP axial ratios. The program may be downloaded
from the EarthRef.org Digital Archive. |
Octubre
2007
THE EARLY ANTHROPOGENIC HYPOTHESIS: CHALLENGES AND RESPONSES
Authors: William F. Ruddiman
Link: Click Here
Abstract
Ruddiman (2003) proposed that late Holocene anthropogenic intervention
caused CH4 and CO2 increases that kept climate from cooling
and that preindustrial pandemics caused CO2 decreases and a
small cooling. Every aspect of this early anthropogenic hypothesis
has been challenged: the timescale, the issue of stage 11 as
a better analog, the ability of human activities to account
|
for
the gas anomalies, and the impact of the pandemics. This review
finds that the late Holocene gas trends are anomalous in all ice
timescales; greenhouse gases decreased during the closest stage
11 insolation analog; disproportionate biomass burning and rice
irrigation can explain the methane anomaly; and pandemics explain
half of the CO2 decrease since 1000 years ago. Only ~25% of the
CO2 anomaly can, however, be explained by carbon from early deforestation.
The remainder must have come from climate system feedbacks, including
a Holocene ocean that remained anomalously warm because of anthropogenic
intervention. |
Octubre 2007
Climate response to projected changes in short-lived species
under an A1B scenario from 2000-2050 in the GISS climate model
Authors: Drew T. Shindell, Greg Faluvegi et al
Link: Click Here
Abstract
We investigate the climate forcing from and response to projected
changes in short-lived species and methane under an A1B scenario
from 2000-2050 in the GISS climate model. We present a meta-analysis
of new simulations of the full evolution of gas and aerosol
species and other existing experiments with variations of the
same model. The comparison highlights the importance of several
physical processes in determining radiative forcing, especially
the effect of climate change on stratosphere-troposphere exchange,
heterogeneous sulfate-nitrate-dust chemistry, and changes in
methane oxidation and
|
natural emissions.
However, the impact of these fairly uncertain physical effects
is substantially less than the difference between alternative
emission scenarios for all short-lived species. The net global
mean annual average direct radiative forcing from the short-lived
species is .02 W/m2 or less in our projections, as substantial
positive ozone forcing is largely offset by negative aerosol direct
forcing. Since aerosol reductions also lead to a reduced indirect
effect, the global mean surface temperature warms by ~0.07°C
by 2030 and ~0.13°C by 2050, adding 19% and 17%, respectively,
to the warming induced by long-lived greenhouse gases. Regional
direct forcings are large, up to 3.8 W/m2. The ensemble-mean climate
response shows little regional correlation with the spatial pattern
of the forcing, however, suggesting that oceanic and atmospheric
mixing generally overwhelms the effect of even large localized
forcings. Exceptions are the polar regions, where ozone and aerosols
may induce substantial seasonal climate changes. |
Octubre
2007
Field relations between the spectral composition of ground
motion and hydrological effects during the 1999 Chi-Chi (Taiwan)
earthquake
Authors: Alexander Wong and Chi-Yuen Wang
Link: Click Here
Abstract
The possibility that the frequency content of ground motion
in earthquakes affects hydrological responses to earthquakes
has yet to be extensively tested in the field. The M w = 7.6
Chi-Chi earthquake provided an opportunity to do so, as widespread
liquefaction and groundwater level changes were recorded by
instrumentation in the
|
Choshui
Alluvial Fan. Analysis of the resulting data shows that ground
motion parameters that measure low-frequency ground motions are
more strongly correlated with coseismic groundwater level change
and the occurrence of liquefaction than parameters that measure
high-frequency motions. Notably, horizontal peak ground acceleration,
a metric often used to enumerate the strength of ground motion
in liquefaction analysis, is weakly correlated with these hydrological
effects. It is unclear from this analysis whether low-frequency
ground motion caused coseismic groundwater level change and liquefaction
or whether the hydrological changes changed the spectral composition
of the observed seismograms during the Chi-Chi earthquake. |
Octubre
2007 | Link: Clic Aquí
Movimiento de la corteza en la zona de Valdivia - 1960:
Entendiendo el ciclo de deformación del terremoto |
Autores: Kelin Wang, Yan Hu et al

Abstract
La deformación temporal observada en los ciclos de terremotos
y la permanente tensión asociada a la convergencia de
placas, nos puede propor-cionar pistas críticas para
entender la geodiná-
mica y el riesgo de terremotos en las zonas de subducción.
En el sector de Valdivia que se vio afectado por el terremoto
Mw 9.5 en el año 1960, nosotros hemos obtenido información
de ambos tipos a partir de lecturas de GPS. Nuestras esti-
maciones de velocidad coinciden con el rango latitudinal completo,
considerando las primeras llegadas. Las nuevas observaciones
revelaron un patrón de oposición del movimiento
entre la costa y los sitios interiores, consistente con las
observaciones septentrionales reportadas previamente. Los datos
apoyan el modelo de prolongada deformación postsísmica
como resultado de relajación tensional viscoelástica
del manto. Nuestras observaciones también propor-
cionan la primera evidencia geodésica sobre el movimiento
diestro del sistema de fallas intra-
volcánico y la consecuente traslación del frente
de la curvatura. El movimiento puede ser modelado con una rapidez
de 6.5 mm/a, incluyendo el 75% del margen paralelo de las placas
de Nazca y Sudamericana, movimiento que disminuye hacia el norte.
Además, las observaciones de GPS indican una disminución
meridional en las velocidades del margen normal del área
costera. Preferimos explicar el fenómeno en términos
de cambios de las propiedades de la superficie de fricción.
Debido a la juventud de la placa subducting y al régimen
termal templado en el sur, el actual sector bloqueado de la
placa puede ser más estrecho. Utilizando un modelo visco-
elástico 3D basado en elementos finitos, noso-
tros demostramos que nuestra explicación es consistente
con las observaciones a primer orden (aunque no es la única
explic. posible).
|
Octubre 2007
Crustal motion in the zone of the 1960 Chile earthquake:
Detangling earthquake-cycle deformation and forearc-sliver translation
Authors: Kelin Wang, Yan Hu et al
Link: Click Here
Abstract
Temporary deformation in great earthquake cycles
and permanent shear deformation associated with oblique plate
convergence both provide critical clues for understanding geodynamics
and earthquake hazard at subduction zones. In the region affected
by the Mw 9.5 great Chile earthquake of 1960, we have obtained
GPS observations that provide information on both types of deformation.
Our velocity solutions for the first time span the entire latitudinal
range of the 1960 earthquake.
The new observations revealed a pattern of opposing (roughly
arc-normal) motion of coastal and inland sites, consistent with
what was reported earlier for the northern part of this region.
This finding supports the model of prolonged postseismic deformation
as a result of viscoelastic stress relaxation in the mantle.
The new observations also provide the first geodetic evidence
for the dextral motion of an intravolcanic arc fault system
and the consequent northward translation of a forearc sliver.
The sliver motion can be modeled using a rate of 6.5 mm/a, accommodating
about 75% of the margin-parallel component of Nazca–South
America relative plate motion, with the rate diminishing to
the north. Furthermore, the new GPS observations show a southward
decrease in margin-normal velocities of the coastal area. We
prefer explaining the southward decrease in terms of changes
in the width or frictional properties of the megathrust seismogenic
zone. Because of the much younger age of the subducting plate
and warmer thermal regime in the south, the currently locked
portion of the plate interface may be narrower.
Using a three-dimensional viscoelastic finite element model
of postseismic and interseismic deformation following the 1960
earthquake, we demonstrate that this explanation, although not
unique, is consistent with the GPS observations to the first
order.
|
Octubre
2007
Morphology of the magnetic field near Titan: Hybrid model
study of the Cassini T9 flyby
Authors: E. Kallio, I. Sillanpää et al
Link: Click Here
Abstract
We study the deformation and morphology of the magnetic field
near Titan by a three-dimensional numerical quasi-neutral hybrid
model (HYB-Titan). We analyze two runs, one in which the sub-rotating
plasma consists of oxygen (O+) ions and protons (H+) and another
with only protons. We
|
find
that both cases result in the generation of Alfvén wing-like
flux tubes. In comparison with the proton-only case, the direction
and magnitude of the magnetic field in the oxygen-rich flow case
are in better agreement with Cassini magnetic field observations
during its Titan T9 flyby on Dec. 26, 2006, suggesting that the
sub-rotating plasma contained heavy ions. The oxygen-rich run
also supports the hypotheses that (1) the sub-rotating flow direction
was offset by ~30°-40° from the ideal rotating flow direction
and that (2) the spacecraft was magnetically connected to the
sunlit side of near-Titan space in the tail when it approached
Titan on the T9 flyby. |
Octubre 2007
Dependence of surface fracture
energy on earthquake size: A derivation from hierarchical self-similar
fault zone geometry
Authors: Kenshiro Otsuki
Link: Click Here
Abstract
Based on the hierarchical self-similar geometry of fault zones
in which fault segments and jogs of different hierarchical ranks
are nested, I
|
formulated surface
fracture energy of jogs and the mean surface fracture energy averaged
over the rupture surface. I incorporated this theory with a new
finding that the distribution of microfractures in fault jogs
is fractal with a universal fractal dimension of 1.56 (2-D measurements),
and derived the result that both fracture energy of jogs and the
mean fracture energy are proportional to the rupture length to
the power of 0.56. |
|
|
Octubre
2007
Imaging mantle discontinuities using least squares Radon transform
Authors: Yuling An, Yu Jeffrey Gu et al
Link: Click Here
Abstract
Differential traveltimes of SS precursors have been widely used
to map large-scale mantle structure and depths of discontinuities.
Measurements are commonly made by stacking individual records
to increase the signal-to-noise ratio pertaining to these mild
reflections. However, ray parameters of the SS precursors are
typically less well constrained and undesired seismic arrivals
with vastly different slownesses (for example, scattered waves)
could potentially contaminate the time domain stacks. To overcome
these pitfalls, we introduce a processing scheme based on the
Radon transform and well-constrained inversions. Our
|
method
is particularly effective in suppressing background noise, constraining
differential ray parameter and traveltime, and detecting weak
reflected or converted phases. We apply the Radon-based method
to delineate the discontinuity structure beneath the northeastern
Pacific Ocean and the northwestern Canada. An elevated 410-km
discontinuity and a thickened transition zone (252 km) are observed
beneath the northern British Columbia, which may be caused by
remnant oceanic lithosphere from the subduction of Kula-Farallon
plate under North America. A thin transition zone is identified
beneath the northern Pacific Ocean and its presence is supported
by a low shear-velocity anomaly from recent tomographic models.
The improved accuracy and resolution using the least squares Radon
transform also offer clear evidence for a regional 520-km discontinuity
and several mild reflectors in the depth ranges of 250-330 km
and 900-1200 km. We do not observe a 220-km discontinuity beneath
the study region. |
Octubre 2007
Quantitative estimates of magnetic
field reconnection properties from electric and magnetic field
measurements
Authors: F. S. Mozer and A. Retinò
Link: Click Here
Abstract
Reconnection occurs in a reconnection magnetic field geometry
when there are positive electric field components tangential
to the magnetopause and a magnetic field component normal to
it. Because these three components are the smallest of the six
electric and magnetic fields, their magnitudes are difficult
to determine because of errors in, or oscillations of, the assumed
constant direction normal to the current sheet. A method is
described for minimizing these
|
errors by appropriate selection of the normal
direction and by analyzing the correlations between the large
normal electric field and the large tangential magnetic field.
The correlation coefficients are equal to ratios of the small
fields, which are combined with the less accurate measurements
of the averages of the small fields to produce best estimates
of the small fields. For more than 120 magnetopause crossings,
about 40% had such correlations that signify static conditions
during those crossings. This method is applied to 22 polar subsolar
magnetopause crossings to show that most were located in the
ion diffusion region, as defined by the change of the total
magnetic field, and that 14 had a large and steady reconnection
rate with a zero parallel electric field. In these events the
reconnection rate decreased with increasing guide magnetic field.
|
Octubre
2007
Hybrid simulation of Titan's magnetic
field signature during the Cassini T9 flyby
Authors: S. Simon, G. Kleindienst et al
Link: Click Here
Abstract
Titan's induced magnetotail is analyzed in terms of a three-dimensional
hybrid model, treating the electrons as a massless, charge-neutralizing
fluid, whereas a completely kinetic approach is used to cover
ion dynamics. Since during the T9 flyby of Cassini, the spacecraft
passed through
|
Titan's
induced magnetotail in a distance of about 4 Titan radii downstream
of the obstacle, this flyby offers a unique chance to study the
magnetic lobe structure in Titan's wake region. The key features
of the measured magnetic field signature have shown to be completely
reproducible in the framework of the hybrid approximation. Besides,
the influence of the ionospheric production rates as well as the
upstream plasma direction is investigated. While changes of the
ion production rate yield only a slight modification of the magnetotail
structure, the magnetic field enhancement in the lobes is strongly
modified by changes of the upstream flow direction. |
Octubre 2007
Postseismic deformation of the
Andaman Islands following the 26 December, 2004 Great Sumatra-Andaman
earthquake
Authors: J. Paul, A. R. Lowry et al
Link: Click Here
Abstract
Two years after the Great Sumatra-Andaman earthquake the 3.1
m WSW coseismic displacement at Port Blair, Andaman Islands,
had increased by 32 cm. Postseismic uplift initially exceeded
1 cm per week and decreased to <1 mm/week. By 2007 points
near Port Blair had
|
risen more than 20 cm,
a 24% reversal of coseismic subsidence. Uplift at eight GPS sites
suggests a gradual eastward shift of the coseismic neutral axis
separating subsidence from uplift. Simulations of the GPS postseismic
displacements as viscoelastic relaxation of coseismic stress change
and as slip on the plate interface indicate that slip down-dip
of the seismic rupture dominates near-field deformation during
the first two years. Postseismic slip beneath the Andaman Islands
released moment equivalent to a magnitude M w = 7.5 earthquake,
and the distribution suggests deep slip in the stable frictional
regime accelerated to catch up to the coseismic rupture. |
Octubre
2007
Detecting tsunami genesis and scales
directly from coastal GPS stations
Author: Y. Tony Song
Link: Click Here
Abstract
Different from the conventional approach to tsunami warnings
that rely on earthquake magnitude estimates, we have found that
coastal GPS stations are able to detect continental slope displacements
of faulting due to big earthquakes, and that the detected seafloor
displacements are able to determine tsunami source energy and
|
scales
instantaneously. This method has successfully replicated three
historical tsunamis caused by the 2004 Sumatra earthquake, the
2005 Nias earthquake, and the 1964 Alaska earthquake, respectively,
and has been compared favorably with the conventional seismic
solutions that usually take hours or days to get through inverting
seismographs. Because many coastal GPS stations are already in
operation for measuring ground motions in real time as often as
once every few seconds, this study suggests a practical way of
identifying tsunamigenic earthquakes for early warnings and reducing
false alarms. |
Octubre 2007
Formation of scoria cone during
explosive eruption at Izu-Oshima volcano, Japan
Authors: K. Mannen and T. Ito et al
Link: Click Here
Abstract
Scoria cones have been regarded to be formed by accumulation
of ballistic bombs ejected by mild eruptions. However, more
recent geological investigations show that some scoria cones
could be formed during explosive eruptions. Here, we demonstrate
how the scoria cone of the 1986 Izu-Oshima eruption was formed
during the explosive
|
eruption. We measured
particle fractionation of the cone and propose a theoretical model
to explain the observation. The model considers lateral transport
of particles by turbulent eddies; particles that reached characteristic
column radius, L, are laterally transported to ωL where
they starts to free-fall. We obtained ω = 1.2 and 2.5 for
larger and smaller particles, respectively, which is consistent
to the observation. In the model, extensive fallout takes place
at the base of the column where it expands rapidly. We suggest
that lateral particle projection and the rapid column expansion
are the key processes to form the cones. |
Octubre
2007
Estimation of the stress field in the region of the 2000
Western Tottori Earthquake: Using numerous aftershock focal mechanisms
Authors: Yohei Yukutake, Yoshihisa Iio et al
Link: Click Here
Abstract
We estimated the regional stress field before and after the
2000 Western Tottori Earthquake (Mw = 6.6) from numerous focal
mechanisms. To constrain the stress field prior to the main
shock, we compared the spatial distribution of the static stress
changes generated by the main shock with the stress field following
the main shock. In the northern and central parts of the aftershock
region, it is inferred that the deviatoric stress magnitude
prior to the main shock was too large to be affected by the
static stress changes. The
|
direction
of the maximum principal stress axis prior to the main shock is
consistent with the tectonic stress field. On the other hand,
in the southern part of the aftershock area, in and around the
region where the main shock slip was large, it is found that the
stress field is spatially heterogeneous after the main shock because
of the large magnitude of the static stress changes. Around the
southern edge of the fault, the spatial distribution of the static
stress changes is consistent with that of the P axis azimuths
of the aftershocks, and it is inferred that the deviatoric stress
magnitude prior to the main shock was small enough to be affected
by the static stress changes (5 MPa). The strength of the preexisting
aftershock fault planes in the southern edge might be exceptionally
weak. It is found that the stress field prior to the main shock
was inhomogeneous on scales smaller than the length of the main
shock fault. |
Octubre
2007
Generation of a ULF wave resonator in the magnetosphere
by neutral gas release
Authors: Manish Mithaiwala, Leonid Rudakov et al
Link: Click Here
Abstract
ULF waves generated by the release of a neutral gas in the equatorial
plane of the Earth's magnetosphere will be trapped between two
turning points forming a resonator. This effectively prolongs
the lifetime of the ensuing turbulence making it more useful
for pitch angle scattering
|
trapped energetic
electrons. In the multispecies plasma environment of the magnetosphere,
the presence of a small amount of helium causes the waves to reflect
when their frequency matches the Buchsbaum frequency as they travel
along field lines away from the equator. For a proton-electron-helium
plasma, the Buchsbaum frequency is near the helium cyclotron frequency.
With typical abundances of helium in the inner radiation belts,
there is virtually no loss of wave energy by cyclotron damping
at the turning points as the waves reflect between conjugate Buchsbaum
points. Wave energy is eventually dissipated by collisional effects
and radial convection. |
Octubre
2007
Groundwater contaminant source identification by a hybrid
heuristic approach
Authors: Hund-Der Yeh, Tung-Hua Chang et al
Link: Click Here
Abstract
This study proposes an approach, SATS-GWT, that combines simulated
annealing (SA), tabu search (TS), and the three-dimensional
groundwater flow and solute transport model (MODFLOW-GWT). Our
approach is used to estimate the source information: source
location, release concentration, and release period. The sampling
concentrations at monitoring points are simulated by the MODFLOW-GWT
with an assumed release concentration and release period at
a known source location. In the source estimation process, the
source location is
|
selected
by TS within the suspected source area, and the trials for release
concentrations and release periods are generated by SA. The MODFLOW-GWT
is employed to compute the simulated concentrations at the monitoring
points with the trial solution. The above mentioned procedures
are repeated until the stopping criterion regarding the differences
of objective function value (OFV) is met. The last trial source
information which yields the best OFV is the final solution. Six
studies on the homogeneous site, two studies on the heterogeneous
site, and one study on the transient flow problem are conducted
in this study. A suggestion regarding the optimal number of monitoring
points and the condition for estimating the source information
is given on the basis of the studies for homogeneous and heterogeneous
aquifers. All results indicate that the proposed SATS-GWT can
give good estimations, even when the sampling concentrations contain
measurement errors. |
Octubre 2007
The study of Earth's magnetism (1269–1950): A foundation
by Peregrinus and subsequent development of geomagnetism and paleomagnetism
Authors: Vincent Courtillot and Jean-Louis Le Mouël
Link: Click Here
Abstract
This paper summarizes the histories of geomagnetism and paleomagnetism
(1269–1950). The role of Peregrinus is emphasized. In
the sixteenth century a debate on local versus global departures
of the field from that of an axial dipole pitted Gilbert against
Le
|
Nautonier. Regular measurements
were undertaken in the seventeenth century. At the turn of the
nineteenth century, de Lamanon, de Rossel, and von Humboldt discovered
the decrease of intensity as one approaches the equator. Around
1850, three figures of rock magnetism were Fournet (remanent and
induced magnetizations), Delesse (remagnetization in a direction
opposite to the original), and Melloni (direction of lava magnetization
acquired at time of cooling). Around 1900, Brunhes discovered
magnetic reversals. In the 1920s, Chevallier produced the first
magnetostratigraphy and hypothesized that poles had undergone
enormous displacements. Matuyama showed that the Earth's field
had reversed before the Pleistocene. Our review ends in the 1940s,
when exponential development of geomagnetism and paleomagnetism
starts. |
Octubre
2007
Uncertainty analysis of earthquake source parameters determined
from InSAR: A simulation study
Authors: John Dawson and Paul Tregoning
Link: Click Here
Abstract
We assessed the accuracy of earthquake source parameters inverted
from simulated Interferometric Synthetic Aperture Radar (InSAR)
data. Using focal mechanisms of Australian earthquakes (1959
to the present), we simulated synthetic two-pass InSAR observations
with realistic spatial noise derived from the characteristics
of actual ERS-2 and ENVISAT InSAR data observed over Australia.
The precision of two-pass satellite SAR interferometry with
ERS-2 and ENVISAT SAR data in the Australian region can approach
±2 mm (1σ) and is routinely at the ±4
mm level. The use of spatially correlated
|
observational
weights has minimal impact on the accuracy of earthquake source
parameters inverted from InSAR data. In most cases single geometry
(i.e., ascending or descending) InSAR observations can be used
to accurately determine earthquake source parameters, although
typically a combined geometry reduces the source parameter uncertainties
by a factor of 1.5. In general, earthquakes of magnitude <4.8
are unlikely to be observable by InSAR although very shallow events
would be detectable. InSAR is insensitive to magnitude 6.2 earthquakes
deeper than 10 km, and magnitude 5.5 deeper than 6 km. For earthquake
magnitudes ≥5.8 (average depth 6.5 km) we could estimate the
epicenter of the rupture with an average accuracy of 0.25 km,
depth to within 0.5 km and the fault orientation to better than
2°. Our findings, based on simulated Australian earthquakes,
are representative of typical intraplate earthquakes and would
be valid in many other regions. To date no actual earthquakes
have been observed by InSAR in Australia. |
Octubre 2007
Neural network emulations for complex multidimensional
geophysical mappings: Applications of neural network techniques
to atmospheric and oceanic satellite retrievals and numerical
modeling
Authors: Vladimir M. Krasnopolsky
Link: Click Here
Abstract
A group of geophysical applications, which from the mathematical
point of view, can be formulated as complex, multidimensional,
nonlinear mappings and which in terms of the neural network
(NN) technique, utilize a particular type of
|
NN, the multilayer perceptron
(MLP), is reviewed in this paper. This type of NN application
covers the majority of NN applications developed in geosciences
like satellite remote sensing, meteorology, oceanography, numerical
weather prediction, and climate studies. The major properties
of the mappings and MLP NNs are formulated and discussed. Three
particular groups of NN applications are presented in this paper
as illustrations: atmospheric and oceanic satellite remote sensing
applications, NN emulations of model physics for developing atmospheric
and oceanic hybrid numerical models, and NN emulations of the
dependencies between model variables for application in data assimilation
systems. |
Octubre
2007
Faults and damage zones in fast-spread crust exposed on
the north wall of the Hess Deep Rift: Conduits and seals in seafloor
hydrothermal systems
Authors: Nicholas W. Hayman and Jeffrey
A. Karson
Link: Click Here
Abstract
The northern escarpments of the Hess Deep Rift provide cross-sectional
views of in situ, aprox 1-Ma-old, upper oceanic crust that underwent
extensive, spreading-related brittle deformation. Most of the
deformation and associated alteration occurred within the locus
of magmatic construction of the East Pacific Rise, in the presence
of high-temperature hydrothermal fluids. Passing laterally from
undeformed host rocks, brittle deformation zones are classified
as (1) damage zones where densely spaced fractures overprint
the primary structure of dikes and lavas, (2) cataclastic zones
where interconnected fractures, comminuted grains, and matrix
minerals define deformational fabrics, and (3) very fine-grained,
gouge-filled fault cores. Relative to
|
the
host rock, damage and cataclastic zones are rich in veins of chlorite
and/or actinolite, and lesser amounts of titanite, epidote, and
quartz. These phases mark relict hydrothermal fluid pathways.
Trace and major element compositions of representative samples
also indicate fault-localized hydrothermal alteration, including
an increase in MgO by several weight percent within cataclastic
and damage zones. In contrast, the fault cores are composed of
very finely comminuted basaltic material and have MgO concentrations
similar to the damage zones. Integrated compositional, textural,
and outcrop-scale structural data inform an evolutionary model
for fault growth from the early, widespread dilational phases
of damage-zone development to more restricted noncoaxial strain
in the cataclastic zones. With continued fault development, gouge
develops and seals the fault cores. While the fault cores are
sealed by gouge, surrounding zones remain conduits to hydrothermal
fluid flow, except where sealed by secondary minerals. Sealed
faults can later be reactivated as conduits with additional increments
of fault slip. The dual behavior of faults as conduits and seals
inevitably leads to compartmentalization of the flow regime in
subaxial and ridge-flank areas. |
Octubre 2007
Focused and diffuse extension in controls of ocean ridge
segmentation in analogue models
Authors: Tatiana Tentler
Link: Click Here
Abstract
This work aims to provide insights into the initiation and evolution
of segments along ocean ridges as they propagate and mature.
It is based on the results of analogue modeling performed in
a centrifuge with applied acceleration of gravity force. The
effect of the varied degree of extension
|
allowed during propagation
of analogue ridge results in segments being focused along the
ridge axes or widely distributed within the entire axial zone.
Three distinct types of segment interaction are documented, controlled
by the relative distribution and orientation of interacting segments
and leading to the development of linear ridges, overlapping spreading
centers, and oblique transform zones. The importance of the initial
distribution of segments in the ridge evolution inferred in the
models suggests its significance for formation of natural ridges
whose mature configurations must be largely inherited from that
of spreading centers formed in the early stages of focused or
diffuse extension. |
Octubre
2007
Lithospheric modification during crustal extension in
the Main Ethiopian Rift
Authors: Tyrone Rooney, Tanya Furman et al
Link: Click Here
Abstract
Quaternary lavas erupted in zones of tectonomagmatic extension
within the Main Ethiopian Rift (MER) preserve details of lithospheric
structure in the East African Rift System. Despite observed
source heterogeneity, basalts, trachybasalts, and basaltic trachyandesites
erupted in the Wonjii Fault Belt (WFB) and the Silti-Debre Zeyit
Fault Zone (SDFZ) form coherent fractionation paths dominated
by variable removal of observed phenocryst phases. Crustal assimilation
is not widespread, though it is observed at the southern end
of the WFB where both fault belts merge; farther north, assimilation
of cumulate phases related to fractional crystallization of
previous magmas is identified. Shallow fractionation conditions
(1 kbar) within
|
the
WFB do not change from north to south. In contrast, lavas erupted
within the contemporaneous SDFZ fractionate at various crustal
depths. These results indicate a better developed magmatic system
beneath the WFB where magmas rose quickly before undergoing more
significant fractionation at near surface levels and a less developed
system beneath the SDFZ. The distribution of magmatism and extant
geophysical data indicate thinned crust and a single rift-centered
zone of magmatic activity northeast of 8°30′N, consistent
with a transitional lithosphere between continental and oceanic
settings. Southwest of 8°30′N, thicker crust and rift-marginal
axes of extension suggest lithosphere with continental affinities.
The WFB is propagating southward in response to extension within
the Red Sea Rift; the northward propagating SDFZ is related to
rifting within the East African Rift System. This region records
the unification of two rift systems, requiring care in interpreting
the MER as simply transitional between continental and oceanic
environments. |
Septiembre 2007
Long-term probabilistic forecast of climate impact on
air quality: Model development and t* distribution
Authors: Shao-Hang Chu
Link: Click Here
Abstract
Models are great tools to test ideas. Their usefulness, however,
depends on their ability to simulate the current reality and
to predict the future. In this study, I have derived a new t*
distribution. I show that a statistical model based on the t*
distribution of station temporal data
|
is capable of predicting
the probability of any future outcome to exceed a specific value
using only the currently available sample statistics assuming
a normal random variable. In an air quality management application
the model has demonstrated categorically an average success rate
of over 80% both in simulating the current ozone nonattainment
areas and in forecasting the rate of future violation of the 8-hour
ozone National Ambient Air Quality Standards in the United States
for up to 12 years. While the predictability of deterministic
climate models is still limited by large uncertainties, the probabilistic
forecast by this model provides a promising alternative in assessing
the climate impact on environment for decades. |
Septiembre
2007
Measurements of 220Rn and 222Rn and CO2 emissions in soil
and fumarole gases on Mt. Etna volcano (Italy): Implications for
gas transport and shallow ground fracture
Authors: S. Giammanco, K. W. W. Sims et
al
Link: Click Here
Abstract
Measurements of 220Rn and 222Rn activity and of CO2 flux in
soil and fumaroles were carried out on Mount Etna volcano in
2005–2006, both in its summit area and along active faults
on its flanks.
|
We
observe an empirical relationship between (220Rn/222Rn) and CO2
efflux. The higher the flux of CO2, the lower the ratio between
220Rn and 222Rn. Deep sources of gas are characterized by high
222Rn activity and high CO2 efflux, whereas shallow sources are
indicated by high 220Rn activity and relatively low CO2 efflux.
Excess 220Rn highlights sites of ongoing shallow rock fracturing
that could be affected by collapse, as in the case of the rim
of an active vent. Depletion both in 220Rn and in CO2 seems to
be representative of residual degassing along recently active
eruptive vents. |
Septiembre 2007
A simplified approach to quantifying predictive and parametric
uncertainty in artificial neural network hydrologic models
Authors: R. K. Srivastav, K. P. Sudheer et al
Link: Click Here
Abstract
One of the principal sources of uncertainty in hydrological
models is the absence of understanding of the complex physical
processes of the hydrological cycle within the system. This
leads to uncertainty in input selection and consequently its
associated parameters, and hence evaluation of uncertainty in
a model becomes important. While there has been considerable
interest in developing methods for uncertainty analysis of artificial
neural network (ANN) models, most of the methods are relatively
complex and/or require assumption about the prior distribution
of the uncertain parameters. This paper presents an effective
and simple way to
|
perform uncertainty
analysis for ANN-based hydrologic model. The method is based on
the concept of bootstrap technique and is demonstrated through
a case study of the Kolar River basin located in India. The method
effectively quantifies uncertainty in the model output and the
parameters arising from variation in input data used for calibration.
In the current study, the uncertainty due to model architecture
and the input vector are not directly considered; they have been
minimized during the model calibration. The results from the case
study suggest that the sampling variability of the training patterns
as well as the initial guess of the parameters of ANN do not have
significant impact on the model performance. However, despite
good generalization properties for the models developed in this
study, most of them fail to capture the hydrograph peak flow characteristics.
The proposed method of uncertainty analysis is very efficient,
can be easily applied to an ANN-based hydrologic model, and clearly
illustrates the strong and weak points of the ANN model developed.
|
Septiembre
2007
Close temporal correspondence between geomagnetic anomalies
and earthquakes during the 2002-2003 eruption of Etna volcano
Authors: G. Currenti, C. Del Negro et al
Link: Click Here
Abstract
The early stages of the 2002-2003 lateral eruption at Mount
Etna were accompanied by slow changes (over some hours) and
some rapid step offsets in the local magnetic field. At five
monitoring locations, the total magnetic field intensity has
been measured using continuously operating Overhauser seismic
events located roughly midway between the magnetic stations.
We quantitatively examine possible geophysical mechanisms, which
could cause the magnetic anomalies. The comparison between magnetic
data, seismicity and surface phenomena implies that piezomagnetic
effects are the primary
|
physical mechanism responsible for the observed magnetic anomalies
although the detailed cause of the rapid high stress change required
is not clear. The modeling of the observed coseismic magnetic
changes in terms of piezomagnetic mechanism provides further evidence
of the complex interaction between volcanic and tectonic processes
during dike propagation along the Northeast Rift. magnetometers
at a sampling rate of 10 s. The very unique aspect of these observations
is the close temporal correspondence between magnetic field offsets
and earthquakes that occurred in the upper northern flank of the
volcano on 27 October 2002 prior to a primary eruption. Rapid
coseismic changes of the magnetic field were clearly identified
for three of the most energetic earthquakes, which were concentrated
along the Northeast Rift at a depth of about 1 km below sea level.
Coseismic magnetic signals, with amplitudes from 0.5 to 2.5 nT,
have been detected for three of the largest |
|
Septiembre 2007
Presentan EIA para explotación del Tatio como alternativa
energética
(Boletín ACHS)
Link: Click Here
II Región.- Se presentó oficialmente
el Estudio de Impacto Ambiental (EIA) del proyecto "Perforación
Geotérmica Profunda El Tatio Fase Uno", a cargo
|
de la Compañía Geotérmica del
Norte S.A, que sometió al Sistema de Evaluación
de Impacto Ambiental (SEIA), ante la Comisión Regional
de Medio Ambiente de la Segunda Región este estudio.
El objetivo del proyecto sería realizar perforaciones
geotérmicas profundas en el sector del ex campamento
CORFO ubicado al sur del campo géiseres del Tatio, con
el propósito de verificar la factibilidad técnica
y económica de generar energía eléctrica
a partir de los recursos geotérmicos existentes. El proyecto
se instalaría a 120 km al este de Calama, aproximadamente
a 100 km al norte de San Pedro de Atacama. El terreno donde
se emplazará pertenece a Bienes Nacionales y es parte
del área concesionada a Geotérmica del Norte,
para explotar recursos geotérmicos. De acuerdo a los
antecedentes, el proyecto considera la perforación de
cuatro pozos de dos mil a dos mil quinientos metros cada
uno. Luego estos serán evaluados para determinar su potencial
geotérmico.
( Sustentable | Ver
Más | Impacto
Ambiental)
|
Septiembre
2007
Effects of mesoscopic-scale fault structure on dynamic earthquake
ruptures: Dynamic formation of geometrical complexity of earthquake
faults
Authors: Ryosuke Ando and Teruo Yamashita
Link: Click Here
Abstract
We develop a new hierarchical earthquake rupture model that
takes into account mesoscopic-scale fault structure; shear branches
nucleated on the main fault are specifically assumed as an example
of mesoscopic-scale fault structure. We numerically investigate
dynamic formation of fault geometry and its effects on dynamic
earthquake rupture process based on this rupture model. As long
as the length of the main fault is below a certain threshold
L m, the growth of these branches is shown to be arrested spontaneously
soon after their nucleation. The spatial distribution of arrested
branches is shown to form a self-
|
similar
geometrical structure. This suggests the existence of a simple
scaling relationship between small and large events as long as
the length of the main fault is below the threshold L m . However,
once the length of the main fault exceeds L m , a limited number
of branches begin unstable growth, and their sizes soon become
comparable to that of the main fault. In other words, mesoscopic-scale
branches are spontaneously transformed into macroscopic-scale
ones. This finding indicates the critical importance of the consideration
of mesoscopic-scale fault structure in understanding rupture dynamics:
Once macroscopic-scale branches are formed, they change the fault
geometry, which will considerably affect the rupture dynamics.
The emergence of macroscopic branches suggests that the above
mentioned simple scaling relation is never valid above a critical
length L m . Our study thus indicates that relationship between
small and large earthquakes is complicated by the spontaneous
transformation of a mesoscopic-scale fault structure into a macroscopic-scale
one. |
Septiembre 2007
PlanetWRF: A general purpose, local to global numerical model
for planetary atmospheric and climate dynamics
Authors: Mark I. Richardson, Anthony D. Toigo et al
Link: Click Here
Abstract
A new planetary atmospheric numerical model, "planetWRF,"
has been developed by modifying the Weather Research and Forecasting
(WRF) model. The model has generalized map projection, multiscale,
and nesting capabilities, blurring the distinction between global
and mesoscale models and enabling investigation of
|
coupling between processes
on all scales, including global. The model can also be run in
one, two, or three dimensions. The conversion of the dynamical
core for global application by altering the map projection grid
and the boundary conditions as well as conversion of the physics
parameterizations and constants for planetary application are
described. Validation of the global dynamical core through use
of standard forcing scenarios is presented. Example results from
a series of simulations for Mars, Titan, and Venus are shown to
demonstrate that the model performs well for a variety of planets
and operating modes (microscale, mesoscale, and global scale).
|
Septiembre
2007
Microbial nanowires: Is the subsurface "hardwired"?
Authors: Dimitrios Ntarlagiannis, Estella A. Atekwana et al
Link: Click Here
Abstract
The Earth's shallow subsurface results from integrated biological,
geochemical, and physical processes. Methods are sought to remotely
assess these interactive processes, especially those catalysed
by micro-organisms. Using saturated sand columns and the metal
reducing bacterium Shewanella oneidensis MR-1, we show that
electrically conductive appendages
|
called
bacterial nanowires are directly associated with electrical potentials.
No significant electrical potentials were detectable in columns
inoculated with mutant strains that produced non-conductive appendages.
Scanning electron microscopy imaging revealed a network of nanowires
linking cells-cells and cells to mineral surfaces, "hardwiring"
the entire length of the column. We hypothesize that the nanowires
serve as conduits for transfer of electrons from bacteria in the
anaerobic part of the column to bacteria at the surface that have
access to oxygen, akin to a biogeobattery. These results advance
understanding of the mechanisms of electron transport in subsurface
environments and of how microorganisms cycle geologic material
and share energy. |
Septiembre 2007
Spin transition and equations of state of (Mg, Fe)O solid
solutions
Authors: Yingwei Fei, Li Zhang et al
Link: Click Here
Abstract
We have performed a series of experiments to investigate the
compositional effect on the compression behavior of (Mg, Fe)O
solid solutions at high pressure. The in-situ synchrotron X-ray
diffraction data revealed abnormal volume contractions at about
40, 60 and 80 GPa for (Mg0.80, Fe0.20)O, (Mg0.61, Fe0.39)O,
|
and (Mg0.42, Fe0.58)O,
respectively. The volume contractions are associated with the
reported electronic transition of high-spin to low-spin in Fe2+,
and caused by the reduction of Fe2+ ionic radius across the transition.
A least-squares fit of the compression data to the Birch-Murnaghan
equation of state yielded bulk modulus K 0 (GPa) = 160 - 10X FeO
for the high-spin (Mg,Fe)O and K 0 = 170(3) GPa for the low-spin
(Mg,Fe)O. The equations of state of (Mg,Fe)O established in this
study are directly applicable to the Earth's lower mantle in composition
and pressure ranges and provide essential data for modeling the
density profile of the lower mantle. |
Septiembre
2007
Ergodic transport through aquifers
of non-Gaussian log conductivity distribution and occurrence of
anomalous behavior
Authors: Aldo Fiori, Igor Jankovic et al
Link: Click Here
Abstract
Three-dimensional advective transport of passive solutes through
isotropic porous formations of stationary non-Gaussian log conductivity
distributions is investigated by using an approximate semianalytical
model, which is compared with accurate numerical simulations.
The study is a continuation of our previous works in which formation
heterogeneity is modeled using spherical nonoverlapping inclusions
and an approximate analytical model was developed. Flow is solved
for average uniform velocity, and transport of an ergodic plume
is quantified by mass flux (traveltime distribution) at a control
plane. The analytical model uses a self-consistent transport
of an ergodic plume is quantified by mass flux (traveltime distribution)
at
|
a
control plane. The analytical model uses a self-consistent argument,
and it is based on the solution for an isolated inclusion submerged
in homogeneous background matrix of effective conductivity. As
demonstrated in the past, this analytical model accurately predicted
the entire distributions of traveltimes in formations of Gaussian
log conductivity distributions, as validated by numerical simulations.
The present study (1) extends the results to formations of non-Gaussian
log conductivity structures (the subordination model), (2) extends
the approximate analytical model to cubical blocks that tessellate
the entire domain, (3) identifies a condition in conductivity
distribution, at the tail of low values, that renders transport
anomalous with macrodispersivity growing without bounds, and (4)
provides links of our work to continuous time random walk (CTRW)
methodology, as applied to subsurface transport. It is found that
a class of CTRW solutions proposed in the past cannot be based
on solution of flow in formations with conductivity distribution
of finite integral scale. |
Agosto 2007
An adaptive neural network embedded genetic algorithm approach
for inverse water quality modeling
Authors: Rui Zou, Wu-Seng Lung et
al
Link: Click Here
Abstract
This paper proposes a neural network (NN)-embedded genetic algorithm
(GA) approach for solving inverse water quality modeling problems
to overcome the computational bottleneck of inverse modeling
by replacing a water quality
|
model with an efficient
NN functional evaluator. An existing one-step, NN-embedded GA
approach is found incapable of solving an inverse water quality
modeling problem because it tends to fail in guiding the global
search process to converge toward the near optima. As a remedy,
an adaptive NN-GA approach is proposed to achieve a gradual convergence
toward the near optima through an iterative network learning method.
The proposed approach is applied to a full-scale, numerical example,
and the result shows that the adaptive NN-GA approach is capable
of obtaining near-optimal solutions for the inverse problem of
a complicated water quality model. |
Agosto
2007
Rheologic controls on slab dynamics
Authors: Magali I. Billen &
Greg Hirth
Link: Click Here
Abstract
Several models have been proposed to relate slab geometry to
parameters such as plate velocity or plate age. However, studies
on the observed relationships between slab geometry and a wide
range of subduction parameters show that there is not a simple
global relationship between slab geometry and any one of these
other subduction parameters for all subduction zones. Numerical
and laboratory models of subduction provide a method to explore
the relative importance of different physical processes in determining
subduction dynamics. Employing 2-D numerical models with a viscosity
structure
|
constrained
by laboratory experiments for the deformation of olivine, we show
that the observed range in slab dip and the observed trends between
slab dip and convergence velocity, subducting plate age, and subduction
duration can be reproduced without trench motion (i.e., slab roll-back)
for locations away from slab edges. Successful models include
a stiff slab that is 100-1000 times more viscous than previous
estimates from models of plate bending, the geoid, and global
plate motions. We find that slab dip in the upper mantle depends
primarily on slab strength and plate boundary coupling, with a
small dependence on subducting plate age. Once the slab sinks
into the lower mantle the primary processes controlling slab evolution
are (1) the ability of the stiff slab to transmit stresses up
dip, (2) resistance to slab descent into the higher-viscosity
lower mantle, and (3) subduction-induced flow in the mantle-wedge
corner. |
Agosto 2007
Infrasonic tremor in the diffraction zone
Authors: David Fee & Milton
Garcés
Link: Click Here
Abstract
Volcanic tremor signals recorded by an infrasound array ~12.5
km from the active vent of Kilauea Volcano show clear diurnal
amplitude variations and originate from a well defined direction
pointing to the active Pu`u `O`o crater complex. Spectral amplitudes
between 0.02-0.3 Hz, a proxy for wind speed, increase when the
wind speeds increase during the daytime.
|
Spectral amplitudes
between 0.5-3 Hz, where infrasonic tremor is concentrated, increase
when the wind speeds decrease at night. Power spectrum analysis
shows that when ambient noise levels drop, infrasonic tremor amplitudes
rise. Array processing confirms this relationship and reveals
tremor detections and amplitudes decrease during the day. These
changes in infrasonic tremor spectral amplitude and coherence
are most likely due to diurnal variations in atmospheric boundary
layers affecting infrasonic propagation within the diffraction
zone. Recognition of these regional atmospheric patterns is necessary
to permit robust infrasonic remote sensing of volcanic processes.
|