Comparison of climate field reconstruction techniques: Application to Europe

This paper presents a comparison of principal component (PC) regression and regularized expectation maximization (RegEM) to reconstruct European summer and winter surface air temperature over the past millennium. Reconstruction is performed within a surrogate climate using the National Center for Atmospheric Research (NCAR) Climate System Model (CSM) 1.4 and the climate model ECHO-G 4, assuming different white and red noise scenarios to define the distortion of pseudoproxy series. We show how sensitivity tests lead to valuable “a priori” information that provides a basis for improving real world proxy reconstructions. Our results emphasize the need to carefully test and evaluate reconstruction techniques with respect to the temporal resolution and the spatial scale they are applied to. Furthermore, we demonstrate that uncertainties inherent to the predictand and predictor data have to be more rigorously taken into account. The comparison of the two statistical techniques, in the specific experimental setting presented here, indicates that more skilful results are achieved with RegEM as low frequency variability is better preserved. We further detect seasonal differences in reconstruction skill for the continental scale, as e.g. the target temperature average is more adequately reconstructed for summer than for winter. For the specific predictor network given in this paper, both techniques underestimate the target temperature variations to an increasing extent as more noise is added to the signal, albeit RegEM less than with PC regression. We conclude that climate field reconstruction techniques can be improved and need to be further optimized in future applications.

Finding Homogeneity in Heterogeneity - A New Approach to Quantifying Landscape Mosaics Developed for the Lao PDR

A key challenge for land change science in general and research on swidden agriculture in particular, is linking land cover information to human–environment interactions over larger spatial areas. In Lao PDR, a country facing rapid and multi-level land change processes, this hinders informed policy- and decision-making. Crucial information on land use types and people involved is still lacking. This article proposes an alternative approach for the description of landscape mosaics. Instead of analyzing local land use combinations, we studied land cover mosaics at a meso-level of spatial scale and interpreted these in terms of human–environmental interactions. These landscape mosaics were then overlaid with population census data. Results showed that swidden agricultural landscapes, involving 17% of the population, dominate 29% of the country, while permanent agricultural landscapes involve 74% of the population in 29% of the territory. Forests still form an important component of these landscape mosaics.

Arctic Winds in the “Twentieth Century Reanalysis”

Climate in the European part of the Arctic underwent a rapid warming between the 1910s and the 1930s. Previous studies have addressed the role of atmospheric circulation in this period based on geopotential height fields because observations of upper-level winds in the Arctic are rare. Here we analyse winds over the Arctic and specifically over Spitsbergen in the “Twentieth Century Reanalyses” (20CR). We compare in situ upper-air wind measurements performed in 1912 and 1913 in Spitsbergen with six-hourly 20CR data. Furthermore, we compare monthly-to-seasonal 20CR winds at 700 hPa over the European Arctic with statistically reconstructed winds at 3 km altitude. Finally, we analyse long-term trends in Arctic winds in 20CR. The general agreement between observed upper-air winds and 20CR on the day-to-day scale is rather poor, which is not surprising given the paucity of observations in the Arctic at that time that constrain 20CR. In contrast, the seasonally averaged winds (which represent a larger spatial scale) in 20CR compare well with statistically reconstructed winds. The analysis of long term near-surface wind time series in 20CR shows arguably artificial trends from 1871 to around the 1950s over sparsely observed regions, particularly oceanic regions. Densely observed regions such as Europe or the USA show no such trends. This analysis shows that great care needs to be taken when working with 20CR in the Arctic and other sparsely observed regions.

The heatwaves in Switzerland in summer 1947

In 1947, Switzerland was affected by a heat period of large spatial and temporal extent and rare occurrence. The heatwaves of 1947 can be compared with the events of 2003 in terms of intensity and duration. The summer of 1947 is studied based on the analysis of MeteoSwiss station data as well as the “Twentieth Century Reanalysis” (20CR) data set. Heatwaves were defined as six consecutive exceedances of the local 90th percentile of temperature. Five different heatwaves were identified which struck Switzerland during the summer of 1947. The most intense heatwave event is analysed in more detail. The meteorological situation was characterized by a high-pressure bridge over Central Europe. Based on a comparison with literature and with observations, the applicability of the 20CR dataset for the meteorological analysis of heatwave events could be demonstrated. The representation of the heat period in summer 1947 in 20CR is satisfactory when compared with station data, albeit with a temperature bias due to differences in topography. Hence, heatwaves cannot be defined using an absolute threshold. We conclude that 20CR is applicable for an overview of the meteorological patterns characterizing a heat wave but may not reproduce local details.

Optimizing the size of the area surveyed for monitoring a Eurasian lynx (Lynx lynx) population in the Swiss Alps by means of photographic capture-recapture

