Double-averaged velocity profiles over fixed dune shapes

Spatially averaged profiles of time averaged velocity, using integrals over thin horizontal slabs (Cartesian double average), are employed in characterizing the flow over fixed dune shapes. For comparison the spatial averaging method of Smith and McLean (1977) that averages along lines at constant distance from the local bed elevation is also investigated. The Cartesian double averaged profiles of the inverse of the velocity shear are nearly constant below the crest elevation, but increase rapidly above the crest level. This results in velocity profiles that increase linearly with distance from the bed below the crest. Above the crest it can be argued that the velocity increases logarithmically, but a power law profile can also be argued. Spatially averaged eddy viscosity profiles are calculated by multiplying the average Reynolds stress by the inverse shear. The resulting profile is more complex than the uniform flow counterpart.     

SWAT.nz: New-Zeland-based “Sand Waves and Turbulence” experimental programme

The SWAT.nz (“New-Zealand-based Sand Waves and Turbulence”) research programme was carried out to advance understanding of subaqueous sand waves. The programme was based around detailed measurements at varying scales of bed morphologies and associated flow fields as sand waves formed from plane-bed conditions and grew to equilibrium. This paper outlines the philosophy and details of the SWAT.nz programme, with the aim of providing insight into experiment and analysis design and methodologies for studies of highly-variable bed surfaces and flows. Example challenges addressed in the SWAT.nz programme include the measurement over large spatial domains of developing flow fields and three-dimensional bed morphology, including flow measurements below roughness (sand-wave) crests, and how to interpret the collected measurements. Insights into sand-wave dynamics that have arisen from the programme are presented to illustrate the values of the SWAT.nz programme and the developed methodologies. Results are presented in terms of mobile-bed processes, and flow-bed interaction and flow processes for fixed-bed roughness and erodible beds, respectively.     

Modelling topoclimatic controls on palaeoglaciers: implications for inferring palaeoclimate from geomorphic evidence

Evidence of small glaciers is often used to infer past atmospheric climate through calculation of steady-state ELAs. However, if topographic niches such as shading or windblown-snow augmented mass-balance then ELAs cannot reflect regional climate and determining the significance of these topoclimates is therefore important. The Brecon Beacons, South Wales, contains upland glacial landforms dating to the Younger Dryas (11,000–10,000 yr BP) when local climate was at the threshold for glaciation. This case study categorises topoclimate using three-dimensional modelling of topography and reconstructed palaeoglaciers from two sites containing mapped moraines whose orientation suggest complex patterns of deglaciation. Ablation season solar radiation is modelled over multiple ice-surfaces as shade from surrounding topography and intensity from ice-surface incidence angle. Snowblow and avalanching potential models are also used and the significance of all topoclimate variables assessed against a mass-balance deficit calculated for each glacier given palaeoclimate models for the region. Results demonstrate that both glaciers were likely to be heavily reliant on topoclimate and that previous studies underestimate the significance of solar radiation. Modelled over multiple ice-surfaces reflecting patterns of recession, results indicate that the distribution of topoclimate variables predicts the style of deglaciation at both sites, possibly explaining the complexity of glacial evidence in this environment.     

A snapshot of environmental iodine and selenium in La Pampa and San Juan provinces of Argentina

Soil and water samples were collected from farmsteads and provincial towns across the provinces of La Pampa and San Juan in Argentina. Inductively coupled plasma mass spectrometry was used for the determination of iodine in water following addition of TMAH to 1% v/v and soils extracted with 5% TMAH. Iodine in agricultural soils was in the range of 1.3–20.9 mg/kg in La Pampa located in central Argentina and 0.1–10.5 mg/kg in San Juan located in the northwest Andean region of Argentina, compared to a worldwide mean of 2.6 mg/kg. Mean selenium concentrations for soils from both provinces were 0.3 mg/kg, compared to a worldwide mean of 0.4 mg/kg. The majority of soils were slightly alkaline at pH 6.7 to 8.8. The organic content of soils in La Pampa was 2.5–5.9% and in San Juan 0.1–2.3%, whilst, mobile water extractable soil-iodine was 1–18% for La Pampa and 2–42% for San Juan. No simple relationship observed for pH and organic content, but mobile iodine (%) was highest when organic content was low, higher for lower total iodine concentrations and generally highest at pH > 7.5. Water drawn for drinking or irrigation of a variety of crops and pasture was found to range from 52 to 395 µg/L iodine and 0.8 to 21.3 µg/L selenium in La Pampa and 16–95 µg/L iodine and 0.6 to 8.2 µg/L selenium in San Juan. The water samples were all slightly alkaline between pH 8 and 10. Water–iodine concentrations were highest at pH 7.8 to 8.8 and in groundwaters positively correlated with conductivity. Raw water entering water treatment works in La Pampa was reduced in iodine content from approximately 50 µg/L in raw water to 1 µg/L in treated drinking water, similar to levels observed in regions experiencing iodine deficiency.     

