Microclimatic characterization of a karstic cave: human impact on microenvironmental parameters of a prehistoric rock art cave (Candamo Cave, northern Spain)

The Candamo Cave contains an important group of paleolithic paintings which have been seriously deteriorated due to mass tourism. In this work, an analysis was carried out of different climatic parameters (CO₂, temperature, humidity, ²²²Rn) during annual cycles with the cave closed to the public and during an experimental period of controlled visits. The effect of visits on the geochemical characteristics of karstic water was also analyzed together with the cave ventilation. The natural variations in the cave air CO₂ were above 3000 ppm, the increase produced through visits was only 100–110 ppm and since the humidity is almost permanently at saturation point, the critical parameter which limits the visitor capacity becomes air temperature. The temperature changes during the annual cycle are of the order of 1  °C in the external part and less than 0.5  °C in the internal part of the cave and a maximum increase of 0.13  °C was observed during the period of the visits. The ²²²Rn and CO₂ concentration minimums in the summer period (July–October) show that this is the most propitious time for visits, since the greatest ventilation is produced in the cave at this time and, therefore, the greatest capacity for recovery. The geochemistry of the water, on the other hand, indicated that this is the period of the year in which processes of wall corrosion can be most easily introduced, although this would be of limited magnitude. The visitor capacity calculated was 29 visitors/day. Received: 29 August 1996 · Accepted: 23 June 1997     

How accurately are climatological characteristics and surface water and energy balances represented for the Colombian Caribbean Catchment Basin?

In Colombia, the access to climate related observational data is restricted and their quantity is limited. But information about the current climate is fundamental for studies on present and future climate changes and their impacts. In this respect, this information is especially important over the Colombian Caribbean Catchment Basin (CCCB) that comprises over 80 % of the population of Colombia and produces about 85 % of its GDP. Consequently, an ensemble of several datasets has been evaluated and compared with respect to their capability to represent the climate over the CCCB. The comparison includes observations, reconstructed data (CPC, Delaware), reanalyses (ERA-40, NCEP/NCAR), and simulated data produced with the regional climate model REMO. The capabilities to represent the average annual state, the seasonal cycle, and the interannual variability are investigated. The analyses focus on surface air temperature and precipitation as well as on surface water and energy balances. On one hand the CCCB characteristics poses some difficulties to the datasets as the CCCB includes a mountainous region with three mountain ranges, where the dynamical core of models and model parameterizations can fail. On the other hand, it has the most dense network of stations, with the longest records, in the country. The results can be summarised as follows: all of the datasets demonstrate a cold bias in the average temperature of CCCB. However, the variability of the average temperature of CCCB is most poorly represented by the NCEP/NCAR dataset. The average precipitation in CCCB is overestimated by all datasets. For the ERA-40, NCEP/NCAR, and REMO datasets, the amplitude of the annual cycle is extremely high. The variability of the average precipitation in CCCB is better represented by the reconstructed data of CPC and Delaware, as well as by NCEP/NCAR. Regarding the capability to represent the spatial behaviour of CCCB, temperature is better represented by Delaware and REMO, while precipitation is better represented by Delaware. Among the three datasets that permit an analysis of surface water and energy balances (REMO, ERA-40, and NCEP/NCAR), REMO best demonstrates the closure property of the surface water balance within the basin, while NCEP/NCAR does not demonstrate this property well. The three datasets represent the energy balance fairly well, although some inconsistencies were found in the individual balance components for NCEP/NCAR.     

Evolution and modulation of tropical heating from the last glacial maximum through the twenty-first century

Twentieth century observations show that during the last 50?years the sea-surface temperature (SST) of the tropical oceans has increased by ~0.5°C and the area of SST >26.5 and 28°C (arbitrarily referred to as the oceanic warm pool: OWP) by 15 and 50% respectively in association with an increase in green house gas concentrations, with non-understood natural variability or a combination of both. Based on CMIP3 projections the OWP is projected to double during twenty-first century in a moderate CO₂ forcing scenario (IPCC A1B scenario). However, during the observational period the area of positive atmospheric heating (referred to as the dynamic warm pool, DWP), has remained constant. The threshold SST (T _H ), which demarks the region of net heating and cooling, has increased from 26.6°C in the 1950s to 27.1°C in the last decade and it is projected to increase to ~28.5°C by 2100. Based on climate model simulations, the area of the DWP is projected to remain constant during the twenty-first century. Analysis of the paleoclimate model intercomparison project (PMIP I and II) simulations for the Last Glacial maximum and the Mid-Holocene periods show a very similar behaviour, with a larger OWP in periods of elevated tropical SST, and an almost constant DWP associated with a varying T _H . The constancy of the DWP area, despite shifts in the background SST, is shown to be the result of a near exact matching between increases in the integrated convective heating within the DWP and the integrated radiative cooling outside the DWP as SST changes. Although the area of the DWP remains constant, the total tropical atmospheric heating is a strong function of the SST. For example the net heating has increased by about 10% from 1950 to 2000 and it is projected to increase by a further 20% by 2100. Such changes must be compensated by a more vigorous atmospheric circulation, with growth in convective heating within the warm pool, and an increase of subsiding air and stability outside the convective warm pool and an increase of vertical shear at the DWP boundaries. This finding is contrary to some conclusions from other studies but in accord with others. We discuss the similarities and differences at length.     

