Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

The use of rogation ceremony records in climatic reconstruction: a case study from Catalonia (Spain)

The use of rogation ceremonies due to environmental causes constitutes an important source of information in paleoclimatic reconstructions. Their specific characteristics and full documental records permit highly reliable series to be reconstructed with daily, monthly, seasonal or annual resolution over periods of several centuries (3–4 centuries in the case of Catalonia). The levels of intensity, reflected in the type of religious ceremony enacted, allows quantification. Comparative analysis is made possible by the similarity of the mechanisms developed in different localities. The use of these series in paleoclimatological studies is a promising line of research, particularly as regards the pro pluvia rogations celebrated in the Mediterranean countries and in South America.     

Daily Meteorological Observations in Cádiz – San Fernando. Analysis of the Documentary Sources and the Instrumental Data Content (1786–1996)

Meteorological observations in the city of Cádiz are acknowledged as having been made from the middle of the 18th century onwards although they were only recorded and preserved in documentary form from 1789 onwards. Data readings were taken at the new Naval Observatory in San Fernando, ten kilometres from Cádiz, from 1797 onwards. Continuous series for temperature and atmospheric pressure at a daily resolution have been compiled and constructed from 1817 until 1996. The data series is composed of thrice daily observations made at Cádiz (1821–1880) by local observers and hourly data from Naval Observatory at San Fernando (1870–1996).     

Short‐ and long‐term evapotranspiration rates at ecological restoration sites along a large river receiving rare flow events

Many large rivers around the world no longer flow to their deltas, due to ever greater water withdrawals and diversions for human needs. However, the importance of riparian ecosystems is drawing increasing recognition, leading to the allocation of environmental flows to restore river processes. Accurate estimates of riparian plant evapotranspiration (ET) are needed to understand how the riverine system responds to these rare events and achieve the goals of environmental flows. In 2014, historic environmental flows were released into the Lower Colorado River at Morelos Dam (Mexico); this once perennial but now dry reach is the final stretch to the mighty Colorado River Delta. One of the primary goals was to supply native vegetation restoration sites along the reach with water to help seedlings establish and boost groundwater levels to foster the planted saplings. Patterns in ET before, during, and after the flows are useful for evaluating whether this goal was met and understanding the role that ET plays in this now ephemeral river system. Here, diurnal fluctuations in groundwater levels and Moderate Resolution Imaging Spectroradiometer (MODIS) data were used to compare estimates of ET specifically at 3 native vegetation restoration sites during 2014 planned flow events, and MODIS data were used to evaluate long‐term (2002–2016) ET responses to restoration efforts at these sites. Overall, ET was generally 0–10 mm d^?1 across sites, and although daily ET values from groundwater data were highly variable, weekly averaged estimates were highly correlated with MODIS‐derived estimates at most sites. The influence of the 2014 flow events was not immediately apparent in the results, although the process of clearing vegetation and planting native vegetation at the restoration sites was clearly visible in the results.     

Carbon fluxes through major phytoplankton groups during the spring bloom and post-bloom in the Northwestern Mediterranean Sea

The carbon flux through major phytoplankton groups, defined by their pigment markers, was estimated in two contrasting conditions of the Northwestern Mediterranean open ocean ecosystem: the spring bloom and post-bloom situations (hereafter Bloom and Post-bloom, respectively). During Bloom, surface chlorophyll a (Chl a) concentration was higher and dominated by diatoms (53% of Chl a), while during Post-bloom Synechococcus (42%) and Prymnesiophyceae (29%) became dominant. The seawater dilution technique, coupled to high pressure liquid chromatography (HPLC) analysis of pigments and flow cytometry (FCM), was used to estimate growth and grazing rates of major phytoplankton groups in surface waters. Estimated growth rates were corrected for photoacclimation based on FCM-detected changes in red fluorescence per cell. Given the 30% average decrease in the pigment content per cell between the beginning and the end of the incubations, overlooking photoacclimation would have resulted in a 0.40 d^?1 underestimation of phytoplankton growth rates. Corrected average growth rates (μ_o) were 0.90±0.20 (SD) and 0.40±0.14 d^?1 for Bloom and Post-bloom phytoplankton, respectively. Diatoms, Cryptophyceae and Synechococcus were identified as fast-growing groups and Prymnesiophyceae and Prasinophyceae as slow-growing groups across Bloom and Post-bloom conditions. The higher growth rate during Bloom was due to dominance of phytoplankton groups with higher growth rates than those dominating in Post-bloom. Average grazing rates (m) were 0.58±0.20 d^?1 (SD) and 0.31±0.07 d^?1. The proportion of phytoplankton growth consumed by microzooplankton grazing (m/μ_o) tended to be lower in Bloom (0.69±0.34) than in Post-bloom (0.80±0.08). The intensity of nutrient limitation experienced by phytoplankton indicated by μ_o/μ_n (where μ_n is the nutrient-amended growth rate), was similar during Bloom (0.78) and Post-bloom (0.73). Primary production from surface water (PP) was estimated with ¹⁴C incubations. A combination of PP and Chl a synthesis rate yielded C/Chl a ratios of 34±21 and 168±75 (g:g) for Bloom and Post-bloom, respectively. Transformation of group-specific Chl a fluxes into carbon equivalents confirmed the dominant role of diatoms during Bloom and Synechococcus and Prymnesiophyceae during Post-bloom.     

Influence of the Madden Julian Oscillation on precipitation and surface air temperature in South America

The regional influence of the Madden–Julian oscillation (MJO) on South America is described. Maps of probability of weekly-averaged rainfall exceeding the upper tercile were computed for all seasons and related statistically with the phase of the MJO as characterized by the Wheeler–Hendon real-time multivariate MJO (RMM) index and with the OLR MJO Index. The accompanying surface air temperature and circulation anomalies were also calculated. The influence of the MJO on regional scales along with their marked seasonal variations was documented. During December–February when the South American monsoon system is active, chances of enhanced rainfall are observed in southeastern South America (SESA) region mainly during RMM phases 3 and 4, accompanied by cold anomalies in the extratropics, while enhanced rainfall in the South Atlantic Convergence Zone (SACZ) region is observed in phases 8 and 1. The SESA (SACZ) signal is characterized by upper-level convergence (divergence) over tropical South America and a cyclonic (anticyclonic) anomaly near the southern tip of the continent. Impacts during March–May are similar, but attenuated in the extratropics. Conversely, in June–November, reduced rainfall and cold anomalies are observed near the coast of the SACZ region during phases 4 and 5, favored by upper-level convergence over tropical South America and an anticyclonic anomaly over southern South America. In September–November, enhanced rainfall and upper-level divergence are observed in the SACZ region during phases 7 and 8. These signals are generated primarily through the propagation of Rossby wave energy generated in the region of anomalous heating associated with the MJO.     

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.