How accurate are volcanic ash simulations of the 2010 Eyjafjallajökull eruption?

In the event of a volcanic eruption the decision to close airspace is based on forecast ash maps, produced using volcanic ash transport and dispersion models. In this paper we quantitatively evaluate the spatial skill of volcanic ash simulations using satellite retrievals of ash from the Eyja allajökull eruption during the period from 7 to 16 May 2010. We find that at the start of this period, 7–10 May, the model (FLEXible PARTicle) has excellent skill and can predict the spatial distribution of the satellite-retrieved ash to within 0.5∘ × 0.5∘ latitude/longitude. However, on 10 May there is a decrease in the spatial accuracy of the model to 2.5∘× 2.5∘ latitude/longitude, and between 11 and 12 May the simulated ash location errors grow rapidly. On 11 May ash is located close to a bifurcation point in the atmosphere, resulting in a rapid divergence in the modeled and satellite ash locations. In general, the model skill reduces as the residence time of ash increases. However, the error growth is not always steady. Rapid increases in error growth are linked to key points in the ash trajectories. Ensemble modeling using perturbed meteorological data would help to represent this uncertainty, and assimilation of satellite ash data would help to reduce uncertainty in volcanic ash forecasts.

Advancing polar prediction capabilities on daily to seasonal time scales

It is argued that existing polar prediction systems do not yet meet users’ needs; and possible ways forward in advancing prediction capacity in polar regions and beyond are outlined. The polar regions have been attracting more and more attention in recent years, fuelled by the perceptible impacts of anthropogenic climate change. Polar climate change provides new opportunities, such as shorter shipping routes between Europe and East Asia, but also new risks such as the potential for industrial accidents or emergencies in ice-covered seas. Here, it is argued that environmental prediction systems for the polar regions are less developed than elsewhere. There are many reasons for this situation, including the polar regions being (historically) lower priority, with less in situ observations, and with numerous local physical processes that are less well-represented by models. By contrasting the relative importance of different physical processes in polar and lower latitudes, the need for a dedicated polar prediction effort is illustrated. Research priorities are identified that will help to advance environmental polar prediction capabilities. Examples include an improvement of the polar observing system; the use of coupled atmosphere-sea ice-ocean models, even for short-term prediction; and insight into polar-lower latitude linkages and their role for forecasting. Given the enormity of some of the challenges ahead, in a harsh and remote environment such as the polar regions, it is argued that rapid progress will only be possible with a coordinated international effort. More specifically, it is proposed to hold a Year of Polar Prediction (YOPP) from mid-2017 to mid-2019 in which the international research and operational forecasting community will work together with stakeholders in a period of intensive observing, modelling, prediction, verification, user-engagement and educational activities.

ENSO related variation of equatorial MRG and Rossby waves and forcing from higher latitudes

The contrasting behaviour of westward-moving mixed Rossby-gravity (WMRG) and the first Rossby (R1) waves in El Niño (EN) and La Niña (LN) seasons is documented with a focus on the Northern Hemisphere winter. The eastward-moving variance in the upper troposphere is dominated by WMRG and R1 structures that appear to be Doppler-shifted by the flow and are referred to as WMRG-E and R1-E. In the East Pacific and Atlantic the years with stronger equatorial westerly winds have the stronger WMRG and WMRG- E. In the East Pacific, R1 is also a maximum in LN. However, R1-E exhibits an eastward-shift between LN and EN. The changes with ENSO phase provide a test-bed for the understanding of these waves. In the East Pacific and Atlantic, the stronger WMRG-E and WMRG with stronger westerlies are in accord with the dispersion relation with simple Doppler-shifting by the zonal flow. The possible existence of free waves can also explain stronger R1 in EN in the Eastern Hemisphere. 1-D free wave propagation theory based on wave activity conservation is also important for R1. However, this theory is unable to explain the amplitude maxima for other waves observed in the strong equatorial westerly regions in the Western Hemisphere, and certainly not their ENSO-related variation. The forcing of equatorial waves by higher latitude wave activity and its variation with ENSO phase is therefore examined. Propagation of extratropical eastward-moving Rossby wave activity through the westerly ducts into the equatorial region where it triggers WMRG-E is favoured in the stronger westerlies, in LN in the East Pacific and EN in the Atlantic. It is also found that WMRG is forced by Southern Hemisphere westward-moving wavetrains arching into the equatorial region where they are reflected. The most significant mechanism for both R1 and R1-E appear to be lateral forcing by subtropical wavetrains.

