P–T evolution of kyanite eclogite from the Pirin Mountains (SW Bulgaria): implications for the Rhodope UHP Metamorphic Complex

A new occurrence of kyanite eclogite in the Pirin Mountains of southwestern Bulgaria within the rocks belonging to the Obidim Unit of the Rhodope Metamorphic Complex is presented. This eclogite provides important information about the peak–pressure conditions despite strong thermal overprint at low pressure. Textural relationships, phase equilibrium modelling and conventional geothermobarometry were used to constrain the metamorphic evolution. Garnet porphyroblasts with inclusions of omphacite (up to 43 mol.% Jd), phengite (up to 3.5 Si p.f.u.), kyanite, polycrystalline quartz, pargasitic amphibole, zoisite and rutile in the Mg‐rich cores (X_Mg = 0.44–0.46) record a prograde increase in P–T conditions from ～2.5 GPa and 650 °C to ～3 GPa and 700–750 °C. Maximum pressure values fall within the stability field of coesite. During exhumation, the peak–pressure assemblage garnet + omphacite + phengite + kyanite was variably overprinted by a lower pressure one forming symplectitic textures, such as diopside + plagioclase after omphacite and biotite + plagioclase after phengite. The development of spinel (X_Mg = 0.4–0.45) + corundum + anorthite assemblage in the kyanite‐bearing domains at ～1.1 GPa and 800–850 °C suggests a thermal overprint in the high‐pressure granulite facies stability field. This thermal event was followed by cooling at ～0.8 GPa under amphibolite facies conditions; retrograde kelyphite texture involving plagioclase and amphibole was developed around garnet. Our results add to the already existing evidence for ultra high pressure (UHP) metamorphism in the Upper Allochthon of the Rhodope Metamorphic Complex as in the Kimi Unit and show that it is more widespread than previously known. Published age data and field structural relations suggest that the Obidim Unit represents Variscan continental crust involved into the Alpine nappe edifice of the Rhodopes and that eclogite facies metamorphism was Palaeozoic, in contrast to the Kimi Unit where age determinations suggest a Jurassic or Cretaceous age for UHP metamorphism. This implies that UHP metamorphism in the Upper Allochthon of the Rhodopes may have occurred twice, during Alpine and pre‐Alpine orogenic events, and that two independent HP/UHP provinces of different age overlap in this area.     

Internal testing of a numerical model of hillslope–channel coupling using laboratory flume experiments

The numerical model COUP 2D simulates the hydrological coupling between hillslopes and the river channel during a rainfall event. In order to test the numerical model, a 1:100 scaled laboratory flume which was modified to incorporate lateral hillslope elements, was used to run a series of experiments in which hillslope angle, channel angle, hillslope discharge and channel discharge were the varying parameters. Overall, there were 18 different experimental configurations with three replicates carried out for each condition, leading to a total of 54 experiments. These conditions were then used to parameterize and run COUP 2D. Internal model outputs of flow depth and flow velocity at four cross‐sections in the channel were compared to the measurements made in the physical model for the same parameter conditions. Statistical comparisons of the measured and modelled data were carried out for each experiment and across all experiments, using two goodness‐of‐fit measures—root mean square error and Nash–Sutcliffe coefficient of efficiency—in order to assess the performance of the model over an entire simulation as well as over all the simulations. The main effects on the goodness‐of‐fit measures for flow depth of each experimental variable, as well as the interactions between variables, were evaluated using statistical modelling. The results show that the model captures flow‐depth variations in response to changing channel and hillslope parameters. Statistical modelling suggests that the main effects on model error are cross‐section position, channel angle and channel discharge. Significant interactions also occur between all the channel variables and between the channel variables and hillslope discharge. The results of the testing procedure have significant implications for the consideration of different model components and for the interaction between data‐ and model evaluation. Copyright © 2007 John Wiley & Sons, Ltd.     

Uncertainty in modelling of faecal coliform overland transport associated with manure application in Maryland

