An Experimental Performance Evaluation of Spatio-Temporal Join Strategies

Many applications capture, or make use of, spatial data that changes over time. This requirement for effective and efficient spatio‐temporal data management has given rise to a range of research activities relating to spatio‐temporal data management. Such work has sought to understand, for example, the requirements of different categories of application, and the modelling facilities that are most effective for these applications. However, at present, there are few systems with fully integrated support for spatio‐temporal data, and thus developers must often construct custom solutions for their applications. Developers of both bespoke solutions and of generic spatio‐temporal platforms will often need to support the fusion of large spatio‐temporal data sets. Supporting such requests in a database setting involves the use of join operations with both spatial and temporal conditions – spatio‐temporal joins. However, there has been little work to date on spatio‐temporal join algorithms or their evaluation. This paper presents an evaluation of several approaches to the implementation of spatio‐temporal joins that build upon widely available indexing techniques. The evaluation explores how several algorithms perform for databases with different spatial and temporal characteristics, with a view to helping developers of generic infrastructures or custom solutions in the selection and development of appropriate spatio‐temporal join strategies.     

Empirical Orthogonal Functions Analysis Applied to the Inverse Problem in Hydrogeology: Evaluation of Uncertainty and Simulation of New Solutions

To fulfil the need to generate more realistic solutions, stochastic inverse simulations in hydrogeology are now constrained on both piezometric head and hydraulic conductivity data. These inverse techniques, often based on geostatistics, allow modifications of an initial solution conditioned only on hydraulic conductivity data to arrive at a final solution that also matches observed heads. By repeating the process as many times as necessary with different initial solutions, one generates an ensemble of final solutions thereby addressing the uncertainty of the inverse problem. This requires a method able to handle the whole ensemble and to work on its relevant characteristics. From this standpoint, the analysis by Empirical Orthogonal Functions (EOF) appears promising. The method builds an orthogonal decomposition of the covariance matrix, calculated over the whole set of solutions, and the areas in space where the first functions have a greater influence corresponding to locations of maximum uncertainty in the solutions. These locations depend both on the hydraulic characteristics of the flow problem and on the spatial distribution of available data. The EOF analysis is used on a synthetic problem that mimics a possible behavior of the Culebra aquifer of the Waste Isolation Pilot Plant (WIPP, New Mexico). The method also allows new solutions to be generated at lower computational cost by a random composition of the functions obtained by the EOF analysis. These new solutions keep the main characteristics of the initial ensemble and because they can be conditioned, they return very good results when they are used to solve the direct problem.     

Analytical Solutions for Water Flow and Solute Transport in the Unsaturated Zone

Analytical solutions for the water flow and solute transport equations in the unsaturated zone are presented. We use the Broadbridge and White nonlinear model to solve the Richards’ equation for vertical flow under a constant infiltration rate. Then we extend the water flow solution and develop an exact parametric solution for the advection-dispersion equation. The method of characteristics is adopted to determine the location of a solute front in the unsaturated zone. The dispersion component is incorporated into the final solution using a singular perturbation method. The formulation of the analytical solutions is simple, and a complete solution is generated without resorting to computationally demanding numerical schemes. Indeed, the simple analytical solutions can be used as tools to verify the accuracy of numerical models of water flow and solute transport. Comparison with a finite-element numerical solution indicates that a good match for the predicted water content is achieved when the mesh grid is one-fourth the capillary length scale of the porous medium. However, when numerically solving the solute transport equation at this level of discretization, numerical dispersion and spatial oscillations were significant.     

A paleointensity study of Cretaceous volcanic rocks from the Western Cordillera,Colombia

