The impact of highly permeable layer on hydraulic system in a coastal aquifer

This study is aimed to understand the hydraulic mechanism of coastal aquifer systems that include highly permeable layers (HPLs). These hydrologic conditions can be found in many volcanic islands that are composed of a series of lava flows discharged into sea or other standing body of water. In the first part, we developed a numerical model based on the geologic and hydrologic data obtained from the eastern Jeju Island, Korea, of which the aquifer contains clinker and hyaloclastite layers. The simulation results reproduced spatial location of fresh‐saline water interface, especially the abrupt decline of interface at the inland part and the thickness variation of transition zone along the cross‐section observed at the eastern Jeju coastal aquifer. We were able to find out that these phenomena are strongly related to the presence of the HPL. In the second part, quantitative analyses were conducted with the use of hypothetical models in order to understand the dynamic characteristics of coastal system that includes HPLs. A series of sensitivity studies were conducted to assess the effect of the horizontal length and vertical depth of HPL on the spatial location of the interface toe and the configuration of transition zone. Various case studies have shown that the seawater intruded into the inland more as the horizontal length of HPL was increased and its vertical depth was decreased. In other simulations including two HPLs, the vertical distance between these two HPLs primarily controlled the flow regime, flux variations, and the configuration of the transition zone. Finally, we performed simulations to evaluate the effect of a rising sea‐level. This study provides more understanding of how the presence of HPL controls the seawater intrusion processes, and the spatial configurations of fresh‐saline water interface at coastal aquifers. Copyright © 2012 John Wiley & Sons, Ltd.     

On seasonal approach to flood frequency modelling. Part II: flood frequency analysis of Polish rivers

The annual peak flow series of Polish rivers are mixtures of summer and winter flows. As Part II of a sequence of two papers, practical aspects of applicability of seasonal approach to flood frequency analysis (FFA) of Polish rivers are discussed. Taking A Two‐Component Extreme Value (TCEV1) model as an example it was shown in the first part that regardless of estimation method, the seasonal approach can give profit in terms of upper quantile estimation accuracy that rises with the return period of the quantile and is the greatest for no seasonal variation. In this part, an assessment of annual maxima (AM) versus seasonal maxima (SM) approach to FFA was carried out with respect to seasonal and annual peak flow series of 38 Polish gauging stations. First, the assumption of mutual independence of the seasonal maxima has been tested. The smoothness of SM and AM empirical probability distribution functions was analysed and compared. The TCEV1 model with seasonally estimated parameters was found to be not appropriate for most Polish data as it considerably underrates the skewness of AM distributions and upper quantile values as well. Consequently, the discrepancies between the SM and AM estimates of TCEV1 are observed. Taking SM and TCEV1 distribution, the dominating season in AM series was confronted with predominant season for extreme floods. The key argument for presumptive superiority of SM approach that SM samples are more statistically homogeneous than AM samples has not been confirmed by the data. An analysis of fitness to SM and AM of Polish datasets made for seven distributions pointed to Pearson (3) distribution as the best for AM and Summer Maxima, whereas it was impossible to select a single best model for winter samples. In the multi‐model approach to FFA, the tree functions, i.e., Pe(3), CD3 and LN3, should be involved for both SM and AM. As the case study, Warsaw gauge on the Vistula River was selected. While most of AM elements are here from winter season, the prevailing majority of extreme annual floods are the summer maxima. The upper quantile estimates got by means of classical annual and two‐season methods happen to be fairly close; what's more they are nearly equal to the quantiles calculated just for the season of dominating extreme floods. Copyright © 2011 John Wiley & Sons, Ltd.     

Response of a tunnel deeply embedded in a viscoelastic medium

This paper presents a three‐dimensional energy‐based solution for the time‐dependent response of a deeply embedded and unsupported semi‐infinite tunnel of circular cross‐section. The tunnel is taken to be excavated quasi‐instantaneously from an infinite rock body that initially exhibits an isotropic stress state and that is made up of a homogeneous, isotropic and viscoelastic material. The viscoelastic behaviour is modelled by means of Burger's model, and the rock is taken to behave volumetrically linear elastic and to exhibit exclusively deviatoric creep. This viscoelastic problem is transformed into the Laplace domain, where it represents a quasi‐elastic problem. The displacement fields in the new solution are taken to be the products of independent functions that vary in the radial and longitudinal directions. The differential equations governing the displacements of the system and appropriate boundary conditions are obtained using the principle of minimum potential energy. The solutions for these governing equations in the Laplace domain are then obtained analytically and numerically using a one‐dimensional finite difference technique. The results are then transformed back into the time domain using an efficient numerical scheme. The accuracy of the new solution is comparable with that of a finite element analysis but requires much less computation effort. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatial interpolation of snow water equivalency using surface observations and remotely sensed images of snow‐covered area