A multi-isotope (δD, δO,Sr/Sr,and δB) approach for identifying saltwater intrusion and resolving groundwater evolution along the Western Caprock Escarpment of the Southern High Plains,New Mexico

Declining water levels in arid and semi-arid regions increase an aquifer’s vulnerability to natural and anthropogenic influences. A multi-isotope (δD, δ¹⁸O, ⁸⁷Sr/⁸⁶Sr, and δ¹¹B) approach was used to resolve the geochemical evolution of groundwater in a declining aquifer in a semi-arid region of the southwestern USA as groundwater composition reacts to source-water mixing, cross-formational flow including saltwater intrusion, water–rock interaction, and likely agricultural recharge. Sub-aquifers or local flow systems are present in the Southern High Plains aquifer along the Western Caprock Escarpment in New Mexico, and the study site’s local flow system contains a Na–Cl, high dissolved-solids groundwater that flows from the escarpment until it mixes with a high quality regional aquifer or regional flow system. The local flow system contains water that is similar in composition to the underlying, upper Dockum Group aquifer. Saltwater found in the upper Dockum Group aquifer likely originates in the adjacent Pecos River Basin and crosses beneath or possibly through the hydrologic divide of the Western Caprock Escarpment. Strontium concentrations of 0.9–31 mg/L and a ⁸⁷Sr/⁸⁶Sr range of 0.70845–0.70906 were sufficient to estimate source-water fractions, mixing patterns, and contributions from chemical weathering through mass balance inverse calculations. Boron concentrations (59–1740 mg/L) and δ¹¹B values (+6.0–+46.0‰) were used to confirm source-water mixing, further evaluate water–rock interaction, and examine the influence of possible agricultural recharge. Alteration of B concentrations and δ¹¹B values in an area of likely agricultural recharge indicated the loss of B and decrease in δ¹¹B values likely from plant uptake, adsorption, and weathering contributions in the soil/vadose zone prior to recharge. The effectiveness of ⁸⁷Sr/⁸⁶Sr and δ¹¹B for resolving the geochemical influences in groundwater in the Southern High Plains along the Western Caprock Escarpment allowed for the reinterpretation of the isotopic composition of water that has been shown to be highly variable in the Southern High Plains. This study shows the utility of a multi-isotope approach for resolving the geochemical evolution of groundwater in an aquifer that has a complex relationship with underlying aquifers and the applicability of these isotopes as indicators of the alteration of source waters from natural or anthropogenic influences.     

Cassini observations of Saturn's inner plasmasphere: Saturn orbit insertion results

We present new and definitive results of Cassini plasma spectrometer (CAPS) data acquired during passage through Saturn's inner plasmasphere by the Cassini spacecraft during the approach phase of the Saturn orbit insertion period. This analysis extends the original analysis of Sittler et al. [2005. Preliminary results on Saturn's inner plasmasphere as observed by Cassini: comparison with Voyager. Geophys. Res. Lett. 32, L14S07, doi:10.1029/2005GL022653] to L∼10 along with also providing a more comprehensive study of the interrelationship of the various fluid parameters. Coincidence data are sub-divided into protons and water group ions. Our revised analysis uses an improved convergence algorithm which provides a more definitive and independent estimate of the spacecraft potential Φ_SC for which we enforce the protons and water group ions to co-move with each other. This has allowed us to include spacecraft charging corrections to our fluid parameter estimations and allow accurate estimations of fluctuations in the fluid parameters for future correlative studies. In the appendix we describe the ion moments algorithm, and minor corrections introduced by not weighting the moments with sinθ term in Sittler et al. [2005] (Correction offset by revisions to instruments geometric factor). Estimates of the spacecraft potential and revised proton densities are presented. Our total ion densities are in close agreement with the electron densities reported by Moncuquet et al. [2005. Quasi-thermal noise spectroscopy in the inner magnetosphere of Saturn with Cassini/RPWS: electron temperatures and density. Geophys. Res. Lett. 32, L20S02, doi:10.1029/2005GL022508] who used upper hybrid resonance (UHR) emission lines observed by the radio and plasma wave science (RPWS) instrument. We show a positive correlation between proton temperature and water group ion temperature. The proton and thermal electron temperatures track each with both having a positive radial gradient. These results are consistent with pickup ion energization via Saturn's rotational electric field. We see evidence for an anti-correlation between radial flow velocity V_R and azimuthal velocity V_φ, which is consistent with the magnetosphere tending to conserve angular momentum. Evidence for MHD waves is also present. We show clear evidence for outward transport of the plasma via flux tube interchange motions with the radial velocity of the flow showing positive radial gradient with functional dependence for 4<L<10 (i.e., if we assume to be diffusive transport then D_LL∼D₀L¹¹ for fixed stochastic time step δt). Previous models with centrifugal transport have used D_LL∼D₀L³ dependence. The radial transport seems to begin at Enceladus’ L shell, L∼4, where we also see a minimum in the W⁺ ion temperature . For the first time, we are measuring the actual flux tube interchange motions in the magnetosphere and how it varies with radial distance. These observations can be used as a constraint with regard to future transport models for Saturn's magnetosphere. Finally, we evaluate the thermodynamic properties of the plasma, which are all consistent with the pickup process being the dominant energy source for the plasma.     

