Bayesian estimation of recent migration rates after a spatial expansion

Approximate Bayesian computation (ABC) is a highly flexible technique that allows the estimation of parameters under demographic models that are too complex to be handled by full-likelihood methods. We assess the utility of this method to estimate the parameters of range expansion in a two-dimensional stepping-stone model, using samples from either a single deme or multiple demes. A minor modification to the ABC procedure is introduced, which leads to an improvement in the accuracy of estimation. The method is then used to estimate the expansion time and migration rates for five natural common vole populations in Switzerland typed for a sex-linked marker and a nuclear marker. Estimates based on both markers suggest that expansion occurred < 10,000 years ago, after the most recent glaciation, and that migration rates are strongly male biased.

Cloud-resolving simulations of deep convection over a heated mountain

Cloud-resolving numerical simulations of airflow over a diurnally heated mountain ridge are conducted to explore the mechanisms and sensitivities of convective initiation under high pressure conditions. The simulations are based on a well-observed convection event from the Convective and Orographically Induced Precipitation Study (COPS) during summer 2007, where an isolated afternoon thunderstorm developed over the Black Forest mountains of central Europe, but they are idealized to facilitate understanding and reduce computational expense. In the conditionally unstable but strongly inhibited flow under consideration, sharp horizontal convergence over the mountain acts to locally weaken the inhibition and moisten the dry midtroposphere through shallow cumulus detrainment. The onset of deep convection occurs not through the deep ascent of a single updraft but rather through a rapid succession of thermals that are vented through the mountain convergence zone into the deepening cloud mass. Emerging thermals rise through the saturated wakes of their predecessors, which diminishes the suppressive effects of entrainment and allows for rapid glaciation above the freezing level as supercooled cloud drops rime onto preexisting ice particles. These effects strongly enhance the midlevel cloud buoyancy and enable rapid ascent to the tropopause. The existence and vigor of the convection is highly sensitive to small changes in background wind speed U0, which controls the strength of the mountain convergence and the ability of midlevel moisture to accumulate above the mountain. Whereas vigorous deep convection develops for U0 = 0 m s−1, deep convection is completely eliminated for U0 = 3 m s−1. Although deep convection is able to develop under intermediate winds (U0 = 1.5 m s−1), its formation is highly sensitive to small-amplitude perturbations in the initial flow.

Early–Middle Holocene vegetation history, climate change and human activities at Lago Riane (Ligurian Apennines, NW Italy)

The radiocarbon-dated palaeoecological study of Lago Riane (Ligurian Apennines, NW Italy) presented here forms part of a wider investigation into the relationships between Holocene vegetation succession, climate change and human activities in the northern Apennines. The record of vegetation history from Lago Riane indicates that, since the end of the last glaciation, climate change and prehistoric human activities, combined with several local factors, have strongly influenced the pattern and timing of natural vegetation succession. The pollen record indicates an important change in vegetation cover at Lago Riane at ~8500–8200 cal. years b.p., coincident with a well-known period of rapid climate change. At ~6100 cal. years b.p., Fagus woodland colonised Lago Riane during a period of climate change and expansion of Late Neolithic human activities in the upland zone of Liguria. A marked decline in Abies woodland, and the expansion of Fagus woodland, at ~4700 cal. years b.p., coincided with further archaeological evidence for pastoralism in the mountains of Liguria during the Copper Age. At ~3900–3600 cal. years b.p. (Early to Middle Bronze Age transition), a temporary expansion of woodland at Lago Riane has been provisionally attributed to a decline in human pressure on the environment during a period of short-term climate change

Dynamic Devensian ice flow in NE England: a sedimentological reconstruction

