Dissolving the mystery of subsurface controls on snowmelt–discharge dynamics in karst mountain watersheds using hydrologic timeseries

Streamflow generation in mountain watersheds is strongly influenced by snow accumulation and melt as well as groundwater connectivity. In mountainous regions with limestone and dolomite geology, bedrock formations can host karst aquifers, which play a significant role in snowmelt–discharge dynamics. However, mapping complex karst features and the resulting surface-groundwater exchanges at large scales remains infeasible. In this study, timeseries analysis of continuous discharge and specific conductance measurements were combined with gridded snowmelt predictions to characterize seasonal streamflow response and evaluate dominant watershed controls across 12 monitoring sites in a karstified 554 km² watershed in northern Utah, USA. Immense surface water hydrologic variability across subcatchments, years and seasons was linked to geologic controls on groundwater dynamics. Unlike many mountain watersheds, the variability between subcatchments could not be well described by typical watershed properties, including elevation or surficial geology. To fill this gap, a conceptual framework was proposed to characterize subsurface controls on snowmelt–discharge dynamics in karst mountain watersheds in terms of conduit flow direction, aquifer storage capacity and connectivity. This framework requires only readily measured surface water and climatic data from nested monitoring sites and was applied to the study watershed to demonstrate its applicability for evaluating dominant controls and climate sensitivity.     

Characteristics of vertical ozone distribution in the lower troposphere in the Yangtze River Delta at Lin’an in the spring of 2001

Deterioration of air quality as a result of rapid eco- nomic development in the Yangtze River Delta (YRD) has attracted much scientific and public attention[1―7]. To explore the air quality of a region, measurement of pollutants in locations reflecting the “background” atmospheric conditions of the region is essential. The Regional Baseline Air Quality Monitoring Station of the Chinese Academy of Meteorological Sciences (CAMS) at Lin’an is located in such an area in Zhe- jiang Pro…     

Sustainable sanitation systems for low income urban areas – A case of the city of Bulawayo,Zimbabwe

Lack of basic sanitation systems threaten environmental and human health in low income urban communities. In 2005, the Government of Zimbabwe carried out a cleanup exercise in urban areas involving the destruction of illegal structures which left many people homeless. As a solution to this problem, the government embarked on an extensive housing construction exercise on unserviced land; the ‘Garikai/Hlalani Kuhle’ development programme. The objective of this paper was to investigate the sanitation status in one such area (Cowdray Park Extension, Bulawayo) and determine a sustainable sanitation system for the improved collection of wastewater from the unserviced low income urban area. The study was carried out between October 2010 and February 2011. The sanitation status as well as the residents’ preferences for improved sanitation and the economic set up of the community for the study area was determined through use of questionnaires to the residents. The local authority was then consulted to recommend sanitation facilities and system for the area that met regulatory requirements. A literature study identified sanitation options that were applicable to low income and high density urban areas. The baseline survey found that 61% of the people in the study area lacked sanitation facilities and practiced open defecation. The majority of the residents (70%) preferred ‘flush and discharge’ system sanitation facilities, which was in line with the local council’s requirements. On-site sanitation options were found not to be feasible as per the council regulations and the findings of the literature study, for areas with a high density of houses. Therefore a sewerage system was designed using the conventional sewerage design approach as well as the simplified sewerage design approach in order to determine the collection system that would best meet the needs of the community. In conclusion, the community was in dire need of a sanitation system and a waterborne offsite sanitation system was found to be a feasible option. The simplified sewerage system was found to deliver the same hydraulic benefits for collecting wastewater as the conventional but was 33% cheaper to construct and hence more affordable for the community. It was recommended that stringent environmental monitoring of the sanitation system be put in place to minimise any potential environmental impacts.     

Benthic life in karst spring—The life cycle and secondary production of benthic macroinvertebrates under the effects of constant water temperature

