A comparative study of the performance of various vertical discretization schemes

Summary Most finite-difference numerical weather prediction models employ vertical discretizations that are staggered, and are low-order (usually second-order) approximations for the important terms such as the derivation of the geopotential from the hydrostatic equation, and the calculation of the vertically integrated divergence. In a sigma-coordinate model the latter is used for computing both the surface pressure change and the vertical velocity. All of the above-mentioned variables can diminish the accuracy of the forecast if they are not calculated accurately, and can have an impact on related quantities such as precipitation.In this study various discretization schemes in the vertical are compared both in theory and in practice. Four different vertical grids are tested: one unstaggered and three staggered (including the widely-used Lorenz grid). The comparison is carried out by assessing the accuracy of the grids using vertical numerics that range from second-order up to sixth-order.The theoretical part of the study examines how faithfully each vertical grid reproduces the vertical modes of the governing equations linearized with a basic state atmosphere. The performance of the grids is evaluated for 2nd, 4th and 6th-order numerical schemes based on Lagrange polynomials, and for a 6th-ordercompact scheme.Our interpretation of the results of the theoretical study is as follows. The most important result is that the order of accuracy employed in the numerics seems to be more significant than the choice of vertical grid. There are differences between the grids at second-order, but these differences effectively vanish as the order of accuracy increases. The sixth-order schemes all produce very accurate results with the grids performing equally well, and with the compact scheme significantly outperforming the Lagrange scheme. A second major result is that for the number of levels typically used in current operational forecast models, second-order schemes (which are used almost universally) all appear to be relatively poor, for other than the lowest modes.The theoretical claims were confirmed in practice using a large number (100) of forecasts with the Australian Bureau of Meteorology Research Centre's operational model. By comparing test model forecasts using the four grids and the different orders of numerics with very high resolution control model forecasts, the results of the theoretical study seem to be corroborated.With 8 Figures     

The effect of surface roughness on the transmission of microwave radiation through a planetary surface

To account for surface roughness, the transmission of microwave radiation through a planetary surface to an observer is treated by a Monte Carlo technique. Sizable effects are found near the limb of the planet, and they should be included in analyses of high-resolution observations and high-precision integrated disk observations.     

Characterization of bacterial community structure in vestimentiferan tubeworm Ridgeia piscesae trophosomes

Terminal‐restriction fragment length polymorphisms (T‐RFLPs) were used to address the question whether there are multiple phylotypes of bacteria within the trophosome of the vestimentiferan tubeworm Ridgeia piscesae in addition to the known endosymbiont. Clone libraries were constructed to aid in the identification of the additional bacterial phylotypes. Individual R. piscesae specimens were collected from Juan de Fuca Ridge, Axial Caldera and Explorer Ridge, Magic Mountain. Clone library analyses revealed only one bacterial phylotype (expected R. piscesae endosymbiont) within the Explorer Ridge trophosomes. However, the Axial Caldera clone library revealed five operational taxonomic units (OTUs). Three of the resulting phylotypes detected (designated RAE) were putative thioautotrophic symbionts within the γ‐Proteobacteria which belonged to the anticipated endosymbiont group (RAE OTU 1), the Maorithyas hadalis symbiont II group (RAE OTU 3), and the Halothiobacillus group (RAE OTU 4). The remaining two phylotypes were most likely opportunistically derived. RAE OTU 5 was an α‐Proteobacterium within the Roseobacter group and RAE OTU 2 was within the Cytophaga–Flavobacterium–Bacteroidetes. T‐RFLP analyses revealed that all 15 trophosomes evaluated contained at least RAE OTUs 1 and 4. Five of the 15 trophosomes contained at least RAE OTUs 1, 3, and 4 and two of the trophosomes contained at least five RAE OTUs. Cluster analysis of the T‐RFLP data revealed three distinct clusters. The number of taxa within the trophosome may be an indication of the general health of the tubeworm. This study strongly suggests that we have discovered and identified more than one bacterial phylotype within the trophosome of vestimentiferan R. piscesae.     

