Bridging the Transition from Mesoscale to Microscale Turbulence in Numerical Weather Prediction Models

With a focus towards developing multiscale capabilities in numerical weather prediction models, the specific problem of the transition from the mesoscale to the microscale is investigated. For that purpose, idealized one-way nested mesoscale to large-eddy simulation (LES) experiments were carried out using the Weather Research and Forecasting model framework. It is demonstrated that switching from one-dimensional turbulent diffusion in the mesoscale model to three-dimensional LES mixing does not necessarily result in an instantaneous development of turbulence in the LES domain. On the contrary, very large fetches are needed for the natural transition to turbulence to occur. The computational burden imposed by these long fetches necessitates the development of methods to accelerate the generation of turbulence on a nested LES domain forced by a smooth mesoscale inflow. To that end, four new methods based upon finite amplitude perturbations of the potential temperature field along the LES inflow boundaries are developed, and investigated under convective conditions. Each method accelerated the development of turbulence within the LES domain, with two of the methods resulting in a rapid generation of production and inertial range energy content associated to microscales that is consistent with non-nested simulations using periodic boundary conditions. The cell perturbation approach, the simplest and most efficient of the best performing methods, was investigated further under neutral and stable conditions. Successful results were obtained in all the regimes, where satisfactory agreement of mean velocity, variances and turbulent fluxes, as well as velocity and temperature spectra, was achieved with reference non-nested simulations. In contrast, the non-perturbed LES solution exhibited important energy deficits associated to a delayed establishment of fully-developed turbulence. The cell perturbation method has negligible computational cost, significantly accelerates the generation of realistic turbulence, and requires minimal parameter tuning, with the necessary information relatable to mean inflow conditions provided by the mesoscale solution.     

Nesting Turbulence in an Offshore Convective Boundary Layer Using Large-Eddy Simulations

The applicability of the one-way nesting technique for numerical simulations of the heterogeneous atmospheric boundary layer using the large-eddy simulation (LES) framework of the Weather Research and Forecasting model is investigated. The focus of this study is on LES of offshore convective boundary layers. Simulations were carried out using two subgrid-scale models (linear and non-linear) with two different closures [diagnostic and prognostic subgrid-scale turbulent kinetic energy (TKE) equations]. We found that the non-linear backscatter and anisotropy model with a prognostic subgrid-scale TKE equation is capable of providing similar results when performing one-way nested LES to a stand-alone domain having the same grid resolution but using periodic lateral boundary conditions. A good agreement is obtained in terms of velocity shear and turbulent fluxes, while velocity variances are overestimated. A streamwise fetch of 14 km is needed following each domain transition in order for the solution to reach quasi-stationary results and for the velocity spectra to generate proper energy content at high wavelengths, however, a pile-up of energy is observed at the low-wavelength portion of the spectrum on the first nested domain. The inclusion of a second nest with higher resolution allows the solution to reach effective grid spacing well within the Kolmogorov inertial subrange of turbulence and develop an appropriate energy cascade that eliminates most of the pile-up of energy at low wavelengths. Consequently, the overestimation of velocity variances is substantially reduced and a considerably better agreement with respect to the stand-alone domain results is achieved.     

Geochemical evidence for sundering of the West Mariana arc in miocene ash from the Parece Vela Basin

Glass and mineral fragments from discrete volcanic ash layers were sampled from DSDP/IPOD Site 450 in the Parece Vela Basin, Philippine Sea and analyzed by electron microprobe. The ashes are interpreted as eruptive products of the adjacent West Mariana arc system between 25 and 14 Ma B.P., and have compositions between basaltic andesite and rhyolite, and rarely, boninite. ‘Continuous’ chemical trends appear to reflect mixing of mafic and silicic magmas. ‘Discontinuous’ trends between these end-members are relatively few, and are consistent with ‘liquid lines’ produced by fractional crystallization. Andesitic tephra become progressively richer in MgO and CaO through the middle Miocene, while boninite appears towards the end of the sequence, between 14 and 15 Ma B.P. Coeval rhyolitic glasses become richer in K₂O and Na₂O, with maximum concentrations at about 15 Ma B.P. Chronologic changes in fractionation type and composition of parent magmas are interpreted to reflect the subaerial volcanic evolution of the West Mariana arc. The appearance of boninite is believed to signal early stages of arc sundering, and corresponds temporally with regional uplift of the sea floor above the carbonate compensation depth, precursor to a new pulse of back-arc spreading.     

A detailed study of the 5-Hz quasi-periodic oscillations in the bright X-ray transient and black hole candidate GRS 1739–278

