Comparison of sediment transport computations using hydrodynamic versus hydrologic models in the Simiyu River in Tanzania

This paper presents the results of a study that compares the sediment routing of the Simiyu River using the hydrologic model, Soil and Water Assessment Tool (SWAT) and the 1D hydrodynamic simulation software for Rivers and Estuaries (SOBEK-RE) model. Routing in SWAT is completed using the simplified Bagnold’s equation and in the SOBEK-RE model is undertaken using the Saint Venant equation. The upstream boundary conditions for the routing modules were derived from the subcatchments sediment yields that were estimated by SWAT using the Modified Universal Soil Loss Equation (MUSLE). The sediment loads extrapolated or interpolated from the sediment rating curve for the catchment outlet were used for calibration and validation purposes. The SWAT model predicted an erosion rate of 2.09 Mt/yr. The total sediment load transported to the main outlet of the catchment simulated by the SWAT and SOBEK-RE models was equal to 2.94 and 2.72 Mt/yr, respectively. Thus the models computed a net erosion in the channels of 0.84 Mt/yr (SWAT) and 0.63 Mt/yr (SOBEK-RE). When comparing the results of the models for the different reaches of the main channel and main tributaries, the models showed different results both in magnitude and in sign (erosion/deposition). However, in a situation where data is scarce (such as grain size, channel geometry), the more complex hydrodynamic model does not necessarily lead to more reliable results.     

Comprehensive hydrologic calibration of SWAT and water balance analysis in mountainous watersheds in northwest China

Model calibration is important for streamflow simulations using distributed hydrological models, especially in highland and cold areas of northwest China with scarce data. Quantitative analysis of water balance based on the accurate simulation is also essential for reasonably planning and managing water resources in these river basins facing a severe water shortage. In this study, a comprehensive method was proposed to calibrate the Soil and Water Assessment Tool (SWAT) model in the Yingluoxia watershed, upstream area of the Heihe River basin; it was based on multi-temporal, multi-variable and multi-site integrated drainage characteristics. Meanwhile a fresh approach of the parameter transferability and model validation was used by applying the set of calibrated parameters in its tributary to other area of the watershed. The results indicated that the method was effective and feasible; the values of Nash–Sutcliffe Efficiency (NSE) and Coefficient of Determination (r²) were greater than 0.81 and as high as 0.94 and the absolute values of the Percent Bias (PBIAS) were less than 2. Based the output of model the water balance in the Yingluoxia watershed was analyzed, that the mean annual precipitation, evapotranspiration, and discharge of the watershed from 1990 to 2000 were 491.8 mm, 334 mm, and 157.8 mm, respectively. The comprehensive calibration method based on multi-temporal, multi-variable and multi-site integrated drainage characteristics can better portray the hydrological processes of watershed and improve the model simulation; and the output of the model then provide a reliable reference for assessing and managing water resource of the watershed.     

Flow path oscillations in transient ground-water simulations of large peatland systems

Transient numerical simulations of the Glacial Lake Agassiz Peatland near the Red Lakes in Northern Minnesota were constructed to evaluate observed reversals in vertical ground-water flow. Seasonal weather changes were introduced to a ground-water flow model by varying evapotranspiration and recharge over time. Vertical hydraulic reversals, driven by changes in recharge and evapotranspiration were produced in the simulated peat layer. These simulations indicate that the high specific storage associated with the peat is an important control on hydraulic reversals. Seasonally driven vertical flow is on the order of centimeters in the deep peat, suggesting that seasonal vertical advective fluxes are not significant and that ground-water flow into the deep peat likely occurs on decadal or longer time scales. Particles tracked within the ground-water flow model oscillate over time, suggesting that seasonal flow reversals will enhance vertical mixing in the peat column. The amplitude of flow path oscillations increased with increasing peat storativity, with amplitudes of about 5 cm occurring when peat specific storativity was set to about 0.05 m⁻¹.     

Late Miocene shortening of the Northern Apennines back-arc

