San Quirce (Palencia,Spain): new chronologies for the Lower to Middle Palaeolithic transition of south-west Europe

San Quirce is an open-air archaeological site situated on a fluvial terrace in the Duero basin (Palencia, northern Iberia). This paper presents new and consistent chronologies obtained for the sedimentary sequence using post-infrared infrared stimulated luminescence (pIR-IR) dating of K-feldspars and single-grain thermally transferred optically stimulated luminescence (TT-OSL) dating of quartz. The new dating results indicate that the sequence is older than ~200 000 years and place San Quirce Level III within marine isotope stages (MIS) 8 and 7, between 274 ± 13 ka and 238 ± 13 ka. The main lithic assemblage at San Quirce comes from Level III. The predominant tool types found in this level are hammerstones, manuports and flakes, with a small proportion of cores and a significant presence of denticulates. Adaptation to local environmental conditions resulted in distinctive cultural habits, which were embedded in the cultural tradition of hominins occupying the site during the final third of the Middle Pleistocene. San Quirce preserves a simple cultural tradition that was employed by local hominins to engage in a diverse array of activities, and highlights the cultural diversity that appears to have been a characteristic feature of the Lower to Middle Palaeolithic transition 300–200 ka.     

Immiscibility between calciocarbonatitic and silicate melts and related wall rock reactions in the upper mantle: a natural case study from Romanian mantle xenoliths

This paper presents the textural, mineralogical and chemical study of veinlets cross-cutting peridotite xenoliths from the lithospheric mantle and brought to the surface by alkaline basalts (Persani Mountains, Romania). The veinlets utilized pre-existing zones of weakness in the host rocks or display a random distribution, lining grain boundaries or cross-cutting any mineral, and always forming an interconnected network. They are filled with carbonate patches included in a silicate matrix. Both products are holocrystalline. Carbonate products have alkali-poor calciocarbonatitic to sövitic compositions, while the silicate matrix composition ranges from monzodioritic to monzonitic and alkali feldspar syenitic, depending on the host-sample, i.e., within a rather alkaline silica-saturated series. The mineral phases present in the silicate matrix (F-apatite, armalcolite, chromite, diopside–enstatite series, plagioclase–sanidine series) are usually present in the carbonate zones, where forsterite is also found. Some minerals cross-cut the interface between both types of zones. Only the matrix is different, feldspathic (oligoclase to sanidine) in the former and pure calcite in the latter. Thus, mineralogical and textural relationships between both products are consistent with an origin with equilibrium liquid immiscibility. Mantle minerals cross-cut by veinlets are sometimes resorbed at grain boundaries, and at the contact of the most alkaline silicate and carbonate melts, subhedral diopside/augite formed at the expense of mantle enstatite or olivine. In terms of mineral chemistry, the compositional variations recorded by vein minerals vary along a continuous trend. They generally superpose to those observed from lherzolites to harzburgites, and exhibit the same range of composition as that observed between rims and cores of mantle minerals cross-cut by veinlets. In detail, the Ca-rich pyroxenes of veinlets are Al-poor and Mg-rich; cpx in the carbonate zones are slightly more Ca-rich than those in the silicate matrix; spinels are relatively Al- and Mg-poor but rather Cr- and Fe-rich. Existence of only one titanium oxide (armalcolite) and various pairs of pyroxenes suggest crystallization temperatures in the range 1100–1200°C and pressures between 10–15 kb. Feldspar compositions in silicate materials, which vary continuously from labradorite to sanidine, are consistent with hypersolvus and dry crystallization conditions. All of these results provide evidence that immiscibility occurred at mantle depth as the liquid was forcibly injected during hydraulic fracturing of the mantle. The compositions of conjugate melts suggest a very large miscibility gap, as expected at high pressure in a dry environment from the experiments of Kjarsgaard and Hamilton [Kjarsgaard, B.A., Hamilton, D.L., 1988. Liquid immiscibility and the origin of alkali-poor carbonatites. Mineral. Mag. 52, 43–55; Kjarsgaard, B.A., Hamilton, D.L., 1989. The genesis of carbonatites by immiscibility. In: Bell, K. (Ed.), Carbonatites: Genesis and Evolution. Unwyn Hyman, London, pp. 388–404.]. The parental melt was carbonate, silica-undersaturated and rich in F, Cl and CO₂. Both immiscible melts were water-undersaturated. The cooling rate until total crystallization in veinlets was very slow, limited and necessarily occurred at mantle depth. Wall rock reactions leading to the formation of Ca-rich pyroxene at the expense of mantle enstatite or olivine occurred only at the contact with somewhat alkali-rich carbonatitic or silicate melts. Calcite, always anhedral, is the last mineral to crystallize. It is a differentiation product formed by magmatic crystallization or wall rock reaction. In some cases, given the rarity of any other minerals, it may be the product of the crystallization of a pure sövite immiscible melt.     

On Monin–Obukhov Scaling in and Above the Atmospheric Surface Layer: The Complexities of Elevated Scintillometer Measurements

