Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

Fingerprinting the 8.2 ka event climate response in a coupled climate model

Using results from coupled climate model simulations of the 8.2 ka climate event that produced a cold period over Greenland in agreement with the reconstructed cooling from ice cores, we investigate the typical pattern of climate anomalies (fingerprint) to provide a framework for the interpretation of global proxy data for the 8.2 ka climate event. For this purpose we developed an analysis method that isolates the forced temperature response and provides information on spatial variations in magnitude, timing and duration that characterise the detectable climate event in proxy archives. Our analysis shows that delays in the temperature response to the freshwater forcing are present, mostly in the order of decades (30 a over central Greenland). The North Atlantic Ocean initially cools in response to the freshwater perturbation, followed in certain parts by a warm response. This delay, occurring more than 200 a after the freshwater pulse, hints at an overshoot in the recovery from the freshwater perturbation. The South Atlantic and the Southern Ocean show a warm response reflecting the bipolar seesaw effect. The duration of the simulated event varies for different areas, and the highest probability of recording the event in proxy archives is in the North Atlantic Ocean area north of 40° N. Our results may facilitate the interpretation of proxy archives recording the 8.2 ka event, as they show that timing and duration cannot be assumed to correspond with the timing and duration of the event as recorded in Greenland ice cores. Copyright © 2011 John Wiley & Sons, Ltd.     

Estimation of daily stream water temperatures with a Bayesian regression approach

Stream water temperature plays a significant role in aquatic ecosystems where it controls many important biological and physical processes. Reliable estimates of water temperature at the daily time step are critical in managing water resources. We developed a parsimonious piecewise Bayesian model for estimating daily stream water temperatures that account for temporal autocorrelation and both linear and nonlinear relationships with air temperature and discharge. The model was tested at 8 climatically different basins of the USA and at 34 sites within the mountainous Boise River Basin (Idaho, USA). The results show that the proposed model is robust with an average root mean square error of 1.25 °C and Nash–Sutcliffe coefficient of 0.92 over a 2‐year period. Our approach can be used to predict historic daily stream water temperatures in any location using observed daily stream temperature and regional air temperature data.     

Watershed structural influences on the distributions of stream network water and solute travel times under baseflow conditions

Watershed structure influences the timing, magnitude, and spatial location of water and solute entry to stream networks. In turn, stream reach transport velocities and stream network geometry (travel distances) further influence the timing of export from watersheds. Here, we examine how watershed and stream network organization can affect travel times of water from delivery to the stream network to arrival at the watershed outlet. We analysed watershed structure and network geometry and quantified the relationship between stream discharge and solute velocity across six study watersheds (11.4 to 62.8 km²) located in the Sawtooth Mountains of central Idaho, USA. Based on these analyses, we developed stream network travel time functions for each watershed. We found that watershed structure, stream network geometry, and the variable magnitude of inputs across the network can have a pronounced affect on water travel distances and velocities within a stream network. Accordingly, a sample taken at the watershed outlet is composed of water and solutes sourced from across the watershed that experienced a range of travel times in the stream network. We suggest that understanding and quantifying stream network travel time distributions are valuable for deconvolving signals observed at watershed outlets into their spatial and temporal sources, and separating terrestrial and in‐channel hydrological, biogeochemical, and ecological influences on in‐stream observations. Copyright © 2016 John Wiley & Sons, Ltd.     

A GLUE‐based uncertainty assessment framework for tritium‐inferred transit time estimations under baseflow conditions

The last decade has seen major technical and scientific improvements in the study of water transfer time through catchments. Nevertheless, it has been argued that most of these developments used conservative tracers that may disregard the oldest component of water transfer, which often has transit times greater than 5 years. Indeed, although the analytical reproducibility of tracers limits the detection of the older flow components associated with the most dampened seasonal fluctuations, this is very rarely taken into account in modelling applications. Tritium is the only environmental tracer at hand to investigate transfer times in the 5‐ to 50‐year range in surface waters, as dissolved gases are not suitable due to the degassing process. Water dating with tritium has often been difficult because of the complex history of its atmospheric concentration, but its current stabilization together with recent analytical improvements open promising perspectives. In this context, the innovative contribution of this study lies in the development of a generalized likelihood uncertainty estimation‐based approach for analysing the uncertainties associated with the modelling of transit time due to both parameter identification and tracer analytical precision issues. A coupled resampling procedure allows assessment of the statistical significance of the transfer time differences found in diverse waters. This approach was developed for tritium and the exponential‐piston model but can be implemented for virtually any tracer and model. Stream baseflow, spring and shallow aquifer waters from the Vallcebre research catchments, analysed for tritium in different years with different analytical precisions, were investigated by using this approach and taking into account other sources of uncertainty. The results showed three groups of waters of different mean transit times, with all the stream baseflow and spring waters older than the 5‐year threshold needing tritium. Low sensitivity of the results to the model structure was also demonstrated. Dual solutions were found for the waters sampled in 2013, but these results may be disambiguated when additional analyses will be made in a few years. Copyright © 2016 John Wiley & Sons, Ltd.     