The inner Northern Apennines (western Tuscany and Tyrrhenian basin) is characterized by a relatively thin continental crust (∼20–25 km), high heat flow (>100 mW m⁻²), and the presence of relevant tectonic elision of stratigraphic sequences, a setting known as Serie Ridotta. These features are normally ascribed to an extensional deformation that affected the back-arc area above the subducting Adria plate since the Early-Middle Miocene (∼16 Ma). However, various geophysical studies image the continental crust to be currently affected by W-dipping thrust faults (and associated basement uplifts) that have not been obliterated by this claimed long-lasting extensional process. These observations raise the question whether the thrusts are older or younger than the continental extension. To address this question we have reprocessed and interpreted the deep seismic reflection profile CROP03/c that crosses the onshore hinterland sector, and investigated the structural setting of some of the Late Miocene-Pliocene hinterland basins (Cinigiano-Baccinello, Siena-Radicofani, Tafone, Albegna and Radicondoli basins) that are situated at the front or in-between the basement uplifts. The analysis of field structures and commercial seismic profiles has allowed the recognition that both substratum and basins’ infill have been intensely shortened. These findings and the architecture of the basins suggest that the latter developed under a contractional regime, which would have started around 8.5 Ma with the onset of the continental sedimentation. This compressive stress state followed an earlier phase of continental extension that presumably started at ∼16 Ma (with the blocking of the Corsica-Sardinia rotation), and thinned both the continental crust and sedimentary cover producing most of the Serie Ridotta. The main phases of basin shortening are bracketed between 7.5 and 3.5 Ma, and thus overlap with the increase in the exhumation rate of the metamorphic cores at ∼6–4 Ma determined through thermochronological data. We therefore propose a correlation between the basin deformation and the activity of the nearby basement thrusts, which would have thus shortened a previously thinned continental crust. This chronology of deformation may suggest a geodynamic model in which the back-arc and hinterland sector of the Northern Apennines was recompressed during Late Miocene-Early Pliocene times. This evolution may be explained through different speculative scenarios involving a blockage of the subduction process, which may vary between end members of complete slab detachment and stalled subduction.     

Streamflow variation due to glacier melting and climate change in upstream Heihe River Basin,Northwest China

Streamflow simulation is often challenging in mountainous watersheds because of incomplete hydrological models, irregular topography, immeasurable snowpack or glacier, and low data resolution. In this study, a semi-distributed conceptual hydrological model (SWAT-Soil Water Assessment Tool) coupled with a glacier melting algorithm was applied to investigate the sensitivity of streamflow to climatic and glacial changes in the upstream Heihe River Basin. The glacier mass balance was calculated at daily time-step using a distributed temperature-index melting and accumulation algorithm embedded in the SWAT model. Specifically, the model was calibrated and validated using daily streamflow data measured at Yingluoxia Hydrological Station and decadal ice volume changes derived from survey maps and remote sensing images between 1960 and 2010. This study highlights the effects of glacier melting on streamflow and their future changes in the mountainous watersheds. We simulate the contribution of glacier melting to streamflow change under different scenarios of climate changes in terms of temperature and precipitation dynamics. The rising temperature positively contributed to streamflow due to the increase of snowmelt and glacier melting. The rising precipitation directly contributes to streamflow and it contributed more to streamflow than the rising temperature. The results show that glacial meltwater has contributed about 3.25 billion m³ to streamflow during 1960–2010. However, the depth of runoff within the watershed increased by about 2.3 mm due to the release of water from glacial storage to supply the intensified evapotranspiration and infiltration. The simulation results indicate that the glacier made about 8.9% contribution to streamflow in 2010. The research approach used in this study is feasible to estimate the glacial contribution to streamflow in other similar mountainous watersheds elsewhere.     

Prediction of ground motion due to the collapse of a large-scale cooling tower under strong earthquakes

