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.

Proper evaluation of the external convective heat transfer for the thermal analysis of cool roofs

Cool materials are characterized by high solar reflectance and high thermal emittance; when applied to the external surface of a roof, they make it possible to limit the amount of solar irradiance absorbed by the roof, and to increase the rate of heat flux emitted by irradiation to the environment, especially during nighttime. However, a roof also releases heat by convection on its external surface; this mechanism is not negligible, and an incorrect evaluation of its entity might introduce significant inaccuracy in the assessment of the thermal performance of a cool roof, in terms of surface temperature and rate of heat flux transferred to the indoors. This issue is particularly relevant in numerical simulations, which are essential in the design stage, therefore it deserves adequate attention. In the present paper, a review of the most common algorithms used for the calculation of the convective heat transfer coefficient due to wind on horizontal building surfaces is presented. Then, with reference to a case study in Italy, the simulated results are compared to the outcomes of a measurement campaign. Hence, the most appropriate algorithms for the convective coefficient are identified, and the errors deriving by an incorrect selection of this coefficient are discussed.

Re-examining the U.K.'s greatest tornado outbreak: forecasting the limited extent of tornadoes along a cold front

On 23 November 1981, a strong cold front swept across the U.K., producing tornadoes from the west to the east coasts. An extensive campaign to collect tornado reports by the Tornado and Storm Research Organisation (TORRO) resulted in 104 reports, the largest U.K. outbreak. The front was simulated with a convection-permitting numerical model down to 200-m horizontal grid spacing to better understand its evolution and meteorological environment. The event was typical of tornadoes in the U.K., with convective available potential energy (CAPE) less than 150 J kg-1, 0-1-km wind shear of 10-20 m s-1, and a narrow cold-frontal rainband forming precipitation cores and gaps. A line of cyclonic absolute vorticity existed along the front, with maxima as large as 0.04 s-1. Some hook-shaped misovortices bore kinematic similarity to supercells. The narrow swath along which the line was tornadic was bounded on the equatorward side by weak vorticity along the line and on the poleward side by zero CAPE, enclosing a region where the environment was otherwise favorable for tornadogenesis. To determine if the 104 tornado reports were plausible, first possible duplicate reports were eliminated, resulting in as few as 58 tornadoes to as many as 90. Second, the number of possible parent misovortices that may have spawned tornadoes is estimated from model output. The number of plausible tornado reports in the 200-m grid-spacing domain was 22 and as many as 44, whereas the model simulation was used to estimate 30 possible parent misovortices within this domain. These results suggest that 90 reports was plausible.

Holocene Climatic Optimum centennial-scale paleoceanography in the NE Aegean (Mediterranean sea)

Combined micropaleontological and geochemical analyses of the high-sedimentation gravity core M-4G provided new centennial-scale paleoceanographic data for sapropel S1 deposition in the NE Aegean Sea during the Holocene Climatic Optimum. Sapropel layer S1a (10.2–8.0 ka) was deposited in dysoxic to oxic bottom waters characterized by a high abundance of benthic foraminiferal species tolerating surface sediment and/or pore water oxygen depletion (e.g., Chilostomella mediterranensis, Globobulimina affinis), and the presence of Uvigerina mediterranea, which thrives in oxic mesotrophic-eutrophic environments. Preservation of organic matter (OM) is inferred based on high organic carbon as well as loliolide and isololiolide contents, while the biomarker record and the abundances of eutrophic planktonic foraminifera document enhanced productivity. High inputs of terrigenous OM are attributed to north Aegean borderland riverine inputs. Both alkenone-based sea surface temperatures (SSTs) and δO18G. bulloides records indicate cooling at 8.2 ka (S1a) and ~7.8 ka (S1 interruption). Sapropelic layer S1b (7.7–6.4 ka) is characterized by rather oxic conditions; abundances of foraminiferal species tolerant to oxygen depletion are very low compared with the U. mediterranea rise. Strongly fluctuating SSTs demonstrate repeated cooling and associated dense water formation, with a major event at 7.4 ka followed by cold spells at 7.0, 6.8, and 6.5 ka. The prominent rise of the carbon preference index within the S1b layer indicates the delivery of less degraded terrestrial OM. The increase of algal biomarkers, labile OM-feeding foraminifera and eutrophic planktonic species pinpoints an enhanced in situ marine productivity, promoted by more efficient vertical convection due to repeated cold events. The associated contributions of labile marine OM along with fresher terrestrial OM inputs after ~7.7 ka imply sources alternative/additional to the north Aegean riverine borderland sources for the influx of organic matter in the south Limnos Basin, plausibly related to the inflow of highly productive Marmara/Black Sea waters.