The Cretaceous Normal Superchron is a period of great interest to investigate global scale variations of the geomagnetic field. Long periods of single polarity are still a matter of debate: up to now there are two contradicting theories, which try to relate geomagnetic field intensity and reversal rate. We aim to shed light on the geomagnetic field strength during the Cretaceous Normal Superchron because data are still scarce and of dissimilar quality. To obtain reliable, absolute paleointensity determinations we investigate volcanic rocks from the Western Cordillera of Colombia. Several age determinations allow relating the samples to an age of about 92.5 Ma. To characterize the samples, we investigate rock magnetic properties and determine the characteristic remanent magnetization behavior. To determine paleointensities, we use a multimethod approach: first, we apply the classic Thellier-Coe protocol, and then, the relatively new multispecimen method. Rock magnetic measurements indicate magnetite as the main ferrimagnetic mineral, a stable magnetization revealed by reversible and nearly reversible thermomagnetic curves, and grain sizes that are either in the pseudosingle domain range or a mixture of single and multidomain grains. Alternating field and thermal demagnetization are rather complex, although we observe a few vector diagrams with a single, essentially uni-vectorial component with a small viscous overprint. Paleointensity determination with the Thellier-Coe protocol was unsuccessful, while with the multispecimen protocol we obtained four successful determinations out of 20. The failure of the Thellier-Coe protocol can be attributed to multidomain grains, which were observed during demagnetization and in rock magnetic experiments, and to the inhomogeneity of the volcanic rocks. Our multispecimen paleointensity determinations support low field strength at around 90 Ma during the Cretaceous Normal Superchron.     

Spatial and temporal changes in reef fish assemblages on disturbed coral reefs,north Pacific coast of Costa Rica

Benthic structure of coral reefs determines the availability of refuges and food sources. Therefore, structural changes caused by natural and anthropogenic disturbances can have negative impacts on reef‐associated communities. During the 1990s, coral reefs from Bahía Culebra were considered among the most diverse ecosystems along the Pacific coast of Costa Rica; however, recently they have undergone severe deterioration as consequence of chronic stressors such as El Niño‐Southern Oscillation and harmful algal blooms. Reef fish populations in this area have also been intensely exploited. This study compared reef fish assemblages during two periods (1995–1996 and 2014–2016), to determine whether they have experienced changes as a result of natural and anthropogenic disturbances. For both periods, benthic composition and reef fish abundance were recorded using underwater visual censuses. Live coral cover (LCC) decreased from 43.09 ± 18.65% in 1995–1996 to 1.25 ± 2.42% in 2014–2016 (U = 36, p < 0.05). Macroalgal cover (%) in 2014–2016 was sixfold higher than mean values reported for the Eastern Tropical Pacific region. Mean (±SD) fish species richness in 1995–1996 (36.67 ± 14.20) was higher than in 2014–2016 (23.00 ± 9.14; U = 20, p < 0.05). Over 40% of reef fish orders observed in 1995–1996 were not detected in the 2014–2016 surveys, including large‐bodied predators. Reduction in abundance of fish predators such as sharks, grunts, and snappers is likely attributed to changes in habitat structure. Herbivorous such as parrotfishes and pufferfishes increased their abundance at sites with low LCC, probably in response to predators decline and increased algal cover. These findings revealed significant degradation and drastic loss of structural complexity in coral reefs from Bahía Culebra, which now are dominated by macroalgae. The large reduction in structural complexity of coral reefs has resulted in the loss of diversity and key ecological roles (e.g., predation and herbivory), thus potentially reducing the resilience of the entire ecosystem.     

3D gravimetric investigation of the Cerro do Jarau structure,Rio Grande do Sul,Brazil

The Cerro do Jarau structure is possibly the third Brazilian basaltic crater formed in continental flood basalt of the Serra Geral Formation, Paraná Basin, a large igneous province (LIP) in southern Brazil. It is a nearly circular landform with a diameter of approximately 13 km that rises 200 m above the plains of the “pampas” in southern Brazil. In this work, Bouguer anomalies were calculated from gravity accelerations measured on the area of this structure. The residual Bouguer map shows a strong positive anomaly trending NE‐SW, located in the northeastern part of the structure, a feature not commonly associated with impact structures. However, the negative anomaly present in its center and the circular positive anomaly surrounding the central portion are typical of impact structures. The residual Bouguer anomaly varies from ?2 mGal to 8 mGal. The positive circular anomaly is not spatially coincident with the rim of the structure. Based on the interpretation of our gravimetric data, the estimated diameter of the structure is 12 km and the central portion has a diameter of approximately 5 km, both slightly smaller than previously suggested. The Bouguer anomaly map was inverted into a 3D density model using a constrained inversion method with a maximum density contrast of 0.5 g cm^?3. This model was interpreted to associate densities with rock types, resulting in a geological model. This geological model is in accordance with the meteorite impact nature of Cerro do Jarau.     

