Non-Darcian flow to a well in an aquifer–aquitard system

In this study, we use a linearization procedure and a finite difference method to solve non-Darcian flow to a well in an aquifer–aquitard system. The leakage effect is considered. Flow in the aquifer is assumed to be non-Darcian and horizontal, whereas flow in the aquitard is assumed to be Darcian and vertical. The Izbash equation [Izbash SV. O filtracii V Kropnozernstom Materiale. USSR: Leningrad; 1931 [in Russian]] is employed to describe the non-Darcian flow. The wellbore storage is also considered in this study. An approximate semi-analytical solution has been obtained by the linearization procedure, and a numerical solution has been obtained by using a finite difference method. The previous solutions for Darcian flow case and non-Darcian flow case without leakage can be described as special cases of the new solutions. The error caused by the linearization procedure has also been analyzed. The relative error caused by the linearization procedure is nearly 100% at early times, and decreases to zero at late times. We have also compared the results in this study with Wen et al. [Wen Z, Huang G, Zhan H. A numerical solution for non-Darcian flow to a well in a confined aquifer using the power law function. J Hydrol, 2008d [in revision]] in which the leakage effect is not considered, and Hantush and Jacob [Hantush MS, Jacob CE. Non-steady radial flow in an infinite leaky aquifer. Trans Am Geophys Union 1955;36(1):95–100] who investigated a similar problem in Darcian flow case. The comparison of this study and Wen et al. (2008d) indicates the dimensionless drawdown in the aquifer with leakage is less than that without leakage, and the leakage has little effect at early times. The comparison between the results of this study and that of Hantush and Jacob (1955) indicates that the dimensionless drawdown in the aquifer for non-Darcian flow is larger at early times and smaller at late times, than their counterparts for Darcian flow. A larger dimensionless non-Darcian conductivity k_D results in a smaller dimensionless drawdown in the aquifer at late times, and leads to a larger dimensionless drawdown in the aquifer at early times. A smaller dimensionless leakage parameter B_D results in a smaller drawdown at late times, and the leakage does not affect the early-time drawdown. The analysis of the dimensionless drawdown inside the well has also been included in this study when the wellbore storage is considered.     

Estimating the role of a ditch network in groundwater recharge in a Mediterranean catchment using a water balance approach

Water balance variables were monitored in a farmed Mediterranean catchment characterized by a dense ditch network to allow for the separate estimation of the diffuse and concentrated recharge terms during flood events. The 27 ha central part of the catchment was equipped with (i) rain gauges, (ii) ditch gauge stations, (iii) piezometers, (iv) neutron probes, and (v) an eddy covariance mast including a 3D sonic anemometer and a fast hygrometer. The water balance was calculated for two autumnal rain and flood events. We also estimated the uncertainty of this approach with Monte Carlo simulations. Results show, that although ditch area represents only 6% of the total study area, concentrated recharge appeared to be the main source of groundwater recharge. Indeed, it was 40–50% of the total groundwater recharge for autumnal events, which are the major annual recharge events. This indicate that both, concentrated and diffuse recharge should be taken into account in any hydrological modeling approach for Mediterranean catchments. This also means that, since they collect overland flow that is often largely contaminated by chemicals, ditches may be a place where groundwater contamination is likely to occur. The uncertainty analysis indicates that recharge estimates based on water balance exhibit large uncertainty ranges. Nevertheless, Monte Carlo simulations showed that concentrated recharge was higher than expected based on their area.     

Improving hydrodynamic modeling of river networks by incorporating data assimilation using a particle filter