The effects of vertical shear and stratification on stationary rossby waves

Abstract

The generation of stationary Rossby waves by sources of potential vorticity in a westerly flow is examined here in the context of a two-layer, quasi-geostrophic, β-plane model. The response in each layer consists of a combination of a barotropic Rossby wave disturbance that extends far downstream of the source, and a baroclinic disturbance which is evanescent or wave-like in character, depending on the shear and degree of stratification. Contributions from each of these modes in each layer are strongly dependent on the basic flows in each layer; the degree of stratification; and the depths of the two layers. The lower layer response is dominated by an evanescent baroclinic mode when the upper layer westerlies are much larger than those in the lower layer. In this case, weak stationary Rossby waves of large wavelengths are confined to the upper layer and the disturbance in the lower layer is confined to the source region.

Increasing the upper layer flow (with the lower layer flow fixed) increases the Rossby wavelength and decreases the amplitude. Decreasing the lower layer flow (with the upper layer flow fixed) decreases the wavelength and increases the amplitude. Stratification increases the contribution from the barotropic wave-like mode and causes the response to be confined to the lower layer.

The finite amplitude response to westerly flow over two sources of potential vorticity is also considered. In this case stationary Rossby waves induced by both sources interact to reinforce or diminish the downstream wave pattern depending on the separation distance of the sources relative to the Rossby wavelength. For fixed separation distance, enhancement of the downstreatm Rossby waves will only occur for a narrow range of flow variables and stratification.     

WASP-1b and WASP-2b: two new transiting exoplanets detected with SuperWASP and SOPHIE

We have detected low-amplitude radial-velocity variations in two stars, USNO-B1.0 1219–0005465 (GSC  02265–00107 = WASP–1  ) and USNO-B1.0 0964–0543604 (GSC  00522–01199 = WASP–2  ). Both stars were identified as being likely host stars of transiting exoplanets in the 2004 SuperWASP wide-field transit survey. Using the newly commissioned radial-velocity spectrograph SOPHIE at the Observatoire de Haute-Provence, we found that both objects exhibit reflex orbital radial-velocity variations with amplitudes characteristic of planetary-mass companions and in-phase with the photometric orbits. Line-bisector studies rule out faint blended binaries as the cause of either the radial-velocity variations or the transits. We perform preliminary spectral analyses of the host stars, which together with their radial-velocity variations and fits to the transit light curves yield estimates of the planetary masses and radii. WASP-1b and WASP-2b have orbital periods of 2.52 and 2.15 d, respectively. Given mass estimates for their F7V and K1V primaries, we derive planet masses 0.80–0.98 and 0.81–0.95 times that of Jupiter, respectively. WASP-1b appears to have an inflated radius of at least 1.33 R _Jup, whereas WASP-2b has a radius in the range 0.65–1.26 R _Jup.     

A SuperWASP search for additional transiting planets in 24 known systems

We present results from a search for additional transiting planets in 24 systems already known to contain a transiting planet. We model the transits due to the known planet in each system and subtract these models from light curves obtained with the SuperWASP (Wide Angle Search for Planets) survey instruments. These residual light curves are then searched for evidence of additional periodic transit events. Although we do not find any evidence for additional planets in any of the planetary systems studied, we are able to characterize our ability to find such planets by means of Monte Carlo simulations. Artificially generated transit signals corresponding to planets with a range of sizes and orbital periods were injected into the SuperWASP photometry and the resulting light curves searched for planets. As a result, the detection efficiency as a function of both the radius and orbital period of any second planet is calculated. We determine that there is a good (>50 per cent) chance of detecting additional, Saturn-sized planets in   P ∼  10 d orbits around planet-hosting stars that have several seasons of SuperWASP photometry. Additionally, we confirm previous evidence of the rotational stellar variability of WASP-10, and refine the period of rotation. We find that the period of the rotation is  11.91 ± 0.05  d, and the false alarm probability for this period is extremely low  (∼10⁻¹³)  .