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

Mayo de 2007
Electrical resistivity image of the Jingsutu Graben at the NE margin of the Ejina Basin (NW China) and implications for the basin development
Authors: M. Becken, S. Hölz et al
Link: Click Here

Abstract
Magnetotelluric and transient electromagnetic data reveal the Graben geometry of a NNE striking depression within crystalline basement at the NE

margin of the Ejina Basin. We interpret this structure as a Late Mesozoic pull-apart basin, which is delineated by NNE-striking normal faults to the east and west. Satellite imagery, SRTM elevation data and previously published geophysical data of our work group suggest that the Graben is part of normal-fault system that continues at least 140 km farther to the south, thus forming the eastern margin of the Ejina Basin. NNE striking normal faults are atypical in a regional context, where NW striking and ENE striking fault systems prevail. Thus, this feature marks a distinct difference of the Ejina Basin when compared to other basin structures in the China-Mongolia border region.

Mayo de 2007
Analysis of Global Positioning System (GPS) radio occultation measurement errors based on Satellite de Aplicaciones Cientificas-C (SAC-C) GPS radio occultation data recorded in open-loop and phase-locked-loop mode
Author: Martin S. Lohmann
Link: Click Here

Abstract
The error characteristics of Global Positioning System (GPS) radio occultation (RO) measurement errors are studied based on Satellite de Aplicaciones Cientificas-C (SAC-C) GPS radio occultation data tracked in both open-loop (OL) and phase-locked-loop (PLL) mode. The error characteristics are derived by applying dynamical error estimation, i.e., without using any external data. The computed error profiles show that the mean measurement errors are the smallest in the height range between about 5-7

km and 20-25 km, about 0.2-1% for bending angles and 0.1-0.2% for refractivity at all latitudes. The largest measurement errors are found in the lower troposphere, where the mean bending angle measurement errors are within the range from 1 to 6%, whereas the mean refractivity measurement errors are within the range from 0.2% to 1%. From the error distributions, it is found that the occultation-to-occultation variability of the measurement errors generally spans one order of magnitude. The bending angle error correlation length is about 1 km and 100 m, at high and low altitudes, respectively, corresponding approximately to the cutoff frequency of the applied noise filters. The widths of the refractivity error autocorrelation functions are notably broader. The variability and the magnitude of the OL measurement errors are larger than for the PLL measurement errors. This is mainly attributed to the ability of OL tracking to track RO signals under atmospheric conditions for which PLL tracking fails.

Mayo de 2007
Viscosity of the Earth's fluid core and torsional oscillations
Authors: Jon Mound, Bruce Buffett et al
Link: Click Here

Abstract
Viscous effects on the temporal behavior and spatial structure of torsional oscillation normal modes of the Earth's fluid core are examined. If the viscosity of the fluid core is equal to that inferred from theoretical and experimental studies of liquid metals, then viscous effects on torsional oscillations are unlikely to be important; however, if the viscosity is as high as allowed from observational constraints, then Ekman layers at

the fluid core boundaries may play an important role in damping torsional oscillations. The observationally inferred decay time of decadal torsional oscillations leads to an upper bound on the viscosity of the fluid core at the core-mantle boundary of the order of 10-2 m2/s. A sufficiently large viscosity would result in torsional oscillation normal modes that undergo pure decay, as opposed to damped oscillatory motion. The viscosity value at which the temporal behavior switches between these regimes depends on the period of the normal mode. Numerical geodynamo models must reach an Ekman number of the order of 10-7 for free-slip boundaries or 10-10 for no-slip boundaries to accurately model the temporal behavior of decadal period torsional oscillations.

Mayo de 2007
Laboratory sandbox validation of transient hydraulic tomography
Authors: X. Liu, W. A. Illmanet al
Link: Click Here

Abstract
Hydraulic tomography is a method that images the hydraulic heterogeneity of the subsurface through the inversion of multiple pumping or cross-hole hydraulic test data. Transient hydraulic tomography is different from steady state hydraulic tomography in that it utilizes transient hydraulic head records to yield the distribution of hydraulic conductivity (K) as well as specific storage (S s ) of an aquifer. In this paper we demonstrate the robustness of transient hydraulic tomography through the use of hydraulic head data obtained from multiple cross-hole pumping tests conducted in a laboratory sandbox with deterministic heterogeneity. We utilize the algorithm developed by Zhu and Yeh (2005) to

conduct the transient inversions and validate the K and S s tomograms using a multimethod and multiscale validation approach previously proposed by Illman et al. (2006). Validation data consist of cross-hole tests not used in the inversion as well as other hydraulic tests that provided local (core, single-hole tests) as well as large-scale (unidirectional flow-through tests) estimates of hydraulic parameters. Results show that the algorithm is able to yield consistent estimates that agree with independently collected local as well as large-scale hydraulic parameter data. In addition, we find that the transient hydraulic tomography requires a fewer number of pumping tests to estimate a similar quality K tomogram when compared with steady state hydraulic tomography, as the former approach utilizes more data from each pumping test. Overall, we find that transient hydraulic tomography is a robust subsurface characterization technique that can delineate the subsurface heterogeneity in both K and S s from multiple pumping or cross-hole hydraulic tests.

Mayo de 2007
Methane thermodynamics in nanoporous ice: A new methane reservoir on Titan
Authors: L. F. Voss, B. F. Henson et al
Link: Click Here

Abstract
A porous, icy regolith has been proposed to house a methane/ethane/nitrogen liquid reservoir on Titan. The thermodynamics of such a mesoporous reservoir with a pore diameter fixed at 2 mm are equivalent to the previously proposed surface ocean. This work shows that if such a porous hydrocarbon/nitrogen reservoir contains a significant volume within pores of nanometer dimension, then the modification of methane thermodynamics through the Kelvin effect significantly changes the reservoir-atmosphere interaction on Titan. We measure the capillary uptake of pure methane liquid into well-characterized ice grown at temperatures and