Verification of the role of the low level jets in Amazon squall lines

In this study we present a climatology of the Amazon squall lines (ASLs), between the years 2000 and 2008, using satellite imagery and European Centre for Medium-Range Weather Forecasts (ECMWF) reanalyses. The ASLs we are interested in are typically formed along the northern coast of Brazil and sometimes propagate for long distances inland. Results show that, on average, an ASL occurs every 2 days. ASLs are more frequent between April and June and less frequent between October and November. The years of 2005 and 2006 showed 25% more cases than the other years. This might be related to an increase of the Atlantic sea surface temperature. Of the total number of ASL cases, 54% propagated less than 170 km, 26% propagated between 170 and 400 km, and 20% propagated more than 400 km. We also studied the occurrence of low level jets (LLJs) associated with the coastal ASLs. Although LLJs are always present in the environment before the formation of the ASL and even on days without ASL cases, important differences were found, mainly related to the LLJ depths. (C) 2010 Elsevier B.V. All rights reserved.

Mechanisms Associated with Large Daily Rainfall Events in Northeast Brazil

The mechanisms resulting in large daily rainfall events in Northeast Brazil are analyzed using data filtering to exclude periods longer than 30 days. Composites of circulation fields that include all independent events do not reveal any obvious forcing mechanisms as multiple patterns contribute to Northeast Brazil precipitation variability. To isolate coherent patterns, subsets of events are selected based on anomalies that precede the Northeast Brazil precipitation events at different locations. The results indicate that at 10 degrees S, 40 degrees W, the area of lowest annual rainfall in Brazil, precipitation occurs mainly in association with trailing midlatitude synoptic wave trains originating in either hemisphere. Closer to the equator at 5 degrees S, 37.5 degrees W, an additional convection precursor is found to the west, with a spatial structure consistent with that of a Kelvin wave. Although these two sites are located within only several hundred kilometers of each other and the midlatitude patterns that induce precipitation appear to be quite similar, the dates on which large precipitation anomalies occur at each location are almost entirely independent, pointing to separate forcing mechanisms.

Cloud streets and land-water interactions in the Amazon

