Streamflow dynamics in a Mediterranean temporary river

Abstract

Mediterranean rivers are characterized by the irregularity of flow, harsh hydrological fluctuations and a profound transformation as the result of human activity. In this study, we investigate the streamflow response of a Mediterranean temporary river in which different groundwater, agriculture and urban contributions play an important role. Streamflow was measured at three nested gauging stations installed along Na Borges River, a lowland agricultural catchment (319 km²) on the island of Mallorca. Based on two hydrological years (2004/05 and 2005/06), potential evapotranspiration and surface water–groundwater interaction generated a succession of four different hydrological periods playing an important role in baseflow dynamics. The runoff coefficients were very low (<5%). At the event scale, groundwater also controlled runoff response, being very different according to hydrogeology, antecedent conditions and human impacts. During dry seasons, wastewater and karstic spring discharges maintain an influent regime into some streams. As a result, intense rainstorms in late summer generated water volumes over the impervious urban surfaces involved, with the result that quickflow was significant because the hydrological pathways were active. Citation Estrany, J., Garcia, C. & Alberich, R. (2010 Estrany, J., Garcia, C. and Batalla, R. J. 2010. Hydrological response of a small Mediterranean agricultural catchment. J. Hydrol., 380(1-2): 180–190. doi:10.1016/j.jhydrol.2009.10.035[Crossref] , [Google Scholar]) Streamflow dynamics in a Mediterranean temporary river. Hydrol. Sci. J. 55(5), 717–736.      

Hydrological impact and potential flooding of convective rain cells in a semi-arid environment

Abstract

t Spatio-temporal storm properties have a large impact on catchment hydrological response. The sensitivity of simulated flash floods to convective rain-cell characteristics is examined for an extreme storm event over a 94 km² semi-arid catchment in southern Israel. High space–time resolution weather radar data were used to derive and model convective rain cells that then served as input into a hydrological model. Based on alterations of location, direction and speed of a major rain cell, identified as the flooding cell for this case, the impacts on catchment rainfall and generated flood were examined. Global sensitivity analysis was applied to identify the most important factors affecting the flash flood peak discharge at the catchment outlet. We found that the flood peak discharge could be increased three-fold by relatively small changes in rain-cell characteristics. We assessed that the maximum flash flood magnitude that this single rain cell can produce is 175 m³/s, and, taking into account the rest of the rain cells, the flash flood peak discharge can reach 260 m³/s.

Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Morin, E. and Yakir, H., 2013. Hydrological impact and potential flooding of convective rain cells in a semi-arid environment. Hydrological Sciences Journal, 59 (7), 1275–1284. http://dx.doi.org/10.1080/02626667.2013.841315     

Catchment-scale vulnerability assessment of groundwater pollution from diffuse sources using the DRASTIC method: a case study

Abstract

The catchment-scale groundwater vulnerability assessment that delineates zones representing different levels of groundwater susceptibility to contaminants from diffuse agricultural sources has become an important element in groundwater pollution prevention for the implementation of the EU Water Framework Directive (WFD). This paper evaluates the DRASTIC method using an ArcGIS platform for assessing groundwater vulnerability in the Upper Bann catchment, Northern Ireland. Groundwater vulnerability maps of both general pollutants and pesticides in the study area were generated by using data on the factors depth to water, net recharge, aquifer media, soil media, topography, impact of vadose zone, and hydraulic conductivity, as defined in DRASTIC. The mountain areas in the study area have “high” (in 4.5% of the study area) or “moderate” (in 25.5%) vulnerability for general pollutants due to high rainfall, net recharge and soil permeability. However, by considering the diffuse agricultural sources, the mountain areas are actually at low groundwater pollution risk. The results of overlaying the maps of land use and the groundwater vulnerability are closer to the reality. This study shows that the DRASTIC method is helpful for guiding the prevention practices of groundwater pollution at the catchment scale in the UK.

Citation Yang, Y. S. & Wang, L. (2010 Yang, Y. S. and Wang, L. 2010. A review of modelling tools for implementation of the EU Water Framework Directive in handling diffuse water pollution. Water Resour. Manage., 24: 1819–1843.  [Google Scholar]) Catchment scale vulnerability assessment of groundwater pollution from diffuse sources using the DRASTIC method: a case study. Hydrol. Sci. J. 55(7), 1206–1216.     