This study reconstructs the depositional environments that accompanied both ice advance and ice retreat of the last British–Irish Ice Sheet in NE England during the Last Glacial Maximum, and proposes three regional ice-flow phases. The Late Devensian (29–22 cal. ka BP) Tyne Gap Ice Stream initially deposited the Blackhall Till Formation during shelf-edge glaciation (Phase I). This subglacial traction till comprises several related facies, including stratified and laminated diamictons, tectonites, and sand and gravel beds deposited both in subglacial canals and in proglacial streams. Eventually, stagnation of the Tyne Gap Ice Stream led to ice-marginal sedimentation in County Durham (Phase II). During the Dimlington Stadial (21 cal. ka BP), the North Sea Lobe advanced towards the coastline of N Norfolk. This resulted initially in sandur deposition (widespread, tabular sand and gravel; the Peterlee Sand and Gravel Formation; Phase II) and ultimately in deposition of the Horden Till Formation (Phase III), a massive subglacial till. As the North Sea Lobe overrode previous formations, it thrusted and stacked sediments in County Durham, and dammed proglacial lakes between the east-coast ice, the Pennine uplands and the remaining Pennine ice. The North Sea Lobe retreated after Heinrich Event 1 (16 ka). This study highlights the complexity of ice flow during the Late Devensian glaciation of NE England, with changing environmental and oceanic conditions forcing a mobile and sensitive ice sheet.

Antarctic Peninsula ice sheet evolution during the Cenozoic Era

The Antarctic Peninsula region is currently undergoing rapid environmental change, resulting in the thinning, acceleration and recession of glaciers and the sequential collapse of ice shelves. It is important to view these changes in the context of long-term palaeoenvironmental complexity and to understand the key processes controlling ice sheet growth and recession. In addition, numerical ice sheet models require detailed geological data for tuning and testing. Therefore, this paper systematically and holistically reviews published geological evidence for Antarctic Peninsula Ice Sheet variability for each key locality throughout the Cenozoic, and brings together the prevailing consensus of the extent, character and behaviour of the glaciations of the Antarctic Peninsula region. Major contributions include a downloadable database of 186 terrestrial and marine calibrated dates; an original reconstruction of the LGM ice sheet; and a new series of isochrones detailing ice sheet retreat following the LGM. Glaciation of Antarctica was initiated around the Eocene/Oligocene transition in East Antarctica. Palaeogene records of Antarctic Peninsula glaciation are primarily restricted to King George Island, where glacigenic sediments provide a record of early East Antarctic glaciations, but with modification of far-travelled erratics by local South Shetland Island ice caps. Evidence for Neogene glaciation is derived primarily from King George Island and James Ross Island, where glaciovolcanic strata indicate that ice thicknesses reached 500–850 m during glacials. This suggests that the Antarctic Peninsula Ice Sheet draped, rather than drowned, the topography. Marine geophysical investigations indicate multiple ice sheet advances during this time. Seismic profiling of continental shelf-slope deposits indicates up to ten large advances of the Antarctic Peninsula Ice Sheet during the Early Pleistocene, when the ice sheet was dominated by 40 kyr cycles. Glacials became more pronounced, reaching the continental shelf edge, and of longer duration during the Middle Pleistocene. During the Late Pleistocene, repeated glacials reached the shelf edge, but ice shelves inhibited iceberg rafting. The Last Glacial Maximum (LGM) occurred at 18 ka BP, after which transitional glaciomarine sediments on the continental shelf indicate ice-sheet retreat. The continental shelf contains large bathymetric troughs, which were repeatedly occupied by large ice streams during Pleistocene glaciations. Retreat after the LGM was episodic in the Weddell Sea, with multiple readvances and changes in ice-flow direction, but rapid in the Bellingshausen Sea. The late Holocene Epoch was characterised by repeated fluctuations in palaeoenvironmental conditions, with associated glacial readvances. However, this has been subsumed by rapid warming and ice-shelf collapse during the twentieth century.

Climate at high-obliquity

