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

Septiembre de 2006
Sulfur's impact on core evolution and magnetic field generation on Ganymede
Authors: Steven A. Hauck II , Jonathan M. Aurnou et al
Link: Click here

Abstract
Analysis of the melting relationships of potential core forming materials in Ganymede indicate that fluid motions, a requirement for a dynamo origin for the satellite's magnetic field, may be driven, in part, either by iron (Fe) "snow" forming below the core-mantle boundary or solid iron sulfide (FeS) floating upward from the deep core. Eutectic melting temperatures and eutectic sulfur contents in the binary Fe-FeS system decrease with increasing pressure within the interval of core pressures on Ganymede (<14 GPa). Comparison of melting temperatures to adiabatic temperature

gradients in the core suggests that solid iron is thermodynamically stable at shallow levels for bulk core compositions more iron-rich than eutectic (i.e., <21 wt % S). Calculations based on high-pressure solid-liquid phase relationships in the Fe-FeS system indicate that iron snow or floatation of solid iron sulfide, depending on whether the core composition is more or less iron-rich than eutectic, is an inevitable consequence of cooling Ganymede's core. These results are robust over a wide range of plausible three-layer internal structures and thermal evolution scenarios. For precipitation regimes that include Fe-snow, we present scaling arguments that give typical Rossby and magnetic Reynolds numbers consistent with dynamo action occurring in Ganymede's core. Furthermore, by applying recently derived scaling relationships relating magnetic field strength to buoyancy flux, we obtain estimates of surface magnetic field strength comparable with observed values.

Septiembre de 2006
Single-crystal elasticity and sound velocities of (Mg0.94Fe0.06)O ferropericlase to 20 GPa
Authors: Jennifer M. Jackson , Stanislav V. Sinogeikin et al
Link: Click here

Abstract
The single-crystal elastic properties of high-spin (Mg0.94Fe0.06)O ferropericlase were measured by Brillouin spectroscopy on a sample compressed to 20 GPa with diamond anvil cells using methanol-ethanol-water as a pressure-transmitting medium. At room pressure, the adiabatic bulk (K 0S) and shear (µ0S) moduli are K 0S = 163 ± 3 GPa and µ0S = 121 ± 2 GPa, in excellent agreement with ultrasonic results from the same bulk sample (Jacobsen et al., 2002). A fit to all our high-pressure Brillouin data using a third-order finite-strain equation of state yields the following pressure derivatives of the adiabatic bulk and shear moduli: K'0S = 3.9 ± 0.2 and µ'0S

= 2.1 ± 0.1. Within the uncertainties, we find that K 0S and K'0S of (Mg0.94Fe0.06)O are unchanged from MgO. However, µ0S and µ'0S of (Mg0.94Fe0.06)O are reduced by 8% and 11%, respectively. The aggregate compressional (V P) and shear (V S) wave velocities are reduced by 4% and 6%, respectively, as compared to MgO. The pressure dependence of the single-crystal elastic moduli and aggregate sound velocities is linear within the investigated pressure range. The elastic anisotropy of (Mg0.94Fe0.06)O is about 10% greater than that of MgO at ambient conditions. At the highest pressure obtained here, the elastic anisotropy of (Mg0.94Fe0.06)O is close to zero. On the basis of our measurements and earlier ultrasonic measurements, we find that the pressure derivatives of shear moduli obtained at room pressure for low iron concentrations (<20 mol% FeO) of high-spin ferropericlase are inconsistent with those inferred from the lower mantle PREM model.

Septiembre de 2006
Critical point theory of earthquakes: Observation of correlated and cooperative behavior on earthquake fault systems
Authors: Chien-chih Chen, John B. Rundle et al
Link: Click here

Abstract
The critical point theory for earthquakes was originally proposed to explain the scaling relations observed in earthquakes, including the Gutenberg-Richter frequency- magnitude relation and the Omori's law for aftershocks. In this model, main shocks, their foreshocks and aftershocks

are all associated with the formation of a correlated, cooperative spatial region with high stress. Until now, only indirect evidence of the existence of these correlated regions has been reported. Here in this paper we present observations and analyses that allow us to directly map the high stress, spatially correlated regions preceding four major earthquakes, i.e. the 1992 Landers (California), 1995 Kobe (Japan), 1999 Chi-Chi (Taiwan) and 1999 Hector Mine (California) earthquakes. We therefore conclude that the locations and extent of large main shocks and their immediate aftershocks can be determined from seismicity data taken prior to the main shocks, and provide additional evidence in support of the critical point theory for earthquakes.

Agosto de 2006
Hydrogeological insights at Stromboli volcano (Italy) from geoelectrical, temperature, and CO2 soil degassing investigations
Authors:A. Finizola , A. Revil et al
Link: Click here

Abstract
Finding the geometry of aquifers in an active volcano is important for evaluating the hazards associated with phreato-magmatic phenomena and incidentally to address the problem of water supply. A combination of electrical resistivity tomography (ERT), self-potential, CO2, and temperature measurements provides insights about the location and pattern of ground water

flow at Stromboli volcano. The measurements were conducted along a NE-SW profile across the island from Scari to Ginostra, crossing the summit (Pizzo) area. ERT data (electrode spacing 20 m, depth of penetration of ~200 m) shows the shallow architecture through the distribution of the resistivities. The hydrothermal system is characterized by low values of the resistivity (<50 O m) while the surrounding rocks are resistive (>2000 O m) except on the North-East flank of the volcano where a cold aquifer is detected at a depth of ~80 m (resistivity in the range 70-300 O m). CO2 and temperature measurements corroborate the delineation of the hydrothermal body in the summit part of the volcano while a negative self-potential anomaly underlines the position of the cold aquifer.

Agosto de 2006
Low Earth orbit satellite navigation errors and vertical total electron content in single-frequency GPS tracking
Authors: Miquel Garcia-Fernàndez and Oliver Montenbruck
Link: Click here