Spatial variability and long-term change in pollen diversity in Nam Co catchment (central Tibetan Plateau): Implications for alpine vegetation restoration from a paleoecological perspective

Pollen diversity offers abundant clues into the floristic diversity and history of vegetation change. Few palynological studies investigated modern pollen diversity or the past floristic diversity on the Tibetan Plateau (TP). Based on modern pollen assemblages from 37 topsoils and 63 surface lake sediments in the Nam Co catchment on the central TP, this study quantitatively explored spatial distribution of modern pollen diversity using Shannon-Wiener index (H) and palynological richness (E(T _n ), n=600). Pollen diversity indices showed spatial variability among vegetation types, reflecting the differences in terrestrial floristic diversity in the lake catchment. Their values were high in the southeastern region of the lake catchment which is covered by alpine steppe, while values were low for alpine meadow and marsh meadow. The pollen diversity in lacustrine pollen assemblage could be an effective proxy to document past floristic diversity. The past floristic diversity in the lake catchment, recovered from a fossil pollen record of NMLC-1, showed a long-term change of ascending overlaid by several rapid diversity changes during the last 8400 years due to the downward shift of altitudinal vegetation belt driven by a general climatic cooling. The results imply that under the environmental challenge of climate warming and vegetation degradation, alpine vegetation restoration in the Nam Co catchment and the central TP should pay attention to altitudinal vegetation belt and zonal vegetation of alpine steppe, and use the long-term change of floristic diversity as a historical analogue.     

Effect of single and multi-site calibration techniques on hydrological model performance,parameter estimation and predictive uncertainty: a case study in the Logone catchment,Lake Chad basin

Understanding hydrological processes at catchment scale through the use of hydrological model parameters is essential for enhancing water resource management. Given the difficulty of using lump parameters to calibrate distributed catchment hydrological models in spatially heterogeneous catchments, a multiple calibration technique was adopted to enhance model calibration in this study. Different calibration techniques were used to calibrate the Soil and Water Assessment Tool (SWAT) model at different locations along the Logone river channel. These were: single-site calibration (SSC); sequential calibration (SC); and simultaneous multi-site calibration (SMSC). Results indicate that it is possible to reveal differences in hydrological behavior between the upstream and downstream parts of the catchment using different parameter values. Using all calibration techniques, model performance indicators were mostly above the minimum threshold of 0.60 and 0.65 for Nash Sutcliff Efficiency (NSE) and coefficient of determination (R ²) respectively, at both daily and monthly time-steps. Model uncertainty analysis showed that more than 60% of observed streamflow values were bracketed within the 95% prediction uncertainty (95PPU) band after calibration and validation. Furthermore, results indicated that the SC technique out-performed the other two methods (SSC and SMSC). It was also observed that although the SMSC technique uses streamflow data from all gauging stations during calibration and validation, thereby taking into account the catchment spatial variability, the choice of each calibration method will depend on the application and spatial scale of implementation of the modelling results in the catchment.     

Catchment environmental change over the 20th Century recorded by sedimentary leaf wax n-alkane δ~(13)C off the Pearl River estuary

The compound-specific stable carbon isotope compositions(δ~(13)C) of leaf wax n-alkanes from two short sediment cores recovered off the Pearl River estuary(PRE) were analyzed to check for their capability of indicating decadal scale catchment environmental change. Sedimentary long-chain n-alkanes exhibited an odd-over-even predominance, with a maximum at n-C_(29) or n-C_(31), indicating their leaf wax origin was from vascular plants. The δ~(13)C values of C_(29) and C_(31) n-alkane in all the sediment samples were in the range of -28.8‰ to -31.2‰, consistent with the C_3 plant-dominated vegetation in the Pearl River catchments. The time series of δ~(13)C records from the two cores were comparable and displayed a decreasing trend from the early 20 th century to the end of the 1970s, followed by a reversal in that change leading to continued increase for ca. 15 years. After being corrected for the effect of atmospheric CO_2 rise and δ~(13)C_(atm) decline, the δ~(13)C_(29) records largely retained their raw changing pattern; the post-1980 increase being more conspicuous. The slightly decreasing trend in corrected δ~(13)C records before around 1980 may have been caused by an increase in precipitation, whereas the subsequent increase of δ~(13)C is likely associated with the observed dry climate and/or intensive anthropogenic deforestation. Our results thus demonstrate that leaf wax n-alkanes buried in the sediments off the PRE may well reflect change in the regional climate and/or human activity in the river catchments over the past century.     