Concerns for microbial safety of surface water facilitate development of predictive models that estimate concentrations and total numbers of pathogen and indicator organisms leaving manure‐fertilized fields in overland flow during runoff events. Spatial variability of bacterial concentrations in applied manure introduces high uncertainty in the model predictions. The objective of this work was to evaluate the uncertainty in model predictions of the manure‐borne bacteria overland transport caused by limited information on the spatial distribution of bacteria in surface‐applied manure. Experiments were carried out at the ARS Beltsville experimental watershed site (OPE3) in Maryland. Dairy bovine manure was applied at a 59·3 t/ha rate on the 3·55 hectare experimental field. Faecal coliform (FC) concentrations in manure measured in 2004, 2005, 2007, and 2009 varied by 4 orders of magnitude each year. Both runoff volume and FC concentrations in runoff water were monitored using a runoff flume equipped with a refrigerated pump sampler. Two runoff events occurred before the manure was incorporated into the soil. A bacteria transport add‐on module simulator of transport with infiltration and runoff (STWIR) was linked with the event‐based kinematic runoff and erosion model (KINEROS2) to simulate convective‐dispersive overland transport, bacteria release from manure, reversible attachment–detachment to soil, and surface straining of infiltrating bacteria. The model was successfully calibrated with the field experiment data. Monte Carlo simulations were carried out to account for the spatial variation in FC in applied manure and uncertainty in the FC distribution in manure caused by the small number of samples. A tenfold and twofold variation in FC concentrations in the runoff were obtained within the 90% probability interval when initial FC spatial distributions in the manure were represented by 5 and 29 samples, respectively. Copyright © 2011 John Wiley & Sons, Ltd.     

Nonlinear bacterial load‐streamflow response on a marine beach

Bacterial concentration (Escherichia coli) is used as the key indicator for marine beach water quality in Hong Kong. For beaches receiving streamflow from unsewered catchments, water quality is mainly affected by local nonpoint source pollution and is highly dependent on the bacterial load contributed from the catchment. As most of these catchments are ungauged, the bacterial load is generally unknown. In this study, streamflow and the associated bacterial load contributed from an unsewered catchment to a marine beach, Big Wave Bay, are simulated using a modelling approach. The physically based distributed hydrological model, MIKE‐SHE, and the empirical watershed water quality model (Hydrological Simulation Program – Fortran) are used to simulate streamflow and daily‐averaged E. coli concentration/load, respectively. The total daily derived loads predicted by the model during calibration (June–July 2007) and validation (July–October 2008) periods agree well with empirical validation data, with a percentage difference of 3 and 2%, respectively. The simulation results show a nonlinear relationship between E. coli load and rainfall/streamflow and reveal a source limiting nature of nonpoint source pollution. The derived load is further used as an independent variable in a multiple linear regression (MLR) model to predict daily beach water quality. When compared with the MLR models based solely on hydrometeorological input variables (e.g. rainfall and salinity), the new model based on bacterial load predicts much more realistic E. coli concentrations during rainstorms. Copyright © 2014 John Wiley & Sons, Ltd.     

Snow accumulation distribution inferred from time‐lapse photography and simple modelling

The spatial and temporal distribution of snow accumulation is complex and significantly influences the hydrological characteristics of mountain catchments. Many snow redistribution processes, such as avalanching, slushflow or wind drift, are controlled by topography, but their modelling remains challenging. In situ measurements of snow accumulation are laborious and generally have a coarse spatial or temporal resolution. In this respect, time‐lapse photography shows itself as a powerful tool for collecting information at relatively low cost and without the need for direct field access. In this paper, the snow accumulation distribution of an Alpine catchment is inferred by adjusting a simple snow accumulation model combined with a temperature index melt model to match the modelled melt‐out pattern evolution to the pattern monitored during an ablation season through terrestrial oblique photography. The comparison of the resulting end‐of‐winter snow water equivalent distribution with direct measurements shows that the achieved accuracy is comparable with that obtained with an inverse distance interpolation of the point measurements. On average over the ablation season, the observed melt‐out pattern can be reproduced correctly in 93% of the area visible from the fixed camera. The relations between inferred snow accumulation distribution and topographic variables indicate large scatter. However, a significant correlation with local slope is found and terrain curvature is detected as a factor limiting the maximal snow accumulation. Copyright © 2010 John Wiley & Sons, Ltd.     

The journal Communications in Numerical Methods in Engineering with Biomedical Applications becomes the International Journal for Numerical Methods in Biomedical Engineering (IJNMBE) from 1st January 2010