Abstract
In the context of space applications, the GPS system is presently a well-established and accepted tracking system. To meet the basic navigation requirements, most satellites in a low Earth orbit are equipped with single-frequency GPS receivers that measure the coarse acquisition code as well as the L1 phase.

However, the resulting kinematic navigation solutions exhibit systematic position errors caused by elevation-dependent ionospheric path delays. In this study a simple analytical model is established, which quantitatively relates the position error to the vertical electron content and the mapping function. This model substantiates the empirical evidence of a mean radial offset that increases in proportion to the total electron content above the satellite. It is furthermore shown that the ratio between this offset and the vertical ionospheric path delay depends on the applied elevation mask angle. Representative ratios of 3-5 are obtained for the mapping function of the Lear ionosphere model and elevation cutoff angles of 10°, 5°, and 0°. This analytical result has further been confirmed by signal simulator tests as well as flight data of the CHAMP satellite.

Agosto de 2006
Global gravity, bathymetry, and the distribution of submarine volcanism through space and time
Authors: A. B. Watts , D. T. Sandwell et al
Link: Click here

Abstract
The seafloor is characterized by numerous seamounts and oceanic islands which are mainly volcanic in origin. Relatively few of these features (<~0.1%), however, have been dated, and so little is known about their tectonic setting. One parameter that is sensitive to whether a seamount formed on, near, or far from a mid-ocean ridge is the elastic thickness, T e , which is a proxy for the long-term strength of the lithosphere. Most previous studies are based on using the bathymetry to calculate the gravity anomaly for different values of T e and then comparing the calculated and observed gravity anomaly. The problem with such an approach is that bathymetry data are usually limited to single-beam echo sounder data acquired along a ship track and these data are too sparse to define seamount shape. We therefore use the satellite-derived gravity anomaly to predict the bathymetry for different values of T e . By comparing the predicted bathymetry to actual shipboard soundings in the vicinity of each locality in the Wessel global seamount database, we have obtained 9758 T e estimates from a wide range of submarine

volcanic features in the Pacific, Indian, and Atlantic oceans. Comparisons where there are previous estimates show that bathymetric prediction is a robust way to estimate T e and its upper and lower bounds. T e at sites where there is both a sample and crustal age show considerable scatter, however, and there is no simple relationship between T e and age. Nevertheless, we are able to tentatively assign a tectonic setting to each T e estimate. The most striking results are in the Pacific Ocean where a broad swath of "on-ridge" volcanism extends from the Foundation seamounts and Ducie Island/Easter Island ridge in the southeast, across the equator, to the Shatsky and Hess rises in the northwest. Interspersed among the on-ridge volcanism are "flank ridge" and "off-ridge" features. The Indian and Atlantic oceans also show a mix of tectonic settings. Off-ridge volcanism dominates in the eastern North Atlantic and northeast Indian oceans, while flank ridge volcanism dominates the northeastern Indian and western south Atlantic oceans. We have been unable to assign the flank ridge and off-ridge estimates an age, but the on-ridge estimates generally reflect, we believe, the age of the underlying oceanic crust. We estimate the volume of on-ridge volcanism to be ~1.1 × 106 km3 which implies a mean seamount addition rate of ~0.007 km3 yr-1. Rates appear to have varied through geological time, reaching their peak during the Late/Early Cretaceous and then declining to the present-day.

Agosto de 2006
Day-to-day variability in equatorial spread F: Is there some physics missing?
Author: Roland T. Tsunoda
Link: Click here

Abstract
Attempts continue to be made, without notable success, to identify the source of day-to-day variability in occurrence of equatorial spread F (ESF). Most seek to uncover a controlling factor in one of the parameters that describe the linear

growth rate of the Rayleigh-Taylor (RT) instability. It is possible, however, that some physics is still missing in that description of ESF generation. Consideration is given here to an F-region response to a large-scale polarization electric field that is generated by a sporadic-E layer instability (Cosgrove and Tsunoda, 2002) and mapped to the bottomside of the F 2 layer. Results indicate that the large-scale wave structure, the most reliable precursor for ESF, may be initiated by this process at the base of the F 2 layer, where plasma drift is westward.

Agosto de 2006
Tsunami travel time prediction using neural networks
Authors: Rahul Barman , B. Prasad Kumar et al
Link: Click here

Abstract
The present work reports the development of a nonlinear technique based on artificial neural network (ANN) for prediction of tsunami travel time in the Indian Ocean. The expected times of arrival

(ETA) computation involved 250 representative coastal stations encompassing 35 countries. A travel time model is developed using ANN approach. The ANN model uses non-linear regression where a Multi-layer Perceptron (MLP) is used to tackle the non-linearity in the computed ETA. The back-propagation feed forward type network is used for training the system using the resilient back-propagation algorithm. The model demonstrates a high degree of correlation, proving its robustness in development of a real-time tsunami warning system for Indian Ocean.
Agosto de 2006
Magnitude estimation using the first three seconds P-wave amplitude in earthquake early warning
Authors: Yih-Min Wu and Li Zhao
Link: Click here

Abstract
Pd is the peak amplitude of displacement in the first three seconds after the arrival of the P wave. We investigated the attenuation of Pd with the hypocentral distance R in southern California as a

function of magnitude M, and obtained the following relationship: log (P d ) = -3.463 + 0.729 × M - 1.374 × log (R) ± 0.305. Given an earthquake location determined by the P-wave arrival times at stations close to the epicenter, this relationship can be used to define a so-called "Pd magnitude" of earthquakes. Our result shows that for earthquakes in southern California the Pd magnitudes agree with the catalog magnitudes with a standard deviation of 0.18 for events less than magnitude 6.5. Therefore, Pd is a robust measurement for estimating the magnitudes of earthquakes and has practical application in earthquake early warning systems.
Agosto de 2006
Distinguishing ionospheric models using Schumann resonance spectra
Authors: Earle R. Williams , Vadim C. Mushtak et al
Link: Click here