partial pressures simulating Titan at present. From these results, a model of a nanoporous reservoir housing the methane/ethane/nitrogen solution is constructed. Modifying a previously published analytical climate model for Titan to incorporate such a reservoir, we report calculations of the coupled reservoir-atmosphere system. This modification leads to a very sensitive dependence of the state of the atmosphere on the total reservoir volume. In particular, assuming a nanoporous volume distribution reported here, calculations show that such a reservoir must comprise at least 1.5% of a 10 km deep ice regolith in order to yield a climate consistent with that currently observed on Titan. In contrast to the mesoporous model, it reduces the sensitivity of the system to temperature, mitigating the runaway greenhouse effect observed in the previous model. Furthermore, the model indicates that in a past epoch a surface ocean subsided into the reservoir as the total surface liquid volume reduced with time.

Mayo de 2007
The 2004 Parkfield earthquake: Test of the electromagnetic precursor hypothesis
Authors: Stephen K. Park and William Dalrymple et al
Link: Click Here

Abstract
A controversy has existed for 30 years concerning the possibility of earthquake prediction using electromagnetic precursors. Long-term electromagnetic monitoring prior to, during, and after the M6.0 earthquake at Parkfield, California, on 28 September 2004 now provides a definitive test of this hypothesis. During the earthquake our instruments recorded clearly documented electrical signals from an earthquake: impulsive changes of up to +2.5 mV at an electrode located 250 m from the rupture zone followed by a

transient decay lasting at least 3 1/4 hours. Similar signals (-1.5 mV) were recorded with the only two M > 5 aftershocks, but their transient decays lasted for only ~17 min. These signals are unambiguously a result of the earthquake and can be useful in studies of fluid flow in faults, a major current topic in earthquake physics. Patterns in the distribution of transient voltages differ from static stress changes, but are consistent with a small coseismic pressure drop on the fault and fluid flow inward. Signals at a distance of 1 km from the fault are less than 1 mV and are statistically no different from zero in our analysis. Transient voltages for the aftershocks have opposite polarity from those recorded with the main shock and are again maximal close to the fault, suggesting subsequent outward flow from the fault. Finally, there is no evidence of any precursory signal, strongly calling into question previously suggested signals preceding smaller or more distant earthquakes.

Abril de 2007
Teoría de Periodo de Recurrencia de Terremotos
Autores: Alexander Saichev y Didier Sornette
Link: Click Aquí

Abstract
La estadística del periodo de recurrencia realizada en grandes áreas ha probado obedecer leyes universales de escala, ya sea que ésta se realice en áreas únicas homógeneas o promediando sobre múltiples regiones. Esas leyes de escala unificadas son caracterizadas por leyes de potencia asintóticas. Por otro lado, Molchan (2005) ha presentado una prueba matemática donde se afirma que si tal ley universal existe, ésta debe ser necesariamente exponencial, lo que está en obvia contradicción con los hechos. Nosotros primeramente generalizamos el argumento de Molchan para mostrar que se puede encontrar una ley aproximadamente unificada y que a la vez es compatible con los hechos cuando incorporamos la Ley de Inducción de Terremotos de Omori-Utsu. A continuación desarrollamos la teoría de la estadística del tiempo entre terremotos en el marco del modelo de inducción sísmica de la Secuencia Post-Terremoto Tipo-Epidémica (ETAS) y mostramos que las observaciones empíricas pueden ser totalmente explicadas. Nuestra expresión teórica ajusta perfectamente la estadística empírica sobre el rango completo de periodos de recurrencia, contabilizado para distintos tipos de regimenes mediante la única utilización de la física de la inducción, la que es cuantificada mediante la Ley de Omori-Utsu. La descripción estadística del periodo de recurrencia sobre múltiples regiones requiere una adicional y sutil derivación matemática a partir de los mapas de la geometría fractal de los epicentros de los terremotos considerando las razones sísmicas promedio en múltiples regiones. Esto nos lleva a predicciones en excelente acuerdo con los datos empíricos para razonables valores de la dimensión fractal d ˜ 1.8, de la razón promedio del agrupamiento n ˜ 0.9 y el exponente de productividad a ˜ 0.9 veces el valor b de la Ley de Gutenberg-Richter.

Theory of earthquake recurrence times
Authors: Alexander Saichev and Didier Sornette
Link: Click Here

Abstract
The statistics of recurrence times in broad areas have been reported to obey universal scaling laws, both for single homogeneous regions and when averaged over multiple regions. These unified scaling laws are characterized by intermediate power law asymptotics. On the other hand, Molchan (2005) has presented a mathematical proof that if such a universal law exists, it is necessarily an exponential, in obvious contradiction with the data. First, we generalize Molchan's argument to show that an approximate unified law can be found which is compatible with the empirical observations when incorporating the impact of the Omori-Utsu law of earthquake triggering. We then develop the theory of the statistics of interevent times in the framework of the Epidemic-Type Aftershock Sequence (ETAS) model of triggered seismicity and show that the empirical observations can be fully explained. Our theoretical expression well fits the empirical statistics over the whole range of recurrence times, accounting for different regimes by using only the physics of triggering quantified by the Omori-Utsu law. The description of the statistics of recurrence times over multiple regions requires an additional subtle statistical derivation that maps the fractal geometry of earthquake epicenters onto the distribution of the average seismic rates in multiple regions. This yields a prediction in excellent agreement with the empirical data for reasonable values of the fractal dimension d ˜ 1.8, the average clustering ratio n ˜ 0.9, and the productivity exponent a ˜ 0.9 times the b value of the Gutenberg-Richter law.

Abril de 2007
Statistics of fracture strength and fluid-induced microseismicity
Authors: Elmar Rothert and Serge A. Shapiro
Link: Click Here

Abstract
In this paper we develop an approach for estimating the strength of rocks by analyzing fluid-induced microseismicity. The strength corresponds to the value of critical pressure in the pore space that must be exceeded in order to activate preexisting fractures, i.e., to trigger earthquakes. We assume that during hydraulic injection experiments in boreholes, microseismicity is mainly triggered by a diffusive process of pore pressure perturbation. An analytical solution can be applied to find time-dependent pore pressure perturbations in rocks caused by fluid injections. Characteristics of the spatiotemporal evolution of microseismic clouds can be then used to estimate minimum and maximum pressures necessary to trigger earthquakes. Moreover, we present a method for