As demand for water continues to escalate in the western Unites States, so does the need for accurate monitoring of the snowpack in mountainous areas. In this study, we describe a simple methodology for generating gridded‐estimates of snow water equivalency (SWE) using both surface observations of SWE and remotely sensed estimates of snow‐covered area (SCA). Multiple regression was used to quantify the relationship between physiographic variables (elevation, slope, aspect, clear‐sky solar radiation, etc.) and SWE as measured at a number of sites in a mountainous basin in south‐central Idaho (Big Wood River Basin). The elevation of the snowline, obtained from the SCA estimates, was used to constrain the predicted SWE values. The results from the analysis are encouraging and compare well to those found in previous studies, which often utilized more sophisticated spatial interpolation techniques. Cross‐validation results indicate that the spatial interpolation method produces accurate SWE estimates [mean R² = 0·82, mean mean absolute error (MAE) = 4·34 cm, mean root mean squared error (RMSE) = 5·29 cm]. The basin examined in this study is typical of many mid‐elevation mountainous basins throughout the western United States, in terms of the distribution of topographic variables, as well as the number and characteristics of sites at which the necessary ground data are available. Thus, there is high potential for this methodology to be successfully applied to other mountainous basins. Copyright © 2010 John Wiley & Sons, Ltd.     

Variation and pattern in the responses of mammal faunas to Late Pleistocene climatic change in southeastern South Australia

We examined mammal occurrence and variability through the Late Pleistocene vertebrate fossil deposit of Grant Hall in Victoria Fossil Cave, Naracoorte, South Australia. To determine long‐term patterns of change, we compared the composition and relative abundance trends of the assemblage with a nearby Middle Pleistocene deposit in Cathedral Cave. Total species richness did not change through the Grant Hall sequence, dated from 93 ± 8 to 70 ± 5 ka. However, species relative abundances varied between ecologically divergent species, and in some cases between species that demonstrate similar environmental preferences. For some species this variation is comparable to that recorded in Cathedral Cave. Of those showing similar trends between the two deposits, the forest inhabitant, Pseudomys fumeus, recorded an 8.6% decline through Grant Hall, coincident with a 9.7% increase in the dry heath/mallee dweller Pseudomys apodemoides. These patterns indicate that climatic transition from relatively warm, moist to cooler, drier conditions impacted some species in similar ways through climatic cycles of the past. However, the majority of the fauna demonstrated complex responses that are individual and variable through time. Statistical tests of species trends from the Grant Hall assemblage caution that large fossil samples are required to validate patterns observed. Copyright © 2012 John Wiley & Sons, Ltd.     

A steepness index unit volume flood hydrograph approach for sub‐daily flow disaggregation

The estimation of sub‐daily flows from daily flood flows is important for many hydrological and hydraulic applications. Flows during flood events often vary significantly within sub‐daily time‐scales, and failure to capture the sub‐daily flood characteristic can result in an underestimation of the instantaneous flood peaks, with possible risk of design failure. It is more common to find a longer record of daily flow series (observed or modelled using daily rainfall series) than sub‐daily flow data. This paper describes a novel approach, known as the steepness index unit volume flood hydrograph approach, for disaggregating daily flood flows into sub‐daily flows that takes advantage of the strong relationship between the standardized instantaneous flood peak and the standardized daily flood hydrograph rising‐limb steepness index. The strength of this relationship, which is considerably stronger than the relationship between the standardized flood peak and the event flood volume, is shown using data from six rivers flowing into the Gippsland Lakes in southeast Australia. The results indicate that the steepness index unit volume flood hydrograph approach can be used to disaggregate modelled daily flood flows satisfactorily, but its reliability is dependent on a model's ability to simulate the standardized daily flood hydrograph rising‐limb steepness index and the event flood volume. Copyright © 2007 John Wiley & Sons, Ltd.     