Because of the thermal stability and the persistence over geological time scales of their aquifers, the permanent karst springs may act as a glacial refuge, as well as serve as “in vivo” laboratories for the testing of various ecological hypotheses. For these reasons, our main hypotheses were: i) to the constant water temperature benthic species will response by specific voltinism and life-cycle length will be different from those that the species usually have in non-spring habitats; ii) under the constant water temperature the particular species and whole community will have higher biomass turnover like in similar non-spring habitats. Investigation was performed in large cave karst spring “Vyvieranie” in the Western Carpathians. In total, 40 species of benthic macroinvertebrates were identified in the studied karst spring. The trajectory of the annual change of the benthic community structure clearly showed a returnable pattern at the species density and biomass level. The constant water temperature throughout the year was reflected by specific life-history trajectories in certain benthic species. While some species had fixed voltinism (e.g. Ephemerella mucronata, Isoperla sudetica), others indicated having a flexible life-cycle pattern. In the species with a flexible life cycle, the trajectory of larval development occurred in two different ways. In the first case, a constant water temperature extended larval development (e.g. Protonemura auberti, Leuctra albida). In the second case, the addition of a winter cohort was recorded (e.g. Gammarus fossarum, Rhyacophila tristis, Protonemura austriaca). The productivity of this benthic community was considerably higher than that previously found in other submountain limestone streams in the Western Carpathians.     

Cold water anomalies in the middle layer of the northern Taiwan Strait in spring—a numerical approach

Ocean Dynamics - Cold water anomalies were found in the middle layer in the northern Taiwan Strait in the spring of 2015. The cold water was located at a depth of approximately...     

Luminescence dating of irrigation systems: Application to a qanat in Aragón,Spain

Optically stimulated luminescence (OSL) techniques have been applied to investigate the potential for dating the deposition of upcast mounds associated with qanat ventilation shafts at the site of a medieval qanat located in Aragón, Spain. Coarse quartz grains, extracted from sediment samples taken from excavated sections of several mounds, possessed sufficiently strong OSL to enable an evaluation of equivalent dose by applying the single aliquot regenerative procedure to small aliquots, each containing an individual bright grain. The OSL dates for both palaeosol and overlying upcast indicate that a chronostratigraphic record has been preserved within the mounds investigated, and micromorphological analysis of thin sections of sediment blocks taken from the mounds is shown to provide an essential means of verifying the characteristics of the strata, in particular, the critical interface of upcast and the ancient ground surface. The earliest OSL dates for basal deposits taken from two separate sections of the same mound are in agreement, placing the mound construction during the first half of the 13th century A.D. However, in two other mounds the OSL dates for the deposition of upcast are internally consistent with the stratigraphy but significantly later, dating to the 16th and 17th centuries A.D. We interpret the differences between the dates for the upcast deposition to be the result of partial erosion of the upper shaft and later repair of the mounds, and this finding underlines the importance of both examining multiple mounds in the same qanat system and the internal structure of each sampled mound. This exploratory work demonstrates the potential for wider application of OSL for dating this important type of subterranean irrigation feature in the study of both the archaeology of human settlement and palaeoenvironmental change in arid regions.     

From rainfall to spring discharge: Coupling conduit flow,subsurface matrix flow and surface flow in karst systems using a discrete–continuum model

Physics-based distributed models for simulating flow in karst systems are generally based on the discrete–continuum approach in which the flow in the three-dimensional fractured limestone matrix continuum is coupled with the flow in discrete one-dimensional conduits. In this study we present a newly designed discrete–continuum model for simulating flow in karst systems. We use a flexible spatial discretization such that complicated conduit networks can be incorporated. Turbulent conduit flow and turbulent surface flow are described by the diffusion wave equation whereas laminar variably saturated flow in the matrix is described by the Richards equation. Transients between free-surface and pressurized conduit flow are handled by changing the capacity term of the conduit flow equation. This new approach has the advantage that the transients in mixed conduit flow regimes can be handled without the Preissmann slot approach. Conduit–matrix coupling is based on the Peaceman’s well-index such that simulated exchange fluxes across the conduit–matrix interface are less sensitive to the spatial discretization. Coupling with the surface flow domain is based on numerical techniques commonly used in surface–subsurface models and storm water drainage models. Robust algorithms are used to simulate the non-linear flow processes in a coupled fashion. The model is verified and illustrated with simulation examples.     

Long-range dependence in the time series of Taiţa River discharges

ABSTRACT

The aim of this article is to offer an insight into the hydrology of southeast Romanian catchments, which have not been studied systematically. The complex behaviour of the Tai?a, one of the most important rivers situated in the northern part of Dobrogea, is demonstrated. Long-range dependence (LRD) in the time series of its streamflow, including the presence of trends, is supported by an estimated Hurst exponent of approximately 0.8, as well as statistical tests.

Editor D. Koutsoyiannis; GUEST editor S. Weijs     

Geometry and dynamics of the freshwater—seawater interface in a coastal aquifer in southeastern Spain

Abstract

