Understanding advances in the simulation of the Madden-Julian Oscillation in a numerical weather prediction model.

The Madden-Julian Oscillation (MJO) is the dominant mode of intraseasonal variability in the Trop- ics. It can be characterised as a planetary-scale coupling between the atmospheric circulation and organised deep convection that propagates east through the equatorial Indo-Pacific region. The MJO interacts with weather and climate systems on a near-global scale and is a crucial source of predictability for weather forecasts on medium to seasonal timescales. Despite its global signifi- cance, accurately representing the MJO in numerical weather prediction (NWP) and climate models remains a challenge. This thesis focuses on the representation of the MJO in the Integrated Forecasting System (IFS) at the European Centre for Medium-Range Weather Forecasting (ECMWF), a state-of-the-art NWP model. Recent modifications to the model physics in Cycle 32r3 (Cy32r3) of the IFS led to ad- vances in the simulation of the MJO; for the first time the observed amplitude of the MJO was maintained throughout the integration period. A set of hindcast experiments, which differ only in their formulation of convection, have been performed between May 2008 and April 2009 to asses the sensitivity of MJO simulation in the IFS to the Cy32r3 convective parameterization. Unique to this thesis is the attribution of the advances in MJO simulation in Cy32r3 to the mod- ified convective parameterization, specifically, the relative-humidity-dependent formulation for or- ganised deep entrainment. Increasing the sensitivity of the deep convection scheme to environmen- tal moisture is shown to modify the relationship between precipitation and moisture in the model. Through dry-air entrainment, convective plumes ascending in low-humidity environments terminate lower in the atmosphere. As a result, there is an increase in the occurrence of cumulus congestus, which acts to moisten the mid-troposphere. Due to the modified precipitation-moisture relationship more moisture is able to build up which effectively preconditions the tropical atmosphere for the transition to deep convection. Results from this thesis suggest that a tropospheric moisture control on convection is key to simulating the interaction between the physics and large-scale circulation associated with the MJO.

On the ability of statistical wind-wave models to capture the variability and long-term trends of the North Atlantic winter wave climate

A dynamical wind-wave climate simulation covering the North Atlantic Ocean and spanning the whole 21st century under the A1B scenario has been compared with a set of statistical projections using atmospheric variables or large scale climate indices as predictors. As a first step, the performance of all statistical models has been evaluated for the present-day climate; namely they have been compared with a dynamical wind-wave hindcast in terms of winter Significant Wave Height (SWH) trends and variance as well as with altimetry data. For the projections, it has been found that statistical models that use wind speed as independent variable predictor are able to capture a larger fraction of the winter SWH inter-annual variability (68% on average) and of the long term changes projected by the dynamical simulation. Conversely, regression models using climate indices, sea level pressure and/or pressure gradient as predictors, account for a smaller SWH variance (from 2.8% to 33%) and do not reproduce the dynamically projected long term trends over the North Atlantic. Investigating the wind-sea and swell components separately, we have found that the combination of two regression models, one for wind-sea waves and another one for the swell component, can improve significantly the wave field projections obtained from single regression models over the North Atlantic.

Decadal predictability of regional scale wind speed and wind energy potentials over Central Europe

Decadal predictions on timescales from one year to one decade are gaining importance since this time frame falls within the planning horizon of politics, economy and society. The present study examines the decadal predictability of regional wind speed and wind energy potentials in three generations of the MiKlip (‘Mittelfristige Klimaprognosen’) decadal prediction system. The system is based on the global Max-Planck-Institute Earth System Model (MPI-ESM), and the three generations differ primarily in the ocean initialisation. Ensembles of uninitialised historical and yearly initialised hindcast experiments are used to assess the forecast skill for 10 m wind speeds and wind energy output (Eout) over Central Europe with lead times from one year to one decade. With this aim, a statistical-dynamical downscaling (SDD) approach is used for the regionalisation. Its added value is evaluated by comparison of skill scores for MPI-ESM large-scale wind speeds and SDD-simulated regional wind speeds. All three MPI-ESM ensemble generations show some forecast skill for annual mean wind speed and Eout over Central Europe on yearly and multi-yearly time scales. This forecast skill is mostly limited to the first years after initialisation. Differences between the three ensemble generations are generally small. The regionalisation preserves and sometimes increases the forecast skills of the global runs but results depend on lead time and ensemble generation. Moreover, regionalisation often improves the ensemble spread. Seasonal Eout skills are generally lower than for annual means. Skill scores are lowest during summer and persist longest in autumn. A large-scale westerly weather type with strong pressure gradients over Central Europe is identified as potential source of the skill for wind energy potentials, showing a similar forecast skill and a high correlation with Eout anomalies. These results are promising towards the establishment of a decadal prediction system for wind energy applications over Central Europe.

