From stable dipolar towards reversing numerical dynamos

We are using a three-dimensional convection-driven numerical dynamo model without hyperdiffusivity to study the characteristic structure and time variability of the magnetic field in dependence of the Rayleigh number (Ra) for values up to 40 times supercritical. We also compare a variety of ways to drive the convection and basically find two dynamo regimes. At low Ra, the magnetic field at the surface of the model is dominated by the non-reversing axial dipole component. At high Ra, the dipole part becomes small in comparison to higher multipole components. At transitional values of Ra, the dynamo vacillates between the dipole-dominated and the multipolar regime, which includes excursions and reversals of the dipole axis. We discuss, in particular, one model of chemically driven convection, where for a suitable value of Ra, the mean dipole moment and the temporal evolution of the magnetic field resemble the known properties of the Earth’s field from paleomagnetic data.     

Seismic Hazard Estimate at the Iberian Peninsula

—?Seismic hazard at the Iberian Peninsula has been evaluated by using a methodology which combines both zonified and non-zonified probabilistic methods. Seismic sources are used when considering zones where certain calculation parameters may be considered homogeneous, as in zonified methods, while, on the other hand, earthquakes are considered wherever it has taken place, as in non-zonified methods. The methodology which is applied in this paper has been originally used to calculate the seismic hazard maps in the United States. In our case, it has been necessary to adapt the method to the specific features of the seismicity in the Iberian Peninsula and its geographical surroundings, not only with respect to its distribution and characteristics, but also with respect to the properties of the seismic catalog used.¶Geographically, the main feature of the result is the fact that it reflects both historical seismicity and current seismic clusters of the region. Despite the smoothing, maps show marked differences between several seismic zones; these differences becoming more noticeable as exposure time increases. Maximum seismic hazard is found to be in the southwestern region of the Peninsula, especially in the area of the Cape St. Vicent, and around Lisbon. The uncertainty of the results, without considering that due to the attenuation laws, as deduced from the other evaluation parameters, is quite stable, being more sensitive to the parameters b and m _max of the Gutenberg-Richter relation.     

A rifted inside corner massif on the Mid-Atlantic Ridge at 5°S

The structure of the Mid-Atlantic Ridge at 5°S was investigated during a recent cruise with the FS Meteor. A major dextral transform fault (hereafter the 5°S FZ) offsets the ridge left-laterally by 80 km. Just south of the transform and to the west of the median valley, the inside corner (IC – the region bounded by the ridge and the active transform) is marked by a major massif, characterized by a corrugated upper surface. Fossil IC massifs can also be identified further to the west. Unusually, a massif almost as high as the IC massif also characterizes the outside corner (OC) south of the inactive fracture zone and to the east of the median valley. This OC massif has axis-parallel dimensions identical to the IC massif and both are bounded on their sides closest to the spreading axis by abrupt, steep slopes. An axial volcanic ridge is well developed in the median valley both south of the IC/OC massifs and in an abandoned rift valley to the east of the OC massif, but is absent along the new ridge-axis segment between the IC and OC massifs. Wide-angle seismic data show that between the massifs, the crust of the median valley thins markedly towards the FZ. These observations are consistent with the formation of the OC massif by the rifting of an IC core complex and the development of a new spreading centre between the IC and OC massifs. The split IC massif presents an opportunity to study the internal structure of the footwall of a detachment fault, from the corrugated fault surface to deeper beneath the fault, without recourse to drilling. Preliminary dredging recovered gabbros from the scarp slope of the rifted IC massif, and serpentinites and gabbros from the intersection of this scarp with the corrugated surface. This is compatible with a concentration of serpentinites along the detachment surface, even where the massif internally is largely plutonic in nature.     

Identification of Mining Blasts at Mid- to Far-regional Distances Using Low Frequency Seismic Signals

—?This paper reports results from two recent monitoring experiments in Wyoming. Broadband seismic recordings of kiloton class delay-fired cast blasts and instantaneous calibration shots in the Black Thunder coal mine were made at four azimuths at ranges from 1° to 2°. The primary focus of this experiment was to observe and to explain low-frequency signals that can be seen at all azimuths and should routinely propagate above noise to mid-regional distances where most events will be recorded by International Monitoring System (IMS) stations.¶The recordings clearly demonstrate that large millisecond delay-fired cast blasts routinely produce seismic signals that have significant spectral modulations below 10?Hz. These modulations are independent of time, the azimuth from the source and the orientation of the sensor. Low-frequency modulations below 5?Hz are seen beyond 9°. The modulations are not due to resonance as they are not produced by the calibration shots. Linear elastic modeling of the blasts that is guided by mine-blast reports fails to reproduce the fine detail of these modulations but clearly indicates that the enhanced “spectral roughness” is due to long interrow delays and source finiteness. The mismatch between the data and the synthetics is likely due to source processes, such as nonlinear interactions between shots, that are poorly understood and to other effects, such as variations of shot time and yield from planned values, that are known to be omnipresent but cannot be described accurately. A variant of the Automated Time-Frequency Discriminant (Hedlin, 1998b), which uses low-frequency spectral modulations, effectively separates these events from the calibration shots.¶The experiment also provided evidence that kiloton class cast blasts consistently yield energetic 2–10 second surface waves. The surface waves are strongly dependent on azimuth but are seen beyond 9°. Physical modeling of these events indicates that the surface waves are due mainly to the extended source duration and to a lesser extent to the slap-down of spalled material. The directionality is largely a path effect. A discriminant that is based on the partitioning of energy between surface and body waves routinely separates these events from the calibration shots.¶The Powder River Basin has essentially no natural seismic activity. How these mining events compare to earthquake observations remains to be determined.     

