Location,location, location: changes in the diversity of animal resources exploited by Tardiglacial humans in northern Spain

Prey diversity of Homo sapiens in north Iberia exhibited shifts associated with climate warming after the Last Glacial Maximum, as deduced from animal remains obtained from Upper Palaeolithic caves in Asturias (northern Spain). Significant association of prey diversity and river proximity was found during adverse climate conditions, indicating opportunistic predation. A high increase of prey diversity was found when small animals (including fish) were taken into account in archaeological records, highlighting the importance of this type of prey in Palaeolithic diets. Copyright © 2009 John Wiley & Sons, Ltd.     

Experimental study of extreme shear stress for shallow flow under simulated rainfall

The raindrop impact and overland flow are two major factors causing soil detachment and particle transportation. In this study, the turbulent characteristics of the shallow rain‐impacted water flow were investigated using a 2‐D fibre‐optic laser Doppler velocimetry (FLDV) and an artificial rainfall simulator. The fluctuating turbulent shear stress was computed using digital data processing techniques. The experimental data showed that the Reynolds shear stress follows a probability distribution with heavy tails. The tail probability increases with an increase of rainfall intensity or raindrop diameter, and it decreases with an increase of Reynolds number. A modified empirical equation was derived using both the raindrop diameter and rainfall intensity as independent variables to provide a better prediction of the Darcy‐Weisbach friction coefficient f under rainfall conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

Bayesian analysis of stage–fall–discharge models for gauging stations affected by variable backwater

Parsimonious stage–fall–discharge rating curve models for gauging stations subject to backwater complications are developed from simple hydraulic theory. The rating curve models are compounded in order to allow for possible shifts in the hydraulics when variable backwater becomes effective. The models provide a prior scientific understanding through the relationship between the rating curve parameters and the hydraulic properties of the channel section under study. This characteristic enables prior distributions for the rating curve parameters to be easily elicited according to site‐specific information and the magnitude of well‐known hydraulic quantities. Posterior results from three Norwegian and one American twin‐gauge stations affected by variable backwater are obtained using Markov chain Monte Carlo simulation techniques. The case studies demonstrate that the proposed Bayesian rating curve assessment is appropriate for developing rating procedures for gauging stations that are subject to variable backwater. Copyright © 2009 John Wiley & Sons, Ltd.     

Simultaneous identification of parameter,initial condition,and boundary condition in groundwater modelling

In this paper, we perform an inverse method to simultaneously estimate aquifer parameters, initial condition, and boundary conditions in groundwater modelling. The parameter estimation is extended to a complete inverse problem that makes the calibrated groundwater flow model more realistic. The adjoint state method, the gradient search method, and the least square error algorithm are combined to build the optimization procedure. Horizontal two‐dimensional groundwater flow in a confined aquifer is exemplified to demonstrate the correlation between unknowns, the contribution of observation, as well as the suitability of applying the inverse method. The correlation analysis shows the connection between storage coefficient and initial condition. Besides, transmissivity and boundary conditions are also highly correlated. More observations at different location and time are necessary to provide sufficient information. A time series of unsteady head is requested for estimation of storage coefficient and initial condition. Observation near boundary is very effective for boundary condition estimation. The observation at pumping well mostly contributes to the estimation of transmissivity. According to all observations, it is possible to identify parameters, initial condition, and boundary condition simultaneously. Furthermore, the results not only illustrate the traditional assumption of known boundary condition but also initial condition, which may cause an incorrect estimation. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical modelling of the potential effects of a dam on a coastal aquifer in S. Spain

This study presents the results of a three‐dimensional variable‐density numerical modelling of the Motril‐Salobreña coastal aquifer and the possible effects of the entry into service in May 2005 of the Rules Dam, located just 17 km from the coast. Present parameters of the Motril‐Salobreña aquifer show that the system's conditions are very similar to a natural regime. The dam will substantially alter aquifer recharge, as the entry flow through the alluvial sediments of the Guadalfeo River will be entirely cut off or drastically reduced. Different scenarios reproducing the possible evolution of the aquifer under operation of the Rules Dam have been modelled. In most cases, results indicate that the conditions of the aquifer would worsen, with a general advance of the freshwater–saltwater interface. The area with most risk of saltwater intrusion is the old mouth of the Guadalfeo River, where the mixing zone could advance 1200 m inland. It is proposed that maintaining a 5–6 Mm³ year^?1 ‘ecological flow’ in the Guadalfeo River could prevent this saline advance. This application demonstrates that variable‐density models are potentially useful tools for estimating the effects of dams on the hydrodynamic and hydrochemical conditions of a coastal aquifer. Copyright © 2009 John Wiley & Sons, Ltd.     

A complex variable solution for stress and displacement field around a lined circular tunnel at great depth