BRAZILIAN OFFSHORE WAVE CLIMATE BASED ON NWW3 REANALYSIS

This paper provides a description of the wave climate off the Brazilian coast based on an eleven-year time series (Jan/1997-Dec/2007) obtained from the NWW3 operational model hindcast reanalysis. Information about wave climate in Brazilian waters is very scarce and mainly based on occasional short-term observations, the present analysis being the first covering such temporal and spatial scales. To define the wave climate, six sectors were defined and analyzed along the Brazilian shelf-break: South (W1), Southeast (W2), Central (W3), East (W4), Northeast (W5) and North (W6). W1, W2 and W3 wave regimes are determined by the South Atlantic High (SAH) and the passage of synoptic cold fronts; W4, W5 and W6 are controlled by the Intertropical Convergence Zone (ITCZ) and its meridional oscillation. The most energetic waves are from the S, generated by the strong winds associated to the passage of cold fronts, which mainly affect the southern region. Wave power presents a decrease in energy levels from south to north, with its annual variation showing that the winter months are the most energetic in W1 to W4, while in W5 and W6 the most energetic conditions occur during the austral summer. The information presented here provides boundary conditions for studies related to coastal processes, fundamental for a better understanding of the Brazilian coastal zone.

Analysis of coprolites from the extinct mountain goat Myotragus balearicus

Humans colonized the Balearic Islands 5-4 ka ago. They arrived in a uniquely adapted ecosystem with the Balearic mountain goat Myotragus balearicus (Bovidae, Antilopinae, Caprini) as the only large mammal. This mammal went extinct rapidly after human arrival. Several hypotheses have been proposed to explain the extinction of M. balearicus. For the present study ancient DNA analysis (Sanger sequencing, Roche-454, Ion Torrent), and pollen and macrofossil analyses were performed on preserved coprolites from M. balearicus, providing information on its diet and paleo-environment. The information retrieved shows that M. balearicus was heavily dependent on the Balearic box species Buxus balearica during at least part of the year, and that it was most probably a browser. Hindcast ecological niche modelling of B. balearica shows that local distribution of this plant species was affected by climate changes. This suggests that the extinction of M. balearicus can be related to the decline and regional extinction of a plant species that formed a major component of its diet. The vegetation change is thought to be caused by increased aridity occurring throughout the Mediterranean. Previous hypotheses relating the extinction of M. balearicus directly to the arrival of humans on the islands must therefore be adjusted. (C) 2013 University of Washington. Published by Elsevier Inn All rights reserved.

Climate-mediated movement of an avian hybrid zone

The interaction between sibling species that share a zone of contact is a multifaceted relationship affected by climate change [ 1, 2 ]. Between sibling species, interactions may occur at whole-organism (direct or indirect competition) or genomic (hybridization and introgression) levels [ 3–5 ]. Tracking hybrid zone movements can provide insights about influences of environmental change on species interactions [ 1 ]. Here, we explore the extent and mechanism of movement of the contact zone between black-capped chickadees (Poecile atricapillus) and Carolina chickadees (Poecile carolinensis) at whole-organism and genomic levels. We find strong evidence that winter temperatures limit the northern extent of P. carolinensis by demonstrating a current-day association between the range limit of this species and minimum winter temperatures. We further show that this temperature limitation has been consistent over time because we are able to accurately hindcast the previous northern range limit under earlier climate conditions. Using genomic data, we confirm northward movement of this contact zone over the past decade and highlight temporally consistent differential—but limited—geographic introgression of alleles. Our results provide an informative example of the influence of climate change on a contact zone between sibling species.