Efficacy of Annual Bacteria Monitoring and Shock Chlorination in Wells Finished in a Floodplain Aquifer

Public health authorities generally recommend annual water-quality monitoring of rural water wells and shock chlorination if coliforms are detected. It is implicitly assumed that shock chlorination is effective in ridding most wells of bacteriological pathogens for months to years. Neither annual monitoring nor shock chlorination was effective in addressing coliform contamination of selected water wells in a small town developed on an alluvial aquifer where septic system effluents are impacting well water quality. Considerable temporal variation in total and fecal coliforms was observed in water wells monitored for a six-month period. Individual wells intermittently met and exceeded the drinking water criteria, indicating annual sampling was insufficient. Shock chlorination of three contaminated wells and their associated distribution systems proved ineffective because colonies apparently originated from outside the wells and reappeared over relatively short time periods (ranging from less than one week up to 21 weeks). The relatively fast and similar rate of recovery of total heterotrophic bacteria suggested they are related to biofilm formation in the wells and not to ground water contamination.     

"THREE DIMENSIONAL MATHEMATICAL MODEL OF SEDIMENT TRANSPORT IN ESTUARINE REGIONS-- A CASE STUDY OF THE HAIHE RIVER MOUTH" Yuchuan Bai,Huanting SHEN and Shixiong HU

The paper presents the 3D finite element simulation of tidal flow and Sediment transport in the estuarine region of the Haihe river. The proposed model adopts sigma-transformation of the hydrodynamic and sediment transport equations. The hydrodynamic and sediment transport models are verified in case of a simple test problem for which analytical solutions are available. Finally the models are applied to muddy Haihe river estuary of North China and it is claimed that hydrodynamic and sediment transport models give a reliable comparison with the observed field data. However, there are certain discrepancies, and some reasonable questions regarding the present state-of-art, in the modeling of three-dimensional multilevel hydrodynamics and sediment transport, which are provided below for answer.     

Lithium and lithium isotope profiles through the upper oceanic crust: a study of seawater-basalt exchange at ODP Sites 504B and 896A

Ocean Drilling Program (ODP) Hole 504B near the Costa Rica Rift is the deepest hole drilled in the ocean crust, penetrating a volcanic section, a transition zone and a sheeted dike complex. The distribution of Li and its isotopes through this 1.8-km section of oceanic crust reflects the varying conditions of seawater alteration with depth. The upper volcanic rocks, altered at low temperatures, are enriched in Li (5.6-27.3 ppm) and have heavier isotopic compositions (δ⁷Li=6.6-20.8‰) relative to fresh mid-ocean ridge basalt (MORB) due to uptake of seawater Li into alteration clays. The Li content and isotopic compositions of the deeper volcanic rocks are similar to MORB, reflecting restricted seawater circulation in this section. The transition zone is a region of mixing of seawater with upwelling hydrothermal fluids and sulfide mineralization. Li enrichment in this zone is accompanied by relatively light isotopic compositions (−0.8-2.1‰) which signify influence of basalt-derived Li during mineralization and alteration. Li decreases with depth to 0.6 ppm in the sheeted dike complex as a result of increasing hydrothermal extraction in the high-temperature reaction zone. Rocks in the dike complex have variable isotopic values that range from −1.7 to 7.9‰, depending on the extent of hydrothermal recrystallization and off-axis low-temperature alteration. Hydrothermally altered rocks are isotopically light because ⁶Li is preferentially retained in greenschist and amphibolite facies minerals. The δ⁷Li values of the highly altered rocks of the dike complex are complementary to those of high-temperature mid-ocean ridge vent fluids and compatible to equilibrium control by the alteration mineral assemblage. The inventory of Li in basement rocks permits a reevaluation of the role of oceanic crust in the budget of Li in the ocean. On balance, the upper 1.8 km of oceanic crusts remains a sink for oceanic Li. The observations at 504B and an estimated flux from the underlying 0.5 km of gabbro suggest that the global hydrothermal flux is at most 8×10⁹ mol/yr, compatible with geophysical thermal models. This work defines the distribution of Li and its isotopes in the upper ocean crust and provides a basis to interpret the contribution of subducted lithosphere to arc magmas and cycling of crustal material in the deep mantle.     

Gravity and magnetic signatures of volcanic plugs related to Deccan volcanism in Saurashtra, India and their physical and geochemical properties