Abstract
A number of models for Schumann resonance (SR) behavior in uniform approximations of the real Earth-ionosphere cavity now populate the literature. These models are treated in terms of variously formulated propagation parameters: as the complex eigenvalue of the propagation problem, as the complex incident angle's sine, as the phase velocity and attenuation rate, or as a

pair of complex characteristic altitudes. This study compares a priori theoretical propagation parameters with the corresponding quantities recovered from calculated Schumann resonance spectra by means of Lorentzian spectral fits. To estimate the ultimate accuracy of the recovery procedure, the influence of source-receiver separation is excluded by assuming a globally uniform distribution of lightning sources. The comparisons show a practically acceptable agreement, within several percent, agreement between recovered and a priori parameters for all models studied. When judged against real Schumann resonance observations, these results shed light on problems with certain models. More importantly, the results reaffirm the ability of procedures based on SR observations to resolve global features of the ionosphere's state and structure.
Agosto de 2006
The apparently isotropic Australian upper mantle
Authors: Maggy Heintz , Brian L. N. Kennett et al
Link: Click here

Abstract
We investigate shear wave splitting measurements performed on two years of data recorded at stations deployed in the TASMAL experiment, a network of 20 broadband seismological stations designed to record data on each side of the controversial Tasman Line in Australia. Whereas a subset of measurements previously performed on one year of data exhibited a curvilinear pattern similar to that of the Tasman

Line, suggesting anisotropy frozen in the lithosphere, considering the whole data set drastically changes the situation: apparent isotropy in the Australian upper mantle is observed at numerous stations. This apparent isotropy together with the EW or NS orientations of the polarization plane of the fast S wave (phi ) observed at some stations is consistent with a two-layer anisotropic system underneath the Australian continent, with a perpendicular orientation of in each layer. From the latest tomographic results, the transition between the Precambrian western and Phanerozoic eastern Australia appears to define blocks of various thickness. Unlike the situation across the TESZ in Europe, these blocks do not seem to be correlated with a different behavior in terms of seismic anisotropy.
Agosto de 2006
Laboratory measurements of P- and S-wave velocities in polycrystalline plagioclase and gabbronorite up to 700 °C and 1 GPa: Implications for the low velocity anomaly in the lower crust
Authors: Yoshio Kono , Masahiro Ishikawa et al
Link: Click here

Abstract
P- (Vp) and S-wave (Vs) velocities of polycrystalline plagioclase (An51±1) and gabbronorite including plagioclase (An49±2) (62.3vol.%) were measured up to 700°C at 0.6, 0.8, and 1.0GPa during heating and cooling. Both

polycrystalline plagioclase and gabbronorite show reversible and discontinuous change in temperature derivative of Vp and Vs, dVp/dT and dVs/dT, respectively, at about 400°C. The dVp/dT and dVs/dT of the polycrystalline plagioclase changes from -0.8 × 10-4 to -3.4 × 10-4 km s-1°C-1 and from -1.1 × 10-4 to -3.3 × 10-4 km s-1°C-1, respectively. The dVp/dT and dVs/dT of the gabbronorite varies from -1.0 × 10-4 to -3.4 × 10-4 km s-1°C-1 and from -0.4 × 10-4 to -3.5 × 10-4 km s-1°C-1, respectively. The reversible and discontinuous changes in dVp/dT and dVs/dT are attributed to a phase transformation of plagioclase. The present data suggest that the sharp decrease in Vp and Vs would cause low velocity anomaly under dry and subsolidus melt-absent conditions in the mid-to-lower crust of relatively high heat flow regions.
Agosto de 2006
Long-range triggered earthquakes that continue after the wave train passes
Author: Emily E. Brodsky
Link: Click here

Abstract
Large earthquakes can trigger distant earthquakes in geothermal areas. Some triggered earthquakes happen while the surface waves pass through a site, but others occur hours or even days later. Does this prolonged seismicity require a special mechanism to store the stress from the seismic waves that differs from ordinary aftershock mechanisms? These questions have driven studies of long-range triggering since the phenomenon's discovery. Here I attempt to answer the questions by examining the statistics

of triggered sequences. Two separate observations are consistent with the prolonged sequences being simply local aftershocks of earthquakes triggered early in the wave train. First, the sequences obey Omori's Law over both short (1 hour) and longer (5 day) time intervals. Secondly, the number of observed triggered earthquakes in the first hour after the wave train can be predicted from the number of earthquakes triggered during the wave train. Even the very vigorous 10-day triggering at Long Valley from the 1992 Landers M w 7.3 earthquakes can be interpreted as the aftershocks of either a local M ˜ 4.1 earthquake or an equivalent combination of several smaller mainshocks. Therefore, long-range triggering does not need to include a mechanism to produce sustained stresses other than the process that generates aftershocks of the earthquakes that occur while the wave train is passing.
Agosto de 2006
Improved modeling of Lunar Prospector neutron spectrometer data: Implications for hydrogen deposits at the lunar poles
Authors: David J. Lawrence , W. C. Feldman et al
Link: Click here

Abstract
New models have been computed for the Lunar Prospector (LP) thermal and epithermal neutron counting rates using the particle transport code MCNPX. This work improves upon previous studies by using one code to model the neutron production, transport, and detection processes, and by examining the sensitivity of epithermal neutrons to elements other than hydrogen. Our modeling results for standard anhydrous lunar