Cloud streets are common feature in the Amazon Basin. They form from the combination of the vertical trade wind stress and moist convection. Here, satellite imagery, data collected during the COBRA-PARA (Caxiuan Observations in the Biosphere, River and Atmosphere of Para) field campaign, and high resolution modeling are used to understand the streets` formation and behavior. The observations show that the streets have an aspect ratio of about 3.5 and they reach their maximum activity around 15:00 UTC when the wind shear is weaker, and the convective boundary layer reaches its maximum height. The simulations reveal that the cloud streets onset is caused by the local circulations and convection produced at the interfaces between forest and rivers of the Amazon. The satellite data and modeling show that the large rivers anchor the cloud streets producing a quasi-stationary horizontal pattern. The streets are associated with horizontal roll vortices parallel to the mean flow that organizes the turbulence causing advection of latent heat flux towards the upward branches. The streets have multiple warm plumes that promote a connection between the rolls. These spatial patterns allow fundamental insights on the interpretation of the Amazon exchanges between surface and atmosphere with important consequences for the climate change understanding.

Interaction of equatorial waves through resonance with the diurnal cycle of tropical heating

In this work, a new theoretical mechanism is presented in which equatorial Rossby and inertio-gravity wave modes may interact with each other through resonance with the diurnal cycle of tropical deep convection. We have adopted the two-layer incompressible equatorial primitive equations forced by a parametric heating that roughly represents deep convection activity in the tropical atmosphere. The heat source was parametrized in the simplest way according to the hypothesis that it is proportional to the lower-troposphere moisture convergence, with the background moisture state function mimicking the structure of the ITCZ. In this context, we have investigated the possibility of resonant interaction between equatorially trapped Rossby and inertio-gravity modes through the diurnal cycle of the background moisture state function. The reduced dynamics of a single resonant duo shows that when this diurnal variation is considered, a Rossby wave mode can undergo significant amplitude modulations when interacting with an inertio-gravity wave mode, which is not possible in the context of the resonant triad non-linear interaction. Therefore, the results suggest that the diurnal variation of the ITCZ can be a possible dynamical mechanism that leads the Rossby waves to be significantly affected by high frequency modes.

The energy cycle and structural evolution of cyclones over southeastern South America in three case studies

In this paper, the Lorenz energy cycle over a limited area was applied for three cyclones with different origins and evolutions, where each of them was formed in an important cyclogenetic region near southeastern South America. The synoptic conditions and energetics were analyzed during each system`s life cycle and showed important relationships between their energy cycle and the evolution of their vertical structure. In the case of the weak baroclinic cyclone which formed on Brazil`s south-southeastern coast, the analysis showed that it originated through a midlevel cutoff low with contribution from barotropic instability. Its evolution would indicate potential transition to a hybrid system if the convective activity were stronger. The system that occurred in the La Plata River mouth had features of an oceanic bomb-type cyclogenesis and showed an important contribution from the available potential energy generation term through the latent heat release by the convection. Meanwhile, the system of the southern Argentina coast presented a classical baroclinic development of extratropical cyclogenesis in the energy cycle, from the wave amplification up to the final occlusion of the associated frontal system. These analyses revealed that the development of some cyclones that occur in eastern South America can present different mechanisms that are not related to the classical extratropical cyclogenesis.

Resonant Wave Interactions in the Presence of a Diurnally Varying Heat Source

Resonant interactions among equatorial waves in the presence of a diurnally varying heat source are studied in the context of the diabatic version of the equatorial beta-plane primitive equations for a motionless, hydrostatic, horizontally homogeneous and stably stratified background atmosphere. The heat source is assumed to be periodic in time and of small amplitude [i.e., O(epsilon)] and is prescribed to roughly represent the typical heating associated with deep convection in the tropical atmosphere. In this context, using the asymptotic method of multiple time scales, the free linear Rossby, Kelvin, mixed Rossby-gravity, and inertio-gravity waves, as well as their vertical structures, are obtained as leading-order solutions. These waves are shown to interact resonantly in a triad configuration at the O(e) approximation, and the dynamics of these interactions have been studied in the presence of the forcing. It is shown that for the planetary-scale wave resonant triads composed of two first baroclinic equatorially trapped waves and one barotropic Rossby mode, the spectrum of the thermal forcing is such that only one of the triad components is resonant with the heat source. As a result, to illustrate the role of the diurnal forcing in these interactions in a simplified fashion, two kinds of triads have been analyzed. The first one refers to triads composed of a k = 0 first baroclinic geostrophic mode, which is resonant with the stationary component of the diurnal heat source, and two dispersive modes, namely, a mixed Rossby-gravity wave and a barotropic Rossby mode. The other class corresponds to triads composed of two first baroclinic inertio-gravity waves in which the highest-frequency wave resonates with a transient harmonic of the forcing. The integration of the asymptotic reduced equations for these selected resonant triads shows that the stationary component of the diurnal heat source acts as an ""accelerator"" for the energy exchanges between the two dispersive waves through the excitation of the catalyst geostrophic mode. On the other hand, since in the second class of triads the mode that resonates with the forcing is the most energetically active member because of the energy constraints imposed by the triad dynamics, the results show that the convective forcing in this case is responsible for a longer time scale modulation in the resonant interactions, generating a period doubling in the energy exchanges. The results suggest that the diurnal variation of tropical convection might play an important role in generating low-frequency fluctuations in the atmospheric circulation through resonant nonlinear interactions.

Cloud condensation nuclei from biomass burning during the Amazonian dry-to-wet transition season

Aircraft measurements of cloud condensation nuclei (CCN) during the Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA) were conducted over the Southwestern Amazon region in September-October 2002, to emphasize the dry-to-wet transition season. The CCN concentrations were measured for values within the range 0.1-1.0% of supersaturation. The CCN concentration inside the boundary layer revealed a general decreasing trend during the transition from the end of the dry season to the onset of the wet season. Clean and polluted areas showed large differences. The differences were not so strong at high levels in the troposphere and there was evidence supporting the semi-direct aerosol effect in suppressing convection through the evaporation of clouds by aerosol absorption. The measurements also showed a diurnal cycle following biomass burning activity. Although biomass burning was the most important source of CCN, it was seen as a source of relatively efficient CCN, since the increase was significant only at high supersaturations.

Precipitation diurnal cycle and summer climatology assessment over South America: An evaluation of Regional Climate Model version 3 simulations

Regional Climate Model version 3 (RegCM3) simulations of 17 summers (1988-2004) over part of South America south of 5 degrees S were evaluated to identify model systematic errors. Model results were compared to different rainfall data sets (Climate Research Unit (CRU), Climate Prediction Center (CPC), Global Precipitation Climatology Project (GPCP), and National Centers for Environmental Prediction (NCEP) reanalysis), including the five summers mean (1998-2002) precipitation diurnal cycle observed by the Tropical Rainfall Measuring Mission (TRMM)-Precipitation Radar (PR). In spite of regional differences, the RegCM3 simulates the main observed aspects of summer climatology associated with the precipitation (northwest-southeast band of South Atlantic Convergence Zone (SACZ)) and air temperature (warmer air in the central part of the continent and colder in eastern Brazil and the Andes Mountains). At a regional scale, the main RegCM3 failures are the underestimation of the precipitation in the northern branch of the SACZ and some unrealistic intense precipitation around the Andes Mountains. However, the RegCM3 seasonal precipitation is closer to the fine-scale analyses (CPC, CRU, and TRMM-PR) than is the NCEP reanalysis, which presents an incorrect north-south orientation of SACZ and an overestimation of its intensity. The precipitation diurnal cycle observed by TRMM-PR shows pronounced contrasts between Tropics and Extratropics and land and ocean, where most of these features are simulated by RegCM3. The major similarities between the simulation and observation, especially the diurnal cycle phase, are found over the continental tropical and subtropical SACZ regions, which present afternoon maximum (1500-1800 UTC) and morning minimum (0900-1200 UTC). More specifically, over the core of SACZ, the phase and amplitude of the simulated precipitation diurnal cycle are very close to the TRMM-PR observations. Although there are amplitude differences, the RegCM3 simulates the observed nighttime rainfall in the eastern Andes Mountains, over the Atlantic Ocean, and also over northern Argentina. The main simulation deficiencies are found in the Atlantic Ocean and near the Andes Mountains. Over the Atlantic Ocean the convective scheme is not triggered; thus the rainfall arises from the grid-scale scheme and therefore differs from the TRMM-PR. Near the Andes, intense (nighttime and daytime) simulated precipitation could be a response of an incorrect circulation and topographic uplift. Finally, it is important to note that unlike most reported bias of global models, RegCM3 does not trigger the moist convection just after sunrise over the southern part of the Amazon.

A numerical investigation of the atmosphere-ocean thermal contrast over the coastal upwelling region of Cabo Frio, Brazil

An one-dimensional atmospheric second order closure model, coupled to an oceanic mixed layer model, is used to investigate the short term variation of the atmospheric and oceanic boundary layers in the coastal upwelling area of Cabo Frio, Brazil (23 degrees S, 42 degrees 08`W). The numerical simulations were carried out to evaluate the impact caused by the thermal contrast between atmosphere and ocean on the vertical extent and other properties of both atmospheric and oceanic boundary layers. The numerical simulations were designed taking as reference the observations carried out during the passage of a cold front that disrupted the upwelling regime in Cabo Frio in July of 1992. The simulations indicated that in 10 hours the mechanical mixing, sustained by a constant background flow of 10 in s(-1), increases the atmospheric boundary layer in 214 in when the atmosphere is initially 2 K warmer than the ocean (positive thermal contrast observed during upwelling regime). For an atmosphere initially -2 K colder than the ocean (negative thermal contrast observed during passage of the cold front), the incipient thermal convection intensifies the mechanical mixing increasing the vertical extent of the atmospheric boundary layer in 360 in. The vertical evolution of the atmospheric boundary layer is consistent with the observations carried out in Cabo Frio during upwelling condition. When the upwelling is disrupted, the discrepancy between the simulated and observed atmospheric boundary layer heights in Cabo Frio during July of 1992 increases considerably. During the period of 10 hours, the simulated oceanic mixed layer deepens 2 in and 5.4 in for positive and negative thermal contrasts of 2 K and -2 K, respectively. In the latter case, the larger vertical extent of the oceanic mixed layer is due to the presence of thermal convection in the atmospheric boundary layer, which in turn is associated to the absence of upwelling caused by the passage of cold fronts in Cabo Frio.

