PWG - Actualidad en Geofísica
HOME | Free E-book | Geophysics

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.

<< Anterior | Home: www.geofisica.cl | Siguiente >>