Sensitivity of the WRF model to PBL parametrisations and nesting techniques: evaluation of wind storms over complex terrain

Simulating surface wind over complex terrain is a challenge in regional climate modelling. Therefore, this study aims at identifying a set-up of the Weather Research and Forecasting Model (WRF) model that minimises system- atic errors of surface winds in hindcast simulations. Major factors of the model configuration are tested to find a suitable set-up: the horizontal resolution, the planetary boundary layer (PBL) parameterisation scheme and the way the WRF is nested to the driving data set. Hence, a number of sensitivity simulations at a spatial resolution of 2 km are carried out and compared to observations. Given the importance of wind storms, the analysis is based on case studies of 24 historical wind storms that caused great economic damage in Switzerland. Each of these events is downscaled using eight different model set-ups, but sharing the same driving data set. The results show that the lack of representation of the unresolved topography leads to a general overestimation of wind speed in WRF. However, this bias can be substantially reduced by using a PBL scheme that explicitly considers the effects of non-resolved topography, which also improves the spatial structure of wind speed over Switzerland. The wind direction, although generally well reproduced, is not very sensitive to the PBL scheme. Further sensitivity tests include four types of nesting methods: nesting only at the boundaries of the outermost domain, analysis nudging, spectral nudging, and the so-called re-forecast method, where the simulation is frequently restarted. These simulations show that restricting the freedom of the model to develop large-scale disturbances slightly increases the temporal agreement with the observations, at the same time that it further reduces the overestimation of wind speed, especially for maximum wind peaks. The model performance is also evaluated in the outermost domains, where the resolution is coarser. The results demonstrate the important role of horizontal resolution, where the step from 6 to 2 km significantly improves model performance. In summary, the combination of a grid size of 2 km, the non-local PBL scheme modified to explicitly account for non-resolved orography, as well as analysis or spectral nudging, is a superior combination when dynamical downscaling is aimed at reproducing real wind fields.

The recent eutrophication of Baldeggersee (Switzerland) as assessed by fossil diatom assemblages

Diatom analyses with an annual resolution were carried out on varves of the hypertrophic Baldeggersee (Central Swiss Plateau) for the timespan ad 1885 to 1993. They reveal seven major changes in the dominant planktonic diatoms. As a result of progressive nutrient enrichment, Baldeggersee changed in the 1910s from a Cyclotella to a Tabellaria fenestrata dominated assemblage, and eventually in the 1950s to a Stephanodiscus parvus dominated diatom assemblage. The timing and direction of diatom-assemblage changes in the varved sediment compare well with sedimentological and limnological observations. Partitioning of the variance in the diatom data revealed that TP is a stronger explanatory variable than temperature for these changes. A diatom-inferred total phosphorus (TP) reconstruction indicates three major steps in eutrophication, occurring at 1909, the mid-1950s and the mid-1970s. Comparison with TP measurements in the water column demonstrates that the diatom-TP inference model used is able to hindcast past TP concentrations reliably. The major steps in eutrophication led to decreases in diatom diversity and also resulted in a progressive increase of calcite grain-size. The lake restoration programme established since 1982 shows no direct impact on the composition of the diatom assemblages. However, the decrease in phosphorus loads since the mid-1970s is reflected in the diatom assemblages and in decreasing diatom-inferred TP concentrations.

Sedimentation at Manop Site H in the eastern equatorial Pacific

Gravity cores recovered from Manganese Nodule Project site H (6°33'N, 92°49'W) show marked downcore variations in the abundance of calcium carbonate, organic carbon, opal, manganese, and other components deposited over the past 400,000 years. Variations in the downcore abundance of organic carbon, which ranges from 0.2 to 1.0%, can be used to hindcast redox conditions in the surface sediments over this time. The results indicate that the depth to the manganese redox boundary varied from about 5 to 25 cm below the seafloor during four major cycles. Downcore variations in solid phase Mn, Ni, and Cu can be produced by such changes in redox conditions. A model which predicts that solid phase Mn can be trapped and buried when the Mn redox boundary migrates rapidly upward is consistent with the observed organic carbon and Mn records and supports the reconstructed redox variations. The history of redox variations at site H can be explained by changes with time in surface water productivity. Major productivity variations at the site occur over 100-kyr cycles, with relatively higher productivity occurring during glacial stages. Thus Quaternary climate changes influence surface water productivity, redox conditions in sediments, and the cycling of transition metals.