The ground motion owing to the collapse of a large-scale cooling tower under strong earthquakes was appropriately predicted using a comprehensive approach. The predicted results can be used for the safety evaluation of nuclear-related facilities adjacent to the cooling tower as well as in the planning of nuclear power plant construction in China. In this study, a cooling tower–soil model was first developed based on a falling weight–soil model, which the authors verified by falling weight tests. Then the collapse process of a cooling tower was simulated, and the collapse-induced ground vibrations were assessed by using the proposed model. Finally, the ground motion, which was a combination of the earthquake-induced ground motion and the collapse-induced ground vibrations, was estimated based on the superposition principle of waves. It was found that the cooling tower may collapse under strong earthquakes with the peak ground accelerations (PGAs) in the range of 0.35–0.45 g in x (EW) and y (NS) directions, respectively. These PGAs are far beyond the PGA range of major earthquakes in the common seismic design in China. The types of the site geologies of towers can significantly affect the collapse-induced ground vibrations. For a typical hard soil consisting of strongly weathered sandy slate, moderate ground vibrations may occur in the considered region. The collapse-induced PGAs were in the range of 0.017–0.046 g for the observed points at distances of 350 m in radial direction. For a rock-like foundation, the collapse-induced radial PGAs may be as high as 0.08 g at distances of 350 m, indicating that the effect of the collapse-induced ground vibrations on the nuclear-related facilities should be seriously assessed in certain scenarios.     

Glacier changes and their impacts on the discharge in the past half-century in Tekes watershed,Central Asia

The glacier is an important and stable water supply in Central Asia. Monitoring the change of glacier and understanding the impacts of glacier change on river discharge are critical to predict the downstream water availability change in future. Glacier changes were discussed and their impacts on river discharge were evaluated by hydrological modeling with a distributed hydrological model SWAT under two land use and land cover scenarios (1970 and 2007) in Tekes watershed, the most important source of water discharge to the Ili River. Compared to the glacier area of 1511 km² in 1970s it decreased by 332 km² in 2007, which resulted in the contribution the discharge from precipitation in the glacier area to the average annual discharge of the watershed changing from 9.8% in the period 1966–1975 to 7.8% in the period 2000–2008. In the month scale, with the decrease of glacier area, the distribution of the contribution of monthly discharge from precipitation in the glacier area to the total of the watershed changed from bimodal pattern to unimodal pattern. By linking a hydrological model to remote sensing image analysis and Chinese glacier inventories to determine glacier area change our approach in quantifying the impacts of glacier changes on hydrology at different scales, will provide quantitative information for stakeholders in making decisions for water resource management.     

The effect of biomanipulation on phosphorus exchange between sediment and water in shallow,tropical Huizhou West Lake,China

Three shallow basins in Huizhou West Lake, China, were compared with respect to phosphorus (P) cycling between sediment and water, binding forms of P in sediment, and macrophyte biomass. The basins had similar sediments and similar depths, but two of the basins were restored by carp fish removal and macrophyte transplantation. These two basins have had clear water, low Chl.a and high macrophyte coverage for seven and ten years, whilst the unrestored control basin had turbid water and higher Chl.a. Judged by diffusive ammonium efflux, sediments in restored basins had higher mineralization rates than the unrestored basin, but the release of total dissolved P were more similar. However, sediments of restored basins released primarily dissolved organic P, while the sediment from the unrestored basin only released dissolved inorganic P. One third of the P release in the unrestored basin occurred from resuspended sediment, while this pathway contributed less than 3% in restored basins where resuspension rates were 10 times lower and the surface sediments affinity for phosphate higher. Besides from the presence of carps in the unrestored basin, the main differences were a large pool of P (700–850 mg P m⁻²) in macrophyte biomass and a smaller pool (∼150 mg m⁻²) as loosely adsorbed P in the sediment of restored basins than in the unrestored (0 in macrophytes and 350 mg P m⁻² as loosely adsorbed). Also, a tendency of higher concentrations of oxidized iron was observed in the surface sediment from restored basins. The study underlines the potential of trophic structure changes to alter internal nutrient cycling in shallow lakes.     

