Paleoflood events recorded by speleothems in caves

Speleothems are usually composed of thin layers of calcite (or aragonite). However, cemented detrital materials interlayered between laminae of speleothemic carbonate have been also observed in many caves. Flowstones comprising discontinuous carbonate layers form due to flowing water films, while flood events introduce fluviokarstic sediments in caves that, on occasion, are recorded as clayey layers inside flowstones and stalagmites. This record provides a potential means of understanding the frequency of palaeofloods using cave records. In this work, we investigate the origin of this type of detrital deposit in El Soplao Cave (Northern Spain). The age of the lowest aragonite layer of a flowstone reveals that the earliest flood period occurred before 500 ka, though most of the flowstone formed between 422 +69/‐43 ka and 400 +66/‐42 ka. This suggests that the cave was periodically affected by palaeoflood events that introduced detrital sediments from the surface as a result of occasional extreme rainfall events, especially at around 400 ka. The mineralogical data enable an evolutionary model for this flowstone to be generated based on the alternation of flood events with laminar flows and carbonate layers precipitation that can be extrapolated to other caves in which detrital sediments inside speleothems have been found. Copyright © 2014 John Wiley & Sons, Ltd.     

Fire severity and surface rock fragments cause patchy distribution of soil water repellency and infiltration rates after burning

Although fire‐induced soil water repellency (SWR) and its effects on soil hydrology and geomorphology have been studied in detail, very few studies have considered the effect of rock fragments resting on the soil surface or partly embedded in soil. In this research, we have studied the effect of rock fragments on the strength and spatial distribution of fire‐induced SWR at different fire severities. A fire‐affected area was selected for this experiment and classified into different zones according to fire severity (unburned, low, moderate and high) and rock fragment cover (low, <20% and high, >60%). During 7 days after fire, SWR and infiltration rates were assessed in the soil surface covered by individual rock fragments and in the midpoint between two adjacent rock fragments (with maximum spacing of 20 cm). SWR increased with fire severity. Rock fragments resting on the soil surface increased the heterogeneity of the spatial distribution of fire‐induced SWR. SWR increased significantly with rock fragment cover in bare areas under moderate and high fire severity, but quantitatively important changes were only observed under high fire severity. In areas with a low rock fragment cover, water repellency from soil surfaces covered by rock fragments increased relative to bare soil surfaces, with increasing SWR. In areas with a high rock fragment cover, SWR increased significantly from non‐covered to covered soil surfaces only after low‐severity burning. Rock fragment cover did not affect infiltration rates, although it decreased significantly in soil surfaces after high‐severity burning in areas under low and high rock fragment cover. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessment of evaporation and water fluxes in a column of dry saline soil subject to different water table levels

Dry saline soils are common in the arid and hyper‐arid basins located in the Chilean Altiplano, where evaporation from shallow groundwater is typically the major component of the water balance. Thus, a good understanding of evaporation processes is necessary for improving water resource planning and management in these regions. In this study, we conducted laboratory experiments with a natural saline soil column to estimate evaporation rates and assess the liquid and water vapor fluxes under different water table levels. Water content, electrical conductivity and temperature at different depths were utilized to assess the liquid and water vapor fluxes in the soil column. We observed movement of water that dissolves salts from the soil and transports them to areas in the column where they accumulate. Isothermal liquid flux was predominant, while thermal and isothermal liquid and thermal water vapor fluxes were negligible, except for deep water table levels where isothermal and thermal water vapor fluxes had similar magnitude but opposite directions. Differences observed in total fluxes for all water table levels were due to different upward and downward fluxes, which depend on changes in water content and temperature within the soil profile. Both the vapor flux magnitude and direction were found to be very sensitive to the choice of empirical parameters used in flux quantification, such as tortuosity and the enhancement factor for local temperature gradients in the air phase within the column. Copyright © 2013 John Wiley & Sons, Ltd.     

Significance tests for covariate-dependent trends in inhomogeneous spatio-temporal point processes

