Estimation of surface roughness (z0) over a stabilizing coastal dune field based on vegetation and topography

Accurate knowledge of the surface roughness and the resultant wind speed are important for many applications, such as climatic models, wind power meteorology, agriculture and erosion hazards, especially on sand dunes in arid and semi‐arid environments, where vegetation cover is scarce. In this study we aimed at quantifying the effects of vegetation cover and topography on surface roughness over a stabilizing dune field on the southern coast of Israel. Forty‐six wind measurements were made at various distances from the coastline, ranging from 10 to 2800 m, and z₀ values were calculated from the wind measurements based on the ratio between the wind gust and the average wind speed. We estimated vegetation cover using the soil adjusted vegetation index (SAVI) from Landsat satellite images for the upwind sector at various lengths, ranging from 15 to 400 m, and based on digital elevation models and differential GPS field measurements we calculated the topographic variable of the relative heights of the stations. z₀ values were positively correlated with the winter SAVI values (r = 0·87 at an upwind length of 200 m) and negatively correlated with the relative height (r = ?0·68 at an upwind length of 200–400 m for the inland dune stations). Using these variables we were able to create a map of estimated z₀ values having an accuracy of over 64%. Such maps provide a better understanding of the spatial variability in both wind speed and sand movement over coastal dune areas. Copyright © 2007 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.     

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.     

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.     

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.     

Cryopreservation of Veliger Larvae of Trumpet Shell,Charonia sauliae: an Essential Preparation to Artificial Propagation

Trumpet shell, Charonia sauliae, is an endangered and valuable species, but its artificial propagation protocol has not been successfully established. To estimate the possibility of cryopreservation for larvae of C. sauliae, which is a potential preparation for its artificial reproduction and further research, in this study a protocol for the cryopreservation of veliger larvae of trumpet shell was optimized. Through a two-step cryopreservation procedure, four kinds of cryoprotectants (ethylene glycol, 1, 2-propanediol,dimethyl sulfoxide and glycerol) were employed at three concentrations (1.0, 1.5 and 2.0 molL-1) respectively and survival rates of larvae were determined after a storage of 1h. The larvae frozen with these four cryoprotectants after 1h storage were cultured, and then survival rates were determined at 24, 72 and 120h after thawing. Dimethyl sulfoxide at a concentration of 1.5 molL-1 showed the best protective effect in all experiments (p<0.05). And survival rates of larvae frozen with dimethyl sulfoxide were determined after 1,7 and 15 d of storage. The survival rates of larvae frozen with 1.5 molL-1 dimethyl sulfoxide after 1h, 1d, 7 d and 15 d of storage were 80.77% ±7.51%, 80.34% ±11.28%, 83.10% ±9.14% and 77.23%±6.22% respectively. No significant differences in survival rates of larvae frozen with dimethyl sulfoxide were observed after various storage periods (p>0.05).     

Early solar system aqueous activity: K isotope evidence from Allende

The alkali element K is moderately volatile and fluid mobile; thus, it can be influenced by both primary processes (evaporation and recondensation) in the solar nebula and secondary processes (thermal and aqueous alteration) in the parent body. Since these primary and secondary processes would induce different isotopic fractionations, K isotopes could become a potential tracer to distinguish them. Using recently developed methods with improved precision (0.05‰, 95% confidence interval), we systematically measured the K isotopic compositions and major/trace elemental compositions of chondritic components (18 chondrules, 3 CAIs, 2 matrices, and 5 bulks) in the carbonaceous chondrite fall Allende. Among all the components analyzed in this study, CAIs, which formed initially under high‐temperature conditions in the solar nebula and were dominated by nominally K‐free refractory minerals, have the highest K₂O content (average 0.53 wt%) and have K isotope compositions most enriched in heavy isotopes (δ⁴¹K: ?0.30 to ?0.25‰). Such an observation is consistent with previous petrologic studies that show CAIs in Allende have undergone alkali enrichment during metasomatism. In contrast, chondrules contain lower K₂O content (0.003–0.17 wt%) and generally lighter K isotope compositions (δ⁴¹K: ?0.87‰ to ?0.24‰). The matrix and bulks are nearly identical in K₂O content and K isotope compositions (0.02–0.05 wt%; δ⁴¹K: ?0.62 to ? 0.46‰), which are, as expected, right in the middle of CAIs and chondrules. This strongly indicates that most of the chondritic components of Allende suffered aqueous alteration and their K isotopic compositions are the ramification of Allende parent‐body processing instead of primary nebular signatures. Nevertheless, we propose the small K isotope fractionations observed (< 1‰) among Allende components are likely similar to the overall range of K isotopic fractionation that occurred in nebular environment. Furthermore, the K isotope compositions seen in the components of Allende in this study are consistent with MC‐ICP‐MS analyses of the components in ordinary chondrites, which also show an absence of large (10‰) isotope fractionations. This is not expected as evaporation experiments in nebular conditions suggest there should be large K isotopic fractionations. Nevertheless, possible nebular processes such as chondrules back exchanging with ambient gas when they formed could explain this lack of large K isotopic variation.     

The evolution of the z distribution of normal neutron stars in the Galaxy

Under the two initial 1‐D one parameter velocity distribution forms (one is normal, the other is exponential), the z direction scale height evolution of normal neutron stars in the Galaxy is studied by numerical simulation. We do statistics for the cases at different time segments, also do statistics for the cumulative cases made of each time segment. The results show in the cumulative cases the evolution curves of the scale heights are smoother than in the each time segment, i.e., the cumulation improve the signal‐to‐noise ratio. Certainly the evolution cases are different at different Galactic disk locations, which also have very large difference from the average cases in the whole disk. In the initial stages of z evolution of normal neutron stars, after the beginning transient states, the cumulative scale heights increase linearly with time, and the cumulative scale height increasing rates have linear relationship with the initial velocity distribution parameters, which have larger fluctuation in the vicinity of the Sun than in the whole disk. We utilize the linear relationship of the cumulative scale height increasing rates vs. the initial velocity distribution parameters in the vicinity of the Sun to make comparison with the observation near the Sun. The results show if there is no magnetic decay, then the deserved initial velocity parameters are obvious lower than the present well known results from some authors; whereas if introducing magnetic decay, for the 1‐D normal case we can make consistence among concerning results using magnetic decay time values which are supported by some authors, while for the 1‐D exponential case the results show the lackness of young pulsar samples in the larger z in the vicinity of the Sun (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)