Communications in Numerical Methods in Engineering, founded by Roland W. Lewis in 1985, will change its title to the ‘International Journal for Numerical Methods in Biomedical Engineering’ and has a revised Aims and Scope. Copyright © 2009 John Wiley & Sons, Ltd.     

Modelling the non‐linear hydrological behaviour of a small Mediterranean forested catchment

A progressive perceptual understanding approach was used to identify a model structure able to represent the non‐linear behaviour of the hydrological cycle in a small intermittent Mediterranean stream. The initial lumped model structure consisting of a series of four connected water tanks (LU3) progressed to a model with five tanks (LU4), and finally to a semidistributed model structure (SD4) in which spatial variability of the evapotranspiration according to the vegetation cover and to the local aspect was considered. In the final model structure, which gave the best fit (Nash–Sutcliffe efficiency index = 0·78), an additional tank representing the riparian zone was included (SD4‐R). Results showed that the abrupt changes of the riparian water table during summer and the formation of a perched water table during the transition from dry to wet conditions were the main mechanisms leading to the non‐linear hydrological behaviour. The transpiration process from the saturated zone and the spatial variability of evapotranspiration resulted in key factors successfully representing the annual water balance. The spatial and temporal validations carried out for each of the four model structures considered in this study supported the hypothesis adopted during the calibration process. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical modelling of the potential effects of a dam on a coastal aquifer in S. Spain

This study presents the results of a three‐dimensional variable‐density numerical modelling of the Motril‐Salobreña coastal aquifer and the possible effects of the entry into service in May 2005 of the Rules Dam, located just 17 km from the coast. Present parameters of the Motril‐Salobreña aquifer show that the system's conditions are very similar to a natural regime. The dam will substantially alter aquifer recharge, as the entry flow through the alluvial sediments of the Guadalfeo River will be entirely cut off or drastically reduced. Different scenarios reproducing the possible evolution of the aquifer under operation of the Rules Dam have been modelled. In most cases, results indicate that the conditions of the aquifer would worsen, with a general advance of the freshwater–saltwater interface. The area with most risk of saltwater intrusion is the old mouth of the Guadalfeo River, where the mixing zone could advance 1200 m inland. It is proposed that maintaining a 5–6 Mm³ year^?1 ‘ecological flow’ in the Guadalfeo River could prevent this saline advance. This application demonstrates that variable‐density models are potentially useful tools for estimating the effects of dams on the hydrodynamic and hydrochemical conditions of a coastal aquifer. Copyright © 2009 John Wiley & Sons, Ltd.     

Remote sensing and hydrological measurement based irrigation performance assessments in the upper Amu Darya Delta,Central Asia

In the Aral Sea Basin, where the Central Asian countries compete for limited water resources, reliable information on the actual water use for eight million ha of irrigated land are rare. In this study, spatially distributed land use data, seasonal actual evapotranspiration, and reference evapotranspiration derived from multitemporal MODIS data were combined with in situ water flow measurements for irrigation performance assessments in the upper Amu Darya Delta. The functioning of the major irrigation and drainage which supplies an agricultural area of 270,000 ha in the Uzbek province Khorezm was analysed using water balancing and adequacy indicators of irrigation water use.An average relative evapotranspiration of 95% indicated fulfilled water demands and partly over-irrigation, whereas values below 75% disclosed inadequate water supply in distant parts of the irrigation system. On the other hand, immense water withdrawals of approximately 24,000 m³ ha⁻¹ recorded at the system boundaries between April and September 2005 clearly exceeded the field water demands for cotton cultivation. Only 46% of the total irrigation amounts were consumed for crop production at field level. Throughout the vegetation period, approximately 58% of the total available water left the region as drainage water. Monthly observations of the depleted fraction and the drainage ratio highlighted drainage problems and rising groundwater levels at regional scale. In the most distant downstream subsystem, a high risk of groundwater and soil salinity during the main irrigation phase was found.A combination of high conveyance losses, hydraulic problems, direct linkages between irrigation and drainage, and low field application efficiencies were identified as major reasons for underperforming irrigation. The findings underlined the necessity of water saving and of reconsidering water distribution in Khorezm. The remote sensing approach was concluded as a reliable data basis for regular performance assessments for all irrigation systems in Central Asia.