Factors controlling sediment yield in a major South American drainage basin: the Magdalena River,Colombia

The Magdalena River, a major fluvial system draining most of the Colombian Andes, has the highest sediment yield of any medium-sized or large river in South America. We examined sediment yield and its response to control variables in the Magdalena drainage basin based on a multi-year dataset of sediment loads from 32 tributary catchments. Various morphometric, hydrologic, and climatic variables were estimated in order to understand and predict the variation in sediment yield. Sediment yield varies from 128 to 2200 t km⁻² yr⁻¹ for catchments ranging from 320 to 59,600 km². The mean sediment yield for 32 sub-basins within the Magdalena basin is ∼690 t km⁻² yr⁻¹. Mean annual runoff is the dominant control and explains 51% of the observed variance in sediment yield. A multiple regression model, including two control variables, runoff and maximum water discharge, explains 58% of the variance. This model is efficient (ME=0.89) and is a valuable tool for predicting total sediment yield from tributary catchments in the Magdalena basin. Multiple correlations for those basins corresponding to the upper Magdalena, middle basin, Eastern Cordillera, and catchment areas greater than 2000 km², explain 75, 77, 89, and 78% of the variance in sediment yield, respectively. Although more variance is explained when dataset are grouped into categories, the models are less efficient (ME<0.72). Within the spatially distributed models, six catchment variables predict sediment yield, including runoff, precipitation, precipitation peakedness, mean elevation, mean water discharge, and relief. These estimators are related to the relative importance of climate and weathering, hillslope erosion, and fluvial transport processes. Time series analysis indicates that significant increases in sediment load have occurred over 68% of the catchment area, while 31% have experienced a decreasing trend in sediment load and thus yield. Land use analysis and increasing sediment load trends indicate that erosion within the catchment has increased over the last 10–20 years.     

The zooplankton communities of small water reservoirs with different trophic conditions in two catchments in western Slovakia

Our study summarizes data from six small water reservoirs in West Slovakia and analyzes the occurrence of zooplankton groups in relation to physico-chemical and catchment variables. The reservoirs are in two different catchments – of the Morava and Váh rivers. A total of 103 species were identified; 64 crustaceans (in both the pelagic and littoral zones) and 39 planktonic Rotifera in the pelagic zones. Significant differences were observed in species richness, abundance and biomass of planktonic crustaceans: 48 species were characteristic of the Váh catchment, while 53 were found in the Morava catchment. The density of zooplankton in the three reservoirs of the Váh River catchment ranged from 102 ind L^?1 to 21,488 ind L^?1 and the zooplankton biomass ranged from 0.12 mg L^?1 to 103.29 mg L^?1. The density of zooplankton in three Morava River catchment reservoirs ranged from 2 ind L^?1 to 3928 ind L^?1 and the zooplankton biomass ranged from 0.1 mg L^?1 to 27.3 mg L^?1. The differences were found to be related to catchment (altitude and catchment affiliation), chemical (BOD5, DO) and biological (Chromophyta, Chlorophyta) factors. Eutrophication of reservoirs in the Váh catchment was mainly due to agriculture and fish management, resulting in high nutrient concentrations. Species richness showed an unimodal response to BOD₅ and N-NH₄ with near optimum low values, 4.6 and 0.19 respectively. The relationship to oxygen content reflects preferences for less eutrophic waters and species richness tended to decrease with increasing DO and to decrease with increasing nutrient content.     

Simulating wind-affected snow accumulations at catchment to basin scales