The contact between freshwater and seawater in coastal aquifers is studied using a relatively simple model for homogeneous aquifers. However, for real aquifers it is not so simple. The desalination plant built to supply water to the city of Almería is situated over the aquifer in the southern part of the River Andarax Delta. Its design capacity is 1100 L s^?1, and it is supplied from boreholes pumping water from beneath the freshwater—seawater contact in this aquifer. Well logs kept over a period of two years have allowed us to accurately define the interface geometry of the freshwater—seawater contact. Lithological data collected from 31 boreholes have also indicated the existence of strata with low hydraulic conductivity, within others of high conductivity. During a simultaneous pumping test of six wells with 690 L s^?1 total discharge, electrical conductivity measurements showed the influx of seawater 6–10 m below sea level and a drawdown of the interface in the piezometers closest to the pumping wells.     

Large-scale aseismic creep in the areas of the strong earthquakes revealed from the GRACE data on the time variations of the Earth’s gravity field

The refinement of the accuracy and resolution of the monthly global gravity field models from the GRACE satellite mission, together with the accumulation of more than a decade-long series of these models, enabled us to reveal the processes that occur in the regions of large (M_w≥8) earthquakes that have not been studied previously. The previous research into the time variations of the gravity field in the regions of the giant earthquakes, such as the seismic catastrophes in Sumatra (2004) and Chile (2010), and the Tohoku mega earthquake in Japan (2011), covered the coseismic gravity jump followed by the long postseismic changes reaching almost the same amplitude. The coseismic gravity jumps resulting from the lower-magnitude events are almost unnoticeable. However, we have established a long steady growth of gravity anomalies after a number of such earthquakes. For instance, in the regions of the subduction earthquakes, the growth of the positive gravity anomaly above the oceanic trench was revealed after two events with magnitudes M_w=8.5 in the Sumatra region (the Nias earthquake of March 2005 and the Bengkulu event of September 2007 near the southern termination of Sumatra Island), after the earthquake with M_w=8.5 on Hokkaido in September 2007, a doublet Simushir earthquake with the magnitudes M_w = 8.3 and 8.1 in the Kuriles in November 2006 and January 2007, and after the earthquake off the Samoa Island in September 2009 (M_w=8.1). The steady changes in the gravity field have also been recorded after the earthquake in the Sichuan region (May 2008, M_w = 8.0) and after the doublet event with magnitudes 8.6 and 8.2, which occurred in the Wharton Basin of the Indian Ocean on April 11, 2012. The detailed analysis of the growth of the positive anomaly in gravity after the Simushir earthquake of November 2006 is presented. The growth started a few months after the event synchronously with the seismic activation on the downdip extension of the coseismically ruptured fault plane zone. The data demonstrating the increasing depth of the aftershocks since March 2007 and the approximately simultaneous change in the direction and average velocity of the horizontal surface displacements at the sites of the regional GPS network indicate that this earthquake induced postseismic displacements in a huge area extending to depths below 100 km. The total displacement since the beginning of the growth of the gravity anomaly up to July 2012 is estimated at 3.0 m in the upper part of the plate’s contact and 1.5 m in the lower part up to a depth of 100 km. With allowance for the size of the region captured by the deformations, the released total energy is equivalent to the earthquake with the magnitude M_w = 8.5. In our opinion, the growth of the gravity anomaly in these regions indicates a large-scale aseismic creep over the areas much more extensive than the source zone of the earthquake. These processes have not been previously revealed by the ground-based techniques. Hence, the time series of the GRACE gravity models are an important source of the new data about the locations and evolution of the locked segments of the subduction zones and their seismic potential.     

An assessment of an oil spill in Gladstone,Australia – Impacts on intertidal areas at one month post-spill

In January 2006, 25 tonnes of heavy fuel oil spilled into the Port of Gladstone in Queensland, Australia, from the breached hull of a bulk carrier ship. Over the following days, approximately 18 tonnes of the oil was recovered, however a certain amount of oil was deposited in the intertidal areas of Port Curtis leaving a highly visible, viscous residue. The objectives of this research were to assess the immediate impacts on the intertidal habitat and to gain baseline information for future comparative assessments.Sediment PAH and metal concentrations, mangrove communities and intertidal macroinvertebrates were assessed within one month post-spill at oil-impacted sites; adjacent sites which were not visibly impacted; and reference sites which were located outside the recorded distribution of the oil spill. Highest PAH concentrations were found at the impacted sites, with concentrations of some PAHs exceeding Australian and New Zealand Sediment Quality Guidelines (ANZECC/ARMCANZ, 2000). These sites contained very few or no crab holes in the high intertidal area, indicating a low crab density in comparison to reference sites. Little immediate impact was evident on the mangrove and macrobenthic communities, however future surveys may show evidence of longer-term impacts on these communities.     