In scintillometry Monin–Obukhov similarity theory (MOST) is used to calculate the surface sensible heat flux from the structure parameter of temperature (C_T²){(C_{T^2})} . In order to prevent saturation a scintillometer can be installed at an elevated level. However, in that case the observation level might be located outside the atmospheric surface layer (ASL) and thus the validity of MOST questioned. Therefore, we examine two concepts to determine the turbulent surface sensible heat flux from the structure parameter at elevated levels with data obtained at 60-m height on the Cabauw tower (the Netherlands). In the first concept (MOSTs) C_T²{C_{T^2}} is still scaled with the surface flux, whereas in the second (MOSTl) C_T²{C_{T^2}} is scaled with the local sensible heat flux. The C_T²{C_{T^2}} obtained from both concepts is compared with direct observations of C_T²{C_{T^2}} using a sonic anemometer/thermometer. In the afternoon (when the measurement height is located within the ASL) both concepts give results that are comparable to the directly observed values of C_T²{C_{T^2}} . In the morning (data outside the ASL), our data do not unequivocally support either of the two concepts. First, the peak in C_T²{C_{T^2}} that occurs when the measurement height is located in the entrainment zone disqualifies the use of MOST. Second, during the morning transition, local scaling shows the correct pattern (zero flux and a minimum in C_T²{C_{T^2}}) but underestimates C_T²{C_{T^2}} by a factor of ten. Third, from the best linear fit a we found that the slope of MOSTl gave better results, whereas the offset is closer to zero for MOSTs. Further, the correlation between the direct observations and MOST-scaled results is low and similar for the two concepts. In the end, we conclude that MOST is not applicable for the morning hours when the observation level is above the ASL.     

A comparison of the climatology of a tropospherestratosphere-mesosphere model with observations