Numerical modeling is a well-recognized method for studying the hydrodynamic processes in river networks. Multi-source measurements also offer abundant information on the patterns and mechanisms within the processes. Therefore, improving hydrodynamic modeling of river networks through the use of data assimilation techniques has become a hot research topic in recent years. The particle filter (PF) is a commonly used data assimilation method and has been proven to be applicable to various nonlinear and non-Gaussian models. In the current study, an improved numerical hydrodynamic model for large-scale river networks is established by incorporating the advanced PF algorithm. Furthermore, the PF method based on the Gaussian likelihood function (GLF) and the method based on the Cauchy likelihood function (CLF) are compared for a complex river network scenario. The feasibility of the PF-based methods was evaluated through application to the Yangtze-Dongting River-lake Network (YDRN) by assimilating water stage data collected at six hydrometric stations during the entire hydrodynamic process in 2003. Additionally, the parameters used in the likelihood function, which affect the assimilation performance, also were explored in the current study. The study results found that the accuracy of the model-derived water stage data was improved when the PF-based methods are utilized, with improvement not only at the data assimilation (calibration) sites but also at three hydrometric stations not used in the data assimilation (i.e., verification sites). The highest average Nash-Sutcliffe Efficiency result for the six assimilation sites were 0.98 while the lowest summed root-mean-square-error result was 1.801 m. The comparison results also indicated that the CLF-based PF outperformed the GLF-based PF when high-accuracy observed data are available. Specifically, the CLF can effectively resolve the filtering failure problem and the dispersion problem of PFs, and further improve the accuracy of the filtering results for a river network scenario. In summary, the CLF-based PF method along with high-accuracy observation data shows promise to provide reliable reference and technical support for hydrodynamic modeling of large-scale river networks.     

The diffuse-to-global spectral irradiance ratio as a cloud-screening technique for radiometric data

Continuous measurements of solar spectral radiation using the multifilter rotating shadow band radiometer (MFRSR-7) are performed at the Actinometric Station of the National Observatory of Athens (ASNOA). The present study utilizes 4 days of continuous observations, from local sunrise to sunset, in order to investigate the daily variation of the radiation components (diffuse and global) as well as their ratio (diffuse-to-global irradiance ratio, DGR) under different atmospheric conditions. DGR has received a great scientific interest, as well as the respective diffuse-to-direct-beam ratio especially for investigating solar irradiance modifications under different atmospheric conditions, aerosol load and optical properties. Apart from this, the present study shows that the DGR can also constitute a powerful tool for cloud screening, i.e. for removing perturbed data due to cloud contamination from automated sun scanning radiometers. The relationship between DGR at a specific wavelength with the respective ratio for the whole MFRSR band (300–1100 nm) is found to exhibit a curvature; this curvature is strongly modified when perturbed irradiance data (possibly caused by clouds) occur. Even though the perturbed data can also be easily identified from the diurnal irradiance variation, the present study is the first to show the effect of perturbed solar spectral data on the DGR.     

Static solution of a crack degenerated from dynamic solution of a propagating crack

ＳｔａｔｉｃｓｏｌｕｔｉｏｎｏｆａｃｒａｃｋｄｅｇｅｎｅｒａｔｅｄｆｒｏｍｄｙｎａｍｉｃｓｏｌｕｔｉｏｎｏｆａｐｒｏｐａｇａｔｉｎｇｃｒａｃｋＳｈｉ－ＹｕＬＩ（李世愚）ａｎｄＹｕｎ－ＴａｉＣＨＥＮ（陈运泰）（ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｐｈｙｓｉｃｓ...     

Application of a multi‐temporal,LiDAR‐derived,digital terrain model in a landslide‐volume estimation

Sediments produced by landslides are crucial in the sediment yield of a catchment, debris flow forecasting, and related hazard assessment. On a regional scale, however, it is difficult and time consuming to measure the volumes of such sediment. This paper uses a LiDAR‐derived digital terrain model (DTM) taken in 2005 and 2010 (at 2 m resolution) to accurately obtain landslide‐induced sediment volumes that resulted from a single catastrophic typhoon event in a heavily forested mountainous area of Taiwan. The landslides induced by Typhoon Morakot are mapped by comparison of 25 cm resolution aerial photographs taken before and after the typhoon in an 83.6 km² study area. Each landslide volume is calculated by subtraction of the 2005 DTM from the 2010 DTM, and the scaling relationship between landslide area and its volume are further regressed. The relationship between volume and area are also determined for all the disturbed areas (V_L = 0.452A_L^1.242) and for the crown areas of the landslides (V_L = 2.510A_L^1.206). The uncertainty in estimated volume caused by use of the LiDAR DTMs is discussed, and the error in absolute volume estimation for landslides with an area >10⁵ m² is within 20%. The volume–area relationship obtained in this study is also validated in 11 small to medium‐sized catchments located outside the study area, and there is good agreement between the calculation from DTMs and the regression formula. By comparison of debris volumes estimated in this study with previous work, it is found that a wider volume variation exists that is directly proportional to the landslide area, especially under a higher scaling exponent. Copyright © 2013 John Wiley & Sons, Ltd.     