reconstruction of the full spectrum of rock strength; that is, we show how to estimate the probability density function of the critical pressure. We verify the approach using numerical data and apply it to real data of injection-induced microseismicity from two Hot Dry Rock tests in crystalline rocks and one hydraulic fracturing experiment in a sedimentary rock. The results show that such an analysis of microseismicity is able to provide us with a completely new feature of natural fractured rocks in situ, namely, the statistics of their strength. We find that very low critical pressures, in the range 0.001-1 MPa, characterize the strength of preexisting cracks for all data sets analyzed. The range of critical pressures is broadly distributed within 3 to 4 orders of magnitude. However, the probability density functions of critical pressure change very quickly from the zero level to significant values and from significant values to the zero level at the lower and upper limits of their nonvanishing value ranges, respectively. The lower bound of critical pressure is possibly defined by the magnitude of tidal-induced stresses, permanently occurring and relaxing in the Earth.

Abril de 2007
Earthquake source characteristics from dynamic rupture with constrained stochastic fault stress
Authors: J. Ripperger and J.-P. Ampuero et al
Link: Click Here

Abstract
One of the challenging tasks in predicting near-source ground motion for future earthquakes is to anticipate the spatiotemporal evolution of the rupture process. The final size of an event but also its temporal properties (propagation velocity, slip velocity) depend on the distribution of shear stress on the fault plane. Though these incipient stresses are not known for future earthquakes, they might be sufficiently well characterized in a stochastic sense. We examine the evolution of dynamic rupture in numerical models of a fault

subjected to heterogeneous stress fields with varying statistical properties. By exploring the parameter space of the stochastic stress characterization for a large number of random realizations we relate generalized properties of the resulting events to the stochastic stress parameters. The nucleation zone of the simulated earthquake ruptures in general has a complex shape, but its average size is found to be independent of the stress field parameterization and is determined only by the material parameters and the friction law. Furthermore, we observe a sharp transition in event size from small to system-wide events, governed mainly by the standard deviation of the stress field. A simplified model based on fracture mechanics is able to explain this transition. Finally, we find that the macroscopic rupture parameters (e.g., moment, moment rate, seismic energy) of our catalog of model quakes are generally consistent with observational data.

Abril de 2007
Seasonal modulation of seismicity in the Himalaya of Nepal
Authors: L. Bollinger, F. Perrier et al
Link: Click Here

Abstract
For the period 1995-2000, the Nepal seismic network recorded 37 ± 8% fewer earthquakes in the summer than in the winter; for local magnitudes ML > 2 to ML > 4 the percentage increases from 31% to 63% respectively. We show the probability of observing this by chance is

less than 1%. We find that most surface loading phenomena are either too small, or have the wrong polarity to enhance winter seismicity. We consider enhanced Coulomb failure caused by a pore-pressure increase at seismogenic depths as a possible mechanism. For this to enhance winter seismicity, however, we find that fluid diffusion following surface hydraulic loading would need to be associated with a six-month phase lag, which we consider to be possible, though unlikely. We favor instead the suppression of summer seismicity caused by stress-loading accompanying monsoon rains in the Ganges and northern India, a mechanism that is discussed in a companion article.

Abril de 2007
Sensitivity of U.S. surface ozone to future emissions and climate changes
Authors: Zhining Tao and Allen Williams et al
Link: Click Here

Abstract
The relative contributions of projected future emissions and climate changes to U.S. surface ozone concentrations are investigated focusing on California, the Midwest, the Northeast, and Texas. By 2050 regional average ozone concentrations

increase by 2-15% under the IPCC SRES A1Fi ("dirty") scenario, and decrease by 4-12% under the B1 (relatively "clean") scenario. However, the magnitudes of ozone changes differ significantly between major metropolitan and rural areas. These ozone changes are dominated by the emissions changes in 61% area of the contiguous U.S. under the B1 scenario, but are largely determined by the projected climate changes in 46% area under the A1Fi scenario. In the ozone responses to climate changes, the biogenic emissions changes contribute strongly over the Northeast, moderately in the Midwest, and negligibly in other regions.

Abril de 2007
Bounds on the viscosity coefficient of continental lithosphere from removal of mantle lithosphere beneath the Altiplano and Eastern Cordillera
Authors: Peter Molnar and Carmala N. Garzione
Link: Click Here

Abstract
The rapid rise of the central Andean plateau between ~10 and 6.8 Ma implies that mantle lithosphere, including eclogitized lower crust, was removed from beneath the region in that time interval; we infer from that removal that the average viscosity coefficient of mantle lithosphere was quite low when removal occurred. Using scaling laws for the growth of perturbations to the thickness of a dense layer over an inviscid substratum (Rayleigh-Taylor instability), we place bounds on the average viscosity coefficient for central Andean lithosphere. When compared with laboratory measurements of flow laws for olivine and eclogite, the allowed range of viscosity coefficients yields bounds on the temperature of ~500-800°C at the Moho beneath this region and suggests that mean stresses across mantle lithosphere during continental deformation are less than ~50 MPa. This range of temperature is comparable with, if a slightly lower than, we might expect for lithosphere approximately doubled in thickness and not yet equilibrated with the doubled crustal radioactivity. The mean deviatoric stress is comparable to that associated with stresses that drive plates and hence shows that lithospheric material is not too strong to prevent removal of its mantle part.

Abril de 2007
Flank instability on Mount Etna: Radon, radar interferometry, and geodetic data from the southwestern boundary of the unstable sector
Authors: Marco Neri, Francesco Guglielmino et al
Link: Click Here

Abstract
Understanding Etnean flank instability is hampered by uncertainties over its western boundary. Accordingly, we combine soil radon emission, interferometric synthetic aperture radar (InSAR), and electronic distance measurement (EDM) data to study the Ragalna fault system (RFS) on the SW flank of the volcano. Valuable synergy developed between our differing techniques, producing consistent results and serving as a model for other studies of partly obscured active faults. The RFS, limited in its surface expression, is revealed as a complex interlinked structure ~14 km long that extends from the edifice base toward the area of summit rifting, possibly linking northeastward to the Pernicana fault system (PFS) to define the unstable sector. Short-term deformation rates on the RFS from InSAR data reach ~7 mm yr-1 in the satellite line of sight on the upslope segment and ~5 mm yr-1 on the prominent central segment. Combining this with EDM data confirms the central segment of the RFS as a dextral transtensive structure, with strike-slip and dip-slip components of ~3.4 and ~3.7 mm yr-1, respectively. We measured thoron (220Rn, half-life 56 s) as well as radon, and probably because of its limited diffusion range, this appears to be a more sensitive but previously unexploited isotope for pinpointing active near-surface faults. Contrasting activity of the PFS and RFS reinforces proposals that the instability they bound is divided into at least three subsectors by intervening faults, while, in section, fault-associated basal detachments also form a nested pattern. Complex temporal and spatial movement interactions are expected between these structural components of the unstable sector.

