How well can we simulate past climates? Evaluating the models using global palaeoenvironmental datasets.

Global syntheses of palaeoenvironmental data are required to test climate models under conditions different from the present. Data sets for this purpose contain data from spatially extensive networks of sites. The data are either directly comparable to model output or readily interpretable in terms of modelled climate variables. Data sets must contain sufficient documentation to distinguish between raw (primary) and interpreted (secondary, tertiary) data, to evaluate the assumptions involved in interpretation of the data, to exercise quality control, and to select data appropriate for specific goals. Four data bases for the Late Quaternary, documenting changes in lake levels since 30 kyr BP (the Global Lake Status Data Base), vegetation distribution at 18 kyr and 6 kyr BP (BIOME 6000), aeolian accumulation rates during the last glacial-interglacial cycle (DIRTMAP), and tropical terrestrial climates at the Last Glacial Maximum (the LGM Tropical Terrestrial Data Synthesis) are summarised. Each has been used to evaluate simulations of Last Glacial Maximum (LGM: 21 calendar kyr BP) and/or mid-Holocene (6 cal. kyr BP) environments. Comparisons have demonstrated that changes in radiative forcing and orography due to orbital and ice-sheet variations explain the first-order, broad-scale (in space and time) features of global climate change since the LGM. However, atmospheric models forced by 6 cal. kyr BP orbital changes with unchanged surface conditions fail to capture quantitative aspects of the observed climate, including the greatly increased magnitude and northward shift of the African monsoon during the early to mid-Holocene. Similarly, comparisons with palaeoenvironmental datasets show that atmospheric models have underestimated the magnitude of cooling and drying of much of the land surface at the LGM. The inclusion of feedbacks due to changes in ocean- and land-surface conditions at both times, and atmospheric dust loading at the LGM, appears to be required in order to produce a better simulation of these past climates. The development of Earth system models incorporating the dynamic interactions among ocean, atmosphere, and vegetation is therefore mandated by Quaternary science results as well as climatological principles. For greatest scientific benefit, this development must be paralleled by continued advances in palaeodata analysis and synthesis, which in turn will help to define questions that call for new focused data collection efforts.

GPS/VRS positioning using atmospheric modeling

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Description of atmospheric conditions at the Pierre Auger Observatory using the Global Data Assimilation System (GDAS)

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malargue and averaged monthly models, the utility of the GDAS data is shown. (C) 2012 Elsevier B.V. All rights reserved.

High Resolution Transmission Spectroscopy as a Diagnostic for Jovian Exoplanet Atmospheres: Constraints from Theoretical Models

We present high resolution transmission spectra of giant planet atmospheres from a coupled 3-D atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9 to 55 day orbital periods around solar-type stars. The results of our 3-D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple 1-D models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blue shifts of up to 3 km s−1, whereas less irradiated planets show almost no net Doppler shifts. Compared to 1-D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheric winds. Finally, high resolution transmission spectra may be useful in studying the atmospheres of exoplanets with optically thick clouds since line cores for very strong transitions should remain optically thick to very high altitude. High resolution transmission spectra are an excellent observational test for the validity of 3-D atmospheric dynamics models, because they provide a direct probe of wind structures and heat circulation. Ground-based exoplanet spectroscopy is currently on the verge of being able to verify some of our modeling predictions, most notably the dependence of SSAS winds on insolation. We caution that interpretation of high resolution transmission spectra based on 1-D atmospheric models may be inadequate, as 3-D atmospheric motions can produce a noticeable effect on the absorption signatures.

Studies on orographic effects in a numerical weather prediction model

Numerical models, used for atmospheric research, weather prediction and climate simulation, describe the state of the atmosphere over the heterogeneous surface of the Earth. Several fundamental properties of atmospheric models depend on orography, i.e. on the average elevation of land over a model area. The higher is the models' resolution, the more the details of orography directly influence the simulated atmospheric processes. This sets new requirements for the accuracy of the model formulations with respect to the spatially varying orography. Orography is always averaged, representing the surface elevation within the horizontal resolution of the model. In order to remove the smallest scales and steepest slopes, the continuous spectrum of orography is normally filtered (truncated) even more, typically beyond a few gridlengths of the model. This means, that in the numerical weather prediction (NWP) models, there will always be subgridscale orography effects, which cannot be explicitly resolved by numerical integration of the basic equations, but require parametrization. In the subgrid-scale, different physical processes contribute in different scales. The parametrized processes interact with the resolved-scale processes and with each other. This study contributes to building of a consistent, scale-dependent system of orography-related parametrizations for the High Resolution Limited Area Model (HIRLAM). The system comprises schemes for handling the effects of mesoscale (MSO) and small-scale (SSO) orographic effects on the simulated flow and a scheme of orographic effects on the surface-level radiation fluxes. Representation of orography, scale-dependencies of the simulated processes and interactions between the parametrized and resolved processes are discussed. From the high-resolution digital elevation data, orographic parameters are derived for both momentum and radiation flux parametrizations. Tools for diagnostics and validation are developed and presented. The parametrization schemes applied, developed and validated in this study, are currently being implemented into the reference version of HIRLAM.