Mean fields from a perpetual January simulation of a GCM extending from the surface to 0.01 hPa (near 80 km) are compared to observations. The zonal mean temperature and wind fields correspond quite well with reality; the low stratosphere, especially in the polar night, is too cold, but warmer than in the original version of the model, with an upper boundary at 25 hPa. Mean fields at standard levels show that the major features of the troposphere are represented by the model, but rather over emphasised; the stratospheric winter polar vortex is too strong, too cold, and too barotropic; it resembles an `undisturbed' January rather than the climatology. Differences in the stationary eddy activity between the extended and orginal versions of the model are noted, and used to explain some differences between the two simulations.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute für Meteorology. Guest Editor for these papers is Dr. L. Dümenil     

Climate change triggered sedimentation and progressive tectonic uplift in a coupled piedmont–axial system: Cuyama Valley,California, USA

Channels on the north‐facing piedmont of the Sierra Madre range in Cuyama Valley, California have alternated between three process regimes during the late Quaternary: (1) vertical incision into piedmont alluvium and older sedimentary deposits; (2) lateral erosion; and (3) sediment accumulation. The state of the piedmont system at a given time has been controlled by upstream sediment flux, regional tectonic uplift and incision of the axial Cuyama River. To better understand the timing and to attempt to interpret causes of past geomorphological processes on the Sierra Madre piedmont, we mapped the surficial geology and dated alluvial deposits using radiocarbon, cosmogenic and optical dating methods. Four primary episodes of sedimentation have occurred since ca. 100 ka, culminating in the most recent period of extensive piedmont sedimentation between 30 and 20 ka. Fill terraces in Cuyama Valley formed by piedmont sediment accumulation followed by vertical incision and lateral erosion are fairly planar and often mantle strath bedrock surfaces. Their vertical spatial arrangement is a record of progressive regional tectonic uplift and concomitant axial Cuyama River channel incision migrating up tributary piedmont channels. Subparallel longitudinal terrace profiles which have a linear age–elevation relationship indicate that multiple episodes of climatically controlled sedimentation overprints ～1 m kyr^?1 of regional uplift affecting the Cuyama River and its tributaries. Sedimentation was probably a result of increased precipitation that caused saturation landsliding in steep catchments. It is possible that increased precipitation during the Last Glacial Maximum was caused by both continental‐scale circulation pattern reorganization and increased Pacific storm frequency and intensity caused by ‘early warming’ of nearby Pacific Ocean surface waters. Older episodes of piedmont sedimentation are difficult to correlate with specific climate regimes, but may correlate with previous periods of increased precipitation. Copyright © 2007 John Wiley & Sons, Ltd.     

Inter-annual variability of wet season freshwater plume extent into the Great Barrier Reef lagoon based on satellite coastal ocean colour observations

Riverine freshwater plumes are the major transport mechanism for nutrients, sediments and pollutants into the Great Barrier Reef (GBR) lagoon and connect the land with the receiving coastal and marine waters. Knowledge of the variability of the freshwater extent into the GBR lagoon is relevant for marine park management to develop strategies for improving ecosystem health and risk assessments. In this study, freshwater extent has been estimated for the entire GBR lagoon area from daily satellite observations of the Moderate Resolution Imaging Spectroradiometer (MODIS) between 2002 and 2010. To enable a reliable mapping of freshwater plumes we applied a physics-based coastal ocean colour algorithm, that simultaneously retrieves chlorophyll-a, non-algal particulate matter and coloured dissolved organic matter (CDOM), from which we used CDOM as a surrogate for salinity (S) for mapping the freshwater extent.     

Scaling Variances of Scalars in a Convective Boundary Layer Under Different Entrainment Regimes

For the presentation and analysis of atmospheric boundary-layer (ABL) data, scales are used to non-dimensionalise the observed quantities and independent variables. Usually, the ABL height, surface sensible heat flux and surface scalar flux are used. This works well, so long as the absolute values of the entrainment ratio for both the scalar and temperature are similar. The entrainment ratio for temperature naturally ranges from −0.4 to −0.1. However, the entrainment ratio for passive scalars can vary widely in magnitude and sign. Then the entrainment flux becomes relevant as well. The only customary scalar scale that takes into account both the surface flux and the entrainment flux is the bulk scalar scale, but this scale is not well-behaved for large negative entrainment ratios and for an entrainment ratio equal to −1. We derive a new scalar scale, using previously published large-eddy simulation results for the convective ABL. The scale is derived under the constraint that scaled scalar variance profiles are similar at those heights where the variance producing mechanisms are identical (i.e., either near the entrainment layer or near the surface). The new scale takes into account that scalar variance in the ABL is not only related to the surface flux of that scalar, but to the scalar entrainment flux as well. Furthermore, it takes into account that the production of variance by the entrainment flux is an order of magnitude larger than the production of variance by the surface flux (per unit flux). Other desirable features of the new scale are that it is always positive (which is relevant when scaling standard deviations) and that the scaled variances are always of order 1–10.     

Simultaneous perturbation stochastic approximation for tidal models

The Dutch continental shelf model (DCSM) is a shallow sea model of entire continental shelf which is used operationally in the Netherlands to forecast the storm surges in the North Sea. The forecasts are necessary to support the decision of the timely closure of the moveable storm surge barriers to protect the land. In this study, an automated model calibration method, simultaneous perturbation stochastic approximation (SPSA) is implemented for tidal calibration of the DCSM. The method uses objective function evaluations to obtain the gradient approximations. The gradient approximation for the central difference method uses only two objective function evaluation independent of the number of parameters being optimized. The calibration parameter in this study is the model bathymetry. A number of calibration experiments is performed. The effectiveness of the algorithm is evaluated in terms of the accuracy of the final results as well as the computational costs required to produce these results. In doing so, comparison is made with a traditional steepest descent method and also with a newly developed proper orthogonal decomposition-based calibration method. The main findings are: (1) The SPSA method gives comparable results to steepest descent method with little computational cost. (2) The SPSA method with little computational cost can be used to estimate large number of parameters.     

Simultaneous calibration of surface flow and baseflow simulations: a revisit of the SWAT model calibration framework

Accurate analysis of water flow pathways from rainfall to streams is critical for simulating water use, climate change impact, and contaminants transport. In this study, we developed a new scheme to simultaneously calibrate surface flow (SF) and baseflow (BF) simulations of soil and water assessment tool (SWAT) by combing evolutionary multi‐objective optimization (EMO) and BF separation techniques. The application of this scheme demonstrated pronounced trade‐off of SWAT's performance on SF and BF simulations. The simulated major water fluxes and storages variables (e.g. soil moisture, evapotranspiration, and groundwater) using the multiple parameters from EMO span wide ranges. Uncertainty analysis was conducted by Bayesian model averaging of the Pareto optimal solutions. The 90% confidence interval (CI) estimated using all streamflows substantially overestimate the uncertainty of low flows on BF days while underestimating the uncertainty of high flows on SF days. Despite using statistical criteria calculated based on streamflow for model selection, it is important to conduct diagnostic analysis of the agreement of SWAT behaviour and actual watershed dynamics. The new calibration technique can serve as a useful tool to explore the trade‐off between SF and BF simulations and provide candidates for further diagnostic assessment and model identification. Copyright © 2011 John Wiley & Sons, Ltd.     

Acoustic tomography in the atmospheric surface layer

Acoustic tomography is presented as a technique for remote monitoring of meteorological quantities. This method and a special algorithm of analysis can directly produce area-averaged values of meteorological parameters. As a result consistent data will be obtained for validation of numerical atmospheric micro-scale models. Such a measuring system can complement conventional point measurements over different surfaces. The procedure of acoustic tomography uses the horizontal propagation of sound waves in the atmospheric surface layer. Therefore, to provide a general overview of sound propagation under various atmospheric conditions a two-dimensional ray-tracing model according to a modified version of Snell’s law is used. The state of the crossed atmosphere can be estimated from measurements of acoustic travel time between sources and receivers at different points. Derivation of area-averaged values of the sound speed and furthermore of air temperature results from the inversion of travel time values for all acoustic paths. Thereby, the applied straight ray two-dimensional tomographic model using SIRT (simultaneous iterative reconstruction technique) is characterised as a method with small computational requirements, satisfactory convergence and stability properties as well as simple handling, especially, during online evaluation.