Comparison of TEX86 and UK'37 temperature proxies in sinking particles in the Cariaco Basin

The Cariaco Basin, a silled, permanently anoxic basin on the continental shelf of Venezuela with a dynamic chemocline (-240-350 m), has been subject of > 20 years of oceanographic observation and sediment trap studies. We evaluated UK'37 and the TEX86 temperature proxies using sinking particles collected in shallow sediment trap samples at 275 m (Trap A) and 455 m (Trap B) (within and below the chemocline). The organic geochemical temperature proxies, UK'37. (based on coccolithophorid alkenone lipids) and TEX86 (based on archaeal glycerol dialkyl glycerol tetraether (GDGT) lipids), use observed relationships between the ratio of specific lipids and measured sea surface temperature to hindcast past sea surface temperatures. In this study, both UK'37 and TEX86 temperature proxies record seasonal temperature variations, including the cooling associated with upwelling events. UK'37-based temperatures are colder than measured sea surface temperatures, and better correlated temperature at the chlorophyll maximum. In sediment trap material collected below the chemocline (Trap B), UK'37 values are higher than those in Trap A. Warmer subchemocline UK'37 based temperatures may be related to autooxidation of sinking particles, either by small amounts of available oxygen or by alternate electron acceptors concentrated in the biologically dynamic chemocline (e.g. intermediate sulfur compounds). The absolute flux weighted TEX86 temperature values measured in sinking particles from Trap A match the measured SST well. The differences in the TEX86 values between Traps A and B are small and reflect less impact of degradation. Overall, the TEX86 temperatures in sinking particles in the Cariaco Basin reflect annual SST.

Simulated climatically disturbed emergence of agricultures in Western Eurasia 8500-3000 BC

Was the spread of agropastoralism from the Fertile Crescent throughout Europe influenced by rapid climatic shifts? We here generate idealized climate events using palaeoclimate records. In a mathematical model of regional sociocultural development, these events disturb the subsistence base of simulated forager and farmer societies. We evaluate the regional simulated transition timings and durations against a published large set of radiocarbon dates for western Eurasia; the model is able to realistically hindcast much of the inhomogeneous space-time evolution of regional Neolithic transitions. Our study shows that the inclusion of climate events improves the simulation of typical lags between cultural complexes, but that the overall difference to a model without climate events is not significant. Climate events may not have been as important for early sociocultural dynamics as endogenous factors.

Simulated agricultural subsistence, timing of transition, and percentage of immigrant farmers

Farming and herding were introduced to Europe from the Near East and Anatolia; there are, however, considerable arguments about the mechanisms of this transition. Were it the people who moved and either outplaced, or admixed with, the indigenous hunter-gatherer groups? Or was it material and information that moved---the Neolithic Package---consisting of domesticated plants and animals and the knowledge of their use? The latter process is commonly referred to as cultural diffusion and the former as demic diffusion. Despite continuous and partly combined efforts by archaeologists, anthropologists, linguists, palaeontologists and geneticists, a final resolution of the debate has not yet been reached. In the present contribution we interpret results from the Global Land Use and technological Evolution Simulator (GLUES). GLUES is a mathematical model for regional sociocultural development, embedded in the geoenvironmental context, during the Holocene. We demonstrate that the model is able to realistically hindcast the expansion speed and the inhomogeneous space-time evolution of the transition to agropastoralism in western Eurasia. In contrast to models that do not resolve endogenous sociocultural dynamics, our model describes and explains how and why the Neolithic advanced in stages. We uncouple the mechanisms of migration and information exchange and also of migration and the spread of agropastoralism. We find that: (1) An indigenous form of agropastoralism could well have arisen in certain Mediterranean landscapes, but not in Northern and Central Europe, where it depended on imported technology and material. (2) Both demic diffusion by migration and cultural diffusion by trade may explain the western European transition equally well. (3) Migrating farmers apparently contribute less than local adopters to the establishment of agropastoralism. Our study thus underlines the importance of adoption of introduced technologies and economies by resident foragers.