By applying the series expansion technique in the complex variable method established by Muskhelishvili, the plane elasticity problem for the stress and displacement field around a lined circular tunnel in conjunction with the consideration of misfit and interaction between the liner and the surrounding geomaterial is dealt with. The tunnel is assumed to be driven in a homogeneous and isotropic geomaterial. The coefficients in the Laurent series expansion of the stress functions are determined. The complex potentials in the liner and the surrounding geomaterial are explicitly derived, respectively. As an example, the case of a lined circular tunnel located in an isotropic initial stress field but subjected to uniform internal pressure is numerically considered. Numerical results indicate that the installation of tunnel liner can reduce the influences of the tunnel excavation on the in situ displacement and stress fields. However, the relative thickness and rigidity of the liner should be in an appropriate range. In addition, the effect of the tunnel excavation upon the displacement field is more significant than that upon the stress field. As far as the stress field in the surrounding geomaterial is concerned, when the ratio between the cover depth of tunnel and the tunnel radius is larger than 5, the results for the stress field in the paper are applicable. When the ratio between the tunnel depth and the tunnel radius is larger than 20, the results are applicable for the displacement field. Copyright © 2008 John Wiley & Sons, Ltd.     

Cratering and modification of wet‐target craters: Projectile impact experiments and field observations of the Lockne marine‐target crater (Sweden)

Abstract— Marine impacts are one category of crater formation in volatile targets. At target water depths exceeding the diameter of the impactor, the zones of vaporization, melting, and excavation of the standard land‐target cratering model develop partially or entirely in the water column. The part of the crater that has a potential of being preserved (seafloor crater) may to a great extent be formed by material emplacement and excavation processes that are very different from land‐target craters. These processes include a high‐energy, water‐jet‐driven excavation flow. At greater water depths, the difference in strength of the target layers causes a concentric crater to evolve. The crater consists of a wide water cavity with a shallow excavation flow along the seabed surrounding a nested, deeper crater in the basement. The modification of the crater is likewise influenced by the water through its forceful resurge to fill the cavity in the water mass and the seafloor. The resurge flow is strongly erosive and incorporates both ejecta and rip‐up material from the seabed surrounding the excavated crater. A combination of field observations and impact experiments has helped us analyze the processes affecting the zone between the basement crater and the maximum extent of the water cavity. The resurge erosion is facilitated by fragmentation of the upper parts of the solid target caused by a) spallation and b) vibrations from the shallow excavation flow and, subsequently, c) the vertical collapse of the water cavity rim wall. In addition, poorly consolidated and saturated sediments may collapse extensively, possibly aided by a violent expansion of the pore water volume when it turns into a spray during passage of the rarefaction wave. This process may also occur at impacts into water‐saturated targets without an upper layer of seawater present. Our results have implications for impacts on both Earth and Mars, and possibly anywhere in the solar system where volatiles exist/have existed in the upper part of the target.     

The Gyldén‐type problem revisited: More refined analytical solutions

We resume and consistently extend our previous researches concerning the Gyldén‐type problem (a two‐body problem with time‐dependent equivalent gravitational parameter). To approach most of the concrete astronomical situations to be modelled in this way, we consider a periodic small perturbation. For the nonresonant case, we present a second‐order analytical solution. For the resonant case, we adopt the most realistic astronomical situation: only one dominant term of the Hamiltonian. In this case we point out a fundamental model of resonance, common to every resonant situation, and, moreover, identical to the first fundamental model of resonance . Considering the simplest model of periodic change of the equivalent gravitational parameter, we .nd that all possible resonances are con.ned to the first fundamental model. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

New observations on high‐pressure phases in a shock melt vein in the Villalbeto de la Peña meteorite: Insights into the shock behavior of diopside

The petrology and mineralogy of shock melt veins in the L6 ordinary chondrite host of Villalbeto de la Peña, a highly shocked, L chondrite polymict breccia, have been investigated in detail using scanning electron microscopy, transmission electron microscopy, Raman spectroscopy, and electron probe microanalysis. Entrained olivine, enstatite, diopside, and plagioclase are transformed into ringwoodite, low‐Ca majorite, high‐Ca majorite, and an assemblage of jadeite‐lingunite, respectively, in several shock melt veins and pockets. We have focused on the shock behavior of diopside in a particularly large shock melt vein (10 mm long and up to 4 mm wide) in order to provide additional insights into its high‐pressure polymorphic phase transformation mechanisms. We report the first evidence of diopside undergoing shock‐induced melting, and the occurrence of natural Ca‐majorite formed by solid‐state transformation from diopside. Magnesiowüstite has also been found as veins injected into diopside in the form of nanocrystalline grains that crystallized from a melt and also occurs interstitially between majorite‐pyrope grains in the melt‐vein matrix. In addition, we have observed compositional zoning in majorite‐pyrope grains in the matrix of the shock‐melt vein, which has not been described previously in any shocked meteorite. Collectively, all these different lines of evidence are suggestive of a major shock event with high cooling rates. The minimum peak shock conditions are difficult to constrain, because of the uncertainties in applying experimentally determined high‐pressure phase equilibria to complex natural systems. However, our results suggest that conditions between 16 and 28 GPa and 2000–2200 °C were reached.