In non-forested mountain regions, wind plays a dominant role in determining snow accumulation and melt patterns. A new, computationally efficient algorithm for distributing the complex and heterogeneous effects of wind on snow distributions was developed. The distribution algorithm uses terrain structure, vegetation, and wind data to adjust commonly available precipitation data to simulate wind-affected accumulations. This research describes model development and application in three research catchments in the Reynolds Creek Experimental Watershed in southwest Idaho, USA. All three catchments feature highly variable snow distributions driven by wind. The algorithm was used to derive model forcings for Isnobal, a mass and energy balance distributed snow model. Development and initial testing took place in the Reynolds Mountain East catchment (0.36 km²) where R² values for the wind-affected snow distributions ranged from 0.50 to 0.67 for four observation periods spanning two years. At the Upper Sheep Creek catchment (0.26 km²) R² values for the wind-affected model were 0.66 and 0.70. These R² values matched or exceeded previously published cross-validation results from regression-based statistical analyses of snow distributions in similar environments. In both catchments the wind-affected model accurately located large drift zones, snow-scoured slopes, and produced melt patterns consistent with observed streamflow. Models that did not account for wind effects produced relatively homogenous SWE distributions, R² values approaching 0.0, and melt patterns inconsistent with observed streamflow. The Dobson Creek (14.0 km²) application incorporated elevation effects into the distribution routine and was conducted over a two-dimensional grid of 6.67 × 10⁵ pixels. Comparisons with satellite-derived snow-covered-area again demonstrated that the model did an excellent job locating regions with wind-affected snow accumulations. This final application demonstrated that the computational efficiency and modest data requirements of this approach are ideally suited for large-scale operational applications.     

Baseflow prediction in a data-scarce catchment with Inselberg topography,Central Mozambique

This study aimed to improve the understanding of hydrological processes in a humid (sub)tropical area in Africa with Inselberg topography. Additionally, the study intended to develop an approach for selective discharge data acquisition to determine water availability for smallholder irrigation in similar data-scarce catchments. During the December 2012–August 2013 field campaign meteorological and river stage data were collected at the Messica catchment in Central Mozambique. The 220 km² catchment has an estimated 1000 ha of irrigated land, developed by smallholder farmers. Baseflow in the perennial tributary streams on the slopes of a meta-sedimentary Inselberg is the source of irrigation water. The baseflow recession curve of one of these tributaries is analysed and the water balance of an average year was determined. Precipitation, potential evapotranspiration, actual evapotranspiration and discharge were estimated to be 1224, 1462, 949 and 266 mm/year respectively. Differential gauging showed that the perennial tributaries gain water; the groundwater contribution increased with approximately 50% over two and a half month relative to the downstream discharge from March to May. In the downstream parts the groundwater contribution per metre stream length is between 30% and 100% higher compared to the upstream parts for two of the tributaries. Nevertheless, due to natural streambed infiltration and irrigation canals, discharge varies over the length of these tributaries. A rainfall–runoff model (HBV) was calibrated using the field data to examine the relation between precipitation characteristics and discharge at the start of the dry season. For precipitation scenarios with low and high intensity precipitation, discharges from June onwards were approximately similar in size according to the calibrated model. This suggest that discharge at the start of the dry season is mainly determined by total precipitation and the timing of precipitation (i.e. early or late in the wet season), not by individual rainfall events or rainfall intensity. It is concluded that the use of selective discharge measurements and low frequency precipitation measurements can effectively be used for water availability assessments in Inselberg catchments. Further research should be conducted to verify the validity of the used techniques in other humid sub-tropical Inselberg areas.     

Estimation of sub-annual inter-catchment groundwater flow using short-term water balance method

Predicting inter-catchment groundwater flow (IGF) is essential because IGF greatly affects stream water discharge and water chemistry. However, methods for estimating sub-annual IGF and clarifying its mechanisms using minimal data are limited. Thus, we quantified the sub-annual IGF and elucidated its driving factors using the short-term water balance method (STWB) for three forest headwater catchments in Japan (named here catchment A, B and As). Our previous study using the chloride mass balance indicated that annual IGF of catchment A (49.0 ha) can be negligible. Therefore, we calculated the daily evapotranspiration (ET) rate using the Priestley–Taylor expression and the 5-year water balance in catchment A (2010–2014). The sub-annual IGF of the three catchments was then calculated by subtracting the ET rate from the difference between rainfall and stream discharge during the sub-annual water balance periods selected using the STWB. The IGF rates of catchment B (7.0 ha), which is adjacent to catchment A, were positive in most cases, indicating that more groundwater flowed out of the catchment than into it, and exhibited positive linear relationships with rainfall and stream discharge. This suggested that as the catchments became wetter, more groundwater flowed out of catchment B. Conversely, the IGF rates of catchment As (5.3 ha), included in catchment A, were negative in most cases, indicating that more groundwater flowed into the catchment than out from it, and exhibited negative linear relationships with rainfall and stream discharge. Given the topography of the catchments studied, infiltration into the bedrock was the probable reason for the IGF outflow from catchment B. We hypothesized that in catchment As, the discrepancy between the actual hydrological boundary and the surface topographic boundary could have caused an IGF inflow. This study provides a useful tool for determining an IGF model structure to be incorporated into rainfall-runoff models.     