Multidecadal climate oscillations detected in a transparency record from a subtropical Florida lake

Synchronous interannual variability in water transparency observed in neighboring lakes has been linked to regional precipitation and resultant runoff of dissolved organic material, but many climate forcings oscillate over time scales longer than most limnological records can detect. A strong relationship (R2 5 0.86) between transparency and the previous two years’ rainfall and lake stage in a 25-yr record from a Florida lake enabled us to hindcast transparency from a longer 75-yr record of rainfall and lake stage. Predictions revealed a ,30-yr cycle in transparency linked to the Atlantic Multidecadal Oscillation (AMO). Transparency was greatest (4–8 m) in the cool phase of the AMO (,1962–1993) associated with below-average rainfall in south Florida and lowest (0.1– 3.0 m) during two warm phases (,1932–1961, 1994–present) associated with above-average, but more variable, annual rainfall. Models that predict effects of large-scale hydrologic restoration projects on solute export from South Florida’s expansive wetlands need to account for recent entry into a warm AMO phase, where teleconnections between the AMO phases and runoff are opposite of those shown for the U.S. interior.

Stable isotopes, radionuclides, and calculated sea surface temperature of two sediment cores in the Sicily Channel

The Eastern Mediterranean Transient (EMT) occurred in the Aegean Sea from 1988 to 1995 and is the most significant intermediate-to-deep Mediterranean overturning perturbation reported by instrumental records. The EMT was likely caused by accumulation of high salinity waters in the Levantine and enhanced heat loss in the Aegean Sea, coupled with surface water freshening in the Sicily Channel. It is still unknown whether similar transients occurred in the past and, if so, what their forcing processes were. In this study, sediments from the Sicily Channel document surface water freshening (SCFR) at 1910±12, 1812±18, 1725±25 and 1580±30 CE. A regional ocean hindcast links SCFR to enhanced deep-water production and in turn to strengthened Mediterranean thermohaline circulation. Independent evidence collected in the Aegean Sea supports this reconstruction, showing that enhanced bottom water ventilation in the Eastern Mediterranean was associated with each SCFR event. Comparison between the records and multi-decadal atmospheric circulation patterns and climatic external forcings indicates that Mediterranean circulation destabilisation occurs during positive North Atlantic Oscillation (NAO) and negative Atlantic Multidecadal Oscillation (AMO) phases, reduced solar activity and strong tropical volcanic eruptions. They may have recurrently produced favourable deep-water formation conditions, both increasing salinity and reducing temperature on multi-decadal time scales.

Spatio-temporal characteristics of Agulhas leakage: a model inter-comparison study

Investigating the variability of Agulhas leakage, the volume transport of water from the Indian Ocean to the South Atlantic Ocean, is highly relevant due to its potential contribution to the Atlantic Meridional Overturning Circulation as well as the global circulation of heat and salt and hence global climate. Quantifying Agulhas leakage is challenging due to the non-linear nature of this process; current observations are insufficient to estimate its variability and ocean models all have biases in this region, even at high resolution . An Eulerian threshold integration method is developed to examine the mechanisms of Agulhas leakage variability in six ocean model simulations of varying resolution. This intercomparison, based on the circulation and thermo- haline structure at the Good Hope line, a transect to the south west of the southern tip of Africa, is used to identify features that are robust regardless of the model used and takes into account the thermohaline biases of each model. When determined by a passive tracer method, 60 % of the magnitude of Agulhas leakage is captured and more than 80 % of its temporal fluctuations, suggesting that the method is appropriate for investigating the variability of Agulhas leakage. In all simulations but one, the major driver of variability is associated with mesoscale features passing through the section. High resolution (<1/10 deg.) hindcast models agree on the temporal (2–4 cycles per year) and spatial (300–500 km) scales of these features corresponding to observed Agulhas Rings. Coarser resolution models (<1/4 deg.) reproduce similar time scale of variability of Agulhas leakage in spite of their difficulties in representing the Agulhas rings properties. A coarser resolution climate model (2 deg.) does not resolve the spatio-temporal mechanism of variability of Agulhas leakage. Hence it is expected to underestimate the contribution of Agulhas Current System to climate variability.