A comparison between large-scale irregularities and scintillations in the polar ionosphere

A comparison tool has been developed by mapping the global GPS total electron content (TEC) and large coverage of ionospheric scintillations together on the geomagnetic latitude/magnetic local time coordinates. Using this tool, a comparison between large-scale ionospheric irregularities and scintillations are pursued during a geomagnetic storm. Irregularities, such as storm enhanced density (SED), middle-latitude trough and polar cap patches, are clearly identified from the TEC maps. At the edges of these irregularities, clear scintillations appeared but their behaviors were different. Phase scintillations (σsub{φ}) were almost always larger than amplitude scintillations (S4) at the edges of these irregularities, associated with bursty flows or flow reversals with large density gradients. An unexpected scintillation feature appeared inside the modeled auroral oval where S4 were much larger than σsub{φ}, most likely caused by particle precipitations around the exiting polar cap patches.

Poorest countries experience earlier anthropogenic emergence of daily temperature extremes

Understanding how the emergence of the anthropogenic warming signal from the noise of internal variability translates to changes in extreme event occurrence is of crucial societal importance. By utilising simulations of cumulative carbon dioxide (CO2) emissions and temperature changes from eleven earth system models, we demonstrate that the inherently lower internal variability found at tropical latitudes results in large increases in the frequency of extreme daily temperatures (exceedances of the 99.9th percentile derived from pre-industrial climate simulations) occurring much earlier than for mid-to-high latitude regions. Most of the world's poorest people live at low latitudes, when considering 2010 GDP-PPP per capita; conversely the wealthiest population quintile disproportionately inhabit more variable mid-latitude climates. Consequently, the fraction of the global population in the lowest socio-economic quintile is exposed to substantially more frequent daily temperature extremes after much lower increases in both mean global warming and cumulative CO2 emissions.

Formulation of a tropical town energy budget (t-TEB) scheme

This work describes the tropical town energy budget (t-TEB) scheme addressed to simulate the diurnal occurrence of the urban heat island (UHI) as observed in the Metropolitan Area of Rio de Janeiro (MARJ; -22A degrees S; -44A degrees W) in Brazil. Reasoning about the tropical urban climate have guided the scheme implementation, starting from the original equations from Masson (Bound-Lay Meteorol 94:357-397, 2000). The modifications include (a) local scaling approaches for obtaining flux-gradient relationships in the roughness sub-layer, (b) the Monin-Obukhov similarity framework in the inertial sub-layer, (c) increasing aerodynamic conductance toward more unstable conditions, and (d) a modified urban subsurface drainage system to transfer the intercepted rainwater by roofs to the roads. Simulations along 2007 for the MARJ are obtained and compared with the climatology. The t-TEB simulation is consistent with the observations, suggesting that the timing and dynamics of the UHI in tropical cities could vary significantly from the familiar patterns observed in mid-latitude cities-with the peak heat island intensity occurring in the morning than at night. The simulations are suggesting that the thermal phase shift of this tropical diurnal UHI is a response of the surface energy budget to the large amount of solar radiation, intense evapotranspiration, and thermal response of the vegetated surfaces over a very humid soil layer.

Climate perspective on the large-scale circulation associated with the transition of the first South Atlantic hurricane