Abril de 2007
Space weather at Venus and its potential consequences for atmosphere evolution
Authors: J. G. Luhmann, W. T. Kasprzak et al
Link: Click Here

Abstract
Space weather storms at the Earth are dominated by the magnetosphere's response to coronal mass ejections, or CMEs, whose disturbances propagate through the solar wind to 1 AU, and to a lesser extent by the pressure ridges associated with the interactions of solar wind streams with different properties. The interplanetary signatures of these events include high solar wind dynamic pressures and high interplanetary magnetic fields, the same parameters that cause compression and/or magnetization of the ionosphere at weakly magnetized Venus. The nature of Venus' response to these events is also

expected to include increased atmosphere erosion by the solar wind interaction, a matter of potential interest for historical extrapolations of atmosphere escape. We consider the possible evidence for enhanced escape during these disturbances in the data from the Pioneer Venus Orbiter (PVO). Available magnetometer and plasma analyzer results are used to identify periods of interest and the interplanetary characteristics of the disturbed periods. It is found that the suprathermal (>36 eV) ion measurements from the Pioneer Venus neutral mass spectrometer (ONMS) provide convincing evidence of the related enhanced escape of mostly O+ suprathermal ions. The results provide the first direct demonstration that space weather has played an important role in Venus (and other weakly magnetized planet) atmosphere loss through time. It is important to pursue study of these effects with the new measurements from Venus Express and, in light of the approaching solar activity cycle rising phase, extremely timely.

Abril de 2007
Hydrometeorological application of a microwave link: 2. Precipitation
Authors: H. Leijnse and R. Uijlenhoet et al
Link: Click Here

Abstract
A method to estimate areal evaporation using a microwave link (radio wave scintillometer) in combination with an energy budget constraint is proposed. This radio wave scintillometry-energy budget method (RWS-EBM) is evaluated for its applicability in different meteorological conditions and for its sensitivity to various variables (the structure parameter of the refractive index of air C n 2, the total available energy R n - G, the wind velocity u, the effective average vegetation height h 0, and the correlation coefficient between the

temperature and humidity fluctuations r TQ ). The method is shown to be best suited for use in wet to moderately dry conditions, where the latent heat flux is at least a third of the total available energy (i.e., Bowen ratio =2). It is important to accurately measure the total available energy and the wind velocity as the RWS-EBM is most sensitive to these variables. The Flevoland field experiment has provided the data, obtained with a 27-GHz radio wave scintillometer (over 2.2 km), a large-aperture scintillometer (also 2.2 km), and four eddy covariance systems, which are used to test the RWS-EBM. Comparing 92 daytime measurements (30-min intervals) of the evaporation estimated using the RWS-EBM to that determined in alternative manners (eddy covariance and two-wavelength scintillometry) leads to the conclusion that the method provides consistent estimates (coefficient of determination r 2 = 0.85 in both cases) under relatively wet conditions.

Abril de 2007
Improved finite-sample Hurst exponent estimates using rescaled range analysis
Author: Khaled H. Hamed
Link: Click Here

Abstract
Rescaled range analysis is one of the classical methods used for detecting and quantifying long-term dependence in time series. However, rescaled range analysis has been shown in several studies to give biased estimates of the Hurst exponent in finite samples. A new estimator based on a modified expression of the expected value of the rescaled range is proposed. A comparison of the modified estimator with other

alternatives shows that the modified estimator offers a great improvement over the classical rescaled range estimator in terms of bias and root mean square error, which makes it comparable to some of the leading estimators. The application of the proposed modified rescaled range estimator to a group of temperature, rainfall, river flow, and tree-ring time series in the Midwest USA demonstrates the extent to which classical rescaled range analysis can give misleading results. Based on a statistical test of significance, the number of time series exhibiting the Hurst effect is reduced from 36 to only 11 out of 56 temperature series, and from 23 to 7 out of 60 rainfall series. On the other hand, the number of time series exhibiting the Hurst effect marginally increased from 8 to 9 out of 49 river flow series, and from 49 to 54 out of 88 tree-ring series.

Abril de 2007
Analysis of glacial earthquakes
Authors: Victor C. Tsai and Göran Ekström
Link: Click Here

Abstract
In 2003, Ekström et al. reported on the detection of a new class of earthquakes that occur in glaciated regions, with the vast majority being in Greenland. The events have a characteristic radiation pattern and lack the high-frequency content typical of tectonic earthquakes. It was proposed that the events correspond to large and sudden sliding motion of glaciers. Here we present an analysis of all 184 such events detected in Greenland between 1993 and 2005. Fitting the teleseismic long-period surface waves to a landslide model of the source, we obtain improved locations, timing, force amplitudes, and force directions. After relocation, the events cluster into seven regions,

all of which correspond to regions of very high ice flow and most of which are named outlet glaciers. These regions are Daugaard Jensen Glacier, Kangerdlugssuaq Glacier, Helheim Glacier, the southeast Greenland glaciers, the northwest Greenland glaciers, Rinks Isbrae, and Jakobshavn Isbrae. Event amplitudes range from 0.1 to 2.0 × 1014 kg m. Force directions are consistent with sliding in the direction of glacial flow over a period of about 50 s. Each region has a different temporal distribution of events. All glaciers are more productive in the summer, but have their peak activity in different months. Over the study period, Kangerdlugssuaq has had a constant number of events each year, whereas Jakobshavn had most events in 1998-1999, and the number of events in Helheim and the northwest Greenland glaciers has increased substantially between 1993 and 2005. The size distribution of events in Kangerdlugssuaq is peaked above the detection threshold, suggesting that glacial earthquakes have a characteristic size.

