An advanced process for evaluating a linear dielectric constant–bulk electrical conductivity model using a capacitance sensor in field conditions

Abstract

The Hilhorst model was used to convert bulk electrical conductivity (σ_b) to pore water electrical conductivity (σ_p) under laboratory conditions by using the linear relationship between the soil dielectric constant (ε_b) and σ_b. In the present study, applying the linear relationship ε_b–σ_b to data obtained from field capacitance sensors resulted in strong positive autocorrelations between the residuals of that regression. We were able to derive an accurate offset of the relationship ε_b–σ_b and to estimate the evolution of σ_p over time by including a stochastic component to the linear model, rearranging it to a time-varying dynamic linear model (DLM), and using Kalman filtering and smoothing. The offset proved to vary for each depth in the same soil profile. A reason for this might be the changes in soil temperature along the soil profile.

Editor D. Koutsoyiannis; Associate editor M.D. Fidelibus     

Coastal dynamics off Northwest Iberia during a stormy winter period

The consequences of a stormy winter period (2009/2010) on the shelf and coastal dynamics off Northwest Iberia are analysed by using model results in combination with the set of available observations in the frame of the Iberian Margin Ocean Observatory (RAIA), a cross-border infrastructure among North Portugal and Galicia (Spain). During the study winter, the frequent arrival of weather fronts forced river plumes to flow along the inner shelf in a fast (>1 m?s^?1) jet-like structure. The buoyant current strongly influenced the outer rías, the name of the estuaries in the region, where a strong decay of surface salinity (<10.5) has been observed. Once the weather front has passed, the wind reversal forced the offshore expansion of river plumes and also the development of a winter upwelling event. Thermohaline patterns in both model and observations revealed an intrusion of warm (>15 °C) and salty (>35.9) waters into the rías associated with the Iberian Poleward Current. Finally, some Lagrangian modelling experiments were performed to analyse the transport ability of the plume and the effect that could have had in the biological material trapped on it. The experiments reveal that an overall northward displacement of surface particles will be expected after several alternate wind events.     

The use of radon as tracer in environmental sciences

Radon can be used as a naturally occurring tracer for environmental processes. By means of grab-sampling or continuous monitoring of radon concentration, it is possible to assess several types of dynamic phenomena in air and water. We present a review of the use of radon and its progeny at the University of Cantabria. Radon can be an atmospheric dynamics indicator related with air mass interchange near land-sea discontinuities as well as for the study of vertical variations of air parameters (average values of different types of stability: 131–580 Bq m^?3). Concerning indoor gas, we present some results obtained at Altamira Cave (Spain): from 222 to 6549 Bq m^?3 (Hall) and from 999 to 6697 Bq m^?3 (Paintings Room). Finally, variations of radon concentration in soil (0.3 to 9.1 kBq m^?3) and underground water (values up to 500 Bq l^?1) provide relevant information about different geophysical phenomena.     

Non-parametric copulas for circular–linear and circular–circular data: an application to wind directions

This paper proposes a nonparametric approach to estimating the dependence relationships between circular variables and other circular or linear variables using copulas. The proposed method is based on the use of Bernstein copulas which are a very flexible class of non-parametric copulas which allows for the approximation of any kind of dependence structure, including non symmetric relationships. In particular, we present a simple procedure to adapt Bernstein copulas to the circular framework and guarantee that the constructed bivariate distributions are strictly continuous. We provide two illustrative case studies, the first on the relation between wind direction and quantity of rainfall in the North of Spain and the second on the dependence between the wind directions in two nearby buoys at the Atlantic ocean.     

A simple model accounting for the uptake,transport, and deposition of wind‐eroded mineral particles in the hyperarid coastal Atacama Desert of northern Chile

As previously observed in marine sediments collected downwind of African or South American continental sources, recent studies of sediment cores collected at the bottom of Mejillones Bay in north Chile (23°S) show a laminated structure in which the amount of particles of aeolian origin and their size create significant differences between the layers. This suggests inter‐annual to inter‐decadal variations in the strength of the local southerly winds responsible for (1) the erosion of the adjacent hyperarid surface of the Mejillones Pampa, and (2) the subsequent transport of the eroded particles towards the bay. A simple model accounting for the vertical uptake, transport, and deposition of the particles initially set into motion by wind at the surface of the pampa is proposed. This model, which could be adapted to other locations, assumes that the initial rate of (vertical) uptake is proportional to the (horizontal) saltation flux quantified by means of White's equation, that particles are lifted to a height (H), increasing with the magnitude of turbulence, and that sedimentation progressively removes the coarsest particles from the air column as it moves towards the bay. In this model, the proportionality constant (A) linking the vertical flux of particles with the horizontal flux, and the injection height (H) control the magnitude and size distribution of the deposition flux in the bay. Their values are determined using the wind speed measured over the pampa and the size distribution of particles collected in sediment traps deployed in the bay as constraints. After calibration, the model is used to assess the sensitivity of the deposition flux to the wind intensity variations. The possibility of performing such quantitative studies is necessary for interpreting precisely the variability of the aeolian material in the sediment cores collected at the bottom of Mejillones Bay. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamics of lava flow fronts, Arenal Volcano, Costa Rica