Long-period irregular pulsations under the conditions of a quiet magnetosphere

Simultaneous observations of high-latitude long-period irregular pulsations at frequencies of 2.0–6.0 mHz (ipcl) and magnetic field disturbances in the solar wind plasma at low geomagnetic activity (Kp ~ 0) have been studied. The 1-s data on the magnetic field registration at Godhavn (GDH) high-latitude observatory and the 1-min data on the solar wind plasma and IMF parameters for 2011–2013 were used in an analysis. Ipcl (irregular pulsations continuous, long), which were observed against a background of the IMF Bz reorientation from northward to southward, have been analyzed. In this case other solar wind plasma and IMF parameters, such as velocity V, density n, solar wind dynamic pressure P = ρV² (ρ is plasma density), and strength magnitude B, were relatively stable. The effect of the IMF Bz variation rate on the ipcl spectral composition and intensity has been studied. It was established that the ipcl spectral density reaches its maximum (~10–20 min) after IMF Bz sign reversal in a predominant number of cases. It was detected that the ipcl average frequency (f) is linearly related to the IMF Bz variation rate (ΔBz/Δt). It was shown that the dependence of f on ΔBz/Δt is controlled by the α = arctan(By/Bx) angle value responsible for the MHD discontinuity type at the front boundary of magnetosphere. The results made it possible to assume that the formation of the observed ipcl spectrum, which is related to the IMF Bz reorientation, is caused by solar wind plasma turbulence, which promotes the development of current sheet instability and surface wave amplification at the magnetopause.     

Solutes in overland flow following fire in eucalyptus and pine forests,northern Portugal

Small bounded plots were used to investigate dissolved nutrient losses in overland flow in eucalyptus and pine forests in northern Portugal following understorey fires. Losses of NON, POCa²⁺, Mg²⁺ and K⁺ were measured over 19 months during the first 2–3 years after fire. Solute losses in overland flow increased after fire in eucalyptus and pine terrain due to the mineralization of litter and vegetation and increased overland flow. The data suggest that (i) elevated losses persist for at least 2 and 3 years at the pine and eucalyptus sites respectively and (ii) soluble losses of POK⁺ in a post‐fire cycle may exceed those adsorbed to eroded material. Losses of POalthough relatively small, are potentially the most detrimental to soil fertility and forest productivity because of the limited opportunity for P replenishment and correspondingly high P losses adsorbed to eroded sediment. Copyright © 2000 John Wiley & Sons, Ltd.     

Variation in hydrologic connectivity as a result of microtopography explained by discharge to catchment size relationship

In hydrological terms, raised bogs are often approximated by simple models as in the acrotelm–catotelm concept. However, raised bogs are often characterized by a pronounced surface topography, causing large changes in connectivity of contributing areas on the bog. In this study, daily regression of measured discharges versus catchment areas is used to quantify the impact of surface topography on catchment connectivity within a raised bog. The resulting coefficient of determination shows the strength of the relationship between the discharge and catchment area over time under different hydrological conditions. Monitoring of discharge, water table, transmissivity, and basic weather data on a raised bog (1.9 km²) in eastern central Estonia took place from May 2008 to June 2010. Contributing areas, calculated based on the outlet's discharge volume (V _Q ) divided by the net precipitation volume ( ), of the outlet containing the central pool‐ridge system varied between 1×10^?3 and 0.7 km², suggesting significant differences in connectivity between hydrological events. Correlation between discharge and theoretical catchment size was high (R ²>0.75) when the water table was close to the surface (less than 5 cm below peat surface), and consequently, transmissivities were also high (up to 1,030m²d^?1), which led to connectivity of local storage elements, such as pools and hollows. However, a water table below this threshold resulted in large parts of the catchment being disconnected. The importance of water table depths on catchment connectivity suggests the need to reconsider the hydrological concept of raised bogs; to incorporate these shallow flow components and better understand residence time and consequently transport of solutes, such as DOC, from patterned peatlands.     