The landfall of Cyclone Catarina on the Brazilian coast in March 2004 became known as the first documented hurricane in the South Atlantic Ocean, promoting a new view oil how large-scale features can contribute to tropical transition. The aim of this paper is to put the large-scale circulation associated with Catarina`s transition in climate perspective. This is discussed in the light of a robust pattern of spatial correlations between thermodynamic and dynamic variables of importance for hurricane formation. A discussion on how transition mechanisms respond to the present-day circulation is presented. These associations help in understanding why Catarina was formed in a region previously thought to be hurricane-free. Catarina developed over a large-scale area of thermodynamically favourable air/sea temperature contrast. This aspect explains the paradox that such a rare system developed when the sea surface temperature was slightly below average. But, although thermodynamics played an important role, it is apparent that Catarina would not have formed without the key dynamic interplay triggered by a high latitude blocking. The blocking was associated with an extreme positive phase of the Southern Annular Mode (SAM) both hemispherically and locally, and the nearby area where Catarina developed is found to be more cyclonic during the positive phase of the SAM. A conceptual model is developed and a `South Atlantic index` is introduced as a useful diagnostic of potential conditions leading to tropical transition in the area, where large-scale indices indicate trends towards more favourable atmospheric conditions for tropical cyclone formation. Copyright (c) 2008 Royal Meteorological Society

Global and diffuse solar irradiances in urban and rural areas in southeast Brazil

The seasonal evolution of daily and hourly values of global and diffuse solar radiation at the surface are compared for the cities of Sao Paulo and Botucatu, both located in Southeast Brazil and representative of urban and rural areas, respectively. The comparisons are based on measurements of global and diffuse solar irradiance carried out at the surface during a six year simultaneous period in these two cities. Despite the similar latitude and altitude, the seasonal evolution of daily values indicate that Sao Paulo receives, during clear sky days, 7.8% less global irradiance in August and 5.1% less in June than Botucatu. On the other hand, Sao Paulo receives, during clear sky days, 3.6% more diffuse irradiance in August and 15.6% more in June than Botucatu. The seasonal variation of the diurnal cycle confirms these differences and indicates that they are more pronounced during the afternoon. The regional differences are related to the distance from the Atlantic Ocean, systematic penetration of the sea breeze and daytime evolution of the particulate matter in Sao Paulo. An important mechanism controlling the spatial distribution of solar radiation, on a regional scale, is the sea breeze penetration in Sao Paulo, bringing moisture and maritime aerosol that in turn further increases the solar radiation scattering due to pollution and further reduces the intensity of the direct component of solar radiation at the surface. Surprisingly, under clear sky conditions the atmospheric attenuation of solar radiation in Botucatu during winter - the biomass burning period due to the sugar cane harvest - is equivalent to that at Sao Paulo City, indicating that the contamination during sugar cane harvest in Southeast Brazil has a large impact in the solar radiation field at the surface.

Anthropogenic heat in the city of So Paulo, Brazil

The main goal of this work is to describe the anthropogenic energy flux (Q (F)) in the city of So Paulo, Brazil. The hourly, monthly, and annual values of the anthropogenic energy flux are estimated using the inventory method, and the contributions of vehicular, stationary, and human metabolism sources from 2004 to 2007 are considered. The vehicular and stationary sources are evaluated using the primary consumption of energy based on fossil fuel, bio fuel, and electricity usage by the population. The diurnal evolution of the anthropogenic energy flux shows three relative maxima, with the largest maxima occurring early in the morning (similar to 19.9 Wm(-2)) and in the late afternoon (similar to 20.3 Wm(-2)). The relative maximum that occurs around noontime (similar to 19.6 Wm(-2)) reflects the diurnal pattern of vehicle traffic that seems to be specific to So Paulo. With respect to diurnal evolution, the energy flux released by vehicular sources (Q (FV)) contributes approximately 50% of the total anthropogenic energy flux. Stationary sources (Q (FS)) and human metabolism (Q (FM)) represent about 41% and 9% of the anthropogenic energy flux, respectively. For 2007, the monthly values of Q (FV), Q (FS), Q (FM), and Q (F) are, respectively, 16.8 +/- 0.25, 14.3 +/- 0.16, 3.5 +/- 0.03, and 34.6 +/- 0.41 MJ m(-2) month(-1). The seasonal evolution monthly values of Q (FV), Q (FS), Q (FM), and Q (F) show a relative minimum during the summer and winter vacations and a systematic and progressive increase associated with the seasonal evolution of the economic activity in So Paulo. The annual evolution of Q (F) indicates that the city of So Paulo released 355.2 MJ m(-2) year(-1) in 2004 and 415.5 MJ m(-2) year(-1) in 2007 in association with an annual rate of increase of 19.6 MJ m(-2) year(-1) (from 2004 to 2006) and 30.5 MJ m(-2) year(-1) (from 2006 to 2007). The anthropogenic energy flux corresponds to about 9% of the net radiation at the surface in the summer and 15% in the winter. The amplitude of seasonal variation of the maximum hourly value of the diurnal variation increases exponentially with latitude.

The spiral structure of the Galaxy revealed by CS sources and evidence for the 4:1 resonance

We present a map of the spiral structure of the Galaxy, as traced by molecular carbon monosulphide (CS) emission associated with IRAS sources which are believed to be compact H II regions. The CS line velocities are used to determine the kinematic distances of the sources in order to investigate their distribution in the galactic plane. This allows us to use 870 objects to trace the arms, a number larger than that of previous studies based on classical H II regions. The distance ambiguity of the kinematic distances, when it exists, is solved by different procedures, including the latitude distribution and an analysis of the longitude-velocity diagram. The study of the spiral structure is complemented with other tracers: open clusters, Cepheids, methanol masers and H II regions. The well-defined spiral arms are seen to be confined inside the corotation radius, as is often the case in spiral galaxies. We identify a square-shaped sub-structure in the CS map with that predicted by stellar orbits at the 4:1 resonance (four epicycle oscillations in one turn around the galactic centre). The sub-structure is found at the expected radius, based on the known pattern rotation speed and epicycle frequency curve. An inner arm presents an end with strong inwards curvature and intense star formation that we tentatively associate with the region where this arm surrounds the extremity of the bar, as seen in many barred galaxies. Finally, a new arm with concave curvature is found in the Sagitta to Cepheus region of the sky. The observed arms are interpreted in terms of perturbations similar to grooves in the gravitational potential of the disc, produced by crowding of stellar orbits.

Precession and nutation in the eta Carinae binary system: evidence from the X-ray light curve

It is believed that eta Carinae is actually a massive binary system, with the wind-wind interaction responsible for the strong X-ray emission. Although the overall shape of the X-ray light curve can be explained by the high eccentricity of the binary orbit, other features like the asymmetry near periastron passage and the short quasi-periodic oscillations seen at those epochs have not yet been accounted for. In this paper we explain these features assuming that the rotation axis of eta Carinae is not perpendicular to the orbital plane of the binary system. As a consequence, the companion star will face eta Carinae on the orbital plane at different latitudes for different orbital phases and, since both the mass-loss rate and the wind velocity are latitude dependent, they would produce the observed asymmetries in the X-ray flux. We were able to reproduce the main features of the X-ray light curve assuming that the rotation axis of eta Carinae forms an angle of 29 degrees +/- 4 degrees with the axis of the binary orbit. We also explained the short quasi-periodic oscillations by assuming nutation of the rotation axis, with an amplitude of about 5 degrees and a period of about 22 days. The nutation parameters, as well as the precession of the apsis, with a period of about 274 years, are consistent with what is expected from the torques induced by the companion star.

Constraining the orbital orientation of eta Carinae from H Paschen lines

During the past decade, several observational and theoretical works have provided evidence of the binary nature of eta Carinae. Nevertheless, there is still no direct determination of the orbital parameters, and the different current models give contradictory results. The orbit is, in general, assumed to coincide with the Homunculus equator although the observations are not conclusive. Among all systems, eta Car has the advantage that it is possible to observe both the direct emission of line transitions in the central source and its reflection by the Homunculus, which is dependent on the orbital inclination. In this work, we studied the orbital phase-dependent hydrogen Paschen spectra reflected by the south-east lobe of the Homunculus to constrain the orbital parameters of eta Car and determine its inclination with respect to the Homunculus axis. Assuming that the emission excess originates in the wind-wind shock region, we were able to model the latitude dependence of the spectral line profiles. For the first time, we were able to estimate the orbital inclination of eta Car with respect to the observer and to the Homunculus axis. The best fit occurs for an orbital inclination to the line of sight of i similar to 60 degrees +/- 10 degrees, and i* similar to 35 degrees +/- 10 degrees with respect to the Homunculus axis, indicating that the angular momenta of the central object and the orbit are not aligned. We were also able to fix the phase angle of conjunction as similar to -40 degrees, showing that periastron passage occurs shortly after conjunction.

THREE-DIMENSIONAL NUMERICAL SIMULATIONS OF MAGNETIZED WINDS OF SOLAR-LIKE STARS

By means of self-consistent three-dimensional magnetohydrodynamics (MHD) numerical simulations, we analyze magnetized solar-like stellar winds and their dependence on the plasma-beta parameter (the ratio between thermal and magnetic energy densities). This is the first study to perform such analysis solving the fully ideal three-dimensional MHD equations. We adopt in our simulations a heating parameter described by gamma, which is responsible for the thermal acceleration of the wind. We analyze winds with polar magnetic field intensities ranging from 1 to 20 G. We show that the wind structure presents characteristics that are similar to the solar coronal wind. The steady-state magnetic field topology for all cases is similar, presenting a configuration of helmet streamer-type, with zones of closed field lines and open field lines coexisting. Higher magnetic field intensities lead to faster and hotter winds. For the maximum magnetic intensity simulated of 20 G and solar coronal base density, the wind velocity reaches values of similar to 1000 km s(-1) at r similar to 20r(0) and a maximum temperature of similar to 6 x 10(6) K at r similar to 6r(0). The increase of the field intensity generates a larger ""dead zone"" in the wind, i.e., the closed loops that inhibit matter to escape from latitudes lower than similar to 45 degrees extend farther away from the star. The Lorentz force leads naturally to a latitude-dependent wind. We show that by increasing the density and maintaining B(0) = 20 G the system recover back to slower and cooler winds. For a fixed gamma, we show that the key parameter in determining the wind velocity profile is the beta-parameter at the coronal base. Therefore, there is a group of magnetized flows that would present the same terminal velocity despite its thermal and magnetic energy densities, as long as the plasma-beta parameter is the same. This degeneracy, however, can be removed if we compare other physical parameters of the wind, such as the mass-loss rate. We analyze the influence of gamma in our results and we show that it is also important in determining the wind structure.

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.

Circadian Pattern of Wheel-Running Activity of a South American Subterranean Rodent (Ctenomys cf knightii)

Circadian rhythms are regarded as essentially ubiquitous features of animal behavior and are thought to confer important adaptive advantages. However, although circadian systems of rodents have been among the most extensively studied, most comparative biology is restricted to a few related species. In this study, the circadian organization of locomotor activity was studied in the subterranean, solitary north Argentinean rodent, Ctenomys knightii. The genus, Ctenomys, commonly known as Tuco-tucos, comprises more than 50 known species over a range that extends from 12S latitude into Patagonia, and includes at least one social species. The genus, therefore, is ideal for comparative and ecological studies of circadian rhythms. Ctenomys knightii is the first of these to be studied for its circadian behavior. All animals were wild caught but adapted quickly to laboratory conditions, with clear and precise activity-rest rhythms in a light-dark (LD) cycle and strongly nocturnal wheel running behavior. In constant dark (DD), the rhythm expression persisted with free-running periods always longer than 24h. Upon reinstatement of the LD cycle, rhythms resynchronized rapidly with large phase advances in 7/8 animals. In constant light (LL), six animals had free-running periods shorter than in DD, and 4/8 showed evidence of splitting. We conclude that under laboratory conditions, in wheel-running cages, this species shows a clear nocturnal rhythmic organization controlled by an endogenous circadian oscillator that is entrained to 24h LD cycles, predominantly by light-induced advances, and shows the same interindividual variable responses to constant light as reported in other non-subterranean species. These data are the first step toward understanding the chronobiology of the largest genus of subterranean rodents.