Transient foundation solution of saturated soil to impulsive concentrated loading

The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures.     

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.     

Distribution of Mantle Up-welling Determined from Plate Motions: A Case for Large-scale Bénard Convection

—?The number and geometric distribution of putative mantle up-welling centers and the associated convection cell boundaries are determined from the lithospheric plate motions as given by the 14 Euler poles of the observational NUVEL model. For an assumed distribution of up-welling centers (called “cell-cores”) the corresponding cell boundaries are constructed by a Voronoi division of the spherical surface; the resulting polygons are called “Bénard cells.” By assuming the flow-kinematics within a cell, the viscous coupling between the flow and the plates is estimated, and the Euler poles for the plates are computed under the assumption of zero-net-torque. The positions of the cell-cores are optimized for the HS2-NUVEL1 Euler poles by a method of successive approximation (“subplex”); convergence to one of many local minima occurred typically after ～20,000 iterations. Cell-cores associated with the fourteen HS2-NUVEL1 Euler poles converge to a relatively small number of locations (8 to 10, depending on interpretation), irrespective of the number of convection cells submitted for optimized distribution (from 6 to 50). These locations are correlated with low seismic propagation velocities in tomography, uniformly occur within hotspot provinces, and may specifically be associated with the Hawaiian, Iceland, Reunion/Kerguelen (Indian Ocean), Easter Island, Melanesia/Society Islands (South Pacific), Azores/Cape Verde/Canary Islands, Tristan da Cunha (South Atlantic), Balleny Islands, and possibly Yellowstone hotspots. It is shown that arbitrary Euler poles cannot occur in association with mantle Bénard convection, irrespective of the number and the distribution of convection cells. Nevertheless, eight of the observational Euler poles – including the five that are accurately determined in HS2-NUVEL1 (Australia, Cocos, Juan de Fuca, Pacific, and Philippine) – are “Bénard-valid” (i.e., can be explained by our Bénard model). Five of the remaining six observational poles must be relocated within their error-ellipses to become Bénard-valid; the Eurasia pole alone appears to be in error by ～115°, and may actually lie near 40°N, 154°E. The collective results strongly suggest Bénard-like mantle convection cells, and that basal shear tractions are the primary factor in determining the directions of the plate motions as given by the Euler poles. The magnitudes of the computed Euler vectors show, however, that basal shear cannot be the exclusive driving force of plate tectonics, and suggest force contributions (of comparable magnitude for perhaps half of the plates) from the lithosphere itself, specifically subducting slab-pull and (continental) collision drag, which are provisionally evaluated. The relationship of the putative mantle Bénard polygons to dynamic chaos and turbulent flow is discussed.     

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.     

Seismic evaluation of existing nuclear facility using ambient vibration test to characterize dynamic behavior of the structure and microtremor measurements to characterize the soil: a case study

The methods used for a building seismic hazard evaluation are presented with the associated results. The goals of the study are (1) to check the soil nature and the existence or not of a possible site effect around the installation and (2) to characterize the dynamic behavior of the building using ambient vibration records.

The results of the soil study with the Nakamura method are very difficult to interpret because they are not stable in space and time. The spectral ratios method has been used with regional earthquake records. The results of the application of this method allowed us to conclude that the installation was free of site effect.

The ambient vibration measurements on the building brought the conclusion to determine the first and second modes of the structure. These results have been used to calibrate numerical model. The modal shapes in plan (high roof) and in elevation (main column) have been evaluated. The damping of the building has been computed using ambient vibration records.     

Probabilistic Estimates of Water Level in the Ural River Mouth Reach at Interaction between the River Runoff and Positive Setups

A probabilistic technique is developed for assessing water level in the mouth reach of the Ural River in the course of interaction between the river runoff and irregular rises in the Caspian Sea level (positive setups) at the predicted background sea levels of –26.5 and –26.0 m abs.     

Mantle Convection Modeling with Viscoelastic/Brittle Lithosphere: Numerical Methodology and Plate Tectonic Modeling

--The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.     

Associations Between Rural Land Uses and Ground Water Quality in the Ogallala Aquifer, Northwest Texas

Recent nitrate, chloride, and bromide concentrations were studied in the Ogallala Aquifer of northwest Texas. The study included 361 wells with a median depth of 92 m in a rural area dominated by agricultural activity and oil and gas production. Only five observations surpassed the 44.3 mg/L standard for nitrate (10 mg/L NO₃-N). Four other observations, and one from the preceding set, exceeded the secondary standard of 250 mg/L for chloride. Maximum concentrations were 91.2 mg/L, 1530 mg/L, and 0.70 mg/L for nitrate, chloride, and bromide, respectively. Chloride/bromide ratios covered a broad range, from 30.4 to 10930, but medians were < 160 for each of two years analyzed. There were statistically significant correlations between nitrate and chloride, and chloride and well depth. Results of this study suggest that agricultural activity has locally impacted ground water in north-west Texas. Regionally, low aquifer recharge rates have curtailed ground water contamination from potentially adverse land uses.     

Optimized Seismic Threshold Monitoring – Part 2: Teleseismic Processing

v--vThis second paper (Part 2) pertaining to optimized site-specific threshold monitoring addresses the application of the method to regions covered by a teleseismic or a combined regional-teleseismic network. In the first paper (Part 1) we developed the method for the general case, and demonstrated its application to an area well-covered by a regional network (the Novaya Zemlya nuclear test site). In the present paper, we apply the method to the Indian and Pakistani nuclear test sites, and show results during the periods of nuclear testing by these two countries in May 1998. Since the coverage by regional stations in these areas is poor, an optimized approach requires the use of selected, high-quality stations at teleseismic distances.¶To optimize the threshold monitoring of these test sites, we use as calibration events either one of the nuclear explosions or a nearby earthquake. From analysis of the calibration events we derive values for array beamforming steering delays, filter bands, short-term averages (STA) lengths, phase travel times (P waves), and amplitude-magnitude relationships for each station. By applying these parameters, we obtain a monitoring capability of both test sites ranging from m_b 2.8-3.0 using teleseismic stations only. When including the nearby Nilore station to monitor the Indian tests, we show that the threshold can be reduced by about 0.4 magnitude units. In particular, we demonstrate that the Indian tests on 13 May, 1998, which were not detected by any known seismic station, must have corresponded to a magnitude (m_b) of less than 2.4.¶We also discuss the effect of a nearby aftershock sequence on the monitoring capability for the Pakistani test sites. Such an aftershock sequence occurred in fact on the day of the last Pakistani test (30 May, 1998), following a large (m_b 5.5) earthquake in Afghanistan located about 1100 km from the test site. We show that the threshold monitoring technique has sufficient resolution to suppress the signals from these interfering aftershocks without significantly affecting the true peak of the nuclear explosion on the threshold trace.     

Precursory seismic crustal deformation in the area of southern Albanides

On the basis of growing evidence thatstrong earthquakes are preceded by a periodof accelerating seismicity of moderatemagnitude earthquakes, an attempt is madeto search for such seismicity pattern in NWAegean area. Accelerating seismic crustaldeformation has been identified in the areaof southern Albanides mountain range(border region between Greece, formerYugoslavia and Albania). Based on certainproperties of this activity and on itssimilarity with accelerating seismicdeformation observed before a strongearthquake which occurred in the sameregion on 26 May 1960 (M = 6.5), we canconclude that a similar earthquake may begenerated in the same region during thenext few years. This conclusion is inagreement with independent results whichhave been derived on the basis of the timepredictable model.