soils show that when hydrogen is not included in a soil, epithermal neutrons are most sensitive to variations in the abundances of Fe, Gd, and Sm, which is consistent with measured epithermal neutron data. We use our current modeling results, in conjunction with known mineral compositions of lunar soils and other lunar global data sets to conclude that the best explanation for a decrease in the counting rate of epithermal neutrons near both lunar poles is the presence of hydrogen. We have further concluded that the average hydrogen abundance near both lunar poles is 100-150 ppm and is likely buried by 10 ± 5 cm of dry lunar soil, a result that is consistent with previous studies. The localized hydrogen abundance for small (<20 km) areas of permanently shaded regions remains highly uncertain and could range from 200 ppm H up to 40 wt% H2O in some isolated regions.
Agosto de 2006
Modeling the isotopic composition of Antarctic snow using backward trajectories: Simulation of snow pit records
Authors: M. M. Helsen , R. S. W. van de Wal et al
Link: Click here

Abstract
The quantitative interpretation of isotope records (d18O, dD, and d excess) in ice cores can benefit from a comparison of observed meteorology with associated isotope variability. For this reason we studied four isotope records from snow pits in western Dronning Maud Land (DML), Antarctica, covering the period 1998-2001. Timing and magnitude of snowfall events on these locations were monitored using sonic height rangers. For the distinguished snowfall events we evaluated the isotopic composition of the moisture during transport by combining backward trajectory calculations with isotopic modeling, using a Rayleigh-type distillation model (MCIM). The initial isotope ratio of the moisture was determined from

monthly mean isotope fields from a general circulation model (ECHAM4). The trajectory analysis showed that the southern Atlantic Ocean is the major moisture source for precipitation in DML. Modeling results along the trajectories revealed that most of the isotopic depletion occurred during the last day of the transport. Finally, a diffusion model was applied to describe the diffusion in the firn layer such that the modeled isotopes could be compared with the observed isotope records. The resulting modeled isotope profiles were mostly in good agreement with the observed seasonal variability in the snow. However, at low temperatures (especially on the Antarctic interior), magnitude of the total distillation was underestimated. Regarding the d excess parameter, our results show a large influence of advection height on the final value of d excess in precipitation. This in turn points to the importance of the vertical structure of d excess over the oceanic source region, which obscures the classical interpretation of this parameter in terms of temperature and relative humidity in the moisture source region.

Agosto de 2006
Monitoring of an infiltration experiment using the self-potential method
Authors: B. Suski , A. Revil et al
Link: Click here

Abstract
An infiltration test was performed from a ditch with the purpose of monitoring the evolution of the piezometric levels using self-potential measurements made at the ground surface. We used a set of 18 piezometers and a network of 41 nonpolarizable (Pb/PbCl2) electrodes. The variations of the self-potential signals are linearly correlated to the piezometric level changes with an apparent voltage coupling coefficient of -5.5

± 0.9 mV m-1. We measured, independently of this infiltration test, the three material properties entering the macroscopic field equations. They are the resistivity distribution of the soil, its mean hydraulic conductivity, and its intrinsic streaming potential coupling coefficient (-5.8 ± 1.1 mV m-1). Then, we modeled numerically the infiltration test and the associated self-potential signals using a two-dimensional finite difference code. The numerical model reproduces fairly well the observed results. This investigation demonstrates the effectiveness of the self-potential method in field conditions to monitor small variations (<0.60 m) of the water table. It offers for the first time a test of the electrokinetic theory in the field with independent evaluation of the material properties entering the field equations.
Agosto de 2006
Initiation propagation and termination of elastodynamic ruptures associated with segmentation of faults and shaking hazard
Author: Bruce E. Shaw
Link: Click here

Abstract
Using a model of a complex fault system, we examine the initiation, propagation, and termination of ruptures and their relationship to fault geometry and shaking hazard. We find concentrations of epicenters near fault step overs and ends; concentrations of terminations near

fault ends; and persistent propagation directivity effects. Taking advantage of long sequences of dynamic events, we directly measure shaking hazards, such as peak ground acceleration exceedance probabilities, without need for additional assumptions. This provides a new tool for exploring shaking hazard from a physics-based perspective, its dependence on various physical parameters, and its correlation with other geological and seismological observables. Using this capability, we find some significant aspects of the shaking hazard can be anticipated by measures of the epicenters. In particular, asymmetries in the relative peak ground motion hazard along the faults appear well correlated with asymmetries in epicentral locations.
Agosto de 2006
High-precision isotopic characterization of USGS reference materials by TIMS and MC-ICP-MS
Authors: Dominique Weis , Bruno Kieffer et al
Link: Click here

Abstract
The Pacific Centre for Isotopic and Geochemical Research (PCIGR) at the University of British Columbia has undertaken a systematic analysis of the isotopic (Sr, Nd, and Pb) compositions and concentrations of a broad compositional range of U.S. Geological Survey (USGS) reference materials, including basalt (BCR-1, 2; BHVO-1, 2), andesite (AGV-1, 2), rhyolite (RGM-1, 2), syenite (STM-1, 2), granodiorite (GSP-2), and granite (G-2, 3). USGS rock reference materials are geochemically well characterized, but there is neither a systematic methodology nor a database for radiogenic isotopic compositions, even for the widely used BCR-1. This investigation represents the first comprehensive, systematic analysis of the isotopic composition and concentration of USGS reference materials and provides an important database for the isotopic community. In addition, the range of equipment at the PCIGR, including a Nu Instruments Plasma MC-ICP-MS, a