Revisiting 2D numerical models for the 19th century outbursts of eta Carinae

We present here new results of two-dimensional hydrodynamical simulations of the eruptive events of the 1840s (the great) and the 1890s (the minor) eruptions suffered by the massive star eta Carinae (Car). The two bipolar nebulae commonly known as the Homunculus and the little Homunculus (LH) were formed from the interaction of these eruptive events with the underlying stellar wind. We assume here an interacting, non-spherical multiple-phase wind scenario to explain the shape and the kinematics of both Homunculi, but adopt a more realistic parametrization of the phases of the wind. During the 1890s eruptive event, the outflow speed decreased for a short period of time. This fact suggests that the LH is formed when the eruption ends, from the impact of the post-outburst eta Car wind (that follows the 1890s event) with the eruptive flow (rather than by the collision of the eruptive flow with the pre-outburst wind, as claimed in previous models; Gonzalez et al.). Our simulations reproduce quite well the shape and the observed expansion speed of the large Homunculus. The LH (which is embedded within the large Homunculus) becomes Rayleigh-Taylor unstable and develop filamentary structures that resemble the spatial features observed in the polar caps. In addition, we find that the interior cavity between the two Homunculi is partially filled by material that is expelled during the decades following the great eruption. This result may be connected with the observed double-shell structure in the polar lobes of the eta Car nebula. Finally, as in previous work, we find the formation of tenuous, equatorial, high-speed features that seem to be related to the observed equatorial skirt of eta Car.