Impact of water withdrawals from groundwater and surface water on continental water storage variations

Humans have strongly impacted the global water cycle, not only water flows but also water storage. We have performed a first global-scale analysis of the impact of water withdrawals on water storage variations, using the global water resources and use model WaterGAP. This required estimation of fractions of total water withdrawals from groundwater, considering five water use sectors. According to our assessment, the source of 35% of the water withdrawn worldwide (4300 km³/year during 1998–2002) is groundwater. Groundwater contributes 42%, 36% and 27% of water used for irrigation, households and manufacturing, respectively, while we assume that only surface water is used for livestock and for cooling of thermal power plants. Consumptive water use was 1400 km³/year during 1998–2002. It is the sum of the net abstraction of 250 km³/year of groundwater (taking into account evapotranspiration and return flows of withdrawn surface water and groundwater) and the net abstraction of 1150 km³/year of surface water. Computed net abstractions indicate, for the first time at the global scale, where and when human water withdrawals decrease or increase groundwater or surface water storage. In regions with extensive surface water irrigation, such as Southern China, net abstractions from groundwater are negative, i.e. groundwater is recharged by irrigation. The opposite is true for areas dominated by groundwater irrigation, such as in the High Plains aquifer of the central USA, where net abstraction of surface water is negative because return flow of withdrawn groundwater recharges the surface water compartments. In intensively irrigated areas, the amplitude of seasonal total water storage variations is generally increased due to human water use; however, in some areas, it is decreased. For the High Plains aquifer and the whole Mississippi basin, modeled groundwater and total water storage variations were compared with estimates of groundwater storage variations based on groundwater table observations, and with estimates of total water storage variations from the GRACE satellites mission. Due to the difficulty in estimating area-averaged seasonal groundwater storage variations from point observations of groundwater levels, it is uncertain whether WaterGAP underestimates actual variations or not. We conclude that WaterGAP possibly overestimates water withdrawals in the High Plains aquifer where impact of human water use on water storage is readily discernible based on WaterGAP calculations and groundwater observations. No final conclusion can be drawn regarding the possibility of monitoring water withdrawals in the High Plains aquifer using GRACE. For the less intensively irrigated Mississippi basin, observed and modeled seasonal groundwater storage reveals a discernible impact of water withdrawals in the basin, but this is not the case for total water storage such that water withdrawals at the scale of the whole Mississippi basin cannot be monitored by GRACE.     

Hydrological effects on gravity and correlations between gravitational variations and level of the Alzette River at the station of Walferdange,Luxembourg

The gravitational effects of water storage variations driven by local precipitation events are modeled for the Walferdange Underground Laboratory for Geodynamics (WULG) in the Grand Duchy of Luxembourg. A modified mass continuity model is implemented, which uses rainfall data from Walferdange as input. In the absence of soil moisture and groundwater level information, the model is empirically parameterized. Model outputs are compared to the gravity time series registered with the Observatory Superconducting Gravimeter CT040 located in the WULG. We find that the model explains 77% of the gravity residuals. In addition, a statistical analysis is carried out to determine the relationship between precipitation, gravity variations and water level changes in the nearby Alzette River. A time delay of 88 ± 34 min between the maximum variation rates of the water level and gravity signal has been calculated. The signals have an admittance of 45 ± 5 cm μGal^?1.     