We present a detailed study of the 5-Hz quasi-periodic oscillation (QPO) recently discovered in the bright X-ray transient and black hole candidate (BHC) GRS     (Borozdin & Trudolyubov) during a Rossi X-ray Timing Explorer observation taken on 1996 March 31. In total 6.6 ksec of on-source data were obtained, divided in two data sets of 3.4 and 3.2 ksec which were separated by ∼2.6 ksec. The 5-Hz QPO was only present during the second data set. The QPO increased in strength from below 2 per cent rms amplitude for photon energies below 4 keV to ∼5 per cent rms amplitude for energies above 10 keV. The soft QPO photons (below 5 keV) lagged the hard ones (above 10 keV) by almost 1.5 rad. Besides the QPO fundamental, its first overtone was detected. The strength of the overtone increased with photon energy (from < 2 per cent rms below 5 keV to ∼8 per cent rms above 10 keV). Although limited statistics did not allow for an accurate determination of the lags of the first overtone, indications are that also for this QPO the soft photons lagged the hard ones. When the 5-Hz QPO was not detected (i.e., during the first part of the observation), a broad noise component was found for photon energies below 10 keV but it became almost a true QPO (with a Q value of ∼1.9) above that energy, with a frequency of ∼3 Hz. Its hard photons preceded the soft ones in a way reminiscent of the 5-Hz QPO, strongly suggesting that both features are physically related. We discuss our finding in the framework of low-frequency QPOs and their properties in BHCs.     

Quantifying the influence of sill intrusion on the thermal evolution of organic‐rich sedimentary rocks in nonvolcanic passive margins: an example from ODP 210‐1276, offshore Newfoundland,Canada

Intrusive magmatism is an integral and understudied component in both volcanic and nonvolcanic passive margins. Here, we investigate the thermal effects of widespread (ca. 20 000 km²) intrusive magmatism on the thermal evolution of organic‐rich sedimentary rocks on the nonvolcanic Newfoundland passive margin. ODP 210‐1276 (45.41°N, 44.79°W) intersects two sills: an older, upper sill and a younger, lower sill that are believed to correspond to the high amplitude ‘U‐reflector’ observed across the Newfoundland Basin. A compilation of previous work collectively provides; (1) emplacement depth constraints, (2) vitrinite reflectance data and (3) ⁴⁰Ar/³⁹Ar dates. Collectively, these data sets provide a unique opportunity to model the conductive cooling of the sills and how they affect thermal maturity of the sedimentary sequence. A finite differences method was used to model the cooling of the sills, with the model outputs then being entered into the EASY%R_o vitrinite reflectance model. The modelled maturation profile for ODP 210‐1276 shows a significant but localized effect on sediment maturity as a result of the intrusions. Our results suggest that even on nonvolcanic margins, intrusive magmatism can significantly influence the thermal evolution in the vicinity of igneous intrusions. In addition, the presence of widespread sills on nonvolcanic passive margins such as offshore Newfoundland may be indicative of regional‐scale thermal perturbations that should be considered in source rock maturation studies.     

Holocene fossil Habrotrocha angusticollis (Bdelloidea: Rotifera) in North America

Holocene fossil shells of the rotifer, Habrotrocha angusticollus (Bdelloidea: Rotifera) is reported from two peatlands in northern Ontario, Canada. H. angusticollis is a common component of the microfauna in Sphagnum peatlands and other wet mossy habitats. Our knowledge of the distribution and ecology of H. angusticollis is limited and this paper is the first detailed report of the shells as fossils in North America. Fossil shells of H. angusticollis may prove to be a valuable paleoecological indicator in peat deposits once Quaternary paleoecologists learn to recognize them and neoecologists extend their surveys to peatlands.     

Numerical simulations of pinhole and single-mode fibre spatial filters for optical interferometers

We use a numerical simulation to investigate the effectiveness of pinhole spatial filters for optical/IR interferometers and to compare them with single-mode optical fibre spatial filters and interferometers without spatial filters. We show that fringe visibility measurements in interferometers containing spatial filters are much less affected by changing seeing conditions than equivalent measurements without spatial filters. This reduces visibility calibration uncertainties, and hence can reduce the need for frequent observations of separate astronomical sources for calibration of visibility measurements. We also show that spatial filters can increase the signal-to-noise ratios (SNRs) of visibility measurements and that pinhole filters give SNRs within 17 per cent of the values obtained with single-mode fibres for aperture diameters up to 3 r ₀. Given the simplicity of the use of pinhole filters we suggest that it represents a competitive, if not optimal, technique for spatial filtering in many current and next generation interferometers.     

Torques and instabilities in intermediate polars

The observed spin periods of the primaries in intermediate polars require the presence of torques additional to the standard expectation of material and magnetic torques from an accretion disc. In an extension of earlier work by Wickramasinghe, Wu and Ferrario, from rates of mass transfer deduced from optical and X-ray fluxes, and assuming surface magnetic fields similar to those observed in the most rapidly rotating low mass main sequence stars, it is shown that the magnetohydrodynamic torque between the partially shielded primary and the secondary is sufficient to account for the observed spin periods. It is further found that there is a range of magnetic moments and mass transfer rates in which synchronized rotation of the primary can occur even though it possesses an accretion disc. This may account for the structures deduced for V795 Her and V2051 Oph.An analysis of the observed outbursts (or lack of) in the intermediate polars shown that wheres many systems are explainable in terms of the standard theory of -discs, in some systems the magnetic structure is supressing expected outbursts and instabilities that are as yet not understood appear in their place.