Thermo Finnigan Triton TIMS, and a Thermo Finnigan Element2 HR-ICP-MS, permits an assessment and comparison of the precision and accuracy of isotopic analyses determined by both the TIMS and MC-ICP-MS methods (e.g., Nd isotopic compositions). For each of the reference materials, 5 to 10 complete replicate analyses provide coherent isotopic results, all with external precision below 30 ppm (2 SD) for Sr and Nd isotopic compositions (27 and 24 ppm for TIMS and MC-ICP-MS, respectively). Our results also show that the first- and second-generation USGS reference materials have homogeneous Sr and Nd isotopic compositions. Nd isotopic compositions by MC-ICP-MS and TIMS agree to within 15 ppm for all reference materials. Interlaboratory MC-ICP-MS comparisons show excellent agreement for Pb isotopic compositions; however, the reproducibility is not as good as for Sr and Nd. A careful, sequential leaching experiment of three first- and second-generation reference materials (BCR, BHVO, AGV) indicates that the heterogeneity in Pb isotopic compositions, and concentrations, could be directly related to contamination by the steel (mortar/pestle) used to process the materials. Contamination also accounts for the high concentrations of certain other trace elements (e.g., Li, Mo, Cd, Sn, Sb, W) in various USGS reference materials.
Agosto de 2006
Magnetohydrodynamic simulation of the interaction between interplanetary strong shock and magnetic cloud and its consequent geoeffectiveness
Authors: Ming Xiong , Huinan Zheng et al
Link: Click here

Abstract
Numerical studies have been performed to interpret the observed "shock overtaking magnetic cloud (MC)" event by a 2.5 dimensional magnetohydrodynamic (MHD) model in the heliospheric meridional plane. Results of an individual MC simulation show that the MC travels with a constant bulk flow speed. The MC is injected with a very strong inherent magnetic field over that in the ambient flow and expands rapidly in size initially. Consequently, the diameter of the MC increases in an asymptotic speed while its angular width contracts gradually. Meanwhile,

simulations of MC-shock interaction are also presented, in which both a typical MC and a strong fast shock emerge from the inner boundary and propagate along the heliospheric equator, separated by an appropriate interval. The results show that the shock first catches up with the preceding MC, then penetrates through the MC, and finally merges with the MC-driven shock into a stronger compound shock. The morphologies of shock front in interplanetary space and MC body behave as a central concave and a smooth arc, respectively. The compression and rotation of the magnetic field serve as an efficient mechanism to cause a large geomagnetic storm. The MC is highly compressed by the overtaking shock. Contrarily, the transport time of the incidental shock influenced by the MC depends on the interval between their commencements. Maximum geoeffectiveness results from when the shock enters the core of preceding MC, which is also substantiated to some extent by a corresponding simplified analytic model. Quantified by the Dst index, the specific result is that the geoeffectiveness of an individual MC is largely enhanced with 80% increment in maximum by an incidental shock.
Agosto de 2006
How probability weighting affects participation in water markets
Authors: Ram Ranjan and Jason F. Shogren.
Link: Click here

Abstract
The behavioral tendency to overestimate probabilities of loss can affect a farmer's participation in water markets. We examine this issue with a theoretical model of a nonexpected utility maximizing farmer who places subjective weights on the actual probabilities of loss of water rights due to market transactions. The farmer bargains over sharing of surpluses with the buyer

of water. The farmer then incorporates the bargaining outcome in his intertemporal expected benefit maximization problem that accounts for the possible loss of water rights due to its sale out of agriculture. Three key results emerge. First, subjective weighting of probabilities leads to discounting of resources when farmers overestimate probabilities of loss. Second, if farmers have idiosyncratic time preferences, total water supply in the market would depend on the level of heterogeneity in the population. Third, considering the case of two farmers, we find that the farmer with lower endowments bears the burden of risk reduction, whereas the one with higher endowments sells more water for profits. As the level of risk increases, however, the relative difference in risk sharing declines.
Agosto de 2006
Mitigating atmospheric noise for InSAR using a high resolution weather model
Authors: J. Foster, B. Brooks et al
Link: Click here

Abstract
A high resolution weather model is used to predict atmospheric delays for the acquisition times of synthetic aperture radar images over Hawaii.

Refraction of the radar by water vapor in the atmosphere in Hawaii leads to apparent ground-motions with wavelengths and magnitudes similar to the actual ground motions generated by tectonic and volcanic processes. We examine the potential for a weather model to help characterize the atmospheric component in InSAR scenes and find that in the best cases it models the observed delays well, reducing the variance at wavelengths of 30 km and greater by ~60%, while even in the worst cases it provides an independent means of quantifying the expected variance in the image due to the atmosphere.
Agosto de 2006
The Mw 7.8, 2001 Kunlunshan earthquake: Extreme rupture speed variability and effect of fault geometry
Authors: D. P. Robinson , C. Brough et al
Link: Click here

Abstract
By analyzing body wave seismograms, we show that the rupture speed on the Main Kunlun Fault during the M w 7.8 2001 Kunlunshan, Tibet, earthquake was highly variable and the rupture process consisted of three stages. In the first stage, the rupture accelerated from rest to an average speed of 3.3 km/s over a distance of 120 km. The rupture then propagated for another 150 km at an apparent rupture speed exceeding the P

wave speed. In the final stage, the earthquake fault bifurcates, and the rupture front slowed down. The region of highest rupture velocity is found to coincide with the region of highest fault slip, has the longest slip duration, and is where off-fault ground cracking is observed in the field. Stress drops are found to be higher in regions of higher rupture speeds. The greatest concentration of aftershocks is located near the fault bifurcation zone and hence coincides with the region of highest fault slip, highest stress drop and highest rupture velocity. The fault width is no more than 10 km in most places and is about 20 km in the region of highest slip. This narrow fault width is attributed to the fact that crust below this depth is sufficiently warm not to permit brittle failure to occur. The remarkable similarity of this earthquake with the 1906 California earthquake, in spite of occurring in very different tectonic regimes, is discussed.
Agosto de 2006
Evidence in support of the climate change-Atlantic hurricane hypothesis
Author: James B. Elsner
Link: Click here

Abstract
The power of Atlantic tropical cyclones is rising rather dramatically and the increase is correlated with an increase in the late summer/early fall sea surface temperature over the North Atlantic. A debate concerns the nature of these increases with some studies attributing them to a natural