The question of climate at high obliquity is raised in the context of both exoplanet studies (e.g. habitability) and paleoclimates studies (evidence for low-latitude glaciation during the Neoproterozoic and the ”Snowball Earth” hypothesis). States of high obliquity, φ, are distinctive in that, for φ ≥54◦, the poles receive more solar radiation in the annual mean than the Equator, opposite to the present day situation. In addition, the seasonal cycle of insolation is extreme, with the poles alternatively “facing” the sun and sheltering in the dark for months. The novelty of our approach is to consider the role of a dynamical ocean in controlling the surface climate at high obliquity, which in turn requires understanding of the surface winds patterns when temperature gradients are reversed. To address these questions, a coupled ocean-atmosphere-sea ice GCM configured on an aquaplanet is employed. Except for the absence of topography and modified obliquity, the set-up is Earth-like. Two large obliquities φ, 54◦ and 90◦, are compared to today’s Earth value, φ=23.5◦. Three key results emerge at high obliquity: 1) despite reversed temper- ature gradients, mid-latitudes surface winds are westerly and trade winds exist at the equator (as for φ=23.5◦) although the westerlies are confined to the summer hemisphere, 2) a habitable planet is possible with mid-latitude temperatures in the range 300-280 K and 3) a stable climate state with an ice cap limited to the equatorial region is unlikely. We clarify the dynamics behind these features (notably by an analysis of the potential vorticity structure and conditions for baroclinic instability of the atmosphere). Interestingly, we find that the absence of a stable partially glaciated state is critically linked to the absence of ocean heat transport during winter, a feature ultimately traced back to the high seasonality of baroclinic instability conditions in the atmosphere.

Late Quaternary lake-level record from northern Eurasia

Lake records from northern Eurasia show regionally coherent patterns of changes during the late Quaternary. Lakes peripheral to the Scandinavian ice sheet were lower than those today but lakes in the Mediterranean zone were high at the glacial maximum, reflecting the dominance of glacial anticyclonic conditions in northern Europe and a southward shift of the Westerlies. The influence of the glacial anticyclonic circulation attenuated through the late glacial period, and the Westerlies gradually shifted northward, such that drier conditions south of the ice sheet were confined to a progressively narrower zone and the Mediterranean became drier. The early Holocene shows a gradual shift to conditions wetter than present in central Asia, associated with the expanded Asian monsoon, and in the Mediterranean, in response to local, monsoon-type circulation. There is no evidence of mid-continental aridity in northern Eurasia during the mid-Holocene. In contrast, the circum-Baltic region was drier, reflecting the increased incidence of blocking anticyclones centered on Scandinavia in summer. There is a gradual transition to modern conditions after ca. 5000 yr B.P. Although these broad-scale patterns are interrupted by shorter term fluctuations, the long-term trends in lake behavior show a clear response to changes in insolation and glaciation.

Recent work at the Lower Palaeolithic site of Corfe Mullen, Dorset, England

Changes in the geological interpretation of the history of the ancient Solent river basin have focused attention on the handaxes discovered in the Corfe Mullen area during quarrying before the Second World War. Recent geological research suggests that the fluvial terrace the handaxes are associated with may pre-date the Anglian glaciation. This is important because it contributes to the question of just when the Solent basin was first occupied by hominins, and how this relates to other areas of possible contemporary pre-Anglian occupation such as the Boxgrove Marine embayment. However, the artefacts were believed to come from the bluff of the river terrace and were thus not in situ. This paper explores that question and re-examines the context from which the handaxes at Corfe Mullen were discovered.

Mixed-phase cloud physics and Southern Ocean cloud feedback in climate models

Increasing optical depth poleward of 45° is a robust response to warming in global climate models. Much of this cloud optical depth increase has been hypothesized to be due to transitions from ice-dominated to liquid-dominated mixed-phase cloud. In this study, the importance of liquid-ice partitioning for the optical depth feedback is quantified for 19 Coupled Model Intercomparison Project Phase 5 models. All models show a monotonic partitioning of ice and liquid as a function of temperature, but the temperature at which ice and liquid are equally mixed (the glaciation temperature) varies by as much as 40 K across models. Models that have a higher glaciation temperature are found to have a smaller climatological liquid water path (LWP) and condensed water path and experience a larger increase in LWP as the climate warms. The ice-liquid partitioning curve of each model may be used to calculate the response of LWP to warming. It is found that the repartitioning between ice and liquid in a warming climate contributes at least 20% to 80% of the increase in LWP as the climate warms, depending on model. Intermodel differences in the climatological partitioning between ice and liquid are estimated to contribute at least 20% to the intermodel spread in the high-latitude LWP response in the mixed-phase region poleward of 45°S. It is hypothesized that a more thorough evaluation and constraint of global climate model mixed-phase cloud parameterizations and validation of the total condensate and ice-liquid apportionment against observations will yield a substantial reduction in model uncertainty in the high-latitude cloud response to warming.

