Long-term deforestation in NW Spain: linking the Holocene fire history to vegetation change and human activities

The Holocene fire regime is thought to have had a key role in deforestation and shrubland expansion in Galicia (NW Spain) but the contribution of past societies to vegetation burning remains poorly understood. This may be, in part, due to the fact that detailed fire records from areas in close proximity to archaeological sites are scarce. To fill this gap, we performed charcoal analysis in five colluvial soils from an archaeological area (Campo Lameiro) and compared the results to earlier studies from this area and palaeo-ecological literature from NW Spain. This analysis allowed for the reconstruction of the vegetation and fire dynamics in the area during the last ca 11 000 yrs. In the Early Holocene, Fabaceae and Betula sp. were dominant in the charcoal record. Quercus sp. started to replace these species around 10 000 cal BP, forming a deciduous forest that prevailed during the Holocene Thermal Maximum until ～5500 cal BP. Following that, several cycles of potentially fire-induced forest regression with subsequent incomplete recovery eventually led to the formation of an open landscape dominated by shrubs (Erica sp. and Fabaceae). Major episodes of forest regression were (1) ～5500–5000 cal BP, which marks the mid-Holocene cooling after the Holocene Thermal Maximum, but also the period during which agropastoral activities in NW Spain became widespread, and (2) ～2000–1500 cal BP, which corresponds roughly to the end of the Roman Warm Period and the transition from the Roman to the Germanic period. The low degree of chronological precision, which is inherent in fire history reconstructions from colluvial soils, made it impossible to distinguish climatic from human-induced fires. Nonetheless, the abundance of synanthropic pollen indicators (e.g. Plantago lanceolata and Urtica dioica) since at least ～6000 cal BP strongly suggests that humans used fire to generate and maintain pasture.     

Linear discriminant analysis to describe the relationship between rainfall and landslides in Bogotá, Colombia

Bogotá is located in the central Andean region of Colombia, which is frequently affected by landslide processes. These processes are mostly triggered during the rainy season in the city. This fact remarks the importance of determining what rain-derived parameters (e.g. intensity, antecedent rain, daily rain) are better related with the occurrence of landslides. For this purpose, the linear discriminant analysis (LDA), a technique derived from multivariate statistics, was used. The application of this type of analysis led to obtain simple mathematical functions that represent the probability of occurrence of landslides in Bogotá. The functions also allow to identify the most relevant variables derived from records of rainfall linked to the generation of landslides. A proof of concept using the proposed methodology was done using historic rainfall data from a 9-km² area of homogenous climatology and geomorphology in the south part of Bogotá. Landslides needed to be grouped for the LDA. Each one of these grouping categories represents landslides that occurred in similar geomorphologic conditions. Another set of events with no landslides was generated synthetically. Results of the proof of concept show that rainfall parameters such as normalized rainfall intensity I _MAP, normalized daily rainfall R _MAP and rainy-days normal RDN have the best statistical correlation with the landslides observed in the zone of analysis.     

Analysis of aquifer vulnerability and water quality using SINTACS and geographic weighted regression

Aquifer vulnerability and water quality were assessed in the Central Valleys of Oaxaca (Mexico) using the SINTACS method, based on a geographic information system. SINTACS layers were prepared using data such as climate (rainfall and temperature), water table, hydraulic conductivity, geology, soil type and topographic model. Maps for water quality index (WQI), contamination index and pollution sources index (PSI) were also obtained by this work. Groundwater quality in the Central Valleys may be affected by two factors, those with an anthropogenic origin and those with natural origin. High vulnerability values are located in the valleys of the basin, where granular sediments are exposed. Low vulnerability values are distributed in the basin??s ranges, where metamorphic rocks are found. Given that many of the zones with the highest groundwater vulnerability values correspond to zones with the greatest PSI values, there is great risk of groundwater contamination for the area of study because external (indicated by PSI) and internal (indicated by SINTACS) factors that cause pollution can be frequently observed in the same place. Geographic weighted regression (GWR) is used to test the dependency between WQI as dependent variable and SINTACS, PSI, Urban localities, Agriculture, Pastures and Rivers as predictors. The results indicate the non-stationary behavior of the dependent variable with respect to the predictors. While the obtained GWR models used to model WQI cannot be used in practical situations to predict the behavior of said variable, they can be used to estimate the degree to which the predictors influence the variable of interest.     

The effect of possible climate change on natural groundwater recharge based on a simple model: a study of four karstic aquifers in SE Spain

The study presented in this paper constitutes an initial approach to the problematic task of evaluating the effects of possible climate change on natural water recharge to aquifers. To estimate such effects, a purpose-designed mathematical model termed Estimation of Recharge in Over-exploited Aquifers (ERAS) has been used. It enables to simulate the monthly water recharge to an aquifer, provided that prior knowledge of the exploitation to which it is subjected and the variation caused by these two actions on the piezometric level of the aquifer is available. The basic data required for its application are: precipitation, temperature, groundwater extraction, stored groundwater surface and storage coefficient. The main advantage presented by this model is its independence of the mechanism by which water is displaced through the ground and within the unsaturated zone. The ERAS code was applied to four over-exploited karstic aquifers in Alto Vinalopó (Alicante, Spain) with the goal of generating a synthesized series of values for natural groundwater recharge in each of the aquifers for the 100 years of the twentieth century. Each series thus obtained after being grouped into decades was subjected to statistical processing, which revealed that in every case a logarithmically decreasing trend was present.     