Abril de 2007
Similar microearthquakes observed in western Nagano, Japan, and implications for rupture mechanics
Authors: Xin Cheng, Fenglin Niu et al
Link: Click Here

Abstract
We have applied a waveform cross correlation technique to study the similarity and the repeatability of more than 21,000 microearthquakes (0 < M < 4.5) in the aftershock zone of the 1984 western Nagano earthquake in central Japan. We find that the seismicity in this particular intraplate fault essentially consists of no repeating earthquakes that occurred on the same patch of the fault in a quasiperiodic manner in the study period between 1995 and 2001. On the other hand, we identify a total of 278 doublets and 62 multiplets (807 events) that occurred consecutively within seconds to days. On the

basis of the relative arrival times of the P and S waves, we have obtained precise relative locations of these consecutive events with an error between several meters to a few tens of meters. There is a clear lower bound on the distances measured between these consecutive events and the lower bound appears to be proportional to the size of the first events. This feature is consistent with what Rubin and Gillard [2000] have observed near the San Juan Bautista section of the San Andreas Fault. Shear stress increases at the edge of an earthquake rupture, and the rupture edge becomes the most likely place where the second events are initiated. The observed minimum distance thus reflects the rupture size of the first events. The minimum distance corresponds to the rupture size calculated from a circular fault model with a stress drop of 10 MPa. We found that using different time windows results in a slight difference in the delay time estimates and the subsequent projection locations, which may reflect the finite size nature of earthquake ruptures.

Abril de 2007
Correlaciones entre terremotos y grandes erupciones volcánicas de lodos
Autores: R. Mellors, D. Kilb et al
Link: Click Aquí

Abstract
Nosotros examinamos la relación entre la potencial inducción de grandes terremotos a partir de las erupciones volcánicas del tipo metano y lodo. Nuestra base de datos consistió en un catálogo de 191 años (1810-2001) abarcando erupciones de 77 volcanes en Azerbaijan, en el Asia Central, los que fueron suplementados con reportes de erupciones volcánicas de lodo en Japón, Rumania, Pakistán y las islas Andaman. Nosotros comparamos la ocurrencia de terremotos históricos regionales (M > 5) con la ocurrencia de erupciones volcánicas de lodo en Azerbaijan y encontramos que el número de pares terremotos/erupciones ocurridos el mismo día es significativamente mayor que lo esperado si las erupciones y terremotos fuesen procesos independientes (los que son modelados por la Ley de Poisson). La correlación temporal entre terremotos y erupciones es más pronunciada para terremotos cercanos (dentro de un radio de ~ 100 Km) que producen intensidades sísmicas iguales o superiores a 6 Mercalli en la ubicación de la erupción. Esto puede significar que la relación magnitud/distancia para inducción de terremotos observada en Azerbaijan es consistente con la relación documentada entre terremotos y erupciones volcánicas de lodo. También encontramos una débil correlación que intensifica el número de erupciones volcánicas de lodo que ocurren dentro de un año después del gran terremoto. La distribución de erupciones anuales se aproxima fuertemente a procesos tipo Poisson. Sin embargo, el tiempo de reposo parece favorecer algún tipo de razón de falla no homogénea, lo que implica que los volcanes necesitan de cierto tiempo para recargarse después de una erupción. La inducción de la erupción parece resultar del paso de una onda sísmica, aunque el mecanismo preciso no es claro.

Abril de 2007
Correlations between earthquakes and large mud volcano eruptions
Authors: R. Mellors, D. Kilb et al
Link: Click Here

Abstract
We examine the potential triggering relationship between large earthquakes and methane mud volcano eruptions. Our data set consists of a 191-year catalog (1810-2001) of eruptions from 77 volcanoes in Azerbaijan, central Asia, supplemented with reports from mud volcano eruptions in Japan, Romania, Pakistan, and the Andaman Islands. We compare the occurrence of historical regional earthquakes (M > 5) with the occurrence of Azerbaijan mud volcano eruptions and find that the number of same-day earthquake/eruption pairs is significantly higher than expected if the eruptions and earthquakes are independent Poisson processes. The temporal correlation between earthquakes and eruptions is most pronounced for nearby earthquakes (within ~100 km) that produce seismic intensities of Mercalli 6 or greater at the location of the mud volcano. This assumed magnitude/distance relationship for triggering observed in the Azerbaijan data is consistent with documented earthquake-induced mud volcano eruptions elsewhere. We also find a weak correlation that heightened numbers of mud volcano eruptions occur within 1 year after large earthquakes. The distribution of yearly eruptions roughly approximates a Poisson process, although the repose times somewhat favor a nonhomogenous failure rate, which implies that the volcanoes require some time after eruption to recharge. The volcanic triggering likely results from some aspect of the seismic wave's passage, but the precise mechanism remains unclear.

Abril de 2007
Magnetic field visualization of magnetic minerals and grain boundary regions using magneto-optical imaging
Authors: Alexander V. Bobyl, Yuri Y. Podladchikov et al
Link: Click Here

Abstract
Magneto-optical imaging based on the Faraday effect has been used to characterize magnetic minerals embedded in a nonmagnetic matrix. We have studied magnetite grains and magnetite-magnetite grain boundary regions in samples of skarns and serpentinites. Distributions of the

remanent magnetic field were measured across at the surface of polished thin sections kept at room temperature. The magneto-optical images resolve directly magnetic structures on length scales ranging from millimeter down to micrometer, thereby revealing the shape and arrangement of the magnetite grains and allow determination of the grain magnetization. We find that (1) for the skarns the intergrain interactions do not affect the magnetic properties of magnetite grains within 0.6-60 µm of each other, while the saturation remanence decreases weakly with increasing grain size from 40 µm to 0.6 mm, and (2) for the serpentinites the magnetic properties of the stripes are size-dependent due to variations in chemical composition.

Abril de 2007
Earthquake swarms driven by aseismic creep in the Salton Trough, California
Authors: R. B. Lohman and J. J. McGuire
Link: Click Here