Smoothing the source time function – a tool to soften inadequate modelling of the medium during inversion of local waveforms?

The inconsistency of records of local earthquakes and synthetic seismograms due to noise contamination of the data and/or to improper Green function because of an inexact structural model often results in high frequency oscillations of the source time function, which prolongs it spuriously. A priori limit on the roughness of the source time function turned out to reduce the spurious oscillations but it keeps its spurious extension unchanged. Thus, it may yield a severe distortion of the seismic moment. However, it reduces the spurious volumetric component of the mechanism resulting from improper modelling of the structure and, thus, may help decide whether the retrieved volumetric component is real or false. The demand for the smoothness is incorporated as a penalty function in the minimization of the least square residuals, which allows us to decide about the degree of the smoothness. The minimization is performed with a simple genetic algorithm, which offers the advantage of a detailed exploration of the model space reducing the danger of being trapped in local minima.     

Development of FLEMOcs – a new model for the estimation of flood losses in the commercial sector

Abstract

The estimation of flood damage is an important component for risk-oriented flood design, risk mapping, financial appraisals and comparative risk analyses. However, research on flood-loss modelling, especially in the commercial sector, has not gained much attention so far. Therefore, extensive data about flood losses were collected for affected companies via telephone surveys after the floods of 2002, 2005 and 2006 in Germany. Potential loss determining factors were analysed. The new Flood Loss Estimation MOdel for the commercial sector (FLEMOcs) was developed on the basis of 642 loss cases. Losses are estimated depending on water depth, sector and company size as well as precaution and contamination. The model can be applied to the micro-scale, i.e. to single production sites as well as to the meso-scale, i.e. land-use units, thus enabling its countrywide application.

Citation Kreibich, H., Seifert, I., Merz, B. & Thieken, A. H. (2010) Development of FLEMOcs – a new model for the estimation of flood losses in the commercial sector. Hydrol. Sci. J. 55(8), 1302–1314.     

Influence of the plasma zonal E × B drift on the electron concentration in the low-latitude ionospheric F region at the minimum of solar activity near the spring equinox

The three-dimensional nonstationary theoretical model of the concentrations and temperatures of electrons and ions in the ionospheric F region and plasmasphere at low and middle latitudes is used to study variations in the concentration NmF2 and height hmF2 of the ionospheric F2 layer under the action of the plasma zonal drift in the direction geomagnetic west-geomagnetic east perpendicularly to the electric E and geomagnetic B fields. The calculated and measured values of NmF2 and hmF2 for 16 ionospheric sounding stations during the quiet geomagnetic period on March 28–29, 1964 at low solar activity are compared. This comparison made it possible to correct the input parameters of the model: [O] from the NRLMSISE-00 model and the meridional component of the neutral wind velocity from the HWW90 model. It is shown that the nighttime NmF2 values decrease up to twice at low solar activity in the low-latitude ionosphere, and the hmF2 values change by up to 16 km, if the plasma zonal E×B drift is not taken into account. Under the daytime conditions, the influence of the plasma zonal E×B drift on NmF2 can be neglected.     

2-D viscoelastic FEM simulation on stress state in the deep part of a subducted slab

Introduction It is found that there are some relationships between the thermal structures of subduction zones and the deep seismicity, while the mechanism relates the thermal structure and the deep seismicity is still unsure (Helffrich, Brodholt, 1991; Furukawa, 1994; Kirby, et al, 1996). From 1980s, geoscientists have constituted a series of numerical simulations on the stress states of subduction slabs. Based on the kinetic computation of Sung and Burns (1976a, b), Goto, et al (1983, 1987…     

2-D elastic FEM simulation on stress state in the deep part of a subducted slab

Introduction Geoscientists concern the explanation of deep focus earthquakes greatly. As a great progress in understanding the formation of deep focus earthquakes, Sung and Burns (1976a, b) proposed that olivine could exist in transition zone in metastable form in the core of some cold subducting slabs, which have got supports from laboratory studies (Rubie, Ross, 1994). Iidaka and Suetsugu (1992) used the seismic data recorded by the dense seismograph networks in Japan to study the travelti…