Integrated Two-Dimensional Modeling of Fluid Flow and Compaction in a Sedimentary Basin

In this paper we employ mixed finite elements and numerically study an integrated two-dimensional model of fluid flow and compaction in a sedimentary basin. This model describes a single phase incompressible flow in a two-dimensional section of a sedimentary basin with vertical compaction. At each time step, an iterative algorithm is used to solve this model. The determination of the grid movement is based on the mass conservation and movement of sediments in the basin, while the mixed method is utilized to solve the fluid flow over the moving grid. Numerical experiments are presented to verify this iterative algorithm and show representative solutions for the model under consideration.     

Wave velocities and attenuation of shaley sandstones as a function of pore pressure and partial saturation

We obtain the wave velocities and quality factors of clay‐bearing sandstones as a function of pore pressure, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the coexistence of two solids (sand grains and clay particles) and a fluid mixture. Additional attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/fluid mixing law that satisfies the Wood and Voigt averages at low and high frequencies, respectively. Pressure effects are accounted for by using an effective stress law. By fitting a permeability model of the Kozeny– Carman type to core data, the model is able to predict wave velocity and attenuation from seismic to ultrasonic frequencies, including the effects of partial saturation. Testing of the model with laboratory data shows good agreement between predictions and measurements.     

Bringing NEOs into focus

Duncan Steel and Mark Bailey report on the well-attended December 2001 G Discussion Meeting whose subject was near-Earth objects.     

The El Salvador earthquakes of January and February 2001: context, characteristics and implications for seismic risk

The small Central American republic of El Salvador has experienced, on average, one destructive earthquake per decade during the last hundred years. The latest events occurred on 13 January and 13 February 2001, with magnitudes M_w 7.7 and 6.6, respectively. The two events, which were of different tectonic origin, follow the patterns of the seismicity of the region although neither event has a known precedent in the earthquake catalogue in terms of size and location. The earthquakes caused damage to thousands of traditionally built houses and triggered hundreds of landslides, which were the main causes of fatalities. The earthquakes have clearly demonstrated trends of increasing seismic risk in El Salvador due to rapid population expansion in areas of high shaking and landslide hazard, exacerbated by deforestation and uncontrolled urbanisation. The institutional mechanisms required for the control of land use and building practice are very weak and present a major obstacle to risk mitigation.     

Relationships between tissue concentrations of chlorinated hydrocarbons (polychlorinated biphenyls and chlorinated pesticides) and antioxidative responses of marine mussels, Perna viridis