Fast or slow melting of the Marinoan snowball Earth? The cap dolostone record

The end of the Neoproterozoic era is punctuated by two global glacial events marked by the presence of glacial deposits overlaid by cap carbonates. Duration of glacial intervals is now consistently constrained to 3-12 million years but the duration of the post-glacial transition is more controversial due to the uncertainty in cap dolostone sedimentation rates. Indeed, the presence of several stratabound magnetic reversals in Brazilian cap dolostones recently questioned the short sedimentation duration (a few thousand years at most) that was initially suggested for these rocks. Here, we present new detailed magnetostratigraphic data of the Mirassol d`Oeste cap dolostones (Mato Grosso, Brazil) and ""bomb-spike"" calibrated AMS (14)C data of microbial mats from the Lagoa Vermelha (Rio de Janeiro, Brazil). We also compile sedimentary, isotopic and microbiological data from post-Marinoan outcrops and/or recent depositional analogues in order to discuss the deposition rate of Marinoan cap dolostones and to infer an estimation of the deglaciation duration in the snowball Earth aftermath. Taken together, the various data point to a sedimentation duration in the range of a few 10(5) years. (C) 2010 Elsevier B.V. All rights reserved.

A carbon isotope challenge to the snowball Earth

The snowball Earth hypothesis postulates that the planet was entirely covered by ice for millions of years in the Neoproterozoic era, in a self-enhanced glaciation caused by the high albedo of the ice-covered planet. In a hard-snowball picture, the subsequent rapid unfreezing resulted from an ultra-greenhouse event attributed to the buildup of volcanic carbon dioxide (CO(2)) during glaciation(1). High partial pressures of atmospheric CO(2) (p(CO2); from 20,000 to 90,000 p. p. m. v.) in the aftermath of the Marinoan glaciation (similar to 635 Myr ago) have been inferred from both boron and triple oxygen isotopes(2,3). These p(CO2) values are 50 to 225 times higher than present-day levels. Here, we re-evaluate these estimates using paired carbon isotopic data for carbonate layers that cap Neoproterozoic glacial deposits and are considered to record post-glacial sea level rise(1). The new data reported here for Brazilian cap carbonates, together with previous ones for time-equivalent units(4-8), provide p(CO2) estimates lower than 3,200 p. p. m. v.-and possibly as low as the current value of similar to 400 p. p. m. v. Our new constraint, and our reinterpretation of the boron and triple oxygen isotope data, provide a completely different picture of the late Neoproterozoic environment, with low atmospheric concentrations of carbon dioxide and oxygen that are inconsistent with a hard-snowball Earth.

PODOCARPUS SP.PL. REFUGIA IN BRAZIL`S ATLANTIC FOREST: ANALYSING THE PAST TO IMPROVE FUTURE PROTECTION

The Atlantic rainforest has the second highest biodiversity in Brazil. It has been shrinking rapidly in area as a result of intensive deforestation, and only 7% of the original cover now remains, as isolated patches or in ecological reserves. In order to obtain new information on the distribution of the Atlantic rainforest during the Quaternary, we examined herbarium data to locate relevant populations and extracted DNA from fresh leaves from 26 populations. The present-day distribution of endemic Podocarpus populations shows that they are widely dispersed across eastern Brazil, and that the expansion of Podocarpus recorded in single Amazonian pollen records may have originated from either western or eastern populations. Genetic analysis enabled us to determine the boundaries of their regional expansion: northern and central populations of P. sellowii appeared between 5 degrees and 15 degrees S some 16,000 years ago; populations of P lambertii or sellowii have appeared between 15 degrees and 23 degrees S at different times since the last glaciation at least; and P lambertii appeared between 23 degrees and 30 degrees S during the recent expansion of Araucaria forests. The combination of botanical, pollen, and molecular analyses proved to be a rapid means of inferring distribution boundaries for sparse populations and their regional evolution within tropical ecosystems. Today the rainforest refugia we identified have become hotspots that are crucial to the survival of the Atlantic forest under unfavourable climatic conditions and, as such, offer the only possible opportunity for this type of forest to expand in the event of future climate change.