Experimental Study on Velocity Profiles with Different Roughness Elements in a Flume

The classical log law for velocity profile is applied to engineering practice. Field observations indicate that the composition of the bed materials obviously influences the shape of vertical velocity distribution. To clearly understand the roughness effect, six types of materials were laid separately at various depths for the investigation of the effects of roughness elements on the vertical velocity distribution. A down-looking 3D acoustic Doppler velocimeter was used to measure the velocity profiles. The experimental results showed that the curve characteristics of velocity profiles are strongly dependent on the roughness scale and related flow parameters. If d/R, Fr, and Re are larger than 0.15, 0.47, and 60 000, respectively, the velocity distribution may resemble an S-shape profile. The inflexion position Z*/H for a given S-shape profile was empirically deduced as Z*/H = ?0.4481d/R + 0.3225. Otherwise, the velocity profile agrees well with the logarithmic law. The findings of this study are useful in engineering practice (i.e., depth-averaged velocity and flow rate estimate). Key words: velocity profiles, acoustic Doppler velocimeter, roughness elements, logarithmic law, S-shape curve.     

Impact of climatic change on river discharge in the driest region of Poland

Abstract

Discharge in the Warta River in Poland has been analysed based on long series of measurements at the Gorzów Wielkopolski gauge station (covering the whole catchment area) and at Poznań (middle and upper catchment area), and the Note? River is characterized by the gauge station at Nowe Drezdenko. The annual mean discharge of the Warta River for the period 1981–2010 was equal to the average value for the last 163 years (209 m³ s^-1), and there was no significant change in comparison with the ratio of runoff in the summer and winter half-years. In the driest region of Poland, the climate has been described on the basis of precipitation and air temperature. The annual mean precipitation for 1981–2010 (544 mm) in the Warta River catchment area was the same as that for the period 1848–2010. The precipitation has been increasing in spring and winter, and decreasing in summer. There is a positive and very significant correlation (r = 0.705) between the annual discharge and annual precipitation totals. The annual mean air temperature has risen by 0.6°C between the periods 1848–1980 and 1981–2010.

Editor D. Koutsoyiannis

Citation Ilnicki, P., Farat, R., Górecki, K., and Lewandowski, P., 2014. Impact of climatic change on river discharge in the driest region of Poland. Hydrological Sciences Journal, 59 (6), 1117–1134. http://dx.doi.org/10.1080/02626667.2013.831979     

The generalized kinetic equation as a model for the nonlinear transfer in third-generation wave models

An alternative model for the nonlinear interaction term S_nl in spectral wave models, the so called generalized kinetic equation (Janssen J Phys Oceanogr 33(4):863–884, 2003; Annenkov and Shrira J Fluid Mech 561:181–207, 2006b; Gramstad and Stiassnie J Fluid Mech 718:280–303, 2013), is discussed and implemented in the third generation wave model WAVEWATCH-III. The generalized kinetic equation includes the effects of near-resonant nonlinear interactions, and is therefore able, in theory, to describe faster nonlinear evolution than the existing forms of S_nl which are based on the standard Hasselmann kinetic equation (Hasselmann J Fluid Mech 12:481–500, 1962). Numerical simulations with WAVEWATCH have been carried out to thoroughly test the performance of the new form of S_nl, and to compare it to the existing models for S_nl in WAVEWATCH; the DIA and WRT. Some differences between the different models for S_nl are observed. As expected, the DIA is shown to perform less well compared to the exact terms in certain situations, in particular for narrow wave spectra. Also for the case of turning wind significant differences between the different models are observed. Nevertheless, different from the case of unidirectional waves where the generalized kinetic equation represents a obvious improvement to the standard forms of S_nl (Gramstad and Stiassnie 2013), the differences seems to be less pronounced for the more realistic cases considered in this paper.     

Contrasted turbulence intensities in the Indonesian Throughflow: a challenge for parameterizing energy dissipation rate