Scouring of a granular bed by dam break: Experimental study and numerical simulation by a VOF-LPT coupling

The first part of the research reported here consists of an experimental campaign to study the scouring of a granular bed(glass beads, sand) induced by a dam break in an open channel. Two configurations are considered: with and without cylinders. In the second part of this study, the volume of fluid method coupled with the shear stress transport turbulent model and the lagrangian particle tracking method is used to simulate the local scour processes. The four-way coupling is realized by consider...     

Multisource hydrologic modeling uncertainty analysis using the IBUNE framework in a humid catchment

Hydrologic cycle is a complex system associated with both certain and uncertain constituents. The propagation of confidence bounds from different uncertainty sources to model output is of great significance for hydrologic modeling. In this paper, we applied the integrated bayesian uncertainty estimator to quantify the effects of parameter, input and model structure uncertainty on hydrologic modeling progressively. Two hydrologic models (Xinanjiang model and TOPMODEL) were applied to a humid catchment under three scenarios. Case I: the shuffled complex evolution metropolis (SCEM-UA) algorithm was conducted to determine the posterior parameter distribution of hydrologic models and analyze the corresponding forecast uncertainty. Case II: input uncertainty was also considered by assuming rain depth bias follows a normal distribution, and integrated with SCEM-UA. Case III: Simulations from two models were combined by the Bayesian model averaging to fully quantify multisource uncertainty effects. Results suggested that, from Case I to II, the containing ratio (percentage of observed streamflow enveloped by 95% confidence interval) obviously increased by an average magnitude of 10% for the study period 2000–2006. Besides, it also found that the width of 95% confidence interval became wider and narrower for Xinanjiang model and TOPMODEL, respectively, from Case I to II. This may indicate that the uncertainty of TOPMODEL results was more remarkable than Xinanjiang model in Case I. By combining results from two models, model structure uncertainty was also considered in Case III. The accuracy of uncertainty bounds further improved with the containing ratio of 95% confidence interval >95%. In addition, the optimized deterministic results from the uncertainty analysis showed that the average Nash–Sutcliffe coefficient increased continually from Case I to II and III (0.82, 0.84 and 0.90, respectively) for the study period. The analysis demonstrated the improvement of modeling accuracy when extra uncertainty sources were also quantified, and this finding also proved the applicability of IBUNE framework in hydrologic modeling.     

Assessment of change in design flood frequency under climate change using a multivariate downscaling model and a precipitation-runoff model

Precipitation and runoff are key elements in the hydrologic cycle because of their important roles in water supply, flood prevention, river restoration, and ecosystem management. Global climate change, widely accepted to be happening, is anticipated to have enormous consequences on future hydrologic patterns. Studies on the potential changes in global, regional, and local hydrologic patterns under global climate change scenarios have been an intense area of research in recent years. The present study contributes to this research topic through evaluation of design flood under climate change. The study utilizes a weather state-based, stochastic multivariate model as a conditional probability model for simulating the precipitation field. An important premise of this study is that large-scale climatic patterns serve as a major driver of persistent year-to-year changes in precipitation probabilities. Since uncertainty estimation in the study of climate change is needed to examine the reliability of the outcomes, this study also applies a Bayesian Markov chain Monte Carlo scheme to the widely used SAC-SMA (Sacramento soil moisture accounting) precipitation-runoff model. A case study is also performed with the Soyang Dam watershed in South Korea as the study basin. Finally, a comprehensive discussion on design flood under climate change is made.     

Note on the flow of a stratified fluid over a stationary obstacle in a channel

Abstract

