Statistical tests for regional seismic phase characterizations

In seismic analysis some assumptions are often made aboutthe data, e.g. stationarity and Gaussianity. This is not obvious for all realseismic data. Here, we use statistical tests for characterization of regionalseismic data. We apply tests for stationarity, symmetry, linearity, andtime-reversibility. In the analysis we use twelve regional seismic events inFennoscandia recorded with the seismic small-aperture arrays NORESS,ARCESS and FINESA at hypocentral distances in the range from 160 to1580 km. For the tests we use the regional phases Pn, secondary P, Sn and Lg-phases and the preceding noise. Two of the eventsare explosions, two are rockbursts and eight are earthquakes. Theperformance and possibilities of using statistical tests based on bispectra, asa complementary tool for conventional analysis of seismic phases isdemonstrated. The preceding noise recorded before the first onset of theP-wave for the twelve events is tested to be spatially stationarybetween each channel within each array and temporal stationary in 21consecutive time windows of 3.2 sec each. Also, the preceding noise issymmetric and linear. The seismic phases defined by the symmetry test asnon-Gaussian (not symmetric) are all linear. This means a linear model canbe used to characterize both the noise and the phases. The first P-phase for the two explosions is characterized as non-Gaussian at thethree arrays. For all 36 possibly first P-phase arrivals at the three arraystations, 23 are non-Gaussian. The second P-phase is non-Gaussian at13 of 36 data records, the S-phase at eleven of 36 and the Lg-phase at nine of 36. For all the four phases more than 32 of possible36 are time-reversible.     

Modeling Lake Ladoga Zoobenthos and Studying Its Role in Phosphorus Exchange at the Water–Bed Interface

A zoobenthos submodel is included in an ecosystem model. The new model is used to assess the role of zoobenthos in the phosphorus exchange at the water–bed interface, in particular, the phosphorus cycle in the water body ecosystem is described more accurately. The calculated zoobenthos distribution over the bed of a lake can be used in modeling the lake ichthyocenoses.     

The BEAM-project: prediction and assessment of mixture toxicities in the aquatic environment

Freshwater and marine ecosystems are exposed to various multi-component mixtures of pollutants. Nevertheless, most ecotoxicological research and chemicals regulation focus on hazard and exposure assessment of individual substances only, the problem of chemical mixtures in the environment is ignored to a large extent. In contrast, the assessment of combination effects has a long tradition in pharmacology, where mixtures of chemicals are specifically designed to develop new products, e.g. human and veterinary drugs or agricultural and non-agricultural pesticides. In this area, two concepts are frequently used and are thought to describe fundamental relationships between single substance and mixture effects: Independent Action (Response Addition) and Concentration Addition. The question, to what extent these concepts may also be applied in an ecotoxicological and regulatory context may be considered a research topic of major importance, as the concepts would allow to make use of already existing single substance toxicity data for the predictive assessment of mixture toxicities. Two critical knowledge gaps are identified: (a) There is a lack of environmental realism, as a huge part of our current knowledge about the applicability of the concepts is restricted to artificial situations with respect to mixture composition or biological effect assessment. (b) The knowledge on what exactly is needed for using the concepts as tools for the predictive mixture toxicity assessment is insufficient. Both gaps seriously hamper the necessary, scientifically sound consideration of mixture toxicities in a regulatory context.In this paper, the two concepts will be briefly introduced, the necessity of considering the toxicities of chemical mixtures in the environment will be demonstrated and the applicability of Independent Action and Concentration Addition as tools for the prediction and assessment of mixture toxicities will be discussed. An overview of the specific aims and approaches of the BEAM project to fill in the identified knowledge gaps is given and first results are outlined.     

Earthen barriers to control lava flows in the 2001 eruption of Mt. Etna

Preceded by four days of intense seismicity and marked ground deformation, a new eruption of Mt. Etna started on 17 July and lasted until 9 August 2001. It produced lava emission and strombolian and phreatomagmatic activity from four different main vents located on a complex fracture system extending from the southeast summit cone for about 4.5 km southwards, from 3000 to 2100 m elevation (a.s.l.). The lava emitted from the lowest vent cut up an important road on the volcano and destroyed other rural roads and a few isolated country houses. Its front descended southwards to about 4 km distance from the villages of Nicolosi and Belpasso. A plan of intervention, including diversion and retaining barriers and possibly lava flow interruption, was prepared but not activated because the flow front stopped as a consequence of a decrease in the effusion rate. Extensive interventions were carried out in order to protect some important tourist facilities of the Sapienza and Mts. Silvestri zones (1900 m elevation) from being destroyed by the lava emitted from vents located at 2700 m and 2550 m elevation. Thirteen earthen barriers (with a maximum length of 370 m, height of 10–12 m, base width of 15 m and volume of 25 000 m³) were built to divert the lava flow away from the facilities towards a path implying considerably less damage. Most of the barriers were oriented diagonally (110–135°) to the direction of the flow. They were made of loose material excavated nearby and worked very nicely, resisting the thrust of the lava without any difficulty. After the interventions carried out on Mt. Etna in 1983 and in 1991–1992, those of 2001 confirm that earthen barriers can be very effective in controlling lava flows.     

Seismic scattering at the top of the mantle Transition Zone