Abstract
In late August 2005, a swarm of more than a thousand earthquakes between magnitudes 1 and 5.1 occurred at the Obsidian Buttes, near the southern San Andreas Fault. This swarm provides the best opportunity to date to assess the mechanisms driving seismic swarms along

transform plate boundaries. The recorded seismicity can only explain 20% of the geodetically observed deformation, implying that shallow, aseismic fault slip was the primary process driving the Obsidian Buttes swarm. Models of earthquake triggering by aseismic creep can explain both the time history of seismic activity associated with the 2005 swarm and the ~1 km/h migration velocity exhibited by this and several other Salton Trough earthquake swarms. A combination of earthquake triggering models and denser geodetic data should enable significant improvements in time-dependent forecasts of seismic hazard in the key days to hours before significant earthquakes in the Salton Trough.

Abril de 2007
Effects of topography on the spin-up of a Venus atmospheric model
Authors: A. Herrnstein and T. E. Dowling
Link: Click Here

Abstract
We study how topography affects the spin-up from rest of a model of the atmosphere of Venus. The simulations are performed with the EPIC model using its isentropic, terrain-following hybrid vertical coordinate, and are forced with the Newtonian-cooling profile used to achieve superrotation in a Venus model with no topography by Lee et al. (2005). We are able to

reproduce their results with our model, which was developed independently and uses a different vertical coordinate. Both groups use a horizontal resolution of 5°, which is dictated by the need for reasonable computer runtime and is not a claim of numerical convergence. We find that the addition of topography substantially changes both the evolution and end state of the model's spin-up: the magnitude of the superrotation is diminished from 55 ms-1 to 35 ms-1, and it reaches steady state faster, in a few years instead of a few decades. A large, stationary eddy associated with Ishtar Terra forms that has a local horizontal temperature anomaly of order 2 K at the 0.7 bar level; such a feature may be observable in high-resolution infrared images.

Abril de 2007
Fractal topography and subsurface water flows from fluvial bedforms to the continental shield
Authors: Anders Wörman, Aaron I. Packman et al
Link: Click Here

Abstract
Surface-subsurface flow interactions are critical to a wide range of geochemical and ecological processes and to the fate of contaminants in freshwater environments. Fractal scaling relationships have been found in distributions of both land surface topography and solute efflux from watersheds, but the linkage between those

observations has not been realized. We show that the fractal nature of the land surface in fluvial and glacial systems produces fractal distributions of recharge, discharge, and associated subsurface flow patterns. Interfacial flux tends to be dominated by small-scale features while the flux through deeper subsurface flow paths tends to be controlled by larger-scale features. This scaling behavior holds at all scales, from small fluvial bedforms (tens of centimeters) to the continental landscape (hundreds of kilometers). The fractal nature of surface-subsurface water fluxes yields a single scale-independent distribution of subsurface water residence times for both near-surface fluvial systems and deeper hydrogeological flows.

Abril de 2007
Evidence that floodwaters filled and overflowed Capri Chasma, Mars
Authors: Neil M. Coleman, Cynthia L. Dinwiddie et al
Link: Click Here

Abstract
Erosional features on a plateau in the Valles Marineris provide evidence that a lake filled Capri Chasma until it overflowed its eastern rim, carving two large spillover channels. The floodwaters surged into the adjacent lowlands of ancestral Ganges Chasma and eastern Eos Chasma. The

channel floor elevation and depth of incision at Daga Vallis reveal that a 1200-m-deep lake water column was drained by the flooding. The width, depth, and steep energy slope (~0.02) of Daga Vallis and the existence of several cataracts, including a 500-m-high dry falls, demonstrate the power of the floodwaters. We estimate a possible peak discharge rate of 1-6 × 108 m3·s-1. The catastrophic flows may have been triggered by the collapse of topographic barriers in eastern Coprates Chasma, providing a gateway for lake waters in the central Valles Marineris to pour eastward toward Capri Chasma and the lowlands beyond. These may have been among the earliest flows in Simud-Tiu Valles, contributing discharge to a possible sea in the northern plains of Mars.

Abril de 2007
Probabilistic envelope curves for design flood estimation at ungauged sites
Author: Attilio Castellarin
Link: Click Aquí

Abstract
Castellarin et al. (2005) introduced the probabilistic regional envelope curves (PRECs) of flood flows and formulated an empirical estimator of the recurrence interval T associated with the PREC for M synthetic cross-correlated and concurrent sequences of annual maxima with length n. In principle, a PREC can be used to estimate the T-year flood (design flood) for any basin in a given region as a function of the catchment area alone. A number of fundamental issues hamper the construction of PRECs and the estimation of their T for historical regional samples of flood flows. Some of these issues are

addressed by this study, which presents an algorithm for the application of the empirical estimator of T to historical annual maximum series of unequal length and assesses how the selection of a particular cross-correlation formula and plotting position affects the accuracy of design flood estimates retrieved from PRECs (PREC flood quantiles). The reliability of the proposed algorithm and the accuracy of the flood quantiles are discussed through a comprehensive cross validation. The results show that the reliability of PREC flood quantiles for ungauged sites is comparable with the reliability of regional estimates produced by the application of the index flood approach. The results of the cross-validation also show that the selection of a particular correlation formula or plotting position may produce biased PREC flood quantiles, and suggest that the accuracy of the approach could be significantly enhanced by relaxing the assumption that the design flood scales with catchment area alone.

Abril de 2007
Using lunar boulders to distinguish primary from distant secondary impact craters
Authors: Gwendolyn D. Bart and H. J. Melosh
Link: Click Aquí

Abstract
A high-resolution study of 18 lunar craters, including both primary and distant secondary craters, shows that the secondary craters produce larger ejecta fragments at a given crater size than do the primary craters. The maximum boulder diameter (B) increases with crater size (D)

according to the power law B = KD 2/3; for primary craters, when B and D are in meters, K is 0.29, whereas for secondary craters, we find that K is 0.46 (60% larger). Next we show that impact fracture theory predicts that secondary craters, because of their lower impact velocity, will produce larger ejecta fragments than primary craters. This result provides an opportunity for distinguishing between primary and secondary craters in high resolution planetary images. The ability to identify distant secondary craters will help constrain primary production rates of small craters and improve surface age determination of small areas based on small crater counts.