Ship-Mounted Real-Time Surface Observational System on board Indian Vessels for Validation and Refinement of Model Forcing Fields

A network of ship-mounted real-time Automatic Weather Stations integrated with Indian geosynchronous satellites Indian National Satellites (INSATs)] 3A and 3C, named Indian National Centre for Ocean Information Services Real-Time Automatic Weather Stations (I-RAWS), is established. The purpose of I-RAWS is to measure the surface meteorological-ocean parameters and transmit the data in real time in order to validate and refine the forcing parameters (obtained from different meteorological agencies) of the Indian Ocean Forecasting System (INDOFOS). Preliminary validation and intercomparison of analyzed products obtained from the National Centre for Medium Range Weather Forecasting and the European Centre for Medium-Range Weather Forecasts using the data collected from I-RAWS were carried out. This I-RAWS was mounted on board oceanographic research vessel Sagar Nidhi during a cruise across three oceanic regimes, namely, the tropical Indian Ocean, the extratropical Indian Ocean, and the Southern Ocean. The results obtained from such a validation and intercomparison, and its implications with special reference to the usage of atmospheric model data for forcing ocean model, are discussed in detail. It is noticed that the performance of analysis products from both atmospheric models is similar and good; however, European Centre for Medium-Range Weather Forecasts air temperature over the extratropical Indian Ocean and wind speed in the Southern Ocean are marginally better.

Investigation of Fundamental Processes Governing Secondary Organic Aerosol Formation in Laboratory Chambers

Our understanding of the processes and mechanisms by which secondary organic aerosol (SOA) is formed is derived from laboratory chamber studies. In the atmosphere, SOA formation is primarily driven by progressive photooxidation of SOA precursors, coupled with their gas-particle partitioning. In the chamber environment, SOA-forming vapors undergo multiple chemical and physical processes that involve production and removal via gas-phase reactions; partitioning onto suspended particles vs. particles deposited on the chamber wall; and direct deposition on the chamber wall. The main focus of this dissertation is to characterize the interactions of organic vapors with suspended particles and the chamber wall and explore how these intertwined processes in laboratory chambers govern SOA formation and evolution.

A Functional Group Oxidation Model (FGOM) that represents SOA formation and evolution in terms of the competition between functionalization and fragmentation, the extent of oxygen atom addition, and the change of volatility, is developed. The FGOM contains a set of parameters that are to be determined by fitting of the model to laboratory chamber data. The sensitivity of the model prediction to variation of the adjustable parameters allows one to assess the relative importance of various pathways involved in SOA formation.

A critical aspect of the environmental chamber is the presence of the wall, which can induce deposition of SOA-forming vapors and promote heterogeneous reactions. An experimental protocol and model framework are first developed to constrain the vapor-wall interactions. By optimal fitting the model predictions to the observed wall-induced decay profiles of 25 oxidized organic compounds, the dominant parameter governing the extent of wall deposition of a compound is identified, i.e., wall accommodation coefficient. By correlating this parameter with the molecular properties of a compound via its volatility, the wall-induced deposition rate of an organic compound can be predicted based on its carbon and oxygen numbers in the molecule.

Heterogeneous transformation of δ-hydroxycarbonyl, a major first-generation product from long-chain alkane photochemistry, is observed on the surface of particles and walls. The uniqueness of this reaction scheme is the production of substituted dihydrofuran, which is highly reactive towards ozone, OH, and NO3, thereby opening a reaction pathway that is not usually accessible to alkanes. A spectrum of highly-oxygenated products with carboxylic acid, ester, and ether functional groups is produced from the substituted dihydrofuran chemistry, thereby affecting the average oxidation state of the alkane-derived SOA.

The vapor wall loss correction is applied to several chamber-derived SOA systems generated from both anthropogenic and biogenic sources. Experimental and modeling approaches are employed to constrain the partitioning behavior of SOA-forming vapors onto suspended particles vs. chamber walls. It is demonstrated that deposition of SOA-forming vapors to the chamber wall during photooxidation experiments can lead to substantial and systematic underestimation of SOA. Therefore, it is likely that a lack of proper accounting for vapor wall losses that suppress chamber-derived SOA yields contribute substantially to the underprediction of ambient SOA concentrations in atmospheric models.

Filtragem adaptativa aplicada à redução de desvios em séries temporais de previsão numérica climática.