Consistent parameter constraints for soil hydraulic functions

Parameters of functions to describe soil hydraulic properties are derived from measurements by means of parameter estimation. Of crucial importance here is the choice of correct constraints in the parameter space. Often, the parameters are mere shape parameters without physical meaning, giving flexibility to the model. A fundamental requirement is that the hydraulic functions are monotonic: the retention function and the conductivity function can only decrease as the capillary suction increases. A stricter physical requirement for the conductivity function is that its decrease with respect to saturation is at least linear. This linear decrease would occur if all pores of a capillary bundle had an equal radius. In the first part of this contribution, we derive constraints for the so-called tortuosity parameter of the Mualem conductivity model, which allow highest possible flexibility on one hand and guarantee physical consistency on the other hand. In combination with the retention functions of Brooks and Corey, van Genuchten, or Durner, such a constraint can be expressed as a function of the pore-size distribution parameters. In the second part, we show that a common modification of retention models, which is applied to reach zero water content at finite suction, can lead to the physically unrealistic case of increasing water content with increasing suction. We propose a solution for this problem by slightly modifying these models and introducing a correct parameter constraint.     

Geoacoustic character,sedimentology and chronology of a cross-shelf Holocene sediment deposit off Cabo Frio,Brazil (southwest Atlantic Ocean)

The Cabo Frio region in the state of Rio de Janeiro, southeast coast of Brazil, is characterized by a local coastal upwelling system and converging littoral sediment transport systems that are deflected offshore at Cabo Frio, as a consequence of which a thick cross-shelf sediment deposit has developed over time. To investigate the evolution of this muddy deposit, geophysical, sedimentological and geochemical data from four sediment cores (3.8–4.1 m in length) recovered in water depths between 88 and 141 m were analyzed. The high-resolution seismic data show variable sediment thicknesses ranging from 1 to 20 m, comprising two sedimentary units separated by a high-impedance layer at a depth of about 10 m below the seafloor at the coring sites. According to the available age datings, the upper sedimentary unit is late Pleistocene to Holocene in age, whereas the lower unit (not dated) must, by implication, be entirely Pleistocene in age. The boomer-seismic reflection signal can be divided into three echo-types, namely transparent (inner shelf), stratified (middle shelf) and reflective (outer shelf), each type seemingly related to the local sediment composition. The upper 4 m of the upper sedimentary unit is dominated by silty sediment on the middle shelf, and by upward-fining sediments (silty sand to sandy silt) on the inner and outer shelf. The downcore trends of P-wave velocity, gamma-ray density and acoustic impedance are largely similar, but generally reversed to those of water and organic carbon contents. Total organic carbon contents increase with decreasing mean grain size, periodic fluctuations suggesting temporal changes in the regional hydrodynamics and primary productivity fuelled by the local upwelling system. The reconstruction of sedimentation rates in the course of the Holocene is based on 35 AMS age datings of organic material recovered from variable downcore depths. These range from a maximum of 13.3 cm/decade near the base of the inner shelf core (7.73–7.70 ka BP) to generally very low values (<0.11 cm/century) over the last thousand years in all cores. Over the last 6 ka there appear to have been three distinct sedimentation peaks, one between 6 and 5 ka BP, another between 4 and 3 ka PB, and one around 1 ka BP. Due to different time intervals between dates, not every peak is equally well resolved in all four cores. Based on the similar sedimentology of the inner and outer shelf cores, an essentially identical sedimentation model is proposed to have been active in both cases, albeit at different times. Thus, already during the last glacial maximum, alongshore sediment transport was deflected offshore by a change in shoreline orientation caused by the Cabo Frio structural high. The source of terrigenous material was probably a barrier-island complex that was subsequently displaced landward in the course of sea-level rise until it stabilized some 6.5 ka BP along the modern coast.