Marine mussels, Perna viridis, were transplanted from a reference site to various polluted sites around Hong Kong. After 30 d of exposure, antioxidative responses in the gills and hepatopancreas and tissue concentrations of chlorinated hydrocarbons [polychlorinated biphenyls (PCBs) and chlorinated pesticides (CPs)] were determined for individual mussels. Glutathione S transferase (GST) and glutathione (GSH) were positively correlated with tissue PCB concentrations. Only one of the enzymatic antioxidants, glutathione peroxidase (GPx), showed significant response to tissue PCB. No significant correlation was found between tissue concentrations of chlorinated hydrocarbons and other enzymatic antioxidants (superoxide dismutase (SOD), catalase (CAT), glutathione reductase (GR) and NADPH DT-diaphorase (DT-d). Oxidative stress, measured as thiobarbituric acid reactive substances, was correlated with chlorinated pesticide concentrations in tissues. This study demonstrated a correlation between GST/GSH and chlorinated hydrocarbons. The apparent lack of correlation between trace organic pollutants and some of the enzymatic antioxidants may be due to the inhibitory effects caused by these chemicals. The above results suggest that more investigations are needed before these enzymes can be used as biomarkers.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.     

Simultaneous radio and optical observations of the mid-latitude atmospheric response to a major geomagnetic storm of 6–8 April 2000

Basic properties of the mid-latitude traveling ionospheric disturbances (TIDs) during the maximum phase of a major magnetic storm of 6–8 April 2000 are shown. Total electron content (TEC) variations were studied by using data from GPS receivers located in Russia and Central Asia. The nightglow response to this storm at mesopause and termospheric altitudes was also measured by optical instruments FENIX located at the observatory of the Institute of Solar-Terrestrial Physics (51.9°N,103.0°E), and MORTI located at the observatory of the Institute of Ionosphere (43.2°N, 77.0°E). Observations of the O (557.7 and 630.0 nm) emissions originating from atmospheric layers centered at altitudes of 90 and 250 km were carried out at Irkutsk and of the O₂(b¹∑_g⁺−X³∑_g⁻) (0-1) emission originating from an atmospheric layer centered at altitude of 94 km was carried out at Almaty. Our radio and optical measurement network observed a storm-induced solitary large-scale wave with duration of 1 h and a wave front width of no less than 5000 km, while it traveled equatorward with a velocity of 200 m/s from 62°N to 38°N geographic latitude. The TEC disturbance, basically displaying an electron content depression in the maximum of the F2 region, reveals a good correlation with growing nightglow emission, the temporal shift between the TEC and emission variation maxima being different for different altitudes. A comparison of the auroral oval parameters with dynamic spectra of TEC variations and optical 630 nm emissions in the frequency range 0.4–4 mHz (250–2500 s periods) showed that as the auroral oval expands into mid-latitudes, also does the region with a developed medium-sale and small-scale TEC structure.     

Riparian alder fens — source or sink for nutrients and dissolved organic carbon? — 2. Major sources and sinks

Natural riparian forest wetlands are known to be effective in their ability to remove nitrate by denitrification and sediments with attached phosphorus via sedimentation. On the other hand, litter input and decomposition is a process of crucial importance in cycling of nitrogen and phosphorus in a forest ecosystem.In this study we investigated the amount of nitrogen and phosphorus entering the alder fen ecosystem through leaf litter and its decomposition and the removal capacity of nitrogen and phosphorus by measuring denitrification and sedimentation in the alder fen.We found an average input of leaf litter during fall 1998 of 226 g m⁻² yr⁻¹ DW with nutrient concentration of 0.17% P and 1.6% N. This means a yearly input of 0.4 g m⁻² yr⁻¹ P and 3.6 g m⁻² yr⁻¹ N. The decomposition of leaf litter using litter bags with small and large mesh size resulted in bags with macroinvertebrates (large mesh size) and without macroinvertebrates (small mesh size). After 57 days the litter bags with macroinvertebrates had a decomposition rate of 79%.Denitrification was measured in May and June of 1997 using the acetylene inhibition technique on intact soil cores and slurry-experiments. The average annual denitrification rate was 0.2 g m⁻² yr⁻¹ N using data from the core experiments. The denitrification rate was higher after addition of nitrate, indicating that denitrification in the riparian alder fen is mainly controlled by nitrate supply.The sedimentation rate in the investigated alder fen ranged from 0.47 kg m⁻² yr⁻¹ DW to 4.46 kg m⁻² yr⁻¹ DW in 1998 depending on the study site and method we used. Sedimentation rates were lower in newly designed plate traps than in cylinder traps. The alder fen also showed lower rates than the adjacent creek Briese. Average phosphorus removal rate was 0.33 g m⁻² yr⁻¹ P.Input sources for the surface water of the alder fen are sediment mineralization and decomposition of leaf litter; output sources are sedimentation and denitrification. This study showed that a nutrient input of 24.58 kg ha⁻¹ yr⁻¹ N, 8.8 kg ha⁻¹ yr⁻¹ P and 419 kg ha⁻¹ yr⁻¹ DOC into the surface water of the alder fen is possible. Alder fens cannot improve water quality of an adjacent river system. This is only true for a nearly pristine alder fen with the hydrology of 10 months flooded conditions and 2 months non-flooding conditions a year.     

Paleolatitudes of the Kerguelen hotspot: new paleomagnetic results and dynamic modeling

The Kerguelen Plateau, a Large Igneous Province in the southern Indian Ocean, was formed as a product of the Kerguelen hotspot in several eruptive phases during the last 120 Myr. We obtained new paleolatitudes for the central and northern Kerguelen Plateau from paleomagnetic investigations on basalts, which were drilled during ODP Leg 183 to the Kerguelen Plateau-Broken Ridge. The paleolatitudes coincide with paleolatitudes from previous investigations at the Kerguelen Plateau and Ninetyeast Ridge (the track of the Kerguelen hotspot) and indicate a difference between paleolatitudes and present position at 49°S of the Kerguelen hotspot. We show that true polar wander, the global motion between the mantle and the rotation axis, cannot explain this difference in latitudes. We present numerical model results of plume conduit motion in a large-scale mantle flow and the resulting surface hotspot motion. A large number of models all predict southward motion between 3° and 10° for the Kerguelen hotspot during the last 100 Myr, which is consistent with our paleomagnetic results.     

How to obtain alert velocity thresholds for large rockslides

A reliable forecast of the failure stage of large rockslides is difficult, because of non-linear time dependency of displacements and seasonal effects. Aim of this paper is to suggest a practical method to prepare alert thresholds for large rockslides, assessing critical values of velocity for carrying out civil protection actions using monitoring data. Adopted data concern the 20 Mm³ Ruinon rockslide (Valfurva, Central Alps, Italy), still evolving and suitable to originate a fast moving rock avalanche. Multitemporal analysis of aerial photos, LIDAR-ALTM laser topography, field survey and geomechanical analyses allowed to infer the rockslide kinematics and better understand data provided by a monitoring network including distometers, extensometers, GPS benchmarks and inclinometers. The analysis of displacement and rainfall data over five years (1997–2001) allowed to recognise three different evolutionary patterns of displacements, showing a continuously increasing rate since 1997. Data representing large-scale behaviour of the rock mass were fitted by power-law curves, according to the “accelerating creep” model by Voight, in order to evaluate a suitable failure time. This was hampered by the large seasonal deviations, which can significantly delay the occurrence of failure. Data were fitted using the Voight’s equation, expressed in terms of displacement, through non-linear estimation techniques, in order to find values of the controlling parameters (A, α and t_f) suitable to represent the mechanical behaviour of the rock mass approaching the failure. This allowed to compute velocity–time theoretical curves and to define different velocity threshold values (pre-alert, alert and emergency) to be used for emergency management.