Hydro-climatic trends in the Abay/Upper Blue Nile basin,Ethiopia

Trends of the three hydro-meteorological variables precipitation, temperature and stream flow, represented by 13, 12, and 9 gauging stations, respectively, within the Abay/Upper Blue Nile basin have been studied to support water management in the region. The Trends were evaluated over different time periods depending on data availability at the stations. The statistical Mann–Kendall and Pettitt tests have been used to assess trends and change points respectively. The tests have been applied to mean annual, monthly, seasonal, 1- and 7-days annual minimum and maximum values for streamflow, while mean annual, monthly and seasonal timescales were applied to meteorological variables. The results are heterogeneous and depict statistically significant increasing/decreasing trends. Besides, it showed significant abrupt change of point upward/downward shift for streamflow and temperature time series. However, precipitation time series did not show any statistically significant trends in mean annual and seasonal scales across the examined stations.Increasing trends in temperature at different weather stations for the mean annual, rainy, dry and small rainy seasons are apparent. The mean temperature at Bahir Dar – typical station in the Lake Tana sub basin, has been increasing at the rate of about 0.5 °C/decade, 0.3 °C/decade in rainy season (June–September), 0.6 °C/decade in small rainy season (March–May), and 0.6 °C/decade in dry season (October–February). Other stations in the Abay/Upper Blue Nile show comparable results. Overall it is found that trends and change point times varied considerably across the stations and catchment to catchment. Identified significant trends can help to make better planning decisions for water management. However, the cause attributes to the observed changes in hydro-meteorological variables need further research. In particular the combined effects of land use/land cover change and climate variability on streamflow of Abay/Blue Nile basin and its tributaries needs to be understood better.     

Changes in diet and niche breadth of a small fish species in response to the flood pulse in a Neotropical floodplain lake

This study tested the hypothesis that the flood pulse affects the diet composition and the niche breadth of Moenkhausia forestii, a small characid fish inhabiting the littoral zone of lakes. To this end, we compared the diet composition (at the population and individual levels) between hydrological periods (high and low water phases) in a floodplain lake of the Upper Paraná River. PERMANOVA revealed differences in the diet between periods (pseudoF_1,38 = 8.5; p < 0.001), with predominant consumption of chironomid larvae and Ephemeroptera (aquatic resources) in the low-water period and an increase in the contribution of terrestrial resources (Hymenoptera, Coleoptera, and Orthoptera) during the high-water period. Based on the PERMDISP results, inter-individual variability in M. forestii diet also differed between periods (F_1,38 = 5.80; p = 0.02), with higher values during the high-water period resulting in a dietary niche expansion. During the low-water period, we observed the dominance of chironomid larvae in the diets of most individuals, resulting in lower inter-individual variability and thus narrower niche breadth. The diet of M. forestii was affected by the flood pulse at both the population and individual levels. The most important difference was found in the origin of food items; during the low-water period, the diet consisted mainly of aquatic resources, and during the high-water period, there was a large contribution of terrestrial resources. This variation is related to the increased availability of allochthonous resources in the high period, when terrestrial areas are flooded by the overflow of the river, thereby increasing the input of resources into the aquatic environment. The increased availability of food resources during this period allowed the expansion of the trophic niche of M. forestii, accompanied by the highest richness (19 items) and the highest evenness of food items. Our findings demonstrated that the flood pulse affected the composition of the M. forestii diet at both the population and individual levels. These results support the importance of the flood pulse, which connects aquatic and terrestrial ecosystems, in providing food resources for fish.     

Optically stimulated luminescence dating of sediments from the Yellow River terraces in Lanzhou: Tectonic and climatic implications