We studied the influence of surveyed area size on density estimates by means of camera-trapping in a low-density felid population (1-2 individuals/100 km(2) ). We applied non-spatial capture-recapture (CR) and spatial CR (SCR) models for Eurasian lynx during winter 2005/2006 in the northwestern Swiss Alps by sampling an area divided into 5 nested plots ranging from 65 to 760 km(2) . CR model density estimates (95% CI) for models M0 and Mh decreased from 2.61 (1.55-3.68) and 3.6 (1.62-5.57) independent lynx/100 km(2) , respectively, in the smallest to 1.20 (1.04-1.35) and 1.26 (0.89-1.63) independent lynx/100 km(2) , respectively, in the largest area surveyed. SCR model density estimates also decreased with increasing sampling area but not significantly. High individual range overlaps in relatively small areas (the edge effect) is the most plausible reason for this positive bias in the CR models. Our results confirm that SCR models are much more robust to changes in trap array size than CR models, thus avoiding overestimation of density in smaller areas. However, when a study is concerned with monitoring population changes, large spatial efforts (area surveyed ≥760 km(2) ) are required to obtain reliable and precise density estimates with these population densities and recapture rates.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Transport of 234U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of U-238, U-234 and Th-230 was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim of U-234 was to assess the in situ, long-term migration behaviour in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both U-238 and Th-230 immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of U-234 from the rock matrix to the pore water by the decay of U-238 will be established. Diffusive redistribution separates U-234 from its immobile parent U-238 resulting in bulk rock U-234/U-238 activity disequilibria. These may provide a means of estimating the mobility of U-234 in the rock if the diffusion rate of U-234 is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by beta-scanning, which shows that the (dominant) clay is richer in both elements. Samples were digested using aqua regia followed by total HF dissolution, yielding two fractions. in all studied samples U was found to be concentrated in the HF digestion fraction. It has a high U/Th ratio and a study by SEM-EDS points to sub-mu m up to several mu m in size zircon grains as the main U-rich phase. This fraction consistently has U-234/U-238 activity ratios below unity. The minute zircon grains constitute the major reservoir of U in the rock and act as constant rate suppliers of U-234 into the rock matrix and the pore water. The aqua regia leach fraction was found to be enriched in Th, and complementary to the HF fraction, having U-234/U-238 activity ratios above unity. It is believed that these U activity ratios reflect the surplus of having U-234 delivered from the zircon grains. Some cm-spaced samples show bulk rock U-234/U-238 activity ratios that are markedly out of equilibrium. In most of them a striking negative correlation between the total U content and the bulk rock U-234/U-238 activity ratios is observed. This is interpreted to indicate net U-234 transfer from regions of higher supply of U-234 towards those of lower supply which is, in most cases, equivalent to transfer from clayey towards carbonate/sandy portions of the rock. In contrast, the 25 cm averaged samples all have uniform bulk rock U-234/U-238 activity ratios in equilibrium, indicating U immobility in the last 1-1.5 Ma on this spatial scale. It is concluded that the small-scale lithological variations which govern U spatial distribution in the Opalinus Clay are the major factor determining U-234 in situ supply rates, regulating its diffusive fluxes and controlling the observed bulk rock U-234/U-238 activity ratios. A simple box-model is presented to simulate the measured bulk rock U-234/U-238 activity ratios and to give an additional insight into the studied system. (C) 2008 Elsevier Ltd. All rights reserved.

Pollen-based quantitative reconstructions of Holocene regional vegetation cover (plant-functional types and land-cover types) in Europe suitable for climate modelling

We present quantitative reconstructions of regional vegetation cover in north-western Europe, western Europe north of the Alps, and eastern Europe for five time windows in the Holocene around 6k, 3k, 0.5k, 0.2k, and 0.05k calendar years before present (bp)] at a 1 degrees x1 degrees spatial scale with the objective of producing vegetation descriptions suitable for climate modelling. The REVEALS model was applied on 636 pollen records from lakes and bogs to reconstruct the past cover of 25 plant taxa grouped into 10 plant-functional types and three land-cover types evergreen trees, summer-green (deciduous) trees, and open land]. The model corrects for some of the biases in pollen percentages by using pollen productivity estimates and fall speeds of pollen, and by applying simple but robust models of pollen dispersal and deposition. The emerging patterns of tree migration and deforestation between 6k bp and modern time in the REVEALS estimates agree with our general understanding of the vegetation history of Europe based on pollen percentages. However, the degree of anthropogenic deforestation (i.e. cover of cultivated and grazing land) at 3k, 0.5k, and 0.2k bp is significantly higher than deduced from pollen percentages. This is also the case at 6k in some parts of Europe, in particular Britain and Ireland. Furthermore, the relationship between summer-green and evergreen trees, and between individual tree taxa, differs significantly when expressed as pollen percentages or as REVEALS estimates of tree cover. For instance, when Pinus is dominant over Picea as pollen percentages, Picea is dominant over Pinus as REVEALS estimates. These differences play a major role in the reconstruction of European landscapes and for the study of land cover-climate interactions, biodiversity and human resources.