Microstructure measurements were performed along two sections through the Halmahera Sea and the Ombai Strait and at a station in the deep Banda Sea. Contrasting dissipation rates (??) and vertical eddy diffusivities (K _z ) were obtained with depth-averaged ranges of \(\sim [9 \times 10^{-10}-10^{-5}]\) W kg^??1 and of \(\sim [1 \times 10^{-5}-2 \times 10^{-3}]\) m² s^??1, respectively. Similarly, turbulence intensity, \(I={\epsilon }/(\nu N^{2})\) with ν the kinematic viscosity and N the buoyancy frequency, was found to vary seven orders of magnitude with values up to \(10^{7}\). These large ranges of variations were correlated with the internal tide energy level, which highlights the contrast between regions close and far from internal tide generations. Finescale parameterizations of ?? induced by the breaking of weakly nonlinear internal waves were only relevant in regions located far from any generation area (“far field”), at the deep Banda Sea station. Closer to generation areas, at the “intermediate field” station of the Halmahera Sea, a modified formulation of MacKinnon and Gregg (2005) was validated for moderately turbulent regimes with 100 < I < 1000. Near generation areas marked by strong turbulent regimes such as “near field” stations within strait and passages, ?? is most adequately inferred from horizontal velocities provided that part of the inertial subrange is resolved, according to Kolmogorov scaling.     

Peculiarities of Rayleigh wave displacement in a transversely isotropic half space

We consider a transversely isotropic medium with vertical axis of symmetry (VTI). Rayleigh wave displacement components in a homogeneous VTI medium contain exp(±kr_jz), where z is the vertical coordinate, k is the wave number, and j?=?1, 2; r_j may be considered as depth-decay factor. In a VTI medium, r_j can be a real or imaginary as in an isotropic medium, or it can be a complex depending on the elastic parameters of the VTI medium; if complex, r₁ and r₂ are complex conjugates. In a homogeneous VTI half space, Rayleigh wave displacement is significantly modified with a phase shift when r_j changes from real to complex with variation of VTI parameters; at the transition, the displacement becomes zero when r₁?=?r₂. In a liquid layer over a VTI half space, Rayleigh waves are dispersive. Here, also Rayleigh wave displacement significantly modified with a phase shift when r_j changes from real to complex with a decrease of period. At very low period, phase velocity of Rayleigh waves becomes less than P-wave velocity in the liquid layer giving rise to Scholte waves (interface waves). The amplitudes of Scholte waves with a VTI half space are found to be significantly larger than those with an isotropic half space.     

Changes to sediment sources following wildfire in a forested upland catchment,southeastern Australia

Few investigations link post‐fire changes to sediment sources and erosion processes with sediment yield response at the catchment scale. This linkage is essential if downstream impacts on sediment transport after fire are to be understood in the context of fire effects across different forest environments. In this study, we quantify changing source contributions to fine sediment (<63 µm) exported from a eucalypt forest catchment (136 ha) burnt by wildfire. The study catchment is one of a pair of research catchments located in the East Kiewa River valley in southeastern Australia that have been the subject of a research program investigating wildfire effects on runoff, erosion, and catchment sediment/nutrient exports. This previous research provided the opportunity to couple insights gained from a range of measurement techniques with the application of fallout radionuclides ¹³⁷Cs and ²¹⁰Pb_ex to trace sediment sources. It was found that hillslope surface erosion dominated exports throughout the 3·5‐year post‐fire measurement period. During this time there was a pronounced decline in the proportional surface contribution from close to 100% in the first six months to 58% in the fourth year after fire. Over the study period, hillslope surface sources accounted for 93% of the fine sediment yield from the burnt catchment. The largest decline in the hillslope contribution occurred between the first and second years after fire, which corresponded with the previously reported large decline in sediment yield, breakdown of water repellency in burnt soils, substantial reduction in hillslope erodibility, and rapid surface vegetation recovery. Coupling the information on sediment sources with hillslope process measurements indicated that only a small proportion of slopes contributed sediment to the catchment outlet, with material derived from near‐channel areas dominating the post‐fire catchment sediment yield response. Copyright © 2011 John Wiley & Sons, Ltd.     

Identification of Significant Flood Areas in Lithuania