Evaluation of SMOS soil moisture retrievals over the central United States for hydro-meteorological application

Soil moisture has been widely recognized as a key variable in hydro-meteorological processes and plays an important role in hydrological modelling. Remote sensing techniques have improved the availability of soil moisture data, however, most previous studies have only focused on the evaluation of retrieved data against point-based observations using only one overpass (i.e., the ascending orbit). Recently, the global Level-3 soil moisture dataset generated from Soil Moisture and Ocean Salinity (SMOS) observations was released by the Barcelona Expert Center. To address the aforementioned issues, this study is particularly focused on a basin scale evaluation in which the soil moisture deficit is derived from a three-layer Xinanjiang model used as a hydrological benchmark for all comparisons. In addition, both ascending and descending overpasses were analyzed for a more comprehensive comparison. It was interesting to find that the SMOS soil moisture accuracy did not improve with time as we would have expected. Furthermore, none of the overpasses provided reliable soil moisture estimates during the frozen season, especially for the ascending orbit. When frozen periods were removed, both overpasses showed significant improvements (i.e., the correlations increased from r = −0.53 to r = −0.65 and from r = −0.62 to r = −0.70 for the ascending and descending overpasses, respectively). In addition, it was noted that the SMOS retrievals from the descending overpass consistently were approximately 11.7% wetter than the ascending retrievals by volume. The overall assessment demonstrated that the descending orbit outperformed the ascending orbit, which was unexpected and enriched our knowledge in this area. Finally, the potential reasons were discussed.     

The Himalayan foreland basin crust and upper mantle

Modeling of multimode surface wave group velocity dispersion data sampling the eastern and the western Ganga basins, reveals a three layer crust with an average V_s of 3.7 km s^?1, draped by ～2.5 km foreland sediments. The Moho is at a depth of 43 ± 2 km and 41 ± 2 km beneath the eastern and the western Ganga basins respectively. Crustal V_p/V_s shows a felsic upper and middle crust beneath the eastern Ganga basin (1.70) compared to a more mafic western Ganga basin crust (1.77). Due to higher radiogenic heat production in felsic than mafic rocks, a lateral thermal heterogeneity will be present in the foreland basin crust. This heterogeneity had been previously observed in the north Indian Shield immediately south of the foreland basin and must also continue northward below the Himalaya. The high heat producing felsic crust, underthrust below the Himalayas could be an important cause for melting of midcrustal rocks and emplacement of leucogranites. This is a plausible explanation for abundance of leucogranites in the east-central Himalaya compared to the west. The uppermost mantle V_s is also significantly lower beneath the eastern Ganga basin (4.30 km s^?1) compared to the west (4.44 km s^?1).     

Hydrologic vulnerability to climate change of the Mandrone glacier (Adamello-Presanella group,Italian Alps)

In order to assess the annual mass balance of the Mandrone glacier in the Central Alps an energy-balance model was applied, supported by snowpack, meteorological and glaciological observations, together with satellite measurements of snow covered areas and albedo. The Physically based Distributed Snow Land and Ice Model (PDSLIM), a distributed multi-layer model for temperate glaciers, which was previously tested on both basin and point scales, was applied.Verification was performed with a network of ablation stakes over two summer periods. Satellite images processed within the Global Land Ice Measurements from Space (GLIMS) project were used to estimate the ice albedo and to verify the position of the simulated transient snowline on specific dates. The energy balance was estimated for the Mandrone and Presena glaciers in the Central Italian Alps. Their modeled balances (−1439 and −1503 mm w.e. year⁻¹, respectively), estimated over a 15 year period, are in good agreement with those obtained with the glaciological method for the Caresèr glacier, a WGMS (World Glacier Monitoring Service) reference located in the nearby Ortles-Cevedale group.Projections according to the regional climate model COSMO-CLM (standing for COnsortium for Small-scale MOdeling model in CLimate Mode) indicate that the Mandrone glacier might not survive the current century and might be halved in size by 2050.     