Crossover from Nonlinearity Controlled to Heterogeneity Controlled Mixing in Two-Phase Porous Media Flows

We use high resolution Monte Carlo simulations to study the dispersive mixing in two-phase, immiscible, porous media flow that results from the interaction of the nonlinearities in the flow equations with geologic heterogeneity. Our numerical experiments show that distinct dispersive regimes occur depending on the relative strength of nonlinearity and heterogeneity. In particular, for a given degree of multiscale heterogeneity, controlled by the Hurst exponent which characterizes the underlying stochastic model for the heterogeneity, linear and nonlinear flows are essentially identical in their degree of dispersion, if the heterogeneity is strong enough. As the heterogeneity weakens, the dispersion rates cross over from those of linear heterogeneous flows to those typical of nonlinear homogeneous flows.     

Morphometric and molecular data on two mitochondrial genes of a newly discovered chimaeran fish (Hydrolagus melanophasma, Chondrichthyes)

Fresh and preserved (type material) specimens of the black ghost chimaera Hydrolagus melanophasma were compared for morphometric characteristics. A molecular comparison was also performed on two mitochondrial gene sequences (12S rRNA and 16S rRNA gene sequences). While significant differences in measurements were found, the differences were not attributable to sexual dimorphism or the quality of the specimens, but to the sample size and the type of statistical tests. The result of the genetic characterization showed that 12S rRNA and 16S rRNA genes represented robust molecular markers that characterized the species.     

Geochemical model of a magmatic–hydrothermal system at the Lastarria volcano, northern Chile

Lastarria volcano (25°10′ S, 68°31′ W; 5,697 m above sea level), located in the Central Andes Volcanic Zone (northern Chile), is characterized by four distinct fumarolic fields with outlet temperatures ranging between 80°C and 408°C as measured between May 2006–March 2008 and April–June 2009. Fumarolic gasses contain significant concentrations of high temperature gas compounds (i.e., SO₂, HCl, HF, H₂, and CO), and isotopic ratios (³He/⁴He, δ¹³C–CO₂, δ¹⁸O–H₂O, and δD–H₂O) diagnostic of magmatic gas sources. Gas equilibria systematics, in both the H₂O-H₂-CO₂-CO-CH₄ and alkane–alkene C₃ system, suggest that Lastarria fumarolic gasses emanate from a superheated vapor that is later cooled and condensed at relatively shallow depths. This two-stage process inhibits the formation of a continuous aquifer (e.g., horizontal liquid layer) at relatively shallow depth. Recent developments in the magmatic gas system may have enhanced the transfer and release of heat causing shallow aquifer vaporization. The consequent pressure increase and aquifer vaporization likely triggered the inflation events beginning in 2003 at the Lastarria volcano.     

Irregular Satellites of Jupiter: a study of the capture direction

The irregular satellites of Jupiter are believed to be captured asteroids or planetesimals. In the present work is studied the direction of capture of these objects as a function of their orbital inclination. We performed numerical simulations of the restricted three-body problem, Sun-Jupiter-particle, taking into account the growth of Jupiter. The integration was made backward in time. Initially, the particles have orbits as satellites of Jupiter, which has its present mass. Then, the system evolved with Jupiter losing mass and the satellites escaping from the planet. The reverse of the escape direction corresponds to the capture direction. The results show that the Lagrangian points L1 and L2 mainly guide the direction of capture. Prograde satellites are captured through these two gates with very narrow amplitude angles. In the case of retrograde satellites, these two gates are wider. The capture region increases as the orbital inclination increases. In the case of planar retrograde satellites the directions of capture cover the whole 360° around Jupiter. We also verified that prograde satellites are captured earlier in actual time than retrograde ones. This paper was presented at the Asteriods, Comets and Meteors meeting held at Búzios, Rio de Janeiro, Brazil in August 2005 and could not be included in the special issue related to that conference.     

Role of time variant and invariant terms on soil moisture spatio-temporal variability in a subhumid coastal wetland

The knowledge of soil moisture spatio-temporal variability is highly relevant for water resources management. This paper reports an analysis of the spatial–temporal variability of soil moisture data for a small to medium-scale soil-sensors network in a coastal wetland of southwestern Spain. Measurements were taken from five sites located in the Doñana National Park over the time-period of one hydrological year from September 2017 to September 2018. The total area of the soil-sensors network shows an extension about 25 × 3 km. Soil moisture data was separated into time invariant (the temporal mean of the whole period at each site) and time-variant terms (the deviations of soil moisture from the mean, or anomalies). The time-invariant component was generally the main contributor to the total spatial variance of soil moisture and it was mostly controlled by the groundwater levels in the area. Nevertheless, the time variant terms have a huge effect on soil moisture variability in very dry states. Characteristic convex time-dependent patterns for this field site were found between spatially averaged soil moisture and its variability. This information could be used for the up and downscaling of soil moisture from satellite data. Those patterns of relation between spatial mean and variability of soil moisture were only affected by heavy rainfalls giving rise to hysteretic behaviour. This study shows that even though groundwater level is a time-variant variable, it significantly affects soil moisture's time-variant but also time-invariant terms due to the different average groundwater level depths at the different sites.