The development of new analytic elements for transient flow and multiaquifer flow

We deal in this paper with an ongoing development of the analytic element method. We present in outline new analytic line elements that are suitable to model general flow fields, i.e., flow fields that possess a continuously varying areal inflow or outflow. These elements are constructed specifically to model the leakage through leaky layers that separate aquifers in leaky systems and to model transient effects. The leakage or release from storage underneath linear features is modeled precisely by the new elements; the singularity in leakage is matched exactly by the approximate solution. Applications are given for a problem involving leakage and for a case of transient flow. We note that the analytic elements can be used also to reproduce the effect of continuously varying aquifer properties, e.g., the hydraulic conductivity or the elevation of the base of the aquifer. In the latter case, the elements would reproduce the rotation of the flow field caused by the variation in properties, rather than the divergence as for the case of leakage.     

The product spectra of gravity and barometric pressure in Europe

A total of 111 000 hourly values of gravity and barometric pressure from stations in Europe is analysed. The data consist of two sets of records from Brussels, an early set of 36 000 h length and a more recent set of 21 000 h length, a set of records from Bad Homburg of 24 000 h length, and a set of records from Strasbourg of 30 000 h length. All of the gravity measurements were made with similar superconducting instruments and the pressure data were recorded simultaneously at each superconducting gravimeter site. The four sets of records have different time bases, and to bring out common features and suppress individual station systematic errors, the product spectrum is introduced. Spectral density estimates are first computed for a common spectral window for each record, and the product spectrum is formed by multiplying individual spectral estimates across records. The cumulative distribution function is found for the product spectrum and confidence intervals are calculated from it by iteration. The product spectrum in gravity reveals a triplet of resonances in the subtidal band which are shown by an automated computer search to be uniquely associated with the translational modes of the solid inner core. The product spectrum in barometric pressure clearly reveals the first 10 solar heating tides in the atmosphere, but otherwise does not show common features with the gravity product spectrum. In particular, the triplet of resonances in the subtidal band of the gravity product spectrum do not show up in the product spectrum of barometric pressure, climinating the atmosphere as their source.     

A geologic study of the Juan del Grijalva landslide,the most important mass movement during the last century in Mexico

On November 4th 2007, along the Grijalva River in the state of Chiapas, Mexico, has occurred one of the largest landslides ever known. This landslide, known as Juan del Grijalva, destroyed the town of the same name, killing 20 people, and moved 55 million cubic meters of rock and debris down slope to completely block the Grijalva River. In order to understand the characteristics and factors that triggered the Juan del Grijalva landslide, geologic studies were conducted at the site. The results indicate that the landslide was composed of a lithologic sequence of thin-bedded shales and thin to medium-thick-bedded sandstones. This was faulted into several blocks dipping in the same sense as the mass movement. The main triggering factor was the increment of the pore pressure into the lithologic unit due to water saturation after 5 days of heavy rain before the incident. According to records from the last century, the Juan del Grijalva mass movement represents one of the largest mass movements recorded all over the world. The risk conditions of the area after the landslide lead to the rapid construction of an artificial channel to drain the accumulating mass of water upstream and therefore prevent a future catastrophic inundation down stream.     

Extending the double-difference location technique to mining applications part I: Numerical study

The location of the seismic event hypocenter is the very first task undertaken when studying any seismological problem. The accuracy of the solution can significantly influence consecutive stages of analysis, so there is a continuous demand for new, more efficient and accurate location algorithms. It is important to recognize that there is no single universal location algorithm which will perform equally well in any situation. The type of seismicity, the geometry of the recording seismic network, the size of the controlled area, tectonic complexity, are the most important factors influencing the performance of location algorithms. In this paper we propose a new location algorithm called the extended double difference (EDD) which combines the insensitivity of the double-difference (DD) algorithm to the velocity structure with the special demands imposed by mining: continuous change of network geometry and a very local recording capability of the network for dominating small induced events. The proposed method provides significantly better estimation of hypocenter depths and origin times compared to the classical and double-difference approaches, the price being greater sensitivity to the velocity structure than the DD approach. The efficiency of both algorithms for the epicentral coordinates is similar.     

Vertical distribution of plankton in the skagerrak in relation to doming of the seasonal thermocline

In the central Skagerrak between Norway and Denmark, the seasonal thermocline tends to be relatively shallow (10 to 15 m) probably due to geostrophic effects. An investigation of phytoplankton distribution in this region in late August 1981 showed the presence of an intense subsurface chlorophyll maximum (up to 30 γ 1⁻ Chl a) at a temperature of <6°C, dominated by coccolithophores and the dinoflagellate, Prorocentrum minimum. Data are presented on the vertical distribution of zooplankton in relation to these phytoplankton populations.     

Marine risk assessment: linear alkylbenzenesulponates (LAS) in the North Sea

The fate of linear alkylbenzenesulponates (LAS) in estuaries and coastal areas of the North Sea has been characterized with simple environmental models. The predicted concentration range in the estuaries around the North Sea (0.9-9 microg LAS l(-1)) was validated by monitoring data (1-9 microg LAS l(-1)). In offshore sites of the North Sea, it is estimated--and experimentally verified for a few sites--that the LAS concentration is below analytical detection limit (i.e., 0.5 microg LAS l(-1)). The effects of LAS on marine organisms have been reviewed. For short-term acute tests, there was no significant difference (p = 0.83) between the mean LC50 values of freshwater and marine organisms (mainly pelagic species tested, 4.1 and 4.3 mg LAS l(-1), respectively). For longer-term chronic tests, it appeared that the sensitivity (mean no-observed effect concentration (NOEC) value) of marine and freshwater organisms (0.3 and 2.3 mg LAS l(-1), respectively) was significantly different pt-test = 0.007). The predicted no-effect-concentrations (PNEC) were 360 and 31 microg LAS l(-1), for freshwater and marine pelagic communities, respectively. Given that the maximum expected estuarine and marine concentrations are 3 to > 30 times lower than the PNEC, the risk of LAS to pelagic organisms in these environments is judged to be low.     

Numerical modeling of tide-surge interaction along Orissa coast of India

The Orissa coast of India is one of the most vulnerable regions of extreme sea levels associated with severe tropical cyclones. There was extensive loss of life and property due to the October 1999 super cyclone, which devastated large part of the Orissa coast. The shallow nature of the head bay, presence of a large number of deltas formed by major rivers of Orissa such as Mahanadi and Dhamra, and high tidal range are responsible for storm surge flooding in the region. Specifically, rising and falling tidal phases influence the height, duration, and arrival time of peak surge along the coast. The objective of the present study is to evaluate the tide-surge interaction during the 1999 Orissa cyclone by using nonlinear vertically integrated numerical models. The pure tidal solution for the head bay region of the Bay of Bengal provides the initial condition for the fine resolution nested grid Orissa model. However, the feedback from the Orissa model does not affect the head bay model as the study provides a one-way interaction. Numerical experiments are performed to study the tide-surge interaction by considering various relative phases of the tidal waves with the surge-wave produced by 1999 Orissa cyclone. The comparison, although utilizing only the limited estimates of tidal data, appears adequate to assert that the principal features are reproduced correctly.