The Bouguer anomaly and the total intensity magnetic maps of Saurashtra have delineated six circular gravity highs and magnetic anomalies of 40-60 mGal (10⁻⁵m/s²) and 800-1000 nT, respectively. Three of them in western Saurashtra coincide with known volcanic plugs associated with Deccan Volcanic Province (DVP), while the other three in SE Saurashtra coincide with rather concealed plugs exposed partially. The DVP represents different phases of eruption during 65.5±2.5 Ma from the Reunion plume. The geochemical data of the exposed rock samples from these plugs exhibit a wide variation in source composition, which varies from ultramafic/mafic to felsic composition of volcanic plugs in western Saurashtra and an alkaline composition for those in SE Saurashtra. Detailed studies of granophyres and alkaline rocks from these volcanic plugs reveal a calc-alkaline differentiation trend and a continental tectonic setting of emplacement. The alkaline plugs of SE Saurashtra are associated with NE-SW oriented structural trends, related to the Gulf of Cambay and the Cambay rift basin along the track of the Reunion plume. This indicates a deeper source for these plugs compared to those in the western part and may represent the primary source magma. The Junagadh plug with well differentiated ring complexes in western Saurashtra shows well defined centers of magnetic anomaly while the magnetic anomalies due to other plugs are diffused though of the same amplitude. This implies that other plugs are also associated with mafic/ultramafic components, which may not be differentiated and may be present at subsurface levels. Paleomagnetic measurements on surface rock samples from DVP in Saurashtra suggest a susceptibility of 5.5×10⁻² SI units with an average Koenigsberger ratio (Q_n) of almost one and average direction of remanent magnetization of D=147.4° and I=+56.1°. The virtual geomagnetic pole (VGP) position computed from the mean direction of magnetization for the volcanic plugs and Deccan basalt of Saurashtra is 30°N and 74°W, which is close to the VGP position corresponding to the early phases of Deccan eruption. Modeling of gravity and magnetic anomalies along two representative profiles across Junagadh and Barda volcanic plugs suggest a bulk density of 2900 and 2880 kg/m³, respectively and susceptibility of 3.14×10⁻² SI units with a Q_n ratio of 0.56 which are within the range of their values obtained from laboratory measurements on exposed rock samples. The same order of gravity and magnetic anomalies observed over the volcanic plugs of Saurashtra indicates almost similar bulk physical properties for them. The inferred directions of magnetization from magnetic anomalies, however, are D=337° and 340° and I=−38° and −50° which represent the bulk direction of magnetization and also indicate a reversal of the magnetic field during the eruption of these plugs. Some of these plugs are associated with seismic activities of magnitude ≤4 at their contacts. Based on this analysis, other circular/semi-circular gravity highs of NW India can be qualitatively attributed to similar subsurface volcanic plugs.     

Induced Microearthquake Patterns in Hydrocarbon and Geothermal Reservoirs: Six Case Studies

—?The injection or production of fluids can induce microseismic events in hydrocarbon and geothermal reservoirs. By deploying sensors downhole, data sets have been collected that consist of a few hundred to well over 10,000 induced events. We find that most induced events cluster into well-defined geometrical patterns. In many cases, we must apply high-precision, relative location techniques to observe these patterns. At three sedimentary sites, thin horizontal strands of activity are commonly found within the location patterns. We believe this reflects fracture containment between stratigraphic layers of differing mechanical properties or states of stress. At a massive carbonate and two crystalline sites, combinations of linear and planar features indicate networks of intersecting fractures and allow us to infer positions of aseismic fractures through their influence on the location patterns. In addition, the fine-scale seismicity patterns often evolve systematically with time. At sedimentary sites, migration of seismicity toward the injection point has been observed and may result from slip-induced stress along fractures that initially have little resolved shear. In such cases, triggering events may be critical to generate high levels of seismic activity. At one crystalline site, the early occurrence of linear features that traverse planes of activity indicate permeable zones and possible flow paths within fractures. We hope the continued development of microseismic techniques and refinement of conceptual models will further increase our understanding of fluid behavior and lead to improved resource management in fractured reservoirs.     

Establishment of a Non-Permanent GPS Network to Monitor the Recent NE-SW Deformation in the Granada Basin (Betic Cordillera,Southern Spain)

The Granada Basin (Central Betic Cordillera), one of the most seismically active areas of the Iberian Peninsula, is currently subjected to NW-SE compression and NE-SW extension. The present day extension is accommodated by normal faults with various orientations but particularly with a NW-SE strike. At the surface, these active NW-SE normal faults are mainly concentrated on the NE part of the Basin. In this part we have selected a 15-km long segment where several active normal faults crop out. Using the marine Tortonian rocks as a reference, we have calculated a minimum extensional rate of 0.15-0.30 mm/year. The observed block rotation, the listric geometry of faults at depth and the distribution of seismicity over the whole Basin, indicate that this rate is a minimum value. In the framework of an interdisciplinary research project a non-permanent GPS-network has been established in the central sector of Betic Cordillera to monitor the crustal deformations. The first two observation campaigns were done in 1999 and 2000.