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.     

A decade of change in the Skidaway River estuary I. Hydrography and nutrients

The Skidaway River estuary is a tidally-dominated subtropical estuary in the southeastern USA surrounded by extensiveSpartina salt marshes. Weekly smapling at high and low tide began in 1986 for hydrography, nutrients, chlorophylla, particulate matter, and microbial and plankton biomass and composition; hydrographic and nutrient data during 1986–1996 are reported here. Salinity varied inversely with river discharge and exhibited variability at all time scales but with no long-term trend. Water temperature typically ranged over 25°C and was without apparent long-term frend. Seasonal cycles in concentrations of NO₃, NH₄, PO₄, Si(OH)₄, and DON were observed, with annual maxima generally occurring in late summer. Superimposed on seasonal cycles, all five nutrients exhibited steady increases in minimum, mean, and maximum concentrations; mean concentrations increased c. 50–150% during the decade. Nutrient concentrations were highly correlated with water temperature over the ten-year period, but weakly related to salinity and discharge. Nutrients were strongly correlated with one another, and the relative ratios among inorganic nutrients showed little long-term trend. Correlations among temperature and nutrient concentrations exhibited considerable inter-annual variability. Major spikes in organic and inorganic nutrient concentrations coincided with significant rainfall events; concentrations increased hyperbolically with rainfall. Although pristine compared to more heavily impacted waterways primarily outside the region, residential development and population density have been increasing rapidly during the past 15–20 years. Land use is apparently altering nutrient loading over the long-term (months-years), and superimposed on this are stochastic meteorological events that accelerate these changes over the short term (days-weeks).     

Kosa (yellow sand) components in precipitation collected at central Japan

In east Asia, acidic gases derived from fossil fuel combustion have increased in the past decades. On the other hand, the Asian dust, also called Kosa (yellow sand) is transported following windstorms from arid lands in the Asian continent. Many researchers have been interested in the reaction between acidic aerosols and Kosa aerosols as well as the long-range transport of these emissions. To investigate the characteristics of chemical components in precipitation on a long-term basis over Japan, precipitation was sequentially collected from April 1984 to March 1997 at Kanazawa located near the coast of the Sea of Japan. Precipitation samples were collected at 1 mm intervals for the first 5 mm rainfall and all volume of rainwater after 6 mm for all precipitation events with an automatic wet only precipitation collector. According to the analyses of precipitation including Kosa aerosols during Kosa periods, the reaction in the air between Kosa and acidic components during the long-range transport was discussed.     

Expansion and contraction of Chinese deserts during the Quaternary

Episodic dune formations during the Quaternary are found in many deserts of China. The causes of desert expansions on different time scales are not the same. Desert extension at about 1.1 and 0.9 Ma ago were the response to the active tectonic movements, whereas the desert evolutions on the ten-thousand years time scale were the response to the orbital scale climatic changes. Spatial scale studies on desert evolution indicate that desert margins shifted greatly during the last glacial maximum (LGM) and the Holocene optimum, its changing from 125°E of the LGM to 105°E of the climatic optimum. Historical desertification in the semiarid China is not a response to climate drought but largely associated with the human impacts (mainly over-cultivation) since about 2300 years ago, which leads to the reworking of the underlying LGM sands.