ESR chronology of alluvial deposits and first human settlements of the Middle Loire Basin (Region Centre,France)

Since 1970, many prehistoric sites have been discovered in the alluvial terraces of several rivers of the Middle Loire River basin (France). During the last decade, the stepped-terraces systems of three of the Loire tributaries, the Creuse, the Cher and the Loir Rivers, were the focus of multidisciplinary studies, involving Quaternary geology, prehistory and geochronology. The geochronological data now available for these systems were mainly obtained by electron spin resonance (ESR) method applied to optically bleached fluvial quartz. These results provide a chronological framework for the evolution of fluvial systems and the human settlements of this area during Lower and Middle Pleistocene. Two phases of human settlement could be distinguished, around 1.1 Ma (Early Palaeolithic – Mode 1) and around 0.7 Ma ago (Acheulean – Mode II).     

Water mass variation in the Mediterranean and Black Seas

The mass-induced sea level variability and the net mass transport between Mediterranean Sea and Black Sea are derived for the interval between August 2002 and July 2008 from satellite-based observations and from model data. We construct in each basin two time series representing the basin mean mass signal in terms of equivalent water height. The first series is obtained from steric-corrected altimetry while the other is deduced from GRACE data corrected for the contamination by continental hydrology. The series show a good agreement in terms of annual and inter-annual signals, which is in line with earlier works, although different model corrections influence the consistency in terms of seasonal signal and trend.In the Mediterranean Sea, we obtain the best agreement using a steric correction from the regional oceanographic model MFSTEP and a continental hydrological leakage correction derived from the global continental hydrological model WaterGAP2. The inter-annual time series show a correlation of 0.85 and a root mean square (RMS) difference of 15 mm. The two estimates have similar accuracy and their annual amplitude and phase agree within 3 mm and 23 days respectively. The GRACE-derived mass-induced sea level variability yields an annual amplitude of 27 ± 5 mm peaking in December and a trend of 5.3 ± 1.9 mm/yr, which deviates within 3 mm/yr from the altimetry-derived estimate.In the Black Sea, the series are less consistent, with lower accuracy of the GRACE-derived estimate, but still show a promising agreement considering the smaller size of the basin. The best agreement is realized choosing the corrections from WaterGAP2 and from the regional oceanographic model NEMO. The inter-annual time series have a correlation and RMS differences of 0.68 and 55 mm, their annual amplitude and phase agree within 4 mm and 6 days respectively. The GRACE-derived seawater mass signal has an annual amplitude of 32 ± 4 mm peaking in April. On inter-annual time scales, the mass-induced sea level variability is stronger than in the Mediterranean Sea, with an increase from 2003 to 2005 followed by a decrease from 2006 to 2008.Based on mass conservation, the mass-induced sea level variations, river runoff and precipitation minus evaporation are combined to derive the strait flows between the basins and with the Atlantic Ocean. At the Gibraltar strait, the net inflow varies annually with an amplitude of 52 ± 10 × 10⁻³ Sv peaking end of September (1 Sv = 10⁶ m³ s⁻¹). The inflow through the Bosphorus strait displays an annual amplitude of 13 ± 3 ×10⁻³ Sv peaking in the middle of March. Additionally, an increase of the Gibraltar net inflow (3.4 ± 0.8 × 10⁻³ Sv/yr) is detected.     

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.     

Imaging Pn velocities beneath the Pannonian basin