联合GPS和GRACE数据分析陕甘宁地区地表垂直形变

采用甘肃省CORS网和中国大陆构造环境监测网络中共48个台站的GPS观测数据，解算得到观测台站的垂直位移，并与GRACE时变重力场Mascon模型解CSR RL05M数据计算得到的垂直形变进行比较，分析区域地表垂直形变特征。结果表明，研究区内台站垂直形变存在局部特征，甘肃庆阳和平凉地区垂直形变与其他地区存在明显差异，相关系数、均方根减少量和周年信号减少量均高于其他地区；扣除趋势项后，观测台站GPS垂直位移与GRACE垂直形变时间序列相关系数均值为0.72，GPS和GRACE周年信号振幅均值分别为6.00 mm和3.70 mm，周年信号减少量和均方根误差减少量均值分别为0.51和0.29；研究区内GPS垂直位移和GRACE垂直形变时间序列一致性较强，GRACE垂直形变能有效解释50%以上的GPS垂直位移周年信号，GPS垂直位移时间序列包含的非构造形变中平均约29%来源于环境负载变化所引起的负荷形变。     

Spatial aggregation of time‐variant stream water ages in urbanizing catchments

We calibrated an integrated flow–tracer model to simulate spatially distributed isotope time series in stream water in a 7.9‐km² catchment with an urban area of 13%. The model used flux tracking to estimate the time‐varying age of stream water at the outlet and both urbanized (1.7 km²) and non‐urban (4.5 km²) sub‐catchments over a 2.5‐year period. This included extended wet and dry spells where precipitation equated to >10‐year return periods. Modelling indicated that stream water draining the most urbanized tributary was youngest with a mean transit time (MTT) of 171 days compared with 456 days in the non‐urban tributary. For the larger catchment, the MTT was 280 days. Here, the response of urban contributing areas dominated smaller and more moderate runoff events, but rural contributions dominated during the wettest periods, giving a bi‐modal distribution of water ages. Whilst the approach needs refining for sub‐daily time steps, it provides a basis for projecting the effects of urbanization on stream water transit times and their spatial aggregation. This offers a novel approach for understanding the cumulative impacts of urbanization on stream water quantity and quality, which can contribute to more sustainable management. Copyright © 2015 John Wiley & Sons, Ltd.     

Analysis of the interannual variability of annual daily extreme water levels in the St Lawrence River and Lake Ontario from 1918 to 2010

We compared the interannual variability of annual daily maximum and minimum extreme water levels in Lake Ontario and the St Lawrence River (Sorel station) from 1918 to 2010, using several statistical tests. The interannual variability of annual daily maximum extreme water levels in Lake Ontario is characterized by a positive long‐term trend showing two shifts in mean (1929–1930 and 1942–1943) and a single shift in variance (in 1958–1959). In contrast, for the St Lawrence River, this interannual variability is characterized by a negative long‐term trend with a single shift in mean, which occurred in 1955–1956. As for annual daily minimum extreme water levels, their interannual variability shows no significant long‐term change in trend. However, for Lake Ontario, the interannual variability of these water levels shows two shifts in mean, which are synchronous with those for maximum water levels, and a single shift in variance, which occurred in 1965–1966. These changes in trend and stationarity (mean and variance) are thought to be due to factors both climatic (the Great Drought of the 1930s) and human (digging of the Seaway and construction of several dams and locks during the 1950s). Despite this change in means and variance, the four series are clearly described by the generalized extreme value distribution. Finally, annual daily maximum and minimum extreme water levels in the St Lawrence and Lake Ontario are negatively correlated with Atlantic multidecadal oscillation over the period from 1918 to 2010. Copyright © 2013 John Wiley & Sons, Ltd.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.