Modeling and inference for spatial and spatio-temporal point processes is an issue that has been broadly investigated in the last years. Application fields such as forestry, epidemiology and ecology have been the main engine driving such raised interest. The inclusion of spatially varying covariates in the models for the intensity function is becoming of particular interest, but little attention has been paid to testing the significance of such covariates. Testing the significance of covariates is important if one seeks to explain which covariates have an effect in the spatial or spatio-temporal distribution of the point pattern observed. We thus provide practical procedures to build statistical tests of significance for covariates that have an effect on the intensity function of a point pattern. Our approximation focuses on the conditional intensity function, by considering nonparametric kernel-based estimators. We calculate thinning probabilities under the conditions of absence and presence of a covariate and compare them through divergence measures. Based on Monte Carlo experiments, we approximate the statistical properties of our tests under a variety of practical scenarios. An application on testing the significance of a covariate in a spatio-temporal data set on wildfires is also developed.     

A smoothed ANOVA model for multivariate ecological regression

Smoothed analysis of variance (SANOVA) has recently been proposed for carrying out disease mapping. The main advantage of this approach is its conceptual simplicity and ease of interpretation. Moreover, it allows us to fix the combination of diseases of particular interest in advance and to make specific inferences about them. In this paper we propose a reformulation of SANOVA in the context of ecological regression studies. This proposal considers the introduction in a non-parametric way of one (or several) covariate(s) into the model, explaining some pre-specified combinations of the outcome variables. In addition, random effects are also incorporated in order to model geographical variation in the combinations of outcome variables not explained by the covariate. Lastly, the model permits the decomposition of the variance in the set of outcome variables into different orthogonal components, quantifying the contribution of every one of them. The proposed model is applied to the geographical analysis of mortality due to malignant stomach neoplasm among women resident in the city of Barcelona (Spain). The available outcome variables are deaths grouped into two time periods, and a socioeconomic deprivation index is included as a covariate. The model has been implemented through INLA, a novel inference tool for Bayesian statistics.     

A physical interpretation of the deterministic fractal–multifractal method as a realization of a generalized multiplicative cascade

In this study, we attempt to offer a solid physical basis for the deterministic fractal–multifractal (FM) approach in geophysics (Puente, Phys Let A 161:441–447, 1992; J Hydrol 187:65–80, 1996). We show how the geometric construction of derived measures, as Platonic projections of fractal interpolating functions transforming multinomial multifractal measures, naturally defines a non-trivial cascade process that may be interpreted as a particular realization of a random multiplicative cascade. In such a light, we argue that the FM approach is as “physical” as any other phenomenological approach based on Richardson’s eddies splitting, which indeed lead to well-accepted models of the intermittencies of nature, as it happens, for instance, when rainfall is interpreted as a quasi-passive tracer in a turbulent flow. Although neither a fractal interpolating function nor the specific multipliers of a random multiplicative cascade can be measured physically, we show how a fractal transformation “cuts through” plausible scenarios to produce a suitable realization that reflects specific arrangements of energies (masses) as seen in nature. This explains why the FM approach properly captures the spectrum of singularities and other statistical features of given data sets. As the FM approach faithfully encodes data sets with compression ratios typically exceeding 100:1, such a property further enhances its “physical simplicity.” We also provide a connection between the FM approach and advection–diffusion processes.     

Prediction of PM10 and SO2 exceedances to control air pollution in the Bay of Algeciras,Spain

In this paper, the authors apply different classification techniques in order to provide 24 h advance forecasts of the daily peaks of SO₂ and PM10 concentrations in the Bay of Algeciras. K-nearest-neighbours, multilayer neural network with backpropagation and support vector machines (SVMs) are the classification methods used. The aim of this research is to obtain a suitable prediction model that would enable us to predict the peaks of pollutant concentrations in critical meteorological situations caused by the widespread existing industry and population in the area. A resampling strategy with twofold crossvalidation has been applied, using different quality indexes to evaluate the performance of the prediction models. SVM models achieved better true positive rate and accuracy (ACC) quality indexes. Results of ACC index value of 0.795 for PM10 and 0.755 for SO₂ showed the ability of the model to predict peaks and non-peaks correctly.     