Gene flow and the maintenance of species boundaries

Hybrid zones are regions where individuals from genetically differentiated populations meet and mate, resulting in at least some offspring of mixed ancestry. Patterns of gene flow (introgression) in hybrid zones vary across the genome, allowing assessment of the role of individual genes or genome regions in reproductive isolation. Here, we document patterns of introgression between two recently diverged species of field crickets. We sampled at a very fine spatial scale and genotyped crickets for 110 highly differentiated single nucleotide polymorphisms (SNPs) identified through transcriptome scans. Using both genomic and geographic cline analysis, we document remarkably abrupt transitions (<100 m) in allele frequencies for 50 loci, despite high levels of gene flow at other loci. These are among the steepest clines documented for any hybridizing taxa. Furthermore, the cricket hybrid zone provides one of the clearest examples of the semi-permeability of species boundaries. Comparisons between data from the fine-scale transect and data (for the same set of markers) from sampling a much larger area in a different region of the cricket hybrid zone reveal consistent patterns of introgression for individual loci. The consistency in patterns of introgression between these two distant and distinct regions of the hybrid zone suggests that strong selection is acting to maintain abrupt discontinuities within the hybrid zone and that genomic regions with restricted introgression likely include genes that contribute to nonecological prezygotic barriers.

Age, Growth and Population Structure of Jumbo Flying Squid, Dosidicus Gigas, Based on Statolith Microstructure Off the Exclusive Economic Zone of Chilean Waters

The jumbo flying squid, Dosidicus gigas, support an important squid fishery off the Exclusive Economic Zone of Chilean waters. However, we only have limited information about their biology. In this study, age, growth and population structure of D. gigas were studied using statoliths from 333 specimens (386 females and 147 males) randomly sampled in the Chinese squid jigging surveys from 2007 to 2008 off the Exclusive Economic Zone of Chile. Mantle lengths (MLs) of the sample ranged from 206 to 702 mm, and their ages were estimated from 150 to 307 days for females and from 127 to 302 days for males. At least two spawning groups were identified, the main spawning peak tended to occur between August and November (austral spring group), and the secondary peak appeared during March to June (austral autumn group). The ML-age relationship was best modelled by a linear function for the austral spring group and a power function for the austral autumn group, and the body weight (BW)-age relationship was best described by an exponential function for both the groups. Instantaneous relative growth rates and absolute growth rates for ML and BW did not differ significantly between the two groups. The growth rate of D. gigas tended to be high at young stages, and then decreased after the sub-adult stage (>180 days old). This study suggests large spatial and temporal variability in key life history parameters of D. gigas, calling for the collection of more data with fine spatial and temporal scales to further improve our understanding of the fishery biology of D. gigas.

Extracting more out of relocation data: building movement models as mixtures of random walks

We present a framework for fitting multiple random walks to animal movement paths consisting of ordered sets of step lengths and turning angles. Each step and turn is assigned to one of a number of random walks, each characteristic of a different behavioral state. Behavioral state assignments may be inferred purely from movement data or may include the habitat type in which the animals are located. Switching between different behavioral states may be modeled explicitly using a state transition matrix estimated directly from data, or switching probabilities may take into account the proximity of animals to landscape features. Model fitting is undertaken within a Bayesian framework using the WinBUGS software. These methods allow for identification of different movement states using several properties of observed paths and lead naturally to the formulation of movement models. Analysis of relocation data from elk released in east-central Ontario, Canada, suggests a biphasic movement behavior: elk are either in an "encamped" state in which step lengths are small and turning angles are high, or in an "exploratory" state, in which daily step lengths are several kilometers and turning angles are small. Animals encamp in open habitat (agricultural fields and opened forest), but the exploratory state is not associated with any particular habitat type.

Daily polynya area and ice production in the North Water polynya region

