High-altitude data assimilation system experiments for the northern summer mesosphere season of 2007

A global numerical weather prediction system is extended to the mesosphere and lower thermosphere (MLT) and used to assimilate high-altitude satellite measurements of temperature, water vapor and ozone from MLS and SABER during May–July 2007. Assimilated temperature and humidity from 100 to 0.001 hPa show minimal biases compared to satellite data and existing analysis fields. Saturation ratios derived diagnostically from these assimilated temperature and water vapor fields at PMC altitudes and latitudes compare well with seasonal variations in PMC frequency measured from the aeronomy of ice in the mesosphere (AIM) satellite. Synoptic maps of these diagnostic saturation ratios correlate geographically with three independent transient mesospheric cloud events observed at midlatitudes by SHIMMER on STPSat-1 and by ground observers during June 2007. Assimilated temperatures and winds reveal broadly realistic amplitudes of the quasi 5-day wave and migrating tides as a function of latitude and height. For example, analyzed winds capture the dominant semidiurnal MLT wind patterns at 55°N in June 2007 measured independently by a meteor radar. The 5-day wave and migrating diurnal tide also modulate water vapor mixing ratios in the polar summer MLT. Possible origins of this variability are discussed.     

Simulating the diurnal cycle of rainfall in global climate models: resolution versus parameterization

The effects of horizontal resolution and the treatment of convection on simulation of the diurnal cycle of precipitation during boreal summer are analyzed in several innovative weather and climate model integrations. The simulations include: season-long integrations of the Non-hydrostatic Icosahedral Atmospheric Model (NICAM) with explicit clouds and convection; year-long integrations of the operational Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts at three resolutions (125, 39 and 16 km); seasonal simulations of the same model at 10 km resolution; and seasonal simulations of the National Center for Atmospheric Research (NCAR) low-resolution climate model with and without an embedded two-dimensional cloud-resolving model in each grid box. NICAM with explicit convection simulates best the phase of the diurnal cycle, as well as many regional features such as rainfall triggered by advancing sea breezes or high topography. However, NICAM greatly overestimates mean rainfall and the magnitude of the diurnal cycle. Introduction of an embedded cloud model within the NCAR model significantly improves global statistics of the seasonal mean and diurnal cycle of rainfall, as well as many regional features. However, errors often remain larger than for the other higher-resolution models. Increasing resolution alone has little impact on the timing of daily rainfall in IFS with parameterized convection, yet the amplitude of the diurnal cycle does improve along with the representation of mean rainfall. Variations during the day in atmospheric prognostic fields appear quite similar among models, suggesting that the distinctive treatments of model physics account for the differences in representing the diurnal cycle of precipitation.     

Lagrangian and Eulerian estimates of circulation in the lee of Kapiti Island,New Zealand

Lagrangian drifters, moored acoustic Doppler current meters and hydrographic observations are combined with wind observations to describe the mean and variable nature of flow around Kapiti Island, New Zealand. Thirteen day-long deployments of up to six Lagrangian drifters show the mean flow is to the southwest, with evidence of stronger flows in the channel separating the island from the mainland, and an island wake in the lee of the island. Vortices in this island wake may be tidally driven. Scaling considerations suggest the flow is strong enough that tidal-generated vortices are shed on each tidal cycle. Both the drifters and mooring data suggest that the d’Urville Current around Kapiti Island has a significant wind-driven component. During north-westerlies, the drifters tend to hug the coast, and south-eastwards flows in the Rauoterangi Channel are accelerated. We suggest the observed correlation is the local expression of a South Taranaki basin scale response to the winds.     

The trace element compositions of S-type granites: evidence for disequilibrium melting and accessory phase entrainment in the source

Within individual plutons, the trace element concentrations in S-type granites generally increase with maficity (total iron and magnesium content and expressed as atomic Fe + Mg in this study); the degree of variability in trace element concentration also expands markedly with the same parameter. The strongly peraluminous, high-level S-type granites of the Peninsular Pluton (Cape Granite Suite, South Africa) are the product of biotite incongruent melting of a metasedimentary source near the base of the crust. Leucogranites within the suite represent close to pure melts from the anatectic source and more mafic varieties represent mixtures of melt and peritectic garnet and ilmenite. Trace elements such as Rb, Ba, Sr and Eu, that are concentrated in reactant minerals in the melting process, show considerable scatter within the granites. This is interpreted to reflect compositional variation in the source. In contrast, elements such as LREE, Zr and Hf, which are concentrated within refractory accessory phases (zircon and monazite), show well-defined negative correlations with increasing SiO₂ and increase linearly with increasing maficity. This is interpreted to reflect coupled co-entrainment of accessory minerals and peritectic phases to the melt: leucocratic rocks cannot have evolved from the more mafic compositions in the suite by a process of fractional crystallisation because in this case they would have inherited the zircon-saturated character of this hypothetical earlier magma. Trace element behaviour of granites from the Peninsular Pluton has been modelled via both equilibrium and disequilibrium trace element melting. In the disequilibrium case, melts are modelled as leaving the source with variable proportions of entrained peritectic phases and accessory minerals, but before the melt has dissolved any accessory minerals. Thus, the trace element signature of the melt is largely inherited from the reactants in the melting reaction, with no contribution from zircon and monazite dissolution. In the equilibrium case, melt leaves the source with entrained crystals, after reaching zircon and monazite saturation. A significant proportion of the rocks of the Peninsular Pluton have trace element concentrations below those predicted by zircon and monazite saturation. In the case of the most leucocratic rocks all compositions are zircon undersaturated; whilst the majority of the most mafic compositions are zircon oversaturated. However, in both cases, zircon is commonly xenocrystic. Thus, the leucocratic rocks represent close to pure melts, which escaped their sources rapidly enough that some very closely match the trace element disequilibrium melting model applied in this study. Zircon dissolution rates allow the residency time for the melt in the source to be conservatively estimated at less than 500 years.     