Ocean–atmosphere coupling and the boreal winter MJO

The influence of ocean–atmosphere coupling on the simulation and prediction of the boreal winter Madden–Julian Oscillation (MJO) is examined using the Seoul National University coupled general circulation model (CGCM) and atmospheric—only model (AGCM). The AGCM is forced with daily SSTs interpolated from pentad mean CGCM SSTs. Forecast skill is examined using serial extended simulations spanning 26 different winter seasons with 30-day forecasts commencing every 5 days providing a total of 598 30-day simulations. By comparing both sets of experiments, which share the same atmospheric components, the influence of coupled ocean–atmosphere processes on the simulation and prediction of MJO can be studied. The mean MJO intensity possesses more realistic amplitude in the CGCM than in AGCM. In general, the ocean–atmosphere coupling acts to improve the simulation of the spatio-temporal evolution of the eastward propagating MJO and the phase relationship between convection (OLR) and SST over the equatorial Indian Ocean and the western Pacific. Both the CGCM and observations exhibit a near-quadrature relationship between OLR and SST, with the former lagging by about two pentads. However, the AGCM shows a less realistic phase relationship. As the initial conditions are the same in both models, the additional forcing by SST anomalies in the CGCM extends the prediction skill beyond that of the AGCM. To test the applicability of the CGCM to real-time prediction, we compute the Real-time Multivariate MJO (RMM) index and compared it with the index computed from observations. RMM1 (RMM2) falls away rapidly to 0.5 after 17–18 (15–16) days in the AGCM and 18–19 (16–17) days in the CGCM. The prediction skill is phase dependent in both the CGCM and AGCM.     

Pleistocene and Holocene aeolian facies along the Huelva coast (southern Spain): climatic and neotectonic implications

The stratigraphic relationships, genesis and chronology, including radiocarbon dating, of the Quaternary sandy deposits forming the El Asperillo cliffs (Huelva) were studied with special emphasis on the influence of neotectonic activity, sea-level changes and climate upon the evolution of the coastal zone. The E-W trending normal fault of Torre del Loro separates two tectonic blocks. The oldest deposits occur in the upthrown block. They are Early to Middle Pleistocene fluviatile deposits, probably Late Pleistocene shallow-marine deposits along an E-W trending shoreline, and Late Pleistocene and Holocene aeolian sands deposited under prevailing southerly winds. Three Pleistocene and Holocene aeolian units accumulated in the downthrown block. Of these, Unit 1, is separated from the overlying Unit 2 by a supersurface that represents the end of the Last Interglacial. Accumulation of Unit 2 took place during the Last Glacial under more arid conditions than Unit 1. The supersurface separating Units 2 and 3 was formed between the Last Glacial maximum at 18 000 ¹⁴C yr BP and ca. 14 000 ¹⁴C yr BP, the latter age corresponding to an acceleration of the rise of sea level. Unit 3 records wet conditions. The supersurface separating Units 3 and 4 fossilised the fault and the two fault blocks. Units 4 (deposited before the 4th millennium BC), 5 (> 2700 ¹⁴C yr BP to 16th century) and 6 (16th century to present) record relatively arid conditions. Prevailing wind directions changed with time from W (Units 2–4) to WSW (Unit 5) and SW (Unit 6).     

High T/P evolution and metamorphic ages of the migmatitic basement of northern Sierras Pampeanas,Argentina: Characterization of a mid-crustal segment of the Famatinian belt

New petrologic, thermobarometric and U-Pb monazite geochronologic information allowed to resolve the metamorphic evolution of a high temperature mid-crustal segment of an ancient subduction-related orogen. The El Portezuelo Metamorphic-Igneous Complex, in the northern Sierras Pampeanas, is mainly composed of migmatites that evolved from amphibolite to granulite metamorphic facies, reaching thermal peak conditions of 670–820 °C and 4.5–5.3 kbar. The petrographic study combined with conventional and pseudosection thermobarometry led to deducing a short prograde metamorphic evolution within migmatite blocks. The garnet-absent migmatites represent amphibolite-facies rocks, whereas the cordierite-garnet-K-feldspar-sillimanite migmatites represent higher metamorphic grade rocks. U-Pb geochronology on monazite grains within leucosome record the time of migmatization between ≈477 and 470 Ma. Thus, the El Portezuelo Metamorphic-Igneous Complex is an example of exhumed Early Ordovician anatectic middle crust of the Famatinian mobile belt. Homogeneous exposure of similar paleo-depths throughout the Famatinian back-arc and isobaric cooling paths suggest slow exhumation and consequent longstanding crustal residence at high temperatures. High thermal gradients uniformly distributed in the Famatinian back-arc can be explained by shallow convection of a low-viscosity asthenosphere promoted by subducting-slab dehydration.     

