Reducing noise associated with the Monte Carlo Independent Column Approximation for weather forecasting models

The Monte Carlo Independent Column Approximation (McICA) is a flexible method for representing subgrid-scale cloud inhomogeneity in radiative transfer schemes. It does, however, introduce conditional random errors but these have been shown to have little effect on climate simulations, where spatial and temporal scales of interest are large enough for effects of noise to be averaged out. This article considers the effect of McICA noise on a numerical weather prediction (NWP) model, where the time and spatial scales of interest are much closer to those at which the errors manifest themselves; this, as we show, means that noise is more significant. We suggest methods for efficiently reducing the magnitude of McICA noise and test these methods in a global NWP version of the UK Met Office Unified Model (MetUM). The resultant errors are put into context by comparison with errors due to the widely used assumption of maximum-random-overlap of plane-parallel homogeneous cloud. For a simple implementation of the McICA scheme, forecasts of near-surface temperature are found to be worse than those obtained using the plane-parallel, maximum-random-overlap representation of clouds. However, by applying the methods suggested in this article, we can reduce noise enough to give forecasts of near-surface temperature that are an improvement on the plane-parallel maximum-random-overlap forecasts. We conclude that the McICA scheme can be used to improve the representation of clouds in NWP models, with the provision that the associated noise is sufficiently small.

Possible impacts of a future grand solar minimum on climate: stratospheric and global circulation changes

It has been suggested that the Sun may evolve into a period of lower activity over the 21st century. This study examines the potential climate impacts of the onset of an extreme ‘Maunder Minimum like’ grand solar minimum using a comprehensive global climate model. Over the second half of the 21st century, the scenario assumes a decrease in total solar irradiance of 0.12% compared to a reference RCP8.5 experiment. The decrease in solar irradiance cools the stratopause (~1 hPa) in the annual and global mean by 1.4 K. The impact on global mean near-surface temperature is small (~−0.1 K), but larger changes in regional climate occur during the stratospheric dynamically active seasons. In Northern hemisphere (NH) winter-time, there is a weakening of the stratospheric westerly jet by up to ~3-4 m s1, with the largest changes occurring in January-February. This is accompanied by a deepening of the Aleutian low at the surface and an increase in blocking over northern Europe and the north Pacific. There is also an equatorward shift in the Southern hemisphere (SH) midlatitude eddy-driven jet in austral spring. The occurrence of an amplified regional response during winter and spring suggests a contribution from a top-down pathway for solar-climate coupling; this is tested using an experiment in which ultraviolet (200–320 nm) radiation is decreased in isolation of other changes. The results show that a large decline in solar activity over the 21st century could have important impacts on the stratosphere and regional surface climate.

Exploring strategies for coupled 4D-Var data assimilation using an idealised atmosphere-ocean model

Operational forecasting centres are currently developing data assimilation systems for coupled atmosphere-ocean models. Strongly coupled assimilation, in which a single assimilation system is applied to a coupled model, presents significant technical and scientific challenges. Hence weakly coupled assimilation systems are being developed as a first step, in which the coupled model is used to compare the current state estimate with observations, but corrections to the atmosphere and ocean initial conditions are then calculated independently. In this paper we provide a comprehensive description of the different coupled assimilation methodologies in the context of four dimensional variational assimilation (4D-Var) and use an idealised framework to assess the expected benefits of moving towards coupled data assimilation. We implement an incremental 4D-Var system within an idealised single column atmosphere-ocean model. The system has the capability to run both strongly and weakly coupled assimilations as well as uncoupled atmosphere or ocean only assimilations, thus allowing a systematic comparison of the different strategies for treating the coupled data assimilation problem. We present results from a series of identical twin experiments devised to investigate the behaviour and sensitivities of the different approaches. Overall, our study demonstrates the potential benefits that may be expected from coupled data assimilation. When compared to uncoupled initialisation, coupled assimilation is able to produce more balanced initial analysis fields, thus reducing initialisation shock and its impact on the subsequent forecast. Single observation experiments demonstrate how coupled assimilation systems are able to pass information between the atmosphere and ocean and therefore use near-surface data to greater effect. We show that much of this benefit may also be gained from a weakly coupled assimilation system, but that this can be sensitive to the parameters used in the assimilation.