climate fluctuation, known as the Atlantic Multidecadal Oscillation (AMO), and others suggesting climate change related to anthropogenic increases in radiative forcing from greenhouse-gases. Here tests for causality using the global mean near-surface air temperature (GT) and Atlantic sea surface temperature (SST) records during the Atlantic hurricane season are applied. Results show that GT is useful in predicting Atlantic SST, but not the other way around. Thus GT "causes" SST providing additional evidence in support of the climate change hypothesis. Results have serious implications for life and property throughout the Caribbean, Mexico, and portions of the United States.
Agosto de 2006
Three years of Atmospheric Infrared Sounder radiometric calibration validation using sea surface temperatures
Authors: H. H. Aumann , Steve Broberg et al
Link: Click here

Abstract
This paper evaluates the absolute accuracy and stability of the radiometric calibration of the Atmospheric Infrared Sounder (AIRS) by analyzing the difference between the brightness temperatures measured at 2616 cm-1 and those calculated at the top of the atmosphere (TOA), using the Real-Time Global Sea Surface Temperature (RTGSST) for cloud-free night tropical oceans between ±30° latitude. The TOA correction is based on radiative transfer. The analysis of the first 3 years of AIRS radiances verifies the absolute calibration at 2616 cm-1 to better than 200 mK, with better than 16 mK/yr stability. The AIRS radiometric calibration uses an internal full aperture wedge blackbody with the

National Institute of Standards and Technology (NIST) traceable prelaunch calibration coefficients. The calibration coefficients have been unchanged since launch. The analysis uses very tight cloud filtering, which selects about 7000 cloud-free tropical ocean spectra per day, about 0.5% of the data. The absolute accuracy and stability of the radiometry demonstrated at 2616 cm-1 are direct consequences of the implementation of AIRS as a thermally controlled, cooled grating-array spectrometer and meticulous attention to details. Comparable radiometric performance is inferred from the AIRS design for all 2378 channels. AIRS performance sets the benchmark for what can be achieved with a state-of-the-art hyperspectral radiometer from polar orbit and what is expected from future hyperspectral sounders. AIRS was launched into a 705 km altitude polar orbit on NASA's Earth Observation System (EOS) Aqua spacecraft on 4 May 2002. AIRS covers the 3.7-15.4 micron region of the thermal infrared spectrum with a spectral resolution of V/delta V= 1200 and has returned 3.7 million spectra of the upwelling radiance each day since the start of routine data gathering in September 2002.
Agosto de 2006
Relationship between geochemical parameters and the occurrence of Dehalococcoides DNA in contaminated aquifers
Authors: Xiaoxia Lu , John T. Wilson et al
Link: Click here

Abstract
Strains of Dehalococcoides are the only microbes known that can completely dechlorinate PCE, TCE, cis-DCE, and vinyl chloride to ethylene. Either naturally occurring strains or bioaugmentation cultures of Dehalococcoides are widely used for in situ bioremediation of contaminated groundwater. Naturally occurring strains have an important role in natural attenuation of PCE, TCE, cis-DCE, and vinyl chloride in groundwater. This study evaluated the relationship between selected biogeochemical parameters and the presence of Dehalococcoides DNA in field-scale plumes. A total of 81 monitoring wells were sampled from 15 groundwater plumes at 10 locations across the United States (one sample per monitoring well). The presence of Dehalococcoides DNA was

determined with an assay based on the polymerase chain reaction (PCR) using DNA primers targeting the 16S rRNA gene of Dehalococcoides. The groundwater samples were also analyzed for concentrations of O2, NO3 -1 plus NO2 -1 - N, CH4, H2, Fe (II), SO4 -2, TOC, Cl-1, and benzene, toluene, ethylbenzene, and xylene (BTEX) compounds and for alkalinity, ORP, electrical conductivity, pH, and temperature. Dehalococcoides DNA was unequivocally detected in 26 wells, most of which exhibited methanogenic conditions. A two-sample Kolmogorov-Smirnov test was used to compare the distribution of each parameter in water where Dehalococcoides DNA was present to the distribution where Dehalococcoides DNA was absent. The only parameters for which the distributions were different at 95% confidence were NO3 -1 plus NO2 -1 - N, CH4, and ORP. Using these three statistically significant geochemical parameters as descriptors, a predictive model for the presence of Dehalococcoides DNA was developed using logistic regression with a binary response. Under conditions where data of direct biochemical assay are not available a calculated probability could be used to properly calibrate computer models of natural attenuation.
Agosto de 2006
Using uplifted Holocene beach berms for paleoseismic analysis on the Santa María Island, south-central Chile
Authors: B. Bookhagen , H. P. Echtler et al
Link: Click here

Abstract
Major earthquakes (M > 8) have repeatedly ruptured the Nazca-South America plate interface of south-central Chile involving meter scale land-level changes. Earthquake recurrence intervals, however, extending beyond limited historical records are virtually unknown, but would provide

crucial data on the tectonic behavior of forearcs. We analyzed the spatiotemporal pattern of Holocene earthquakes on Santa María Island (SMI; 37°S), located 20 km off the Chilean coast and approximately 70 km east of the trench. SMI hosts a minimum of 21 uplifted beach berms, of which a subset were dated to calculate a mean uplift rate of 2.3 ± 0.2 m/ky and a tilting rate of 0.022 ± 0.002 °/ky. The inferred recurrence interval of strandline-forming earthquakes is ~180 years. Combining coseismic uplift and aseismic subsidence during an earthquake cycle, the net gain in strandline elevation in this environment is ~ 0.4 m per event.

Agosto de 2006
Modeling the geomagnetic effects caused by the solar eclipse of 11 August 1999
Authors: J. J. Curto , B. Heilig at al
Link: Click here

Abstract
The solar eclipse of 11 August 1999 was total along a belt crossing Central Europe, where there