The optical counterparts to Be/X-ray binaries in the Magellanic Clouds

The fields of eight X-ray sources in the Magellanic Clouds believed to be Be/X-ray binaries have been searched for possible Be-star counterparts. BVR _c and H α CCD imaging was employed to identify early-type emission stars through colour indices and H α fluxes. Spectroscopy of five sources confirms the presence of H α emission in each case. Based on the positional coincidence of emission-line objects with the X-ray sources, we identify Be-star counterparts to the ROSAT sources RX J0032.9-7348, RX J0049.1-7250, RX J0054.9-7226 and RX J0101.0-7206, and to the recently discovered ASCA source AX J0051-722. We confirm the Be star nature of the counterpart to the HEAO1 source H0544-66. In the field of the ROSAT source RX J0051.8-7231 we find that there are three possible counterparts, each showing evidence for H α emission. We find a close double in the error circle of the EXOSAT source EXO 0531.1-6609, each component of which could be a Be star associated with the X-ray source.     

Variations in the broad-band spectra of BL Lac objects: millimetre observations of an X-ray-selected sample

Observations at millimetre wavelengths are presented for a representative sample of 22 X-ray-selected BL Lac objects (XBLs). This sample comprises 19 high-energy cut-off BL Lac objects (HBLs), 1 low-energy cut-off BL Lac object (LBL) and 2 'intermediate' sources. Data for LBLs, which are mostly radio-selected BL Lac objects (RBLs), are taken from the literature. It is shown that the radio–millimetre spectral indices of HBLs     are slightly steeper than those of the LBLs     . A correlation exists between α _5–230 and 230 GHz luminosity. While this correlation could be an artefact of comparing two populations of BL Lac objects with intrinsically different radio properties, it is also consistent with the predictions of existing unified schemes that relate BL Lac objects to Fanaroff–Riley class I radio galaxies.
The HBLs have significantly flatter submillimetre–X-ray spectral indices     than the LBLs     although the two intermediate sources also have intermediate values of α _230–X∼−0.9. It is argued that this difference cannot be explained entirely by the viewing-angle hypothesis and requires a difference in physical-source parameters. The α _230–X values for the HBLs are close to the canonical value found for large samples of radio sources and thus suggest that synchrotron radiation is the mechanism that produces the X-ray emission. As suggested by Padovani & Giommi, the inverse-Compton mechanism is likely to dominate in the LBLs requiring the synchrotron spectra of these sources to steepen or cut off at lower frequencies than those of the HBLs.     

Study of near-surface models for large-eddy simulations of a neutrally stratified atmospheric boundary layer

In large-eddy simulations (LES) of the atmospheric boundary layer (ABL), near-surface models are often used to supplement subgrid-scale (SGS) turbulent stresses when a major fraction of the energetic scales within the surface layer cannot be resolved with the temporal and spatial resolution at hand. In this study, we investigate the performance of both dynamic and non-dynamic eddy viscosity models coupled with near-surface models in simulations of a neutrally stratified ABL. Two near-surface models that are commonly used in LES of the atmospheric boundary layer are considered. Additionally, a hybrid Reynolds- averaged/LES eddy viscosity model is presented, which uses Prandtl’s mixing length model in the vicinity of the surface, and blends in with the dynamic Smagorinsky model away from the surface. Present simulations show that significant portions of the modelled turbulent stresses are generated by the near-surface models, and they play a dominant role in capturing the expected logarithmic wind profile. Visualizations of the instantaneous vorticity field reveal that flow structures in the vicinity of the surface depend on the choice of the near-surface model. Among the three near-surface models studied, the hybrid eddy viscosity model gives the closest agreement with the logarithmic wind profile in the surface layer. It is also observed that high levels of resolved turbulence stresses can be maintained with the so-called canopy stress model while producing good agreement with the logarithmic wind profile.     

Generation, mobilization and crystallization of impact-induced alkali-rich melts in granitic target rocks: Evidence from the Araguainha impact structure, central Brazil

This paper provides important insights into the generation, extraction and crystallization of clast-laden impact melt rocks from the Araguainha impact structure, central Brazil. Despite the mixed nature of the Araguainha target rocks (comprising a 2 km thick sequence of sedimentary rocks and underlying granitic basement), the exposed melt bodies are characterised by an alkali-rich granitic matrix embedding mineral and rock fragments derived only from the target granite. The melt rocks occur in the form of a massive impact melt sheet overlying the eroded central uplift structure, and as melt veins in the granite of the core of the central uplift. Bulk-rock major and trace element data (including platinum group elements) indicate that the precursor melts were generated locally, principally by partial melting of the target granite, without any contribution from the sedimentary sequence or the projectile. The dense network of melt veins was formed in isolation, by selective melting of plagioclase and alkali feldspar within the granite target. Plagioclase and alkali feldspar melted discretely and congruently, producing domains in the matrix of the melt veins, which closely match the stoichiometry of these minerals. The compositionally discrete initial melt phases migrated through a dense network of microfractures before being assembled into larger melt veins. Freezing of the melt veins was substantially fast, and the melt components were quenched in the form of alkali-feldspar and plagioclase schlieren in the matrix of the melt veins. The overlying impact melt rock is, in contrast, characterised by a granophyric matrix consisting of albite, sanidine, quartz, biotite and chlorite. In this case, melt components appear to have been more mobile and to have mixed completely to form a granitic parental melt. We relate the melting of the minerals to post-shock temperatures that exceeded the melting point of feldspars.