Processing of zero-derived words in English: an fMRI investigation

Derivational morphological processes allow us to create new words (e.g. punish (V) to noun (N) punishment) from base forms. The number of steps from the basic units to derived words often varies (e.g., nationalitybridge-V) i.e., zero-derivation ( Aronoff, 1980). We compared the processing of one-step (soakingstructure, our findings suggest that morphological processing is based on underlying morphological complexity, independent of overt affixation. This study is the first to provide evidence for the processing of zero derivation, and demonstrates that morphological processing cannot be reduced to surface form-based segmentation.

Geotechnical systems that evolve with ecological processes

Geotechnical systems, such as landfills, mine tailings storage facilities (TSFs), slopes, and levees, are required to perform safely throughout their service life, which can span from decades for levees to “in perpetuity” for TSFs. The conventional design practice by geotechnical engineers for these systems utilizes the as-built material properties to predict its performance throughout the required service life. The implicit assumption in this design methodology is that the soil properties are stable through time. This is counter to long-term field observations of these systems, particularly where ecological processes such as plant, animal, biological, and geochemical activity are present. Plant roots can densify soil and/or increase hydraulic conductivity, burrowing animals can increase seepage, biological activity can strengthen soil, geochemical processes can increase stiffness, etc. The engineering soil properties naturally change as a stable ecological system is gradually established following initial construction, and these changes alter system performance. This paper presents an integrated perspective and new approach to this issue, considering ecological, geotechnical, and mining demands and constraints. A series of data sets and case histories are utilized to examine these issues and to propose a more integrated design approach, and consideration is given to future opportunities to manage engineered landscapes as ecological systems. We conclude that soil scientists and restoration ecologists must be engaged in initial project design and geotechnical engineers must be active in long-term management during the facility’s service life. For near-surface geotechnical structures in particular, this requires an interdisciplinary perspective and the embracing of soil as a living ecological system rather than an inert construction material.

Mechanisms underlying temperature extremes in Iberia: a Lagrangian perspective

The mechanisms underlying the occurrence of temperature extremes in Iberia are analysed considering a Lagrangian perspective of the atmospheric flow, using 6-hourly ERA-Interim reanalysis data for the years 1979–2012. Daily 2-m minimum temperatures below the 1st percentile and 2-m maximum temperatures above the 99th percentile at each grid point over Iberia are selected separately for winter and summer. Four categories of extremes are analysed using 10-d backward trajectories initialized at the extreme temperature grid points close to the surface: winter cold (WCE) and warm extremes (WWE), and summer cold (SCE) and warm extremes (SWE). Air masses leading to temperature extremes are first transported from the North Atlantic towards Europe for all categories. While there is a clear relation to large-scale circulation patterns in winter, the Iberian thermal low is important in summer. Along the trajectories, air mass characteristics are significantly modified through adiabatic warming (air parcel descent), upper-air radiative cooling and near-surface warming (surface heat fluxes and radiation). High residence times over continental areas, such as over northern-central Europe for WCE and, to a lesser extent, over Iberia for SWE, significantly enhance these air mass modifications. Near-surface diabatic warming is particularly striking for SWE. WCE and SWE are responsible for the most extreme conditions in a given year. For WWE and SCE, strong temperature advection associated with important meridional air mass transports are the main driving mechanisms, accompanied by comparatively minor changes in the air mass properties. These results permit a better understanding of mechanisms leading to temperature extremes in Iberia.

The impact of oceanic heat transport on the atmospheric circulation

A general circulation model of intermediate complexity with an idealized Earth-like aquaplanet setup is used to study the impact of changes in the oceanic heat transport on the global atmospheric circulation. Focus is on the atmospheric mean meridional circulation and global thermodynamic properties. The atmosphere counterbalances to a large extent the imposed changes in the oceanic heat transport, but, nonetheless, significant modifications to the atmospheric general circulation are found. Increasing the strength of the oceanic heat transport up to 2.5 PW leads to an increase in the global mean near-surface temperature and to a decrease in its equator-to-pole gradient. For stronger transports, the gradient is reduced further, but the global mean remains approximately constant. This is linked to a cooling and a reversal of the temperature gradient in the tropics. Additionally, a stronger oceanic heat transport leads to a decline in the intensity and a poleward shift of the maxima of both the Hadley and Ferrel cells. Changes in zonal mean diabatic heating and friction impact the properties of the Hadley cell, while the behavior of the Ferrel cell is mostly controlled by friction. The efficiency of the climate machine, the intensity of the Lorenz energy cycle and the material entropy production of the system decline with increased oceanic heat transport. This suggests that the climate system becomes less efficient and turns into a state of reduced entropy production as the enhanced oceanic transport performs a stronger large-scale mixing between geophysical fluids with different temperatures, thus reducing the available energy in the climate system and bringing it closer to a state of thermal equilibrium.

Scale interactions between the MJO and the western maritime continent

State-of-the-art regional climate model simulations that are able to resolve key mesoscale circulations are used, for the first time, to understand the interaction between the large-scale convective environment of the MJO and processes governing the strong diurnal cycle over the islands of the Maritime Continent (MC). Convection is sustained in the late afternoon just inland of the coasts due to sea breeze convergence. Previous work has shown that the variability in MC rainfall associated with the MJO is manifested in changes to this diurnal cycle; land-based rainfall peaks before the active convective envelope of the MJO reaches the MC, whereas oceanic rainfall rates peak whilst the active envelope resides over the region. The model simulations show that the main controls on oceanic MC rainfall in the early active MJO phases are the large-scale environment and atmospheric stability, followed by high oceanic latent heat flux forced by high near-surface winds in the later active MJO phases. Over land, rainfall peaks before the main convective envelope arrives (in agreement with observations), even though the large-scale convective environment is only moderately favourable for convection. The causes of this early rainfall peak are convective triggers from land-sea breeze circulations that are strong due to high surface insolation and surface heating. During the peak MJO phases cloud cover increases and surface insolation decreases, which weakens the strength of the mesoscale circulations and reduces land-based rainfall, even though the large-scale environment remains favourable for convection at this time. Hence, scale interactions are an essential part of the MJO transition across the MC.

Still searching for graves: an analytical strategy for interpreting geophysical data used in the search for “unmarked” graves

Searching for and mapping the physical extent of unmarked graves using geophysical techniques has proven difficult in many cases. The success of individual geophysical techniques for detecting graves depends on a site-by-site basis. Significantly, detection of graves often results from measured contrasts that are linked to the background soils rather than the type of archaeological feature associated with the grave. It is evident that investigation of buried remains should be considered within a 3D space as the variation in burial environment can be extremely varied through the grave. Within this paper, we demonstrate the need for a multi-method survey strategy to investigate unmarked graves, as applied at a “planned” but unmarked pauper’s cemetery. The outcome from this case study provides new insights into the strategy that is required at such sites. Perhaps the most significant conclusion is that unmarked graves are best understood in terms of characterization rather than identification. In this paper, we argue for a methodological approach that, while following the current trends to use multiple techniques, is fundamentally dependent on a structured approach to the analysis of the data. The ramifications of this case study illustrate the necessity of an integrated strategy to provide a more holistic understanding of unmarked graves that may help aid in management of these unseen but important aspects of our heritage. It is concluded that the search for graves is still a current debate and one that will be solved by methodological rather than technique-based arguments.

How do green roofs mitigate urban thermal stress under heat waves?

As the climate warms, heat waves (HW) are projected to be more intense and to last longer, with serious implications for public health. Urban residents face higher health risks because urban heat islands (UHIs) exacerbate HW conditions. One strategy to mitigate negative impacts of urban thermal stress is the installation of green roofs (GRs) given their evaporative cooling effect. However, the effectiveness of GRs and the mechanisms by which they have an effect at the scale of entire cities are still largely unknown. The Greater Beijing Region (GBR) is modeled for a HW scenario with the Weather Research and Forecasting (WRF) model coupled with a state-of-the-art urban canopy model (PUCM) to examine the effectiveness of GRs. The results suggest GR would decrease near-surface air temperature (ΔT2max = 2.5 K) and wind speed (ΔUV10max = 1.0 m s-1) but increase atmospheric humidity (ΔQ2max = 1.3 g kg-1). GRs are simulated to lessen the overall thermal stress as indicated by apparent temperature (ΔAT2max = 1.7 °C). The modifications by GRs scale almost linearly with the fraction of the surface they cover. Investigation of the surface-atmosphere interactions indicate that GRs with plentiful soil moisture dissipate more of the surface energy as latent heat flux and subsequently inhibit the development of the daytime planetary boundary layer (PBL). This causes the atmospheric heating through entrainment at the PBL top to be decreased. Additionally, urban GRs modify regional circulation regimes leading to decreased advective heating under HW.

Seasonal drought stress in the Amazon: Reconciling models and observations

The Amazon Basin is crucial to global circulatory and carbon patterns due to the large areal extent and large flux magnitude. Biogeophysical models have had difficulty reproducing the annual cycle of net ecosystem exchange (NEE) of carbon in some regions of the Amazon, generally simulating uptake during the wet season and efflux during seasonal drought. In reality, the opposite occurs. Observational and modeling studies have identified several mechanisms that explain the observed annual cycle, including: (1) deep soil columns that can store large water amount, (2) the ability of deep roots to access moisture at depth when near-surface soil dries during annual drought, (3) movement of water in the soil via hydraulic redistribution, allowing for more efficient uptake of water during the wet season, and moistening of near-surface soil during the annual drought, and (4) photosynthetic response to elevated light levels as cloudiness decreases during the dry season. We incorporate these mechanisms into the third version of the Simple Biosphere model (SiB3) both singly and collectively, and confront the results with observations. For the forest to maintain function through seasonal drought, there must be sufficient water storage in the soil to sustain transpiration through the dry season in addition to the ability of the roots to access the stored water. We find that individually, none of these mechanisms by themselves produces a simulation of the annual cycle of NEE that matches the observed. When these mechanisms are combined into the model, NEE follows the general trend of the observations, showing efflux during the wet season and uptake during seasonal drought.

Scale architecture of the palmar and plantar epidermis of Polychrus acutirostris Spix, 1825 (Iguania, Polychrotidae) and its relationship to arboreal locomotion

In tetrapod squamates, the diversity of micro-ornamentations of the epidermis of the contact areas of hands and feet is generally associated with constraints and modalities related to locomotion. Polychrus acutirostris is a medium-sized lizard that occurs in open heterogeneous habitats in South America, such as the cerrados, caatingas, and fallow lands. It progresses slowly on branches of various diameters in its arboreal environment. It can also move more rapidly on the ground. The hands and feet are prehensile and may be considered an adaptation for grasping and climbing. Epidermal surfaces from the palmar and plantar areas of the hands and feet of P. acutirostris were prepared for SEM examination, and studied at various magnifications. They show three major levels of complexity: (1) scale types, organized in gradients of size and imbrication, (2) scalar ornamentations, organized by increasing complexity and polarity, and (3) presence of Oberhautchen showing typically iguanian honeycomb micro-ornamentations. The shape and surface structure of the scales with their pattern of micro-ornamental peaks, which improve grip, and the grasping hands and feet indicate that P. acutirostris is morpho-functionally specialized for arboreality. (C) 2009 Elsevier GmbH. All rights reserved.

Chemical bonding and composition of silicon nitride films prepared by inductively coupled plasma chemical vapor deposition

Thin silicon nitride films were prepared at 350 degrees C by inductively coupled plasma chemical vapor deposition on Si(100) substrates under different NH(3)/SiH(4) or N(2)/SiH(4) gas mixture. The chemical composition and bonding structure of the deposited films were investigated as a function of the process parameters, such as the gas flow ratio NH(3)/SiH(4) or N(2)/SiH(4) and the RF power, using X-ray photoelectron spectroscopy (XPS). The gas flow ratio was 1.4, 4.3, 7.2 or 9.5 and the RF power, 50 or 100 W. Decomposition results of Si 2p XPS spectra indicated the presence of bulk Si, under-stoichiometric nitride, stoichiometric nitride Si(3)N(4), oxynitride SiN(x)O(y), and stoichiometric oxide SiO(2), and the amounts of these compounds were strongly influenced by the two process parameters. These results were consistent with those obtained from N 1s XPS spectra. The chemical composition ratio N/Si in the film increased with increasing the gas flow ratio until the gas flow ratio reached 4.3, reflecting the high reactivity of nitrogen, and stayed almost constant for further increase in gas flow ratio, the excess nitrogen being rejected from the growing film. A considerable and unexpected incorporation of contaminant oxygen and carbon into the depositing film was observed and attributed to their high chemical reactivity. (C) 2010 Elsevier B.V. All rights reserved.

Thermoluminescence and synchrotron radiation studies on the persistent luminescence of BaAl(2)O(4)/Eu(2+),Dy(3+)

The persistent luminescence materials, barium aluminates doped with Eu(2+) and Dy(3+) (BaAl(2)O(4): Eu(2+),Dy(3+)), were prepared with the combustion synthesis at temperatures between 400 and 600 degrees C as well as with the solid state reaction at 1500 degrees C. The concentrations of Eu(2+)/Dy(3+) (in mol% of the Ba amount) ranged from 0.1/0.1 to 1.0/3.0. The electronic and defect energy level structures were studied with thermoluminescence (TL) and synchrotron radiation (SR) spectroscopies: UV-VUV excitation and emission, as well as with X-ray absorption near-edge structure (XANES) methods. Theoretical calculations using the density functional theory (DFT) were carried out in order to compare with the experimental data. (C) 2010 Elsevier Inc. All rights reserved.

Intralamellar structural modifications related to the proton exchanging in K(4)Nb(6)O(17) layered phase

Basic structural aspects about the layered hexaniobate of K(4)Nb(6)O(17) composition and its proton-exchanged form were investigated mainly by spectroscopic techniques. Raman spectra of hydrous K(4)Nb(6)O(17) and H(2)K(2)Nb(6)O(17)center dot H(2)O show significant modifications in the 950-800 cm(-1) region (Nb-O stretching mode of highly distorted NbO(6) octahedra). The band at 900 cm(-1) shifts to 940 cm(-1) after the replacement of K(+) ion by proton. Raman spectra of the original materials and the related deuterated samples are similar suggesting that no isotopic effect occurs. Major modifications were observed when H(2)K(2)Nb(6)O(17) was dehydrated: the relative intensity of the band at 940 cm(-1) decreases and new bands seems to be present at about 860-890 cm(-1). The H(+) ions should be shielded by the hydration sphere what preclude the interaction with the layers. Removing the water molecules, H(+) ions can establish a strong interaction with oxygen atoms, decreasing the bond order of Nb-O linkage. X-ray absorption near edge structure studies performed at Nb K-edge indicate that the niobium coordination number and oxidation state remain identical after the replacement of potassium by proton. From the refinement of the fine structure, it appears that the Nb-Nb coordination shell is divided into two main contributions of about 0.33 and 0.39 nm, and interestingly the population, i.e., the number of backscattering atoms is inversed between the two hexaniobate materials. 2009 Elsevier Ltd. All rights reserved.