Esta dissertação apresenta resultados da aplicação de filtros adaptativos, utilizando os algoritmos NLMS (Normalized Least Mean Square) e RLS (Recursive Least Square), para a redução de desvios em previsões climáticas. As discrepâncias existentes entre o estado real da atmosfera e o previsto por um modelo numérico tendem a aumentar ao longo do período de integração. O modelo atmosférico Eta é utilizado operacionalmente para previsão numérica no CPTEC/INPE e como outros modelos atmosféricos, apresenta imprecisão nas previsões climáticas. Existem pesquisas que visam introduzir melhorias no modelo atmosférico Eta e outras que avaliam as previsões e identificam os erros do modelo para que seus produtos sejam utilizados de forma adequada. Dessa forma, neste trabalho pretende-se filtrar os dados provenientes do modelo Eta e ajustá-los, de modo a minimizar os erros entre os resultados fornecidos pelo modelo Eta e as reanálises do NCEP. Assim, empregamos técnicas de processamento digital de sinais e imagens com o intuito de reduzir os erros das previsões climáticas do modelo Eta. Os filtros adaptativos nesta dissertação ajustarão as séries ao longo do tempo de previsão. Para treinar os filtros foram utilizadas técnicas de agrupamento de regiões, como por exemplo o algoritmo de clusterização k-means, de modo a selecionar séries climáticas que apresentem comportamentos semelhantes entre si. As variáveis climáticas estudadas são o vento meridional e a altura geopotencial na região coberta pelo modelo de previsão atmosférica Eta com resolução de 40 km, a um nível de pressão de 250 hPa. Por fim, os resultados obtidos mostram que o filtro com 4 coeficientes, adaptado pelo algoritmo RLS em conjunto com o critério de seleção de regiões por meio do algoritmo k-means apresenta o melhor desempenho ao reduzir o erro médio e a dispersão do erro, tanto para a variável vento meridional quanto para a variável altura geopotencial.

WASP-3b: a strongly irradiated transiting gas-giant planet

We report the discovery of WASP-3b, the third transiting exoplanet to be discovered by the WASP and SOPHIE collaboration. WASP-3b transits its host star USNO-B1.01256-0285133 every 1.846834 +/- 0.000002 d. Our high-precision radial velocity measurements present a variation with amplitude characteristic of a planetary-mass companion and in phase with the light curve. Adaptive optics imaging shows no evidence for nearby stellar companions, and line-bisector analysis excludes faint, unresolved binarity and stellar activity as the cause of the radial velocity variations. We make a preliminary spectroscopic analysis of the host star and find it to have Teff = 6400 +/- 100K and log g = 4.25 +/- 0.05 which suggests it is most likely an unevolved main-sequence star of spectral type F7-8V. Our simultaneous modelling of the transit photometry and reflex motion of the host leads us to derive a mass of 1.76+0.08-0.14 MJ and radius 1.31+0.07-0.14 RJ for WASP-3b. The proximity and relative temperature of the host star suggests that WASP-3b is one of the hottest exoplanets known, and thus has the potential to place stringent constraints on exoplanet atmospheric models.

Modelling the H Lyman lines in evolved late-type stars

The importance of partial redistribution (PRD) in the modelling of the Lyman a and Lyman ß emission lines of hydrogen in stellar atmospheres is examined using simple atmospheric models of a range of late-type stars. These models represent the subgiant Procyon (F5 IV-V), and the two giants ß Gem (K0 III) and a Tau (K5 III). These stars are selected to span a wide range of surface gravities: 1.25 <log g <4.00. The calculations are performed using the computer code MULTI with the modifications made by Hubeny & Lites. It is found that PRD effects are highly significant, both in the direct prediction of the Lyman line profiles and in the application of hydrostatic equilibrium to calculate the atmospheric electron density in static atmospheric models.

Modelling the chromosphere and transition region of ε Eri (K2 V)

Measurements of ultraviolet line fluxes from Space Telescope Imaging Spectrograph and Far-Ultraviolet Spectroscopic Explorer spectra of the K2-dwarf e Eri are reported. These are used to develop new emission measure distributions and semi-empirical atmospheric models for the chromosphere and lower transition region of the star. These models are the most detailed constructed to date for a main-sequence star other than the Sun. New ionization balance calculations, which account for the effect of finite density on dielectronic recombination rates, are presented for carbon, nitrogen, oxygen and silicon. The results of these calculations are significantly different from the standard Arnaud & Rothenflug ion balance, particularly for alkali-like ions. The new atmospheric models are used to place constraints on possible first ionization potential (FIP)-related abundance variations in the lower atmosphere and to discuss limitations of single-component models for the interpretation of certain optically thick line fluxes. © 2005 RAS.

Ragweed Pollen Source Inventory for France - the Second Largest Centre of Ambrosia in Europe