Variability of extreme events in the Colombian Pacific and Caribbean catchment basins

This paper analyses the behavior of extreme events of surface precipitation and temperature inside the Pacific and Caribbean Catchment Basins in Colombia using several datasets such as observations, reconstructed data, NCEP-NCAR and ERA-40 reanalyses and data from the regional model REMO. We use an extreme value method that selects the time series excesses over a nonstationary threshold and adjusts them to a generalized Pareto distribution. The goodness of fit is evaluated through a test that includes the Cramer–von Mises, Kolmogorov–Smirnov and Anderson–Darling statistics and the p values generated by parametric bootstrap resampling. The test not only evaluates the goodness of fit but also the threshold choice. The parameters are presented in maps that allow recognition of the features of the extreme behaviour inside the catchment basins, and differences and similarities between them. Maps of return periods for the maximum extreme events are also presented. A strong influence of the El Niño–Southern oscillation on the extreme events of both temperature and precipitation is found in the two catchment basins.     

Changes in cloudiness over the Amazon rainforests during the last two decades: diagnostic and potential causes

This study shows a decrease of seasonal mean convection, cloudiness and an increase of surface shortwave down-welling radiation during 1984?C2007 over the Amazon rainforests based on the analysis of satellite-retrieved clouds and surface radiative flux data. These changes are consistent with an increase in surface temperature, increased atmospheric stability, and reduction of moisture transport to the Amazon based on in situ surface and upper air meteorological data and reanalysis data. These changes appear to link to the expansion of the western Pacific warm pool during the December?CFebruary season, to the positive phase of the Atlantic Multidecadal Oscillation and increase of SST over the eastern Pacific SST during the March?CMay season, and to an increase of the tropical Atlantic meridional SST gradient and an expansion of the western Pacific warm pool during September?CNovember season. The resultant increase of surface solar radiation during all but the dry season in the Amazon could contribute to the observed increases in rainforest growth during recent decades.     

Groundwater optimization model for sustainable management of the Valley of Puebla aquifer, Mexico

The Valley of Puebla aquifer (VPA), at the central region of Mexico, is subject to intensive exploitation to satisfy the urban and industrial demand in the region. As a result of this increased exploitation, a number of state and federal agencies in charge of water management are concerned about the problems associated with the aquifer (decline of groundwater table, deterioration in water quality, poor well productivity and increased pumping and water treatment costs). This study presents a groundwater management model that combines “MODFLOW” simulation with optimization tools “MODRSP”. This simulation–optimization model for groundwater evaluates a complex range of management options to identify the strategies that best fit the objectives for allocating resources in the VPA. Four hypothetical scenarios were defined to analyze the response of the hydrogeological system for future pumping schemes. Based on the simulation of flow with the MODFLOW program, promising results for the implementation of the optimization of water quantity were found in scenarios 3 and 4. However, upon comparison and analysis of the feasibility of recovery of the piezometric level (considering the policy of gradual reductions of pumping), scenario 4 was selected for optimization purposes. The response functions of scenario 4 were then obtained and optimized, establishing an extraction rate of 204.92 millions of m³/year (Mm³/year). The reduction in groundwater extraction will be possible by substituting the volume removed by 35 wells (that should be discontinued) by the same volume of water from another source.     

Bounding surface‐based modeling of compacted silty sand exhibiting suction dependent postpeak strain softening

This article focuses on modeling the strain hardening‐softening response of statically compacted silty sand as observed from a comprehensive series of suction‐controlled, consolidated‐drained triaxial tests accomplished in a fully automated, double‐walled triaxial test system via the axis‐translation technique. The constitutive model used in this work is based on the theory of Bounding Surface (BS) plasticity and is formulated within a critical state framework. The essential BS model parameters are calibrated using the full set of triaxial test results and then used for predictions of compacted silty sand response at matric suction states varying from 50 to 750 kPa. Complementary simulations using the Barcelona Basic Model have also been included, alongside BS model predictions, in order to get further enlightening insights into some of the main limitations and challenges facing both frameworks within the context of the experimental evidence resulting from the present research effort. In general, irrespective of the value of matric suction applied, the Barcelona Basic Model performs relatively well in predicting response at peak and critical state failure under low net confining pressure while the Bounding Surface Model performs relatively well under high net confining pressures.