Neoproterozoic glacial dynamics revealed by provenance of diamictites of the Bebedouro Formation, Sao Francisco Craton, Central Eastern Brazil

One of the main questions on Neoproterozoic geology regards the extent and dynamics of the glacial systems that are recorded in all continents. We present evidence for short transport distances and localized sediment sources for the Bebedouro Formation, which records Neoproterozoic glaciomarine sedimentation in the central-eastern Sao Francisco Craton (SFC), Brazil. New data are presented on clast composition, based on point counting in thin section and SHRIMP dating of pebbles and detrital zircon. Cluster analysis of clast compositional data revealed a pronounced spatial variability of clast composition on diamictite indicating the presence of individual glaciers or ice streams feeding the basin. Detrital zircon ages reveal distinct populations of Archean and Palaeoproterozoic age. The youngest detrital zircon dated at 874 +/- 9 Ma constrains the maximum depositional age of these diamictites. We interpret the provenance of the glacial diamictites to be restricted to sources inside the SFC, suggesting deposition in an environment similar to ice streams from modern, high latitude glaciers.

Neopaleozoic striated surfaces in the granitic basement and clast pavement remains of Salto, SP

Near Guarau Ceramic, localized southwest of Salto city in the State of Sao Paulo, two granite outcrops, distant some tens of meters from each other, display Neopaleozoic striated surfaces. These surfaces are in contact with diamictites from the Itarare Subgroup. The striae correspond to sub parallel grooves with millimetric spacing and depth, oriented about N48E and dipping 12 degrees to 42 degrees towards SE. Observed features and association with diamictites indicate an origin by glacial abrasion due to ice movement from southeast towards northwest. About 1.8 km east of Salto, unconsolidated material containing flat-iron-shaped and striated clasts was found on top of granite outcrops, interpreted as clasts pavement remains.

An early Permian brachiopod–gastropod fauna from the Calytrix Formation, Barbwire Terrace, Canning Basin, Western Australia

A small brachiopod–gastropod fauna from a core close to the base of the Calytrix Formation within the Grant Group includes the brachiopods Altiplecus decipiens (Hosking), Myodelthyrium dickinsi (Thomas), Brachythyrinella narsarhensis (Reed), Neochonetes (Sommeriella) obrieni Archbold, Tivertonia barbwirensis sp. nov. and the gastropod Peruvispira canningensis sp. nov. The fauna has affinities with that of the late Sakmarian‒early Artinskian Nura Nura Member directly overlying the Grant Group in other parts of the basin but, as with all lower Cisuralian (and Pennsylvanian) glacial strata in Western Australia, its precise age remains poorly constrained, especially in terms of correlation to international stages. Although the Calytrix fauna lies within the Pseudoreticulatispora confluens Palynozone, the only real constraint on its age (and that of the associated glacially influenced strata) is from Sakmarian (Sterlitamakian) and stratigraphically younger faunas. A brief review of radiometric ages from correlative strata elsewhere in Gondwana shows that those ages need to be updated. The presence of Asselian strata and the position of the Carboniferous‒Permian boundary remain unclear in Western Australia.Arturo César Taboada [ataboada@unpata.edu.ar], CONICET-Laboratorio de Investigaciones en Evolución y Biodiversidad (LIEB), Facultad de Ciencias Naturales, Sede Esquel, Universidad Nacional de la Patagonia ‘San Juan Bosco’, Edificio de Aulas, Ruta Nacional 259, km. 16,5, Esquel U9200, Chubut, Argentina; Arthur Mory [arthur.mory@dmp.wa.gov.au], Geological Survey of Western Australia, 100 Plain Street, East Perth, WA 6004, School of Earth and Environment, The University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia; Guang R. Shi [grshi@deakin.edu.au], School of Life and Environmental Sciences, Deakin University, Melbourne Burwood Campus, 221 Burwood Highway, Burwood, Victoria 3125, Australia; David W. Haig [david.haig@uwa.edu.au], School of Earth and Environment (M004), The University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia; María Karina Pinilla [mkpinilla@fcnym.unlp.edu.ar], División Paleozoología Invertebrados, Museo de Ciencias Naturales de La Plata, Paseo del Bosque s/n, 1900 La Plata, Buenos Aires, Argentina.