Abril de 2007
A theoretical model of the inner proton radiation belt
Authors: R. S. Selesnick, M. D. Looper et al
Link: Click Aquí

Abstract
A theoretical calculation provides inner radiation belt proton intensities as a function of time and of the three adiabatic invariants, M, K, and L, in the kinetic energy range from ~10 MeV to ~4 GeV and the L range from 1.1 to 2.4. Long residence times for trapped protons of up to several thousand years require similarly long input time series for the geomagnetic field, solar activity, and solar proton fluences. Additional inputs include galactic cosmic ray spectra, nuclear scattering cross sections, and the neutral and plasma densities in the atmosphere, ionosphere, and plasmasphere. Trapped proton sources are cosmic ray albedo

neutron decay (CRAND), calculated from a Monte Carlo particle transport simulation, and solar proton injection using a derived empirical injection efficiency that is ~10-4 at 10 MeV. Radial diffusion provides inward transport of injected solar protons. Calculated intensities at energies 100 MeV and for L 1.3 are dominated by solar protons, CRAND being the dominant source otherwise. Losses are by ionization of the neutral atmosphere, energy transfer to plasma electrons, and inelastic nuclear scattering. Numerical trajectory tracing determines trapping limits and drift shell averages of the albedo neutron intensity and of neutral and plasma densities for loss rate calculations. Geomagnetic secular variations cause adiabatic energy and drift shell changes. Intensities are greater than they would be in a constant geomagnetic field by factors up to ~10, a result of long proton residence times and the presently decreasing geomagnetic dipole moment.

Abril de 2007
Interpreting in situ pressure and hydraulic properties with borehole penetrometers in ocean drilling: DVTPP and Piezoprobe deployments at southern Hydrate Ridge, offshore Oregon
Authors: H. Long, P. B. Flemings et al
Link: Click Aquí

Abstract
Two borehole penetrometers, Fugro-McClelland's Piezoprobe and the Ocean Drilling Program's (ODP) DVTPP, were deployed 50 m below seafloor at Site 1244 on ODP Leg 204 to measure formation pressure at southern Hydrate Ridge, offshore Oregon. Pore pressure is interpreted to be hydrostatic and the sediment's coefficient of consolidation is interpreted to lie between 6.92 and 7.8 × 10-7 m2/s, which is in approximate

agreement with laboratory measurements. The Piezoprobe pressure reaches 90% of dissipation 14 times sooner than the DVTPP. The observed and modeled pore pressure responses illustrate how penetrometer geometry impacts our ability to interpret in situ properties and demonstrate under what conditions these tools can be effectively used. Because of its narrow tip, the Piezoprobe disturbs a narrower zone than the DVTPP does. This generates a narrower zone of pressure increase around the piezoprobe, which dissipates much faster than the DVTPP. As consolidation proceeds, pressure dissipation of the Piezoprobe is retarded and forms a "bench," or flat spot, on the dissipation curve. Owing to its distinct two-radius geometry, it is possible to apply a consistent method to estimate in situ pressure from partial dissipation record based on the position of the "bench."

Abril de 2007
Groundwater noble gas, age, and temperature signatures in an Alpine watershed: Valuable tools in conceptual model development
Authors: Andrew H. Manning, Jonathan Saul Caine et al
Link: Click Here

Abstract
Bedrock groundwater in alpine watersheds is poorly understood, mainly because of a scarcity of wells in alpine settings. Groundwater noble gas, age, and temperature data were collected from springs and wells with depths of 3-342 m in Handcart Gulch, an alpine watershed in Colorado. Temperature profiles indicate active groundwater circulation to a maximum depth (aquifer thickness) of about 200 m, or about 150 m below the water table. Dissolved noble gas data show unusually high excess air concentrations (>0.02 cm3 STP/g, ?Ne > 170%) in the bedrock,

consistent with unusually large seasonal water table fluctuations (up to 50 m) observed in the upper part of the watershed. Apparent 3H/3He ages are positively correlated with sample depth and excess air concentrations. Integrated samples were collected from artesian bedrock wells near the trunk stream and are assumed to approximate flow-weighted samples reflecting bedrock aquifer mean residence times. Exponential mean ages for these integrated samples are remarkably consistent along the stream, four of five being from 8 to 11 years. The tracer data in combination with other hydrologic and geologic data support a relatively simple conceptual model of groundwater flow in the watershed in which (1) permeability is primarily a function of depth; (2) water table fluctuations increase with distance from the stream; and (3) recharge, aquifer thickness, and porosity are relatively uniform throughout the watershed in spite of the geological complexity of the Proterozoic crystalline rocks that underlie it.

Marzo de 2007
Chemical composition of Earth's primitive mantle and its variance: 1. Method and results
Authors: Tanya Lyubetskaya and Jun Korenaga
Link: Click Aquí

Abstract
We present a new statistical method to construct a model for the chemical composition of Earth's primitive mantle along with its variance. Earth's primitive mantle is located on the melting trend exhibited by the global compilation of mantle peridotites, using cosmochemical constraints on the relative abundances of refractory lithophile elements (RLE). This so-called pyrolite approach involves the least amount of assumptions, thereby being probably most satisfactory compared to other approaches. Its previous implementations, however, suffer from questionable statistical treatment of noisy geochemical data, leaving the uncertainty of model composition poorly

quantified. In order to properly take into account how scatters in peridotite data affect this geochemical inference, we combine the following statistical techniques: (1) modeling a nonlinear melting trend in the multidimensional compositional space through the principal component analysis, (2) determining the primitive mantle composition on the melting trend by simultaneously imposing all of cosmochemical constraints with least squares, and (3) mapping scatters in original data into the variance of the final model through the bootstrap resampling technique. Whereas our model is similar to previous models in terms of Mg, Si, and Fe abundances, the RLE contents are at ~2.16 ± 0.37 times the CI chondrite concentration, which is lower than most of previous estimates. The new model is depleted by >20% in a number of incompatible elements including heat-producing elements, U, Th, and K, and this depleted nature is further amplified (up to 60%) in terms of predicted composition for the present-day mantle.