The North Water (NOW) Polynya is a regularly-forming area of open-water and thin-ice, located between northwestern Greenland and Ellesmere Island (Canada) at the northern tip of Baffin Bay. Due to its large spatial extent, it is of high importance for a variety of physical and biological processes, especially in wintertime. Here, we present a long-term remote sensing study for the winter seasons 1978/1979 to 2014/2015. Polynya characteristics are inferred from (1) sea ice concentrations and brightness temperatures from passive microwave satellite sensors (Advanced Microwave Scanning Radiometer (AMSR-E and AMSR2), Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave Imager/Sounder (SSM/I-SSMIS)) and (2) thin-ice thickness distributions, which are calculated using MODIS ice-surface temperatures and European Center for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis data in a 1D thermodynamic energy-balance model. Daily ice production rates are retrieved for each winter season from 2002/2003 to 2014/2015, assuming that all heat loss at the ice surface is balanced by ice growth. Two different cloud-cover correction schemes are applied on daily polynya area and ice production values to account for cloud gaps in the MODIS composites. Our results indicate that the NOW polynya experienced significant seasonal changes over the last three decades considering the overall frequency of polynya occurrences, as well as their spatial extent. In the 1980s, there were prolonged periods of a more or less closed ice cover in northern Baffin Bay in winter. This changed towards an average opening on more than 85% of the days between November and March during the last decade. Noticeably, the sea ice cover in the NOW polynya region shows signs of a later-appearing fall freeze-up, starting in the late 1990s. Different methods to obtain daily polynya area using passive microwave AMSR-E/AMSR2 data and SSM/I-SSMIS data were applied. A comparison with MODIS data (thin-ice thickness < 20 cm) shows that the wintertime polynya area estimates derived by MODIS are about 30 to 40% higher than those derived using the polynya signature simulation method (PSSM) with AMSR-E data. In turn, the difference in polynya area between PSSM and a sea ice concentration (SIC) threshold of 70% is fairly low (approximately 10%) when applied to AMSR-E data. For the coarse-resolution SSM/I-SSMIS data, this difference is much larger, particularly in November and December. Instead of a sea ice concentration threshold, the PSSM method should be used for SSM/I-SSMIS data. Depending on the type of cloud-cover correction, the calculated ice production based on MODIS data reaches an average value of 264.4 ± 65.1 km**3 to 275.7 ± 67.4 km**3 (2002/2003 to 2014/2015) and shows a high interannual variability. Our achieved long-term results underline the major importance of the NOW polynya considering its influence on Arctic ice production and associated atmosphere/ocean processes.

Reefs shift from net accretion to net erosion along a natural environmental gradient

Coral reefs persist in an accretion-erosion balance and ocean acidification resulting from anthropogenic CO2 emissions threatens to shift this balance in favor of net reef erosion. Corals and calcifying algae, largely responsible for reef accretion, are vulnerable to environmental changes associated with ocean acidification, but the direct effects of lower pH on reef erosion has received less attention, particularly in the context of known drivers of bioerosion and natural variability. This study examines the balance between reef accretion and erosion along a well-characterized natural environmental gradient in Kane'ohe Bay, Hawai'i using experimental blocks of coral skeleton. Comparing before and after micro-computed tomography (µCT) scans to quantify net accretion and erosion, we show that, at the small spatial scale of this study (tens of meters), pH was a better predictor of the accretion-erosion balance than environmental drivers suggested by prior studies, including resource availability, temperature, distance from shore, or depth. In addition, this study highlights the fine-scale variation of pH in coastal systems and the importance of microhabitat variation for reef accretion and erosion processes. We demonstrate significant changes in both the mean and variance of pH on the order of meters, providing a local perspective on global increases in pCO2. Our findings suggest that increases in reef erosion, combined with expected decreases in calcification, will accelerate the shift of coral reefs to an erosion-dominated system in a high-CO2 world. This shift will make reefs increasingly susceptible to storm damage and sea-level rise, threatening the maintenance of the ecosystem services that coral reefs provide.

Indias urbanas en el Buenos Aires tardocolonial. Familia y grupos domésticos

Nos proponemos analizar la integración de las indias a la vida urbana de Buenos Aires en el período colonial tardío. Este estudio se circunscribe al análisis de la conformación de los grupos domésticos. Consideramos que ellos son de una importancia fundamental ya que constituyeron los ámbitos en los cuales se transmitió y transformó el orden social. Para introducirnos a esta cuestión tomaremos en primera instancia los empadronamientos de la ciudad de Buenos Aires de los años 1744 y 1778. Dejamos para más adelante el abordaje de esta cuestión a partir de otras fuentes y con una escala de análisis reducida

Indias urbanas en el Buenos Aires tardocolonial. Familia y grupos domésticos

Nos proponemos analizar la integración de las indias a la vida urbana de Buenos Aires en el período colonial tardío. Este estudio se circunscribe al análisis de la conformación de los grupos domésticos. Consideramos que ellos son de una importancia fundamental ya que constituyeron los ámbitos en los cuales se transmitió y transformó el orden social. Para introducirnos a esta cuestión tomaremos en primera instancia los empadronamientos de la ciudad de Buenos Aires de los años 1744 y 1778. Dejamos para más adelante el abordaje de esta cuestión a partir de otras fuentes y con una escala de análisis reducida