is a high density of magnetic observatories. We studied the transient geomagnetic effects produced by the temporary cutoff of the ionizing solar radiations during the eclipse. It is the first case in which the phenomenon has been analysed in a multisite context. For observatories along the umbra path, a model based on the equivalent electric current schema is proposed.
Agosto de 2006
A nonparametric stochastic downscaling framework for daily rainfall at multiple locations
Authors: R. Mehrotra and Ashish Sharma
Link: Click here

Abstract
Use of General Circulation Models (GCMs) for climate change impact assessment is often limited by their incapability at representing local features and dynamics at spatial scales finer than the in-built GCM grid scale. This has led to the development of downscaling techniques for transfer of coarse GCM simulated weather output to finer spatial resolutions. This paper presents a nonparametric stochastic spatial downscaling framework for multisite daily rainfall occurrence and amount. At site rainfall occurrences are downscaled using a nonparametric nonhomogeneous hidden Markov model (NNHMM) that represents spatial dependence across the rainfall occurrence field using a dynamic weather state indicative of the centroid and average wetness fraction of the rainfall

occurrence field. The rainfall amounts on the wet days are downscaled using a nonparametric kernel density approach that accommodates variations in the rainfall downscaling model at individual locations. Spatial dependence in the rainfall amounts is simulated by driving each of the single-site amounts model with spatially correlated random numbers. The proposed framework is applied for downscaling of rainfall at a network of 30 rain gauge stations around Sydney in Australia, and its performance is evaluated. The analyses of the results show that the logic of providing separate treatments for rainfall occurrence and amounts at individual locations imparts considerable accuracy in the representation of characteristics of interest in hydrologic studies. These characteristics include representation of rainfall spell patterns, spatial distribution of the rainfall occurrence and amount fields, representation of low and high rainfall extremes at individual stations and across the field, as well as common indicators of water balance and variability that are of importance in a catchment scale water balance simulation.
Agosto de 2006
Configuration of the Indian Moho beneath the NW Himalaya and Ladakh
Authors: S. S. Rai, K. Priestley et al
Link: Click here

Abstract
Teleseismic receiver function analysis of seismograms recorded on a ~700 km long profile of 17 broadband seismographs traversing the NW Himalaya shows a progressive northward deepening of the Indian Moho from ~40 km

beneath Delhi south of the Himalayan foredeep to ~75 km beneath Taksha at the Karakoram Fault. Similar studies by Wittlinger et al. (2004) to the north of the Karakoram Fault show that the Moho continues to deepen to ~90 km beneath western Tibet before shallowing substantially to 50-60 km at the Altyn Tagh Fault. The continuity of the Indian Moho imaged in the receiver functions reported here, along with those of Wittlinger et al. (2004), suggest that in this part of the Himalayan orogen the Indian plate may penetrate as far as the Bangong Suture, and possibly as far north as the Altyn Tagh.
Agosto de 2006
A new approach to directly determine the secular variation from magnetic satellite observations
Authors: M. Mandea and N. Olsen
Link: Click here

Abstract
Observatory monthly means provide an excellent opportunity to study the temporal changes of the geomagnetic field at a given location. Unfortunately, determination of the global pattern of changes using observatory data is hampered by their uneven distribution. Satellite data provide excellent global coverage, but the spacecraft

movement makes direct comparisons of satellite and observatory data difficult. To investigate short-period secular variation in observatory and satellite data, we developed an approach to extract satellite monthly means for "virtual observatories" at 400 km altitude, using CHAMP magnetic measurements. Comparison of these virtual observatory monthly means with the corresponding ground values shows a remarkably well-correlated signal at time-scales of months to years, which is beyond the temporal resolution limit of recent global models. Here, we describe this newly developed approach, its validation, and discuss how it can be used to understand short-period changes of the recent geomagnetic field.
Agosto de 2006
Rapid estimation of first-order rupture characteristics for large earthquakes using surface waves: 2004 Sumatra-Andaman earthquake
Authors: Charles J. Ammon , Aaron A. Velasco et al
Link: Click here

Abstract
Broadband surface waves from large earthquakes can be rapidly processed to estimate seismic moment, faulting duration, and general slip distribution, supplying important information for hazard response efforts, including tsunami warnings. Deconvolution of surface-wave Green's

functions removes propagation effects, yielding azimuthally varying effective source time functions. Even a single fault-perpendicular time function can provide robust first-order estimates of the temporal history of seismic radiation and associated slip distribution. Simplified two-dimensional finite-fault modeling using a few surface-wave time functions can reliably characterize the smooth component of the overall rupture process. These procedures are demonstrated using 13 Rayleigh wave observations for the 26 December 2004 Sumatra-Andaman earthquake, yielding satisfactory agreement with models based on more complete seismic data sets. Depending on source and station locations, stable slip estimates can be obtained within 15-60 minutes of the onset of a large earthquake.
Agosto de 2006
Seismic velocity and attenuation structures in the top 400 km of the Earth's inner core along equatorial paths
Authors: Wen-che Yu , Lianxing Wen et al
Link: Click here

Abstract
We study seismic velocity and attenuation structures in the top 400 km of the Earth's inner core based on modeling of differential traveltimes, amplitude ratios, and waveforms of the PKiKP-PKIKP phases observed at the epicentral distance range of 120°-141° and the PKPbc-PKIKP phases observed at the distance range of 146°-160° along equatorial paths. Our data are selected from the seismograms recorded in the Global Seismographic Network from 1990 to 2001 and many regional seismic networks. The observed PKiKP-PKIKP and PKPbc-PKIKP phases exhibit distinctive "east-west" hemispheric patterns: (1) At the distance ranges of 131°-141° and 146°-151°, PKIKP phases arrive about 0.3 s earlier than the theoretical arrivals based on the Preliminary Reference Earth Model (PREM) for the PKIKP phases sampling the "eastern hemisphere" (40°E-180°E) of the inner core and about 0.4 s later for those sampling the "western hemisphere" (180°W-40°E). At the distance range of 151°-160°, PKIKP phases arrive about 0.7 s earlier than the predicted arrivals based on PREM for those sampling the eastern hemisphere and about 0.1 s later for those sampling the western hemisphere.

