A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

Groundwater resources in the Jabal Al Hass region,northwest Syria: an assessment of past use and future potential

In many cases, the development of groundwater resources to boost agricultural production in dry areas has led to a continuous decline in groundwater levels; this has called into question the sustainability of such exploitation. In developing countries, limited budgets and scarce hydrological data often do not allow groundwater resources to be assessed through groundwater modeling. A case study is presented of a low-cost water-balance approach to groundwater resource assessments in a 1,550 km² semi-arid region in northwestern Syria. The past development of irrigated agriculture and its effect on the groundwater system were studied by analysis of Landsat images and long-term groundwater level changes, respectively. All components of the groundwater balance were determined. Groundwater recharge was estimated using the chloride mass balance method. Over the past three decades, groundwater levels have declined, on average, 23 m, coinciding with a two-fold increase in the groundwater-irrigated area. Groundwater resources are currently depleted by a value that lies between 9.5×10⁶ and 118×10⁶ m³ year^?1, which is larger than can be compensated for by a future decrease in natural discharge or changes in boundary conditions. However, groundwater resources are likely to be sufficient to supply domestic and livestock needs in the area.     

High heat flow anomalies on an old oceanic plate observed seaward of the Japan Trench

Thirty-three new measurements on the seaward slope and outer rise of the Japan Trench along a parallel of 38°45′N revealed the existence of high heat flow anomalies on the subducting Pacific plate, where the seafloor age is about 135 m.y.. The most prominent anomaly with the highest value of 114 mW/m² is associated with a small mound on the outer rise, which was reported to be a kind of mud volcano. On the seaward slope of the trench, heat flow is variable: high (70–90 mW/m²) at some locations and normal for the seafloor age (about 50 mW/m²) at others. The spatial variation of heat flow may be related to development of normal faults and horst/graben structures due to bending of the Pacific plate before subduction, with fluid flow along the fault zones enhancing the vertical heat transfer. Possible heat sources of the high heat flow anomalies are intra-plate volcanism in the last several million years like that discovered recently on the Pacific plate east of the Japan Trench.     

Giant ancient landslide in the Alma water gap (Crimean Mountains,Ukraine): notes to the predisposition,structure, and chronology

Large-scale ancient landslides of the area of more than 5 km² and volume exceeding 200 × 10⁶ m³ are characteristic features of the valleys incised in the northern periphery of the Crimean Mountains (Ukraine). The largely affected area is located in the outermost cuesta range of the Crimean Mountains which consists of rigid Sarmatian limestones overlying weak Middle Miocene and Upper Palaeogene deposits. A giant landslide arose in the Alma water gap as a reflection of several coincident preparatory factors such as suitable bedrock stratification, smectite-rich bedrock exposed to swelling activity, presence of faults parallel to the valley trend, and river capture event which preceded the landslide event. The occurrence of such ancient megaslides is particularly interesting in the area which is characterized by low precipitation (<500 mm/year) and weak contemporary seismicity. It probably reflects a more dynamic environment in humid phases of the Holocene; however, seismic triggering along the Mesozoic suture zone cannot be rejected. Compressional features such as gravitational folds in the central and distal parts of the landslide, which probably correlate with the whole landslide genesis or its significant reactivation, arose, according to the radiocarbon dating, during the Holocene climatic optimum in the Atlantic period. The slope deformation has been relatively quiescent since that time, except minor historic reactivization which took place in the frontal part of the landslide. We suppose that the studied landslide could be classified as a transitional type of slope deformation with some signs of spreading and translational block slides.     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

Evolution of cool skin and direct air-sea gas transfer coefficient during daytime

Absorption of solar radiation within the thermal molecular sublayer of the ocean can modify the temperature difference across the cool skin as well as the air-sea gas transfer. Our model of renewal type is based on the assumption that the thermal and diffusive molecular sublayers below the ocean surface undergo cyclic growth and destruction, the heat and gas transfer between the successive burst events are performed by molecular diffusion. The model has been upgraded to include heating due to solar radiation. The renewal time is parameterized as a function of the surface Richardson number and the Keulegan number. A Rayleigh number criterion characterizes the convective instability of the cool skin under solar heating. Under low wind speed conditions, the solar heating can damp the convective instability, strongly increasing the renewal time and correspondingly decreasing the interfacial gas exchange. In the ocean, an additional convective instability caused by salinity flux due to evaporation becomes of importance in such cases. The new parameterization is compared with the cool skin data obtained in the western equatorial Pacific during the Tropical Ocean Global Atmosphere Coupled Ocean Atmosphere Response Experiment in February 1993. In combination with a model of the diurnal thermocline it describes main features of the field data both in nighttime and daytime. Under low wind speed conditions (< 5 m s^-1) diurnal variations of the sea surface temperature due to the formation of a diurnal thermocline were substantially larger than those across the cool skin. Under wind speeds > 5 m s^-1, diurnal variations of the surface temperature due to the variations of the thermal molecular sublayer become more important.     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Simco,a simulation software for CONICA,the coude nir camera for the very large telescope

CONICA is an acronym for COudé Near Infrared Camera. It is one of the four currently planned Infrared instruments for ESO's Very Large Telescope (VLT) in Chile. This multimode instrument is to be installed at the Coudé-focus (of Unit Telescope no. 1), where adaptive optics and speckle interferometry will also be available. High angular resolution imaging (to the diffraction limit) will be possible in the 1–5 m range, as well as spectroscopy with low- and medium resolution and polarimetry by means of Wollaston prism and wiregrid analysers. Various softwares are developed for this instrument. One of them: the Simulation Software has a threefold aim: 1. provide the user with feedback information on his/her choices of observational parameters. This is be achieved by displaying the calculated performance and throughput of the combined Source-Atmosphere-Telescope-Camera-Detector system in various formats, such as images, tables, isophotical images ... for point as well as for extended sources (also annular, double, etc.); 2. verify if his expectations are realistic, before actually using CONICA itself. 3. give feedback on the design to the developers of the instrument.On leave from his original institute, now at ESTEC in the ISO-project