The Pannonian depression is an extensional back-arc basin in central Europe and is an integral part of the Alpine–Carpathian orogenic mountain belts. It can be characterized by thinned lower crust, shallow Moho discontinuity, high surface heat flow and Moho temperature, implying recent active tectonic processes. Imaging the velocity structure of the upper mantle may help us to better understand the structure and formation of the Pannonian region.In this paper, P_n traveltimes from regional earthquakes are used to tomographically image the lateral velocity variations in the uppermost mantle beneath the Pannonian basin. The set of linear tomographic equations, built up of the time term equation for each source–receiver pair, is solved by a truncated singular value decomposition algorithm. The explicit computation of the generalized inverse of the tomographic equations makes it possible to deduce both the resolution matrix and the model covariance matrix, allowing us to estimate the resolution and reliability of the solution.The mean compressional wave velocity in the uppermost mantle beneath the Pannonian basin is 7.9 km/s, substantially lower than the average continental P_n velocity of 8.1 km/s. It is mostly due to the high Moho temperature having values on average 400–500 °C more than those in the surrounding areas. The velocity anomalies range from −0.3 to 0.3 km/s relative to the mean velocity of 7.9 km/s. Due to high Moho temperature, below the North Hungarian range low (7.6–7.7 km/s) velocities can be found. High-velocity anomalies of around 8.1 km/s can be detected along the W-SW boundaries of Hungary and at the junction of the Pannonian basin and the Southern Carpathians. The Great Hungarian Plain shows average (7.9 km/s) P_n velocities.     

Resistivity imaging of Pleistocene alluvial aquifers in a contractional tectonic setting: A case history from the Po plain (Northern Italy)

In this work we present the hydrogeophysical imaging of a key sector of the Quaternary Po foreland basin (northern Italy), focussing on the reconstruction of clastic aquifers and aquitards in a complex tectono-sedimentary subsurface architecture. The study area includes the relic reliefs of Casalpusterlengo and Zorlesco, two smooth morphological features involving uplifted and gently folded Pleistocene marine to alluvial sediments, plausibly linked to the buried Northern Apennines thrust and fold belt. The geophysical data include 35 Direct Current Vertical Electrical Soundings collected over a 37 km² wide area, acquired with Schlumberger array and maximum half-spacing of 500 m. 1-D resistivity-depth profiles were computed for each VES. An integrated hydrostratigraphic approach was applied, to constrain the interpretation of the geophysical data along several cross-sections, including the comparison of resistivity soundings to stratigraphic logs, borehole electric logs and the pore-water properties.The resistivity interfaces, traceable with the same laterally continuous vertical polarity, were used to develop an electrostratigraphic model in order to portray the stacking of electrostratigraphic units down to 200 m below ground surface. Their vertical associations show a general upward increase of electrical resistivity. This assemblage mimics the regional coarsening upwards depositional trend, from the conductive units of the Plio-Pleistocene marine-to-transitional depositional systems to the resistive units of the Middle–Late Pleistocene fluvial and alluvial plain depositional systems. Middle Pleistocene depositional systems host an alternation of North-dipping, high-to-intermediate permeability aquifer systems (70–180 Ωm, thickness of 5–70 m) separated by low permeability aquitards (20–50 Ωm, thickness up to 40 m). These units pinch out against the Casalpusterlengo and Zorlesco relic reliefs, where they cover the uplifted and folded regional aquitard (20–50 Ωm) formed by Pliocene-Lower Pleistocene clays to sandy silts with gravel lenses in agreement with borehole data. In the deepest part of the local stratigraphy, a broad low-resistivity anomaly (< 10 Ωm) was clearly mapped through the study area. By comparison with electrical borehole logs in deep oil-wells, it could be interpreted as the fresh–saltwater interface due to the presence of connate waters and brines hosted by the marine-to-transitional shales.     

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.     

Thermal effects of Caledonian foreland basin formation,based on fission track analyses applied on basement rocks in central Sweden