The stratigraphic chronology of Yellow River terraces was investigated and studied in Lanzhou Basin, western Chinese Loess Plateau. The optically stimulated luminescence (OSL) dating results show that terraces T₁, T₂ and T₃ formed at 8 ka, 20 ka and 70 ka, respectively. Lateral accretion of the riverbed facies gravel sediments occurred during interglacial periods while vertical aggradations of the terrace sediments deposited predominantly under cold and dry glacial period. A thick layer of aeolian loess with a basal age about 35 ka indicates a remarkable drop of air temperature and a dry, cold climate. The temporal correlation between terrace formation and tectonic movement has not yet been established in this research, but the stratigraphic chronology of the terrace sections provides the timing of the terrace formation, the incision rate of the Yellow River, and the slip rate of the fault horizon.     

Solar activity imprints in tree ring width from Chile (1610–1991)

We have investigated the solar activity signal in tree ring data from two locations in Chile. The tree ring time series extended over a period of ∼400 yr. Spectral and wavelet analysis techniques were employed. We have found evidence for the presence of the solar activity Schwabe (∼11 yr), Hale (∼22 yr), fourth-harmonic of the 208-yr Suess cycle (∼52 yr) and Gleissberg (∼80 yr) cycles. The Gleissberg cycle of tree ring data is in anti-phase with solar activity. Wavelet and cross-wavelet techniques revealed that the periods found are intermittent, possibly because solar activity signals observed in tree rings are mostly due to solar influence on local climate (rainfall, temperature, and cloud cover) where trees grow up. Further, cross-wavelet analysis between sunspot and tree ring time series showed that the cross power around the 11 yr solar cycle is more significant during periods of high solar activity (grand maximum) than during periods of low solar activity (grand minimum). As Glaciar Pio XI is practically at the Pacific Ocean level, the tree-ring response may be stronger due to the heating of the Pacific Ocean water following an increase of the solar radiation incidence rather than at the higher altitudes of Osorno region.     

Seasonal and spatial variation in suspended matter,organic carbon,nitrogen, and nutrient concentrations of the Senegal River in West Africa

The Senegal River is of intermediate size accommodating at present about 3.5 million inhabitants in its catchment. Its upstream tributaries flow through different climatic zones from the wet tropics in the source area in Guinea to the dry Sahel region at the border between Senegal and Mauritania. Total suspended matter, particulate and dissolved organic carbon and nitrogen as well as nutrient concentrations were determined during the dry and wet seasons at 19 locations from the up- to downstream river basin. The aims of the study were to evaluate the degree of human interference, to determine the dissolved and particulate river discharges into the coastal sea and to supply data to validate model results. Statistical analyses showed that samples from the wet and dry season are significantly different in composition and that the upstream tributaries differ mainly in their silicate and suspended matter contents. Nutrient concentrations are relatively low in the river basin, indicating low human impact. Increasing nitrate concentrations, however, show the growing agriculture in the irrigated downstream areas. Particulate organic matter is dominated by C4 plants during the wet season and by aquatic plankton during the dry season. The total suspended matter (TSM) discharge at the main gauging station Bakel was about 1.93 Tg yr⁻¹ which is in the range of the only available literature data from the 1980s. The calculated annual discharges of particulate organic carbon (POC), dissolved organic carbon (DOC) and dissolved organic nitrogen (DON) are 55.8 Gg yr⁻¹, 54.1 Gg yr⁻¹, and 5.3 Gg yr⁻¹, respectively. These first estimates from the Senegal River need to be verified by further studies.     

Non-steady long-term uplift rates and Pleistocene marine terrace development along the Andean margin of Chile (31°S) inferred from 10Be dating