A new method is introduced to produce a uniform stratified flow over a stationary obstacle in an open channel. The flow is achieved by discharging the flow from the channel through a sink. The details of the sink are unimportant. The flow speed is limited only by the sink capacity. Selective withdrawal at lower densimetric Froude numbers is effectively eliminated through the use of a contraction. The standing, free-surface, long wave arising from the initiation of the flow is also eliminated by the contraction. Experiments are conducted for flow over a sphere for a range of Reynolds numbers from O(10²) to O(10³) and a range of Richardson numbers from O(10^?1) to O(10). Dye and neutrally buoyant droplets are used for quantitative analysis of the wake structure. The wake is also probed by a hot-film anemometer. The frequency of vortex shedding is obtained. Comparison with data from towed experiments is also presented.     

Reduction of local scour at river conf luences using a collar

River confluences (RCs) are important features within river systems where the three dimensional (3D) flow structures and the downstream mixing of flows can cause deep scour holes. Despite this, few methods have been proposed to control scouring at RCs. In this study, application of a collar was experimentally examined for local scour control at the point where two rivers flow together. In parti-cular, experimental tests were done with and without collar application at three different locations. The results reveal that the scour depth is directly proportional to the discharge ratio, i.e. the ratio of lateral discharge to that in the channel downstream of the confluence, and the densimetric Froude number (Frg). In addition, installation of a collar at RCs can decrease the scour depth up to 100%, thus completely avoiding the scour process. The results also show that by increasing the Frg the optimal installation location for a collar changes and moves towards the river bed level. Using a collar can also reduce the height of the point bar formed downstream of the confluence. The outcomes of the study allow deri-vation of an equation for predicting scour depth when a collar is applied as a countermeasure. The analysis of this equation shows that the estimates are mostly affected by the Frg.     

Experimental study on lag time for a small watershed

This paper is focused on the lag time of flow confluence. Lag is described empirically by the precipitation factor in hydrological modelling, but the traditional way to establish the relationship between lag and intensity is not very satisfactory in arid and semi‐arid regions. A total 215 rainfall–runoff experiments were conducted to reveal the effects of net rain duration on lag and intensity. The results show that a correlation between lag and rain intensity exists under certain conditions. There is a critical value in net rain intensity (0·8–1 mm min^?1) in the correlation curve of lag and net rain intensity in the given experimental watershed in the laboratory. Similarly, a critical value of net rain duration also exists. This value is the total confluence time. The features of lag time change significantly if intensity or duration is less or greater than the critical value. The study also explores the joint effects of net rain intensity and duration on lag. The formula established among lag, intensity and duration resulted in a better fit. Therefore, the two‐parameter (intensity and duration) empirical formula for lag is better than the traditional single‐parameter (intensity) method. This two‐parameter correlation formula can also be applied to a temporally and spatially uneven runoff processes. A typical field watershed is selected to test the results of the experiments. Copyright © 2006 John Wiley & Sons, Ltd.     

A framework approach to address the trend and causes of flood stage change in a river reach downstream of a dam influenced by tributaries

The evaluation of the trend of flood stage changes in alluvial rivers downstream of dams is important for flood management. However, the flood stage associated with a given discharge generally is nonstationary in river reaches with multiple tributaries. This is not only because of the dam-induced shifting in the cross-sectional area and/or channel roughness but also because of the backwater induced by high flows from the tributaries. To determine the total trend of the flood stage and quantify the separate contributions of hydrological and geomorphic effects, the current study proposed a framework approach consisting of hydrological analysis and multiscenario numerical modeling. By this means, the trend in the flood stage could be distinguished from the stage oscillation driven by varying factors, including extreme hydrologic events. The effects of chronic changes, including channel incision and flow resistance increase, also were quantitatively separated. This framework was applied to the Chenglingji–Datong (CD) reach downstream of the Three Gorges Dam (TGD) in the Yangtze River, China. The results indicated that the effect of the roughness increase counterbalanced the effect of channel incision when the flow discharge was beyond the bankfull level. The backwater effect induced by tributary inflow was the major cause of the flood stage rise in recent years. The method presented in the current study provides a useful tool for managers and engineers to obtain better insight into the driving mechanisms of flood stage changes in river reaches that are downstream of dams. These findings indicate that the flood stage may not decline or may even occasionally increase, although the cross-sectional area was enlarged by channel incision. Special attention should be given to the flood risk situation in the study reach after the TGD began operation.     

The effects of a copper-based antifouling paint on mortality and enzymatic activity of a non-target marine organism