We document strong seismic scattering from around the top of the mantle Transition Zone in all available high resolution explosion seismic profiles from Siberia and North America. This seismic reflectivity from around the 410 km discontinuity indicates the presence of pronounced heterogeneity in the depth interval between 320 and 450 km in the Earth’s mantle. We model the seismic observations by heterogeneity in the form of random seismic scatterers with typical scale lengths of kilometre size (10-40 km by 2-10 km) in a 100-140 km thick depth interval. The observed heterogeneity may be explained by changes in the depths to the α-β-γ spinel transformations caused by an unexpectedly high iron content at the top of the mantle Transition Zone. The phase transformation of pyroxenes into the garnet mineral majorite probably also contributes to the reflectivity, mainly below a depth of 400 km, whereas we find it unlikely that the presence of water or partial melt is the main cause of the observed strong seismic reflectivity. Subducted oceanic slabs that equilibrated at the top of the Transition Zone may also contribute to the observed reflectivity. If this is the main cause of the reflectivity, a substantial amount of young oceanic lithosphere has been subducted under Siberia and North America during their geologic evolution. Subducted slabs may have initiated metamorphic reactions in the original mantle rocks.     

Contemporary crustal deformation of the Chinese continent and tectonic block model

We obtain the preliminary result of crustal deformation velocity field for the Chinese continent by analyzing GPS data from the Crustal Motion Observation Network of China (CMONOC), particularly the data from the regional networks of CMONOC observed in 1999 and 2001. We delineate 9 technically active blocks and 2 broadly distributed deformation zones out of a dense GPS velocity field, and derive block motion Euler poles for the blocks and their relative motion rates. Our result reveals that there are 3 categories of deformation patterns in the Chinese continent. The first category, associated with the interior of the Tibetan Plateau and the Tianshan orogenic belt, shows broadly distributed deformation within the regions. The third category, associated with the Tarim Basin and the region east of the north-south seismic belt of China, shows block-like motion, with deformation accommodated along the block boundaries only. The second category, mainly associated with the borderland of the Tibetan Plateau, such as the Qaidam, Qilian, Xining (in eastern Qinghai), and the Diamond-shaped (in western Sichuan and Yunnan) blocks, has the deformation pattern between the first and the third, i.e. these regions appear to deform block-like, but with smaller sizes and less strength for the blocks. Based on the analysis of the lithospheric structures and the deformation patterns of the regions above, we come to the inference that the deformation modes of the Chinese continental crust are mainly controlled by the crustal structure. The crust of the eastern China and the Tarim Basin is mechanically strong, and its deformation takes the form of relative motion between rigid blocks. On the other hand, the northward indentation of the Indian plate into the Asia continent has created the uplift of the Tibetan Plateau and the Tianshan Mountains, thickened their crust, and raised the temperature in the crust. The lower crust thus has become ductile, evidenced in low seismic velocity and high electric conductivity observed. The brittle part of the crust, driven by the visco-plastic flow of the lower crust, deforms extensively at all scales. The regions of the second category located at the borderland of the Tibetan Plateau are at the transition zone between the regions of the first and the third categories in terms of the crustal structure. Driven by the lateral boundary forces, their deformation style is also between the two, in the form of block motion and deformation with smaller blocks and less internal strength.     

Fine structure of Pn velocity beneath Sichuan-Yunnan region

We use 23298 Pn arrival-time data from Chinese national and provincial earthquake bulletins to invert fine structure of Pn velocity and anisotropy at the top of the mantle beneath the Sichuan-Yunnan and its adjacent region. The results suggest that the Pn velocity in this region shows significant lateral variation; the Pn velocity varies from 7.7 to 8.3 km/s. The Pn-velocity variation correlates well with the tectonic activity and heat flow of the region. Low Pn velocity is observed in southwest Yunnan, Tengchong volcano area, and the Panxi tectonic area. These areas have very active seismicity and tectonic activity with high surface heat flow. On the other hand, high Pn velocity is observed in some stable regions, such as the central region of the Yangtze Platform; the most pronounced high velocity area is located in the Sichuan Basin, south of Chengdu. Pn anisotropy shows a complex pattern of regional deformation. The Pn fast direction shows a prominent clockwise rotation pattern from east of the Tibetan block to the Sichuan-Yunnan diamond block to southwest Yunnan, which may be related to southeastward escape of the Tibetan Plateau material due to the collision of the Indian Plate to the Eurasia Plate. Thus there appears to be strong correlation between the crustal deformation and the upper mantle structure in the region. The delay times of events and stations show that the crust thickness decreases from the Tibetan Plateau to eastern China, which is consistent with the results from deep seismic sounding.     

An experimental study of nominally symmetric and asymmetric structures isolated with the FPS

This experimental investigation deals with the earthquake behaviour of a nominally symmetric and a mass‐asymmetric three‐storey structural model isolated with the frictional pendulum system (FPS). Both accidental and natural torsion are evaluated in the structure by using recorded accelerations in all building floors and measured deformations at the isolation level. A 3D‐shaking table was used to subject the model to five different ground motions, including impulsive as well as far‐field subduction‐zone type earthquakes. Results show that the analytical predictions of the earthquake behaviour of the isolated structure, as obtained from a physical model of the FPS, are in close agreement with the true complex inelastic measured behaviour of the FPS. Besides, experimental results also validate previous observations about the importance of accounting for the variability of the normal loads in modelling the earthquake behaviour of FPS isolators. Measured torsional deformation amplifications at the base of the building vary, in the mean, from 2.5% to 6% for the symmetric and asymmetric structural configurations, respectively. In relation to the fixed base structure, the reduction factors for the base shear of the isolated structure are, in the mean, about 3.9 for both configurations. Finally, it is concluded that the FPS is capable of controlling the lateral–torsional motions of mass‐asymmetric structures quite effectively by aligning the centre of mass of the superstructure with the centre of pendular and frictional resistance of the isolation system. Copyright © 2003 John Wiley & Sons, Ltd.     

On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.