Arenal Volcano has effused basaltic andesite lava flows nearly continuously since September, 1968. The two different kinds of material in flows, lava and lava debris, have different rheologic properties and dynamic behavior. Flow morphology depends on the relationship between the amount and distribution of the lava and the debris, and to a lesser extent the ground morphology.Two main units characterize the flows: the channel zone and the frontal zone. The channel zone consists of two different units, the levées and the channel proper. A velocity profile in the channel shows a maximum value at the plug where the rate of shear is zero, and a velocity gradient increasing outward until, at the levées, the velocity becomes zero. Cooling produces a marked temperature gradient in the flow, leading to the formation of debris by brittle fracture when a critical value of shear rate to viscosity is reached. When the lava supply ceases, much of this debris and part of the lava is left behind after the flow nucleus drains out, forming a collapsed channel.Processes at the frontal zone include levée formation, debris formation, the change in shape of the front, and the choice of the flow path. These processes are controlled primarily by the rheological properties of the lava.Frontal zone dynamics can be understood by fixing the flow front as the point of reference. The lava flows through the channel into the front where it flows out into the levées, thereby increasing the length of the channel and permitting the front to advance. The front shows a relationship of critical height to the yield strength (τ₀) surface tension, and slope; its continued movement is activated by the pressure of the advancing lava in the channel behind. For an ideal flow (isothermal, homogeneous, and isotropic) the ratio of the section of channel proper to the section of levées is calculated and the distance the front will have moved at any time t_x can be determined once the amount of lava available to the front is known. Assuming that the velocity function of the front {G(t)} during the collapsing stage is proportional to the entrance pressure of the lava at the channel-front boundary, an exponential decrease of velocity through time is predicted, which shows good agreement with actual frontal velocity measurements taken on two flows. Local variations in slope have a secondary effect on frontal velocities.Under conditions of constant volume the frontal zone can be considered as a machine that consumes energy brought in by the lava to perform work (front advancement). While the front will use its potential energy to run the process, the velocity at which it occurs is controlled by the activation energy that enters the system as the kinetic energy of the lava flowing into the front. A relation for the energy contribution due to frontal acceleration is also derived. Finally the entrance pressure, that permits the front to deform, is calculated. Its small value confirms that the lava behaves very much like a Bingham plastic.     

Two-stage aquifer pumping subject to slow desorption and persistent sources

This work focuses on improving pump-and-treat remediation by optimizing a two-stage operational scheme to reduce volumes extracted when confronted with nonequilibrium desorption, low-permeability units, and continuous contaminant sources such as non-aqueous phase liquids (NAPL). Q1 and Q2 are the initial short-term high pumping rate and later long-term low pumping rate, respectively. A two-dimensional ground water flow and transport management model was used to test the proposed strategy for plumes developed from finite (NAPL-free) and continuous (NAPL-driven) contaminant sources in homogeneous and nonhomogeneous (zoned) aquifers. Remediation scenarios were simulated over durations of 2000, 6000, and 15,000 d to determine (1) the optimal time to switch from a preset Q1 to Q2 and (2) the value of Q2. The problem was constrained by mass removal requirements, maximum allowable downgradient concentrations, and practical bounds on Q2. Q1 was fixed at preset values 50% to 200% higher than the single-stage pumping rates (i.e., steady pumping rates during entire remediation period) necessary to achieve a desired cleanup level and capture the plume. Results for the NAPL-free homogeneous case under nonequilibrium desorption conditions achieved the same level of cleanup as single-stage pumping, while reducing extracted volumes by up to 36%. Comparable savings were obtained with NAPL-driven sources only when the source concentration was reduced by at least 2 orders of magnitude. For the zoned aquifer, the proposed strategy provided volume savings of up to 24% under NAPL-free and reduced source conditions.     

Runup uncertainty on planar beaches

Parameterization of wave runup is of paramount importance for an assessment of coastal hazards. Parametric models employ wave (e.g., H_s and L_p) and beach (i.e., β) parameters to estimate extreme runup (e.g., R_2%). Thus, recent studies have been devoted to improving such parameterizations by including additional information regarding wave forcing or beach morphology features. However, the effects of intra-wave dynamics, related to the random nature of the wave transformation process, on runup statistics have not been incorporated. This work employs a phase- and depth- resolving model, based on the Reynolds-averaged Navier-Stokes equations, to investigate different sources of variability associated with runup on planar beaches. The numerical model is validated with laboratory runup data. Subsequently, the role of both aleatory uncertainty and other known sources of runup variability (i.e., frequency spreading and bed roughness) is investigated. Model results show that aleatory uncertainty can be more important than the contributions from other sources of variability such as the bed roughness and frequency spreading. Ensemble results are employed to develop a new parametric model which uses the Hunt (J Waterw Port Coastal Ocean Eng 85:123–152, 1959) scaling parameter \(\beta \left (H_{s}L_{p}\right )^{1/2}\).