Simulation of probability distributions commonly used in hydrological frequency analysis

Random variable simulation has been applied to many applications in hydrological modelling, flood risk analysis, environmental impact assessment, etc. However, computer codes for simulation of distributions commonly used in hydrological frequency analysis are not available in most software libraries. This paper presents a frequency‐factor‐based method for random number generation of five distributions (normal, log–normal, extreme‐value type I, Pearson type III and log‐Pearson type III) commonly used in hydrological frequency analysis. The proposed method is shown to produce random numbers of desired distributions through three means of validation: (1) graphical comparison of cumulative distribution functions (CDFs) and empirical CDFs derived from generated data; (2) properties of estimated parameters; (3) type I error of goodness‐of‐fit test. An advantage of the method is that it does not require CDF inversion, and frequency factors of the five commonly used distributions involves only the standard normal deviate. Copyright © 2006 John Wiley & Sons, Ltd.     

Clay mineral,geochemical and Sr–Nd isotopic fingerprinting of sediments in the Murray–Darling fluvial system,southeast Australia

The Willyama Supergroup of the Broken Hill region in southern Australia consists of supracrustal sedimentary and magmatic rocks, formed between 1810 and 1600 Ma. A statistical analysis of nearly 2000 SHRIMP U–Pb zircon spot ages, compiled from published and unpublished sources, provides evidence for three distinct tectonostratigraphic successions and four magmatic events during this interval. Succession 1 includes Redan Geophysical Zone gneisses and the lower part of the Thackaringa Group (Cues Formation). These rocks were deposited after 1810 Ma and host granite sills of the first magmatic event (1710–1700 Ma). Succession 2 includes the upper Thackaringa Group (Himalaya Formation), the Broken Hill Group and the Sundown Group and was deposited between 1710 and 1660 Ma. These rocks all contain detrital zircons from the first magmatic event (1710–1700 Ma) and in some cases from the second magmatic event (1690–1680 Ma). The second magmatic event (1690–1680 Ma) was bimodal, resulted from crustal extension, and was coeval with deposition of the Broken Hill Group and deepening of the basin. With this event a mafic sill swarm focused in the Broken Hill Domain. Mafic sills lack any trace of inheritance, unlike the granitoids that commonly contain inherited zircons typical of the supracrustal sediments. Succession 3, the Paragon Group and equivalents were deposited after 1660 Ma, but before a regional metamorphic event at 1600 Ma. Metamorphism was closely followed by inversion of the succession into a fold‐and‐thrust belt, accompanied by a fourth late to post‐orogenic magmatic event (ca 1580 Ma) characterised by granite intrusion and regional acid volcanism (the local equivalents of the Gawler Range Volcanics in South Australia).     

A virtual experiment technique on the elementary behaviour of granular materials with discrete element method

Multi‐scale investigations aided by the discrete element method (DEM) play a vital role for current state‐of‐the‐art research on the elementary behaviour of granular materials. Similar to laboratory tests, there are three important aspects to be considered carefully, which are the proper stress/strain definition and measurement, the application of target loading paths and the designed experiment setup, to be addressed in the present paper. Considering the volume sensitive characteristics of granular materials, in the proposed technique, the deformation of the tested specimen is controlled and measured by deformation gradient tensor involving both the undeformed configuration and the current configuration. Definitions of Biot strain and Cauchy stress are adopted. The expressions of them in terms of contact forces and particle displacements, respectively, are derived. The boundary of the tested specimen consists of rigid massless planar units. It is suggested that the representative element uses a convex polyhedral (polygonal) shape to minimize possible boundary arching effects. General loading paths are described by directly specifying the changes in the stress/strain invariants or directions. Loading can be applied in the strain‐controlled mode by specifying the translations and rotations of the boundary units, or in the stress‐controlled mode by using a servo‐control mechanism, or in the combination of the two methods to realize mixed boundary conditions. Taking the simulation results as the natural consequences originated from a complex system, virtual experiments provide particle‐scale information database to conduct multi‐scale investigations for better understanding in granular material behaviours and possible development of the constitutive theories provided the qualitative similarity between the simulation results from virtual experiments and observations on real material behaviour. Copyright © 2011 John Wiley & Sons, Ltd.