(2) Amplitude ratios of the PKIKP/PKiKP phases at the distance range of 131°-141° and of the PKIKP/PKPbc phases at the distance range of 146°-151° are, in general, smaller for the PKIKP phases sampling the eastern hemisphere than for those sampling the western hemisphere. At distances greater than 151°, the PKIKP/PKPbc amplitude ratios become indistinguishable for the two hemispheres. These observations can be best explained by two different types of seismic velocity and attenuation models along equatorial paths, one for each hemisphere, in the top 400 km of the inner core. For the eastern hemisphere, the velocity structure has a velocity increase of 0.748 km/s across the inner core boundary (ICB), a small velocity gradient of 0.0042 (km/s)/100 km in the top 235 km, followed by a steeper velocity gradient of 0.1 (km/s)/100 km extending from 235 km to 375 km, and a velocity gradient of 0.01 (km/s)/100 km in the deeper portion of the inner core; the attenuation structure has an average Q value of 300 in the top 300 km and an average Q value of 600 in the deeper portion of the inner core. For the western hemisphere, the velocity structure has a velocity increase of 0.645 km/s across the ICB and a velocity gradient of 0.049 (km/s)/100 km in the top 375 km; the attenuation structure has an average Q value of 600 in the top 375 km of the inner core. Our results suggest that the inner core hemispheric variations in velocity extend deeper than 375 km below the ICB and the top 235 km of the inner core in the eastern hemisphere is anomalous compared to the rest of the inner core in having a small velocity gradient, high velocity, and high attenuation.
Agosto de 2006
Characterization of permeability anisotropy using wavelet analysis
Authors:Roseanna M. Neupauer , Kaye L. Powell et al
Link: Click here

Abstract
Wavelet analysis is an image analysis technique that can extract local information at multiple scales. Because of this capability, wavelet analysis can be used to identify dominant scales in statistically heterogeneous and anisotropic random fields. We develop and test a wavelet

analysis method for identifying dominant scales and orientations in permeability fields and for identifying boundaries between regions with different dominant orientations. We evaluate three different wavelets (fully anisotropic Morlet wavelet, Mexican hat wavelet, and Cauchy wavelet) and show that the Morlet wavelet is the most effective of these three wavelets in identifying dominant orientations. We also investigate the use of several different quantitative wavelet measures in identifying dominant scales and orientations in permeability fields. The technique is demonstrated using both a synthetic data set with known characteristics and a laboratory-collected permeability data set from Massillon sandstone
Agosto de 2006
Application of Hopfield neural network for extracting Doppler spectrum from ocean echo
Authors: Renzhuo Gui , Zijie Yang et al
Link: Click here

Abstract
This paper proposes the method of a Hopfield-type neural network (HNN) for extracting Doppler spectrum from ocean echo. First, it introduces the basic principle of HNN for optimized processing.

Second, expanding the principle of utilizing autoregression (AR) to estimate frequency spectrum, we point out how to apply HNN in spectrum estimation. Last, the three methods are utilized to process actual data, that is, the conventional fast Fourier transform method, modern spectrum estimation-AR method, and the spectrum estimation method based on HNN. The results obtained by the three methods prove that the spectrum estimation method based on HNN is feasible for extracting the Doppler spectrum from ocean echo.
Julio de 2006
P and S velocity structure of the upper mantle beneath the Transantarctic Mountains, East Antarctic craton, and Ross Sea from travel time tomography
Authors: Timothy Watson , Andrew Nyblade et al
Link: Click here

Abstract
P and S wave travel times from teleseismic earthquakes recorded by the Transantarctic Mountains Seismic Experiment (TAMSEIS) have been used to tomographically image upper mantle structure beneath portions of the Transantarctic Mountains (TAM), the East Antarctic (EA) craton, and the West An tarctic rift system (WARS) in the vicinity of Ross Island, Antarctica. The TAM form a major tectonic boundary that divides the stable EA craton and the tectonically active WARS. Relative arrival times were determined using a multichannel cross-correlation technique on teleseismic P and S phases from earthquakes with mb = 5.5. 3934 P waves were used from 322 events, and 2244 S waves were used from 168 events. Relative travel time residuals were inverted for upper mantle structure using VanDecar's method. The P wave

tomography model reveals a low-velocity anomaly in the upper mantle of approximately dVp = -1 to -1.5% in the vicinity of Ross Island extending laterally 50 to 100 km beneath the TAM from the coast, placing the contact between regions of fast and slow velocities well inland from the coast beneath the TAM. The magnitude of the low-velocity anomaly in the P wave model appears to diminish beneath the TAM to the north and south of Ross Island. The depth extent of the low-velocity anomaly is not well constrained, but it probably is confined to depths above ~200 km. The S wave model, within resolution limits, is consistent with the P wave model. The low-velocity anomaly within the upper mantle can be attributed to a 200-300 K thermal anomaly, consistent with estimates obtained from seismic attenuation measurements. The presence of a thermal anomaly of this magnitude supports models invoking a thermal buoyancy contribution to flexurally driven TAM uplift, at least in the Ross Island region of the TAM. Because the magnitude of the anomaly to the north and south of Ross Island may diminish, the thermal contribution to the uplift of the TAM could be variable along strike, with the largest contribution in the Ross Island region. The tomography results reveal faster than average velocities beneath East Antarctica, as expected for cratonic upper mantle.
.

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

geophysics job course geophysics training education geophysics geophysics monitoring geology geophysics borehole geophysics geophysics mining exploration geophysics environmental geophysics geophysics introduction engineering geophysics