Antifouling paints are used on a wide range of underwater structures in order to protect them from the development of fouling organisms. The leaching of the toxic substances from the matrix of the paint causes toxic effects not only to the fouling organisms but also on other "non-target" biota. The present study addresses the impact of the antifouling paint Flexgard VI-II on brine shrimp nauplii selected as convenient test organisms. The surface to volume (S/V) concept developed by Persoone and Castritsi-Catharios (1989) was used to determine S/V-LC50s for the test biota exposed to PVC test panels of 400-1000 mm2 surface coated with the antifouling paint in test vessels containing 20 ml seawater. Total ATPase and Mg2+-ATPase were also analyzed for coated surface areas inducing less than 50% mortality in the brine shrimp nauplii. The calculated S/V-LC50 (24 h) was 24.6 mm2/ml, which shows the high toxic character of the antifouling paint. Decreased enzymatic activities were noted in the brine shrimp nauplii exposed to test panels of 50 and 100 mm2 in 20 ml seawater. The present study indicates that the "surface to volume" concept is an interesting methodology that can be applied with both lethal and sublethal effect criteria for the determination of toxic stress from leaches of painted surfaces.     

Inferring catchment flow path responses using a data-driven model: an exploratory study based on a generalized additive model

Identifying physical catchment processes from streamflow data, such as quick- and slow-flow paths, remains challenging. This study is designed to explore whether a flexible nonparametric regression model (generalized additive model, GAM) can be used to infer different flow paths. This assumes that the data relationship in data-driven models is also a reflection of catchment physical processes. The GAM, using time-lagged flow covariates, was fitted to synthetic rainfall–runoff data simulated using simple linear reservoirs. Partial plots of the time-lagged covariates show that the model could differentiate simple and more complex flow paths in simulated synthetic data with short and long memory systems and varying between dry and wet climates. Further analysis of data from real catchments showed that the model could differentiate catchments dominated by slow flow and by quick flow. Therefore, this study indicates that GAM can be used to identify catchment storages and delay processes from streamflow data.     

Hydrodynamic pressure on arch dams—by a mapping finite element method

A novel Finite Element Method (FEM) is proposed for the analysis of the uncoupled hydrodynamic pressures generated on arch dams due to a steady-state ground acceleration. In this method the equation governing hydrodynamic pressures and also the prescribed boundary conditions are all transformed from the Cartesian space to a logarithmically condensed cylindrical polar space; in this process the physical configuration of the reservoir-dam is also mapped into an ‘image’ domain. The transformed governing equation is then solved in the image domain, subject to the transformed boundary conditions, using standard finite elements. Because physical dimensions are logarithmically condensed in the image space, the proposed method is particularly suitable for dealing with large or very large aspect-ratio reservoir-dam systems, economically and efficiently. The high degree of accuracy which the proposed method is capable of, and also the simple way in which it can be applied to complex reservoir-dam shapes, have been demonstrated by means of examples. The method has also been applied to study the uncoupled hydrodynamic pressures on the upstream face of a cylindrical arch dam, generated by a steady-state vertical ground acceleration.     

Estuarine shoreline processes in a dynamic low-energy system

Estuarine shorelines are often classified as low-energy coasts and are, therefore, expected to undergo little variation. Port Stephens (SE Australia) is a ria-like drowned river valley microtidal estuary located on a wave-dominated coast. The outer part of the estuary is tide-dominated and has a large shallow flood-tide delta, which is also affected by waves. The northern (predominantly low-energy) shoreline of outer Port Stephens is a continuous stretch of sand comprising areas of high mobility and areas of relative stability terminating in a western extending sand spit. This paper investigates the effects of periodic high-energy conditions during which waves penetrate into the estuary by analysing two types of storms, low to moderate (more frequent type) storms and severe to extreme (low frequency) storms. It is established that low to moderate storms cause generalised erosion over the northern shoreline. On the contrary, severe to extreme storms, while causing erosion on parts of the beach, can transport new sediment across the flood-tide delta and deposit it to build a mobile shore attached sandwave. Long-term (decadal) trends identified in the study area are in agreement with short- and medium-term results. Moving into the estuary are four complementary zones of sediment transport which include: (1) sandwave formation and westward migration; (2) a relatively stable area between the sandwave and an erosion zone; (3) an erosion zone undergoing shoreline retreat and finally (4) a depositional terminus causing westward extension of the sand spit.     