During preliminary flood risk assessment in Lithuania 54 significant flood areas (SFA) were identified. The detailed flood hazard and risk maps were prepared for these areas in 2014. European Union Floods Directive does not indicate the concrete criteria for SFA delineation. The uncertainty analysis shows that the total length of SFA is not very sensitive to used methodology. In some rivers the uncertainties of 100 year flood peek discharge (Q_1%) were large, but the variation of SFA boundary location was relatively small due to properties of hydrological network. The catchment area and Q_1% change rapidly near the junction with large tributaries, so the boundaries of SFA are usually attached to these junctions. The formal criteria are mostly used to evaluate the possibility of significant floods, but the delineation of SFA is usually based on subjective decision.     

Decaying grid turbulence experiments in a stratified fluid: flow measurements and statistics

Laboratory experiments of decaying grid stratified turbulence were performed in a two-layer fluid and varying the stratification intensity. Turbulence was generated by towing an array of cylinders in a square vessel and the grid was moved at a constant velocity along the total vertical extent of the tank. In order to investigate the influence of the stratification intensity on the turbulence decay, both 2C-PIV and stereo PIV were used to provide time resolved velocity fields in the horizontal plane and the out-of-plane velocity. As expected, a faster decay of the turbulence level along the vertical axis and the collapse in a quasi-horizontal motion increased with the buoyancy frequency, N. In order to characterise the decay process we investigated the time evolution of the vortex statistics, the turbulence scales and the kinetic energy and enstrophy of the horizontal flow. The exponents recovered in the corresponding scaling laws were compared with the theoretical predictions and with reference values obtained in previous experimental studies. Both the spectral analysis and the evolution of characteristic length scales indicate that, in the examined range of N, the dynamics is substantially independent of the stratification intensity. The results obtained were explained in terms of the scaling analysis of decaying turbulence in strongly stratified fluids introduced by Brethouwer et al. (J Fluid Mech 585:343–368.  https://doi.org/10.1017/S0022112007006854, 2007).     

Spatio-temporal streamflow generation in a small,steep headwater catchment in western Japan

Headwaters contribute a substantial part of the flow in river networks. However, spatial variations of streamflow generation processes in steep headwaters have not been well studied. In this study, we examined the spatio-temporal variation of streamflow generation processes in a steep 2.98-ha headwater catchment. The time when baseflow of the upstream section exceeded that downstream was coincident with the time when the riparian groundwater switched from downwelling to upwelling. This suggests that upwelling of the riparian groundwater increased considerably in the upstream section during the wet period, producing a shift in the relative size of baseflow between the upstream and downstream sections. The timing of fluctuations among hillslope soil moisture, hillslope groundwater and streamflow reveals that the hillslope contributed to storm flow, but this contribution was limited to the wet period. Overall, these results suggest that streamflow generation has strong spatial variations, even in small, steep headwater catchments.

EDITOR A. Castellarin ASSOCIATE EDITOR X. Chen     

Effects of temporally correlated infiltration on water flow in an unsaturated–saturated system

Effects of temporally correlated infiltration on water flow in an unsaturated–saturated system were investigated. Both white noise and exponentially correlated infiltration processes were considered. The moment equations of the pressure head (ψ) were solved numerically to obtain the variance and autocorrelation functions of ψ at 14 observation points. Monte Carlo simulations were conducted to verify the numerical results and to estimate the power spectrum of ψ (S _ψψ ). It was found that as the water flows through the system, the variance of the ψ (\( \sigma_{\psi }^{2} \)) were damped by the system: the deeper in the system, the smaller the \( \sigma_{\psi }^{2} \), and the larger the correlation timescale of the infiltration process (λ _I ), the larger the \( \sigma_{\psi }^{2} \). The unsaturated–saturated system gradually filters out the short-term fluctuations of ψ and the damping effect is most significant in the upper part of the system. The fluctuations of ψ is non-stationary at early time and becomes stationary as time progresses: the larger the value of λ _I , the longer the non-stationary period. The correlation timescale of the ψ (λ _ψ ) increases with depth and approaches a constant value at depth: the larger the value of λ _I , the larger the value of λ _ψ . The results of the estimated S _ψψ is consistent with those of the variance and autocorrelation function.