Increasing evidence from fission track studies in Sweden indicate that large parts of the Fennoscandian Shield have been affected by a large-scale thermotectonic event in the Palaeozoic. In this study the results of 17 apatite fission track analyses from central Sweden are presented collected along three NW–SE transects trending from the Bothnian Sea to the Caledonides. On the Bothnian coast samples have been collected directly from the Sub-Cambrian Peneplain. The sedimentary cover protecting this surface until recently is responsible for the thermal increase detected through apatite fission track (FT) thermochronology.The apatite FT ages range between 516 ± 46 Ma (±1σ) on the Bothnian coast around sea level to 191 ± 11 Ma in the Caledonides (～500–1500 m.a.s.l.). The mean track lengths vary from 11.3 ± 2.2 μm (±1σ) in the east to 14.2 ± 2.8 μm in the west, indicating a longer stay in the PAZ in the east, versus a continuous cooling pattern in the west. This pattern in combination with other geological constraints indicates that the crystalline basement rocks near the Caledonian deformation front in the west experienced higher temperatures after the formation of the Sub-Cambrian Peneplain followed by denudation, compared with the basement rocks in the east near the Bothnian coast.The apatite FT data near the Caledonian deformation front indicates prevailing temperatures of more than 110 ± 10 °C prior to the Mid Palaeozoic, causing a resetting of the apatite fission track clock. The temperatures were progressively lower away from the deformation front. Apatite fission track analysis of samples collected from the Sub-Cambrian Peneplain along the Bothnian coast indicate maximum temperatures of 90 ± 15 °C during Late Silurian–Early Devonian time. This heating event is argued to be the result of burial beneath a developing foreland basin in front of the Caledonian orogeny. Assuming a geothermal gradient of 20 °C/km, this temperature increase can be converted to a total burial of the samples. The resulting geometry of this basin can be described as an asymmetrical basin at least 3.5 km deep in the vicinity of the Caledonian deformation front decreasing to about 2.5 km on the Bothnian coast, continuing further onto Finland. The width of this basin was in thus in the order of 600 km. Whether this was formed completely synorogenic or partly synorogenic, broadening after cessation of the orogeny, could not be revealed.The Late Palaeozoic and Mesozoic thermal evolution of this area is related to the extensional tectonics in the North Atlantic Domain.     

The impact of subsurface conceptualization on land energy fluxes

There is a significant body of work demonstrating the importance of hydrologic control on land energy feedbacks. Yet, quantitative data on aquifer conductivity can be difficult to assemble. Furthermore, how subsurface uncertainty propagates into land-surface processes is not well understood. This study analyzes the impact of aquifer characterization on land energy fluxes, using a coupled hydrology–land-surface model. Four gridded subsurface conductivity fields are developed for the Upper Klamath basin using two data sources and different levels of imposed heterogeneity. Each model is forced with the same transient, observed meteorology for 3 years prior to the final year presented here. Results are analyzed to quantify the impact of subsurface heterogeneity on groundwater surface water interactions and spatial patterns in hydrologic variables. Analysis shows that heterogeneity does not fundamentally alter the connection between groundwater and land surface processes. However, differences between scenarios impact the extent and location of the critical zone.     

Active segment of the 12 November 2003 Mw 5.6 earthquake at Salsipuedes oceanic basin,Gulf of California. Mexico

We analyzed records of eight seismic stations of the autonomous broadband seismograph network of a joint project between Utrecht University (the Netherlands), California Institute of Technology, and Centro de Investigación Científica y de Estudios Superiores de Ensenada (CICESE). These stations recorded the Mw 5.6 earthquake that occurred on 12 November 2003 at Salsipuedes basin in the middle of the Gulf of California 2 km west of the island Angel de la Guarda. This event was located at 29.16º N and 113.37º W, 30 km northeast of Bahia de los Angeles. A foreshock and hundreds of aftershocks were recorded in the 48 hours after its origin time. With the location of 29 earthquakes we identified the active segment, perpendicular to the main transform fault NW–SE of Canal de Ballenas, representing the transtensional boundary between the Pacific and North American plates. The direction of the active fault described is consistent with the normal fault mechanism reported by the National Earthquake Information Center (strike=39º, dip=34º, slip=–44º).From the duration magnitude of 456 aftershocks, we calculated a b-value of 1.14±0.28; furthermore, we calculated a seismic moment of (3.5 ±3.3) X10¹⁷Nm, a source radius of 3.7 ± 2.63 km, and a static stress drop of 3.94 ± 1.15 MPa (39.4 ± 11.5 bar.).