Marzo de 2007
First-principles study of density, viscosity, and diffusion coefficients of liquid MgSiO3 at conditions of the Earth's deep mantle
Authors: J. T. K. Wan, T. S. Duffy et al
Link: Click Here

Abstract
Constant-pressure constant-temperature ab initio molecular dynamics simulations at high temperatures have been used to study MgSiO3, the major constituent of the Earth's lower mantle to conditions of the Earth's core-mantle boundary. The calculated equilibrium volumes and densities are compared with simulations using an orthorhombic perovskite configuration under the same conditions. For molten MgSiO3, we have determined the diffusion coefficients and shear viscosities at different thermodynamic conditions. Our results provide new constraints on the properties of molten MgSiO3 at conditions near

the core-mantle boundary. The volume of the liquid is greater than that of the solid throughout the pressure-temperature conditions examined, and the volume change on fusion ranges from 5% at 88 GPa and 3500 K to 2.9% at 120 GPa and 5000 K. Existing experimental constraints on solid-liquid partition coefficients for Fe suggest that Fe is preferentially partitioned into the liquid. Such enrichment of Fe increases the density of the liquid, thus allowing the possibility of negatively buoyant melts from (Mg,Fe)SiO3 perovskite compositions at deep lower mantle conditions for plausible values of solid-liquid partition coefficients for Fe. At 120 GPa and 4500-5000 K, the diffusion coefficient of liquid MgSiO3 is 2-3 × 10-5 cm 2/s and the diffusion rates of the different chemical species are similar. The shear viscosity is estimated using Zwanzig's formula to be 19-31 cP under these conditions. On the basis of our calculated diffusivities, MgSiO3 is above the glass transition temperature at 120 GPa and 4500 K.

Marzo de 2007
El Próximo Gran Terremoto
Author: Robert McCaffrey
Link: Click Aquí

Abstract
El geólogo Robert McCaffrey, del Rensselaer Polytechnic Institute, ha advertido en la revista Science de los riesgos del próximo Gran Terremoto. Tal como explica al respecto un comunicado del citado instituto, estos episodios geológicos están asociados a las así llamadas zonas de subducción, lugares donde se encuentran las placas tectónicas y donde una placa se desliza sobre la otra. Estos fenómenos ocurren cada 200 o 1.000 años. El problema es que, aunque los sismólogos intentan anticipar la evolución previsible de estas zonas de subducción, con la finalidad de prevenir el próximo Gran Terremoto, la realidad es que están muy lejos de conseguirlo: el terremoto Sumatra-Andaman, de diciembre de 2004 (más de 250.000 muertos), estaba entre las zonas con más bajas posibilidades de ocurrencia establecidas por los sismólogos. Eso quiere decir que cualquier zona de subducción puede darnos la sorpresa en cualquier momento, ya sea en el cinturón de fuego del Pacífico, en el Mar de China Meridional o en la placa de Nazca, bajo la cordillera de los Andes (...) McCafrey plantea la necesidad de elaborar políticas de prevención en todas ellas porque muchas están densamente pobladas y los impactos previsibles pueden ser catastróficos.

Eduardo Martínez - www.tendencias21.net
Tema relacionado:
Cómo predecir terremotos - 1
Cómo predecir terremotos - 2

The Next Great Earthquake
Robert McCaffrey

The 2004 Sumatra-Andaman earthquake occurred at a surprising time and place; the lessons learned may help costal communities in the future.
The author is in the Department of Earth and Environmental Sciences, Rensselaer Polytechnic Institute, Troy, NY 12180, USA. E-mail: mccafr@rpi.edu

COMMUNIQUÉ
The 2004 Sumatra-Andaman earthquake and resulting tsunami are now infamous for the damage they caused, but at the time many scientists believed this area was unlikely to create a quake of such magnitude. In the March 23 issue of the journal Science, a geophysicist from Rensselaer Polytechnic Institute urges the public and policy makers to consider all subduction-type tectonic boundaries to be “locked, loaded, and dangerous.” | See More

Marzo de 2007
A comparison of annual vertical crustal displacements from GPS and Gravity Recovery and Climate Experiment (GRACE) over Europe
Authors: T. van Dam , J. Wahr et al
Link: Click here

Abstract
We compare approximately 3 years of GPS height residuals (with respect to the International Terrestrial Reference Frame) with predictions of vertical surface displacements derived from the Gravity Recovery and Climate Experiment (GRACE) gravity fields for stations in Europe. An annual signal fit to the residual monthly heights, corrected for atmospheric pressure and barotropic ocean loading effects, should primarily represent surface displacements due to

long-wavelength variations in water storage. A comparison of the annual height signal from GPS and GRACE over Europe indicates that at most sites, the annual signals do not agree in amplitude or phase. We find that unlike the annual signal predicted from GRACE, the annual signal in the GPS heights is not coherent over the region, displaying significant variability from site to site. Confidence in the GRACE data and the unlikely possibility of large-amplitude small-scale features in the load field not captured by the GRACE data leads us to conclude that some of the discrepancy between the GPS and GRACE observations is due to technique errors in the GPS data processing. This is evidenced by the fact that the disagreement between GPS and GRACE is largest at coastal sites, where mismodeling of the semidiurnal ocean tidal loading signal can result in spurious annual signals.

Marzo de 2007
Effect of shock on the magnetic properties of pyrrhotite, the Martian crust, and meteorites
Authors: Karin L. Louzada, Sarah T. Stewart et al
Link: Click here

Abstract
We performed planar shock recovery experiments on natural pyrrhotite at pressures up to 6.9 GPa. We find that high-field isothermal remanent magnetization in pyrrhotite is demagnetized up to 90% by shock due to preferential removal of low coercivity components of magnetization. Contrary

to static experiments, we do not observe complete demagnetization. Post shock permanent changes in magnetic properties include increasing saturation isothermal remanent magnetization, bulk coercivity and low-temperature memory, and changes in squareness of hysteresis. These changes are consistent with an increase in the volume fraction of single domain grains. The lack of magnetic anomalies over large Martian impact basins is not expected to be solely due to shock demagnetization of the crust. We find that pyrrhotite-bearing rocks and meteorites can retain records of Martian magnetic fields even if shocked to pressures approaching 7 GPa. However, some paleointensity techniques may underestimate this field.

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