Open to closed system transition traced through the TDIC isotopic signature at the aquifer recharge stage, implications for groundwater C dating

¹⁴C dating models are limited when considering recent groundwater for which the carbon isotopic signature of the total dissolved inorganic carbon (TDIC) is mainly acquired in the unsaturated zone. Reducing the uncertainties of dating thus implies a better identification of the processes controlling the carbon isotopic composition of the TDIC during groundwater recharge. Geochemical interactions between gas, water and carbonates in the unsaturated zone were investigated for two aquifers (the carbonate-free Fontainebleau sands and carbonate-bearing Astian sands, France) in order to identify the respective roles of CO₂ and carbonates on the carbon isotopic signatures of the TDIC; this analysis is usually approached using open or closed system terms. Under fully open system conditions, the seasonality of the ¹³C values in the soil CO₂ can lead to important uncertainties regarding the so-called “initial ¹⁴C activity” used in ¹⁴C correction models. In a carbonate-bearing unsaturated zone such as in the Astian aquifer, we show that an approach based on fully open or closed system conditions is not appropriate. Although the chemical saturation between water and calcite occurs rapidly within the first metre of the unsaturated zone, the carbon isotopic contents (δ¹³C) of the CO₂ and the TDIC evolve downward, impacted by the dissolution-precipitation of the carbonates. In this study, we propose a numerical approach to describe this evolution. The δ¹³C and the A¹⁴C (radiocarbon activity) of the TDIC at the base of the carbonate-bearing unsaturated zone depends on (i) the δ¹³C and the A¹⁴C of the TDIC in the soil determined by the soil CO_2, (ii) the water’s residence time in the unsaturated zone and (iii) the carbonate precipitation-dissolution fluxes. In this type of situation, the carbonate δ¹³C-A¹⁴C evolutions indicate the presence of secondary calcite and permit the calculation of its accretion flux, equal to . More generally, for other sites under temperate climate and with similar properties to the Astian sands site, this approach allows for a reliable determination of the carbon isotopic composition at the base of the unsaturated zone as the indispensable “input function” data of the carbon cycle into the aquifer.     

Vegetation indices changes in the cloud forest of La Gomera Island (Canary Islands) and their hydrological implications

Vegetation indices derived from remote sensing data still remain to be used for analysing the relationship between climatic factors and vegetation seasonal phenology in middle latitudes with subtropical conditions forests such as the Canarian laurel forest. The Garajonay National Park, located in the La Gomera Island, protects one of the best preserved examples of the Macaronesian laurel forest, due to the cloud banks produced by trade winds, with rainfall and temperature field data showing a clear Mediterranean climatic pattern. We have analysed seasonal vegetation indices trend for different types of forest inside the Garajonay National Park using normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI) products derived from moderate resolution imaging spectrometer (MODIS) Aqua data for two hydrological years (October 2003 to September 2005) in relationship with the existing field climatic data: rainfall, net fog water and temperature. Maximum annual EVI maps show the highest vegetation indices in the laurel forest of La Gomera that occur during the dry season, mainly in late spring to early summer, with EVI temporal profiles showing that valley‐bottom laurel forest areas have the most clear seasonal trend. Difference maps of EVI values between months with the lowest and highest rainfall of each hydrological year clearly confirm the highest photosynthetic activity in the laurel forest during the dry season. In addition, these forests show a significative temporal correlation between EVI values and the temperature in the forest (p < 0·001). Our results prove the absence of summer drought stress in the laurel forest implying that the fog drip income is high enough to maintain enough soil moisture to allow the forest fully transpire when temperatures are higher. As the laurel forest of La Gomera occurs in the main recharge area of the island's aquifer system, our analysis of EVI data suggests that fog drip constitutes a key hydrological factor. Copyright © 2010 John Wiley & Sons, Ltd.