Pleistocene uplift of the Chilean coast is recorded by the formation of wave-cut platforms resulting from marine erosion during sea-level highstands. In the Altos de Talinay area (~ 31°S), we have identified a sequence of 5 wave-cut platforms. Using in situ produced ¹⁰Be exposure ages we show that these platforms were formed during interglacial periods at 6, 122, 232, 321 and 690 ka. These ages correspond to marine isotopic stages (MIS) or substages (MISS) 1, 5e, 7e, 9c and 17. Shoreline angle elevations used in conjunction with our chronology of wave-cut platform formation, illustrate that surface uplift rates vary from 103 ± 69 mm/ka between 122 and 6 ka, to 1158 ± 416 mm/ka between 321 and 232 ka. The absence of preserved platforms related to the MIS 11, 13 and 15 highstands likely reflects slow uplift rates during these times. We suggest that since 700 ka, the Altos de Talinay area was predominantly uplifted during 2 short periods following MIS 17 and MISS 9c. This episodic uplift of the Chilean coast in the Pleistocene may result from subduction related processes, such as pulses of tectonic accretion at the base of the forearc wedge.     

Scattering coefficients in the mantle revealed from the seismogram envelope analysis based on the multiple isotropic scattering model

At longer periods, scattered ScS waves sometimes dominate over coda waves at large lapse times. Examining recordings of seismic envelopes at 9 IRIS seismic stations of regional earthquakes with focal depths deeper than 150 km in periods from 1 to 20 s for a wide lapse time range up to 2000 s, we found significant frequency dependence. The coda decay gradient at short periods is steeper than that at longer periods; however, the change of coda gradient associated with the ScS arrival becomes distinct as the period becomes longer. In particular, a clear offset of coda amplitude appears in central Asia for 10 and 15 s period bands. The multiple isotropic scattering process of S-waves in the heterogeneous mantle can be simply simulated by using the Monte Carlo simulation method based on the radiative transfer theory in scattering media. Assuming a two-plane-layer attenuation structure and smoothed velocity model of the PREM, we estimated the average total scattering coefficients of S-waves such as 7.52 × 10^− 4∼1.32 × 10^− 3 km^− 1 and 2.08 × 10^− 4∼6.23 × 10^− 4 km^− 1 at 4 s, and 4.51 × 10^− 4∼7.37 × 10^− 4 km^− 1 and 2.80 × 10^− 5∼2.71 × 10^− 4 km^− 1 at 10 s, for the lithosphere and the upper mantle and for the lower mantle, respectively. Our results indicate that scattering occurs mostly in the lithosphere and the upper mantle and support that medium heterogeneity spreads over the whole mantle though its scattering power is small. Strong scattering occurs beneath central Asia and Papua New Guinea, whereas the scattering beneath Italy and regions of east Russia is much weaker. The numerical calculation enables us to confirm that much stronger scattering than intrinsic attenuation causes the offset behavior with coda decay gradient change after the ScS arrival for 4 and 10 s period bands in some regions.     

Winter variability of physical processes and sediment-transport events on the Eel River shelf,northern California

A 5-yr data set of near-bed current and suspended-sediment concentration measured within 2 m of the seabed in 60-m water depth has been analyzed to evaluate the interannual variability of physical processes and sediment transport events on the Eel River continental shelf, northern California. This data set encompasses a wide range of shelf conditions with winter events characterized as: Major Flood (1996/97), strong El Niño (1997/98), strong La Niña (1998/99), and Major Storm (1999/00). Data were collected at a site located 25 km north of the Eel River mouth, on the landward edge of the mid-shelf mud deposit. During the winter months sediment resuspension is forced primarily by near-bed oscillatory flows, and sediment transport occurs both as suspended load and as gravity-driven (fluid-mud) flows. Winter conditions that caused periods of increased sediment transport existed on average for 142 d yr⁻¹ over the total record, ranging between 89 d in the Major Flood year (1996/97) and 171 d in the La Niña year (1998/99). Hourly averaged values of significant wave height varied between 0.5 and 10.7 m and hourly averaged values of near-bed orbital velocities ranged between 0 and 125 cm s⁻¹. During the five winters, sediment threshold conditions were exceeded an average of 35% of the time, ranging from 19% in the Major Flood year (1996/97) to 52% in the La Niña year (1998/99). Mean concentration of suspended sediment, measured at 30 cmab, ranged from values close to 0–8 g l⁻¹. Among winters, major sediment flux events exhibited different patterns due to varying combinations of physical processes including river floods, waves, and shelf circulation. Within winters, the major period of sediment flux varied from a 3-d fluid mud event (Major Flood winter) to a 50-d period of persistent southerlies (El Niño winter) and a winter of continuous storm cycles (La Niña winter). Winter-averaged suspended-sediment concentration appeared to vary in response to river discharge, while total sediment flux responded to storm intensity. The net sediment flux appeared to depend on timing of river discharge and shelf conditions. On the Eel River shelf, the mid-shelf mud deposit apparently is not emplaced by deposition from the river plume, but by secondary processes from the inner shelf including off-shelf transport of sediment suspensions and gravity-driven fluid-mud flows. Thus, these inner-shelf processes redistribute sediment supplied by the Eel River (a point source) making the inner shelf a line source of sediment that forms and nourishes the mid-shelf deposit. Large-scale shelf circulation patterns and interannual variability of the physical forcing are also important in determining the locus of the mid-shelf deposit, and both are influenced by climate variations. Post-depositional alteration of the deposit also depends on the subsequent shelf conditions following major floods.     