Fully probabilistic seismic risk assessment considering local site effects for the portfolio of buildings in Medellín,Colombia

A fully probabilistic seismic risk analysis using a comprehensive approach is conducted for Medellin, the second largest city of Colombia, using a building by building database constructed and complemented from aerial images, considering characteristics such as building use categories, socio-economic levels and replacement values. The seismic hazard used for the analysis corresponds to the most updated study available in the country with the same model that was included in the national building code maps definition. Spectral transfer functions are determined for each of the seismic microzonation zones in order to take into account the dynamic soil response and amplification effects in the risk analysis. Several building types are defined for the city and individual vulnerability functions are assigned to each of them. Risk results are presented in the state of the art metrics such as the loss exceedance curve, probable maximum losses for different return periods, average annual losses and risk maps. The obtained results can be classified by use and socio-economic sectors as well as by structural systems that may help the stakeholders to identify where the risk concentrates.     

Modeling flocculation in a hypertidal estuary

When fine particles are involved, cohesive properties of sediment can result in flocculation and significantly complicate sediment process studies. We combine data from field observations and state-of-the-art modeling to investigate and predict flocculation processes within a hypertidal estuary. The study site is the Welsh Channel located at the entrance of the Dee Estuary in Liverpool Bay. Field data consist of measurements from a fixed site deployment during 12–22 February 2008. Grain size, suspended sediment volume concentration, and current velocity were obtained hourly from moored instruments at 1.5 m above bed. Near-bottom water samples taken every hour from a research vessel are used to convert volume concentrations to mass concentrations for the moored measurements. We use the hydrodynamic model Proudman Oceanographic Laboratory Coastal Ocean Modelling System (POLCOMS) coupled with the turbulence model General Ocean Turbulence Model (GOTM) and a sediment module to obtain three-dimensional distributions of suspended particulate matter (SPM). Flocculation is identified by changes in grain size. Small flocs were found during flood and ebb periods—and correlate with strong currents—due to breakup, while coarse flocs were present during slack waters because of aggregation. A fractal number of 2.4 is found for the study site. Turbulent stresses and particle settling velocities are estimated and are found to be related via an exponential function. The result is a simple semiempirical formulation for the fall velocity of the particles solely depending on turbulent stresses. The formula is implemented in the full three-dimensional model to represent changes in particle size due to flocculation processes. Predictions from the model are in agreement with observations for both settling velocity and SPM. The SPM fortnight variability was reproduced by the model and the concentration peaks are almost in phase with those from field data.     

Methods for estimating the velocities of the Brazil Current in the pre-salt reservoir area off southeast Brazil (23∘ S–26∘ S)

The Brazil Current (BC) is likely the least observed and investigated subtropical western boundary current in the world. This study proposes a simple and systematic methodology to estimate quasi-synoptic cross-sectional speeds of the BC within the Santos Basin (23^° S–26^° S) based on the dynamic method using several combinations of data: Conductivity, temperature, and depth (CTD), temperature profiles, CTD and vessel-mounted Acoustic Doppler Current Profiler (VMADCP), and temperature profiles and VMADCP. All of the geostrophic estimates agree well with lowered Acoustic Doppler Current Profiler (LADCP) velocity observations and yield volume transports of -5.56 ±1.31 and 2.50 ±1.01 Sv for the BC and the Intermediate Western Boundary Current (IWBC), respectively. The LADCP data revealed that the BC flows southwestward and is ～100 km wide, 500 m deep, and has a volume transport of approximately -5.75 ±1.53 Sv and a maximum speed of 0.59 m s^?1. Underneath the BC, the IWBC flows northeastward and has a vertical extent of approximately 1,300 m, a width of ～60 km, a maximum velocity of ～0.22 m s?¹, and a volume transport of 4.11 ± 2.01 Sv. Our analysis indicates that in the absence of the observed velocities, the isopycnal (σ ₀) of 26.82 kg m^?3 (～500 dbar) is an adequate level of no motion for use in geostrophic calculations. Additionally, a simple linear relationship between the temperature and the specific volume anomaly can be used for a reliable first estimate of the BC-IWBC system in temperature-only transects.