(40)Ar/(39)Ar ages and Sr-Nd-Pb-Os geochemistry of CAMP tholeiites from Western Maranhao basin (NE Brazil)

The Central Atlantic Magmatic Province (CAMP), emplaced at the Triassic-Jurassic (T-J) boundary (-200 Ma), is among the largest igneous provinces on Earth. The Maranhao basin in NE Brazil is located around 700 km inland and 2000 km from the site of the earliest Pangea disruption. The CAMP tholeiites occur only in the western part of the basin and have been described as low and high-Ti. Here we document the occurrence of two sub-groups among the high-Ti tholeiites in the Western Maranhao basin. The major and trace elements and the Sr-Nd-Pb isotopic ratios define three chemical groups corresponding to the low-Ti (TiO(2)< 1.3 wt.%), high-Ti (TiO(2)-2.0 wt.%) and evolved high-Ti (TiO(2 >)3 wt.%) western Maranhao basin tholeiites (WMBT). The new (40)Ar/(39)Ar plateau ages obtained on plagioclase separates for high-Ti (199.7 +/- 2.4 Ma) and evolved high-Ti WMBT (197.2 +/- 0.5 Ma and 198.2 +/- 0.6 Ma) are indistinguishable and identical to those of previously analyzed low-Ti WMBT (198.5 +/- 0.8 Ma) and to the mean (40)Ar/(39)Ar age of the CAMP (199 +/- 2.4 Ma). We also present the first Re-Os isotopic data for CAMP basalts. The low and high-Ti samples display mantle-like initial ((187)Os/(188)Os)(i) ranging from 0.1267 to 0.1299, while the evolved high-Ti samples are more radiogenic (((187)Os/ (188)Os)(i) up to 0.184) We propose that the high-Ti WMBT were derived from the sub-lithospheric asthenosphere, and contaminated during ascent by interaction with the subcontinental lithospheric mantle (SCLM). The evolved high-Ti WMBT were derived from the same asthenospheric source but experienced crustal contamination. The chemical characteristics of the low-Ti group can be explained by partial melting of the most fertile portions of the SCLM metasomatized during paleo-subduction. Alternatively, the low-Ti WMBT could be derived from the sub-lithospheric asthenosphere but the resulting melts may have undergone contamination by the SCLM. The occurrences of high-Ti basalts are apparently not restricted to the area of initial continental disruption which may bring into question previous interpretations such as those relating high-Ti CAMP magmatism to the initiation of Atlantic ridge spreading or as the expression of a deep mantle plume. We propose that the CAMP magmatism in the Maranhao basin may be attributed to local hotter mantle conditions due to the combined effects of edge-driven convection and large-scale mantle warming under the Pangea supercontinent. The involvement of a mantle-plume with asthenosphere-like isotopic characteristics cannot be ruled out either as one of the main source components of the WMBT or as a heat supplier. (C) 2010 Elsevier BM. All rights reserved.