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.     

Estimating hydraulic conductivity using grain-size analyses, aquifer tests, and numerical modeling in a riverside alluvial system in South Korea

Hydraulic conductivity (K) for an alluvial system in a riverbank filtration area in Changwon City, South Korea, has been studied using grain-size distribution, pumping and slug tests, and numerical modeling. The alluvial system is composed of layers: upper fine sand, medium sand, lower fine sand, and a highly conductive sand/gravel layer at the base. The geometric mean of K for the sand/gravel layer (9.89?×?10^?4 m s^?1), as determined by grain-size analyses, was 3.33 times greater than the geometric mean obtained from pumping tests (2.97?×?10^?4 m s^?1). The geometric mean of K estimates obtained from slug tests (3.08?×?10^?6 m s^?1) was one to two orders of magnitude lower than that from pumping tests and grain-size analyses. K estimates derived from a numerical model were compared to those derived from the grain-size methods, slug tests and pumping tests in order to determine the degree of deviation from the numerical model. It is considered that the K estimates determined by the slug tests resemble the uppermost part of the alluvial deposit, whereas the K estimates obtained by grain-size analyses and pumping tests are similar to those from the numerical model for the sand/gravel layer of the riverside alluvial system.     

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.     

The contact binary system YY Eri

Relatively new and unanalysed photometric data-sets of the contact binary system YY Eri are presented. The light curves have been analysed using information limit optimization techniques, and the Binary Maker program of Bradstreet (1992).Comparison of the results allows insight into determinacy questions affecting the W UMa type of light curve, which relates to our general understanding of contact binaries. A limiting contact configuration cannot be ruled out on the basis of empirical data-analysis alone.Very careful analysis of the photospheric flux distribution over the surface of the Sun may help establish appropriate values of the gravity brightening parameter for cool dwarf stars. At present, however, empirical photometric information on contact binaries is non-discriminatory: i.e. alternative cosmogonies can find alternative support from available evidence. The implication is then for more and better observational data to allow better real independent parameter determination.     

Deformation-induced garnet zoning

Compositional zoning patterns in garnet porphyroblasts from kyanite-bearing samples of the Devonian Littleton Formation, north-central Massachusetts, reveal complex patterns of growth that are related to multiple deformation and metamorphic events. Garnet porphyroblasts exhibit asymmetrical and irregular zoning patterns in X_Mn, X_Ca and Fe/(Fe + Mg). Zoning reversals in Mn and Fe/(Fe + Mg) and patch distribution in Ca appear to occur around the boundaries of the textural zones. Also, the compositions of the garnet at the textural boundaries are variable for all traverses. These observations suggest that the garnet zoning was not only modified from diffusion processes, but was also influenced by pre-existing microfabrics through the effects of preferential dissolution and resorption in partial disequilibrium. Relationships between chemical and textural truncations indicate that the zoning patterns of garnet were strongly modified from preferential dissolution and precipitation during the development of successive foliations that occurred in zones of high strain/stress (cleavage seams) and zones of low strain/stress, respectively.