含起伏地形地层中点源激发的震电响应研究

地震诱导电磁现象是国内外地学领域十分关注的前沿问题,前人对地震波和电磁场耦合波场的认识主要是基于规则模型获得的.为研究含起伏地表和地下界面的地层中震电波场激发、传播特性,本文采用有限元软件COMSOL Multiphysics模拟点震源激发的电磁场.首先给出频率域二维SHTE模式震电耦合方程组,然后利用COMSOL软件建立计算模型,并求解出点力源激发震电波场的频率域响应,最后利用FFT变换得到地震波场和电磁场的时间域波形.模拟结果表明,震电波场中存在三种类型的电磁信号,第一种是震源直接激发的电磁波;第二种是地震波在分界面处激发的电磁波(包括自由表面、地下不同介质分界面);第三种是伴随地震波的同震信号,前两种电磁波比地震波更早到达远处观测台站,对地震预警有重要意义.此外,研究还发现:当地震波传播至地表并沿着地表传播时,在地表附近空气层中同样记录到了伴随地震波传播的电磁扰动信号,该信号与相同水平源距条件下、地下观测点接收到的电磁信号相同,这与前人的一些观测结果相符.本文研究结果为今后地震电磁信号的解释提供了理论证据.     

Extension of observed flood series by combining a distributed hydro-meteorological model and a copula-based model

Long flood series are required to accurately estimate flood quantiles associated with high return periods, in order to design and assess the risk in hydraulic structures such as dams. However, observed flood series are commonly short. Flood series can be extended through hydro-meteorological modelling, yet the computational effort can be very demanding in case of a distributed model with a short time step is considered to obtain an accurate flood hydrograph characterisation. Statistical models can also be used, where the copula approach is spreading for performing multivariate flood frequency analyses. Nevertheless, the selection of the copula to characterise the dependence structure of short data series involves a large uncertainty. In the present study, a methodology to extend flood series by combining both approaches is introduced. First, the minimum number of flood hydrographs required to be simulated by a spatially distributed hydro-meteorological model is identified in terms of the uncertainty of quantile estimates obtained by both copula and marginal distributions. Second, a large synthetic sample is generated by a bivariate copula-based model, reducing the computation time required by the hydro-meteorological model. The hydro-meteorological modelling chain consists of the RainSim stochastic rainfall generator and the Real-time Interactive Basin Simulator (RIBS) rainfall-runoff model. The proposed procedure is applied to a case study in Spain. As a result, a large synthetic sample of peak-volume pairs is stochastically generated, keeping the statistical properties of the simulated series generated by the hydro-meteorological model. This method reduces the computation time consumed. The extended sample, consisting of the joint simulated and synthetic sample, can be used for improving flood risk assessment studies.     

Estuarine Habitat Assessment for Construction of a Submarine Transmission Line

The present paper describes a submarine survey using the acoustic discrimination system QTC VIEW (Series V) as an exploratory tool to adjust final route alignment of a new pipeline. By using acoustic sound survey as an exploratory tool described in this paper to adjust final route alignment of a new pipeline to minimize the environmental impact caused and ultimately to avoid any mitigation measures. The transmission pipeline extended from the shore line of Abu-Qir Bay, on the Mediterranean Sea in Egypt, out to 70 nautical miles at sea (60 m water depth). Four main surface sediment types were defined in the study area, namely fine sand, silty sand, silt and clay. Results of the acoustic classification revealed four acoustic classes. The first acoustic class corresponded to fine sand, absence of shell debris and very poor habitats characteristics. The second acoustic class is predominant in the study area and corresponds to the region occupied by silt. It is also characterized by intermediate diversity of macrobenthic invertebrate community which is mainly characterized by polychaeta. The third acoustic class is characterized by silt to silty clay. It is characterized by a high diversity of macrobenthic invertebrate community which is mainly polychaeta with an intermediate diversity of gastropoda and bivalvia. The final acoustic class is characterized by clay and high occurrence of shell debris of gastropoda, bivalvia and polychaeta.