Seismotectonic significance of the 29 January 1989 Etne earthquake, SW Norway

The integrated analysis of geological, seismological and field observations with lineament data derived from satellite images allows the identification of a possible seismogenic fault zone for an earthquake which occurred near Etne in southwestern Norway, on 29 February 1989. The hypocentre of the earthquake was located at the mid-crust at a depth of 13.8±0.9 km which is typical of small intraplate earthquakes. The Etne earthquake occurred as a result of normal faulting with a dextral strike-slip component on a NW–SE trending fault. Available geological and lineament data indicate correlation of the inferred seismogenic fault with the NW–SE trending Etne fault zone. An aeromagnetic anomaly related to the Etne fault zone forms a regional feature intersecting both Precambrian basement and allochthonous Caledonian rocks. Based on these associations the occurrence of the Etne event is ascribed to the reactivation of a zone of weakness along the Etne fault zone. Slope-instabilities developed in the superficial deposits during the Etne event demonstrate the existence of potentially hazardous secondary-effects of such earthquakes even in low seismicity areas such as southwestern Norway.     

Inter-annual variability in larval supply to populations of three invertebrate taxa in the northern California Current

We investigated sources of inter-annual variability in larval supply to crab and sea urchin populations at Bodega Head and Point Reyes in northern California. During the spring and summer upwelling seasons of the years 1992 through 1997 we monitored the weekly settlement rates of nine species of crabs and two species of sea urchins. As observed in previous studies, daily values of alongshore windstress, temperature and salinity provided evidence for the poleward flow of relatively warm, low salinity water from south of Point Reyes, an apparent retention zone, during upwelling relaxation events. In years dominated by these events (1992, 1993, 1995 and 1996) we observed that alongshore windstress, temperature and salinity were coherent and temperature was significantly correlated with cancrid crab settlement. During these years the magnitude of cancrid crab settlement and the fraction of cancrid crabs relative to other crab species settling were high. Over four years of concurrent sampling there was consistently greater cancrid crab settlement at the Point Reyes site, within the retention zone, than at Bodega Head. Settlement of non-cancrid crabs (porcellanids, grapsids, pagurids and majids) was not as closely linked to intra-annual patterns of upwelling and relaxation, possibly due to the shorter seasonal availability of larvae allowing for the influence of fewer relaxation events. Settlement of this group among years was positively correlated with environmental indicators of strong seasonal upwelling; high salinity, Bakun upwelling index and low temperature. Sea urchin settlement events were observed in June and July of 1992, 1994 and 1997 during warming periods when salinity and temperature were increasing and alongshore windstress was low. Across the six years of the study, we found that cancrid crab larvae had a more even seasonal availability than larvae of non-cancrid species, which settled in greatest numbers during the early portion of the upwelling season. Sea urchins settled in greatest numbers during the later part of the upwelling season. Together these patterns demonstrate the taxon-specific way that inter-annual variability in larval supply is forced by the coincidence of larval availability with favorable physical transport mechanisms.     

Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

Mass-dependent fractionation of quadruple stable sulfur isotope system as a new tracer of sulfur biogeochemical cycles

Sulfur isotope studies of post-Archean terrestrial materials have focused on the ratio ³⁴S/³²S because additional isotopes, ³³S and ³⁶S, were thought to carry little information beyond the well-known mass-dependent relationship among multiple-isotope ratios. We report high-precision analyses of Δ³³S and Δ³⁶S values, defined as deviations of ³³S and ³⁶S from ideal mass-dependent relationships, for international reference materials and sedimentary sulfides of Phanerozoic age by using a fluorination technique with a dual-inlet isotope ratio mass spectrometer. Measured variations in Δ³³S and Δ³⁶S are explained as resulting from processes involve branching reactions (two or more reservoirs formed) or mixing. Irreversible processes in closed systems (Rayleigh distillation) amplify the isotope effect. We outline how this new isotope proxy can be used to gain new insights into fundamental aspects of the sulfur biogeochemical cycle, including additional constraints on seawater sulfate budget and processes in sedimentary sulfide formation. The isotope systematics discussed here cannot explain the much larger variation of Δ³³S and Δ³⁶S observed in Archean rock records. Furthermore, Phanerozoic samples we have studied show a characteristic Δ³³S and Δ³⁶S relationship that differs from those measured in Archean rocks and laboratory photolysis experiments. Thus, high precision analysis of Δ³³S and Δ³⁶S can be used to distinguish small non-zero Δ³³S and Δ³⁶S produced by mass-dependent processes from those produced by mass-independent processes in Archean rocks and extraterrestrial materials.     

中国农村剩余劳动力转移的趋向和措施

在对中国农村剩余劳动力现状描述的基础上 ,对剩余劳动力形成的原因进行了分析与阐述 ;结合目前形势 ,指出未来农村剩余劳动力的转移前景不容乐观 ;并在此基础上阐述了中国剩余劳动力转移途径和主要措施     

Adaptive Multiscale Permeability Estimation

With multiscale permeability estimation one does not select parameterization prior to the estimation. Instead, one performs a hierarchical search for the right parameterization while solving a sequence of estimation problems with an increasing parameterization dimension. In some previous works on the subject, the same refinement is applied all over the porous medium. This may lead to over-parameterization, and subsequently, to unrealistic permeability estimates and excessive computational work. With adaptive multiscale permeability estimation, the new parameterization at an arbitrary stage in the estimation sequence is such that new degrees of freedom are not necessarily introduced all over the porous medium. The aim is to introduce new degrees of freedom only where it is warranted by the data. In this paper, we introduce a novel adaptive multiscale estimation. The approach is used to estimate absolute permeability from two-phase pressure data in several numerical examples.     

Constraint energy minimizing generalized multiscale finite element method in the mixed formulation

This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high-conductivity channels (which connect the boundaries of a coarse block). Next, we solve local problems to construct multiscale basis functions for the velocity field. These local problems are formulated in the oversampled domain, taking into account some constraints with respect to auxiliary spaces. The latter allows fast spatial decay of local solutions and, thus, allows taking smaller oversampled regions. The number of basis functions depends on small eigenvalues of the local spectral problems. Moreover, multiscale pressure basis functions are needed in constructing the velocity space. Our multiscale spaces have a minimal dimension, which is needed to avoid contrast dependence in the convergence. The method’s convergence requires an oversampling of several layers. We present an analysis of our approach. Our numerical results confirm that the convergence rate is first order with respect to the mesh size and independent of the contrast.     

The Ensemble Kalman Filter: theoretical formulation and practical implementation

The purpose of this paper is to provide a comprehensive presentation and interpretation of the Ensemble Kalman Filter (EnKF) and its numerical implementation. The EnKF has a large user group, and numerous publications have discussed applications and theoretical aspects of it. This paper reviews the important results from these studies and also presents new ideas and alternative interpretations which further explain the success of the EnKF. In addition to providing the theoretical framework needed for using the EnKF, there is also a focus on the algorithmic formulation and optimal numerical implementation. A program listing is given for some of the key subroutines. The paper also touches upon specific issues such as the use of nonlinear measurements, in situ profiles of temperature and salinity, and data which are available with high frequency in time. An ensemble based optimal interpolation (EnOI) scheme is presented as a cost-effective approach which may serve as an alternative to the EnKF in some applications. A fairly extensive discussion is devoted to the use of time correlated model errors and the estimation of model bias.Responsible Editor: Jörg-Olaf Wolff