France, in particular the Rhône-Alpes region, is one of the three main centres of ragweed (Ambrosia) in Europe. The aim of this study is to develop a gridded ragweed pollen source inventory for all of France that can be used in assessments, eradication plans and by atmospheric models for describing concentrations of airborne ragweed pollen. The inventory combines information about spatial variations in annual Ambrosia pollen counts, knowledge of ragweed ecology, detailed land cover information and a Digital Elevation Model. The ragweed inventory consists of a local infection level on a scale of 0–100% (where 100% is the highest plant abundance per area in the studied region) and a European infection level between 0% and 100% (where 100% relates to the highest identified plant abundance in Europe using the same methodology) that has been distributed onto the EMEP grid with 5 km × 5 km resolution. The results of this analysis showed that some of the highest mean annual ragweed pollen concentrations were recorded at Roussillon in the Rhône-Valley. This is reflected by the inventory, where the European infection level has been estimated to reach 67.70% of the most infected areas in Europe i.e. Kecskemét in central Hungary. The inventory shows that the Rhône Valley is the most heavily infected part of France. Central France is also infected, but northern and western parts of France are much less infected. The inventory can be entered into atmospheric transport models, in combination with other components such as a phenological model and a model for daily pollen release, in order to simulate the dispersion of ragweed pollen within France as well as potential long-distance transport from France to other European countries.

Improvement in the Accuracy of Back Trajectories Using WRF to Identify Pollen Sources in Southern Iberian Peninsula

Airborne pollen transport at micro-, meso-gamma and meso-beta scales must be studied by atmospheric models, having special relevance in complex terrain. In these cases, the accuracy of these models is mainly determined by the spatial resolution of the underlying meteorological dataset. This work examines how meteorological datasets determine the results obtained from atmospheric transport models used to describe pollen transport in the atmosphere. We investigate the effect of the spatial resolution when computing backward trajectories with the HYSPLIT model. We have used meteorological datasets from the WRF model with 27, 9 and 3 km resolutions and from the GDAS files with 1 ° resolution. This work allows characterizing atmospheric transport of Olea pollen in a region with complex flows. The results show that the complex terrain affects the trajectories and this effect varies with the different meteorological datasets. Overall, the change from GDAS to WRF-ARW inputs improves the analyses with the HYSPLIT model, thereby increasing the understanding the pollen episode. The results indicate that a spatial resolution of at least 9 km is needed to simulate atmospheric flows that are considerable affected by the relief of the landscape. The results suggest that the appropriate meteorological files should be considered when atmospheric models are used to characterize the atmospheric transport of pollen on micro-, meso-gamma and meso-beta scales. Furthermore, at these scales, the results are believed to be generally applicable for related areas such as the description of atmospheric transport of radionuclides or in the definition of nuclear-radioactivity emergency preparedness.

An Analysis of Local Wind and Air Mass Directions and Their Impact on Cladosporium Distribution Using HYSPLIT and Circular Statistics

Ecological studies that examine species-environment relationships are often limited to several meteorological parameters, i.e. mean air temperature, relative humidity, precipitation, vapour pressure deficit and solar radiation. The impact of local wind, its speed and direction are less commonly investigated in aerobiological surveys mainly due to difficulties related to the employment of specific analytical tools and interpretation of their outputs. Identification of inoculum sources of economically important plant pathogens, as well as highly allergenic bioaerosols like Cladosporium species, has not been yet explored with remote sensing data and atmospheric models such as Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT). We, therefore, performed an analysis of 24 h intra-diurnal cycle of Cladosporium spp. spores from an urban site in connection with both the local wind direction and overall air mass direction computed by HYSPLIT. The observational method was a volumetric air sampler of the Hirst design with 1 h time resolution and corresponding optical detection of fungal spores with light microscopy. The atmospheric modelling was done using the on-line data set from GDAS with 1° resolution and circular statistical methods. Our results showed stronger, statistically significant correlation (p ≤ 0.05) between high Cladosporium spp. spore concentration and air mass direction compared to the local wind direction. This suggested that a large fraction of the investigated fungal spores had a regional origin and must be located more than a few kilometers away from the sampling point.

Future Air Quality Related Health Effects in Europe and the Nordic Region—Sensitivity to Changes in Climate, Anthropogenic Emissions, Demography and Building Stock

Future changes in population exposures to ambient air pollution are inherently linked with long-term trends in outdoor air quality, but also with changes in the building stock. Moreover, the burden of disease is further driven by the ageing of the European populations. This study aims to assess the impact of changes in climate, emissions, building stocks and population on air pollution related human health impacts across Europe in the future. Therefore an integrated assessment model combining atmospheric models and health impacts has been setup for projections of the future developments in air pollution related premature mortality. The focus is here on the regional scale impacts of exposure to surface ozone (O3), Secondary Inorganic Aerosols (SIA) and primary particulate matter (PPM).