Quantitative assessment of the impact of climate variability and human activities on runoff changes: a case study in four catchments of the Haihe River basin,China

Quantitative evaluation of the effect of climate variability and human activities on runoff is of great importance for water resources planning and management in terms of maintaining the ecosystem integrity and sustaining the society development. In this paper, hydro‐climatic data from four catchments (i.e. Luanhe River catchment, Chaohe River catchment, Hutuo River catchment and Zhanghe River catchment) in the Haihe River basin from 1957 to 2000 were used to quantitatively attribute the hydrological response (i.e. runoff) to climate change and human activities separately. To separate the attributes, the temporal trends of annual precipitation, potential evapotranspiration (PET) and runoff during 1957–2000 were first explored by the Mann–Kendall test. Despite that only Hutuo River catchment was dominated by a significant negative trend in annual precipitation, all four catchments presented significant negative trend in annual runoff varying from ?0.859 (Chaohe River) to ?1.996 mm a^?1 (Zhanghe River). Change points in 1977 and 1979 are detected by precipitation–runoff double cumulative curves method and Pettitt's test for Zhanghe River and the other three rivers, respectively, and are adopted to divide data set into two study periods as the pre‐change period and post‐change period. Three methods including hydrological model method, hydrological sensitivity analysis method and climate elasticity method were calibrated with the hydro‐climatic data during the pre‐change period. Then, hydrological runoff response to climate variability and human activities was quantitatively evaluated with the help of the three methods and based on the assumption that climate and human activities are the only drivers for streamflow and are independent of each other. Similar estimates of anthropogenic and climatic effects on runoff for catchments considered can be obtained from the three methods. We found that human activities were the main driving factors for the decline in annual runoff in Luanhe River catchment, Chaohe River catchment and Zhanghe River catchment, accounting for over 50% of runoff reduction. However, climate variability should be responsible for the decrease in annual runoff in the Hutuo River catchment. Copyright © 2012 John Wiley & Sons, Ltd.     

Selected trace metals (As,Cd and Hg) distribution and contamination in the coastal wetland sediment of the northern Beibu Gulf,South China Sea

Contamination with As, Cd and Hg, their spatial and temporal distribution are reported from the coastal wetland sediments of the northern Beibu Gulf, South China Sea. The content of As, Cd, Hg and TOC in surface sediments is 8.1 ± 5.8 μg g^?1, 0.08 ± 0.14 μg g^?1, 0.034 ± 0.028 μg g^?1 and 0.45 ± 0.39%, respectively. The mean sedimentation rates are 0.93–1.37 cm year^?1 during 1920s to 2008 determined by ²¹⁰Pb and ¹³⁷Cs dating in three cores. The vertical profiles of As, Cd and Hg content in the cores retrieved from Qin and Nanliu River estuaries show increasing trends during 1985–2008 due to anthropogenic impact caused by local economic development. Locally the surface sediments have potential ecological risk of As to benthos according to the NOAA sediment quality guidelines.     

Modern sediment accumulation on the Po shelf,Adriatic Sea

Sediment accumulation over the past century on the continental shelf near the Po delta varies with distance from the most active distributary channels. Near the Pila and Goro distributaries, sediment accumulation is rapid (1–4 cm yr⁻¹) and occurs in pulses. In these areas, the seabed is dominated by physical sedimentary structures that can be related to flood sedimentation. Between the two distributaries and in the southern portion of the dispersal system, sediment accumulation is slower (rates reach a minimum of 0.23 cm yr⁻¹ at ∼50 km from the Pila mouth) and steady-state, reflecting more continuous dispersal of sediment during non-flood periods. Sedimentary strata in these locations are composed of finer (clayey silt), mottled sediment. The similarity in the spatial distribution of long-term (100-yr) sediment accumulation to deposition resulting from the 2000 flood event suggests that the Po shelf is flood-dominated.About half of the sediment delivered by the Po River on a 100-yr time scale can be accounted for in the seabed deposit within ∼50 km of the Pila mouth. The remaining sediment is likely transported southward by the prevailing circulation, and this sediment coalesces with inputs from the Apennine Rivers.     

Late Holocene intermediate water variability in the northeastern Atlantic as recorded by deep-sea corals

The Nd isotopic composition of the aragonite skeleton of fossil deep-sea corals (Lophelia pertusa, Madrepora oculata and Desmophyllum dianthus) located in the northeastern Atlantic at water depths between 635 and 1300 m was investigated to reconstruct changes in the Atlantic mid-depth gyre circulation during the past millennium. The coral εNd values varied systematically from ? 11.8 to ? 14.4 during the past 1500 years, reflecting variations in seawater εNd and thus water mass provenance. Low εNd values (εNd = ? 14) occurred during the warm Medieval Climatic Anomaly (MCA) (between 1000 AD and 1250 AD) and during the most recent period (1950 AD to 2000 AD), interrupted by a period of significantly higher εNd values (~?12.5) during the Little Ice Age (LIA) (between 1350 AD and 1850 AD). One long-lived branching coral even recorded an abrupt systematic rise from low to high εNd values around 1250 AD over the course of its 10-year growth period.These variations are interpreted to result from variable contributions of the subpolar and subtropical Atlantic intermediate water masses, which today are characterized by εNd values of ? 15 and ~?11, respectively. The low εNd values observed during the warm MCA and during recent times imply a strong eastward extension of the mid-depth subpolar gyre (SPG) induced by a dominant positive phase of the North Atlantic oscillation (NAO). During the LIA, water from the subtropical gyre (STG) and potentially from the Mediterranean Sea Water (MSW) propagated further northward, as indicated by the higher coral εNd values. This pattern suggests a negative mean state of the NAO during the LIA, with weaker and more southerly located Westerlies and a westward contraction of the SPG. Variations in the contributions of the two gyres imply changes in the heat and salt budgets at intermediate depths during the past millennia that may have contributed to changes in the properties of North Atlantic inflow into the Nordic Seas and thus deep-water formation.     

Trends in the levels of Escherichia coli in commercially harvested bivalve shellfish from England and Wales, 1999–2008

Temporal trends in Escherichia coli concentrations in bivalve shellfish were examined using data collected from 57 production areas around the coast of England and Wales during 1999–2008. Downward trends were detected in annual geometric means of E. coli in shellfish from 12% of the sampling points. The percentage of class B areas (E. coli ⩽ 4600/100 g shellfish in 90% of samples) increased from 69% to 86% during the 10-year period. The improvement in the microbial quality of shellfish is associated with sewerage improvement schemes largely implemented during 2000–2005. Upward trends were detected in 9% of the points. The causes of these increases are not known. It is recommended that quantitative sanitary profiling of shellfish waters and cost–benefit appraisal over long-term planning horizons are considered as part of sewerage investment programmes under the Water Framework Directive. This would allow greater scope to secure protection and improvement of shellfish water quality.