QUAGMIRE v1.3: a quasi-geostrophic model for investigating rotating fluids experiments

QUAGMIRE is a quasi-geostrophic numerical model for performing fast, high-resolution simulations of multi-layer rotating annulus laboratory experiments on a desktop personal computer. The model uses a hybrid finite-difference/spectral approach to numerically integrate the coupled nonlinear partial differential equations of motion in cylindrical geometry in each layer. Version 1.3 implements the special case of two fluid layers of equal resting depths. The flow is forced either by a differentially rotating lid, or by relaxation to specified streamfunction or potential vorticity fields, or both. Dissipation is achieved through Ekman layer pumping and suction at the horizontal boundaries, including the internal interface. The effects of weak interfacial tension are included, as well as the linear topographic beta-effect and the quadratic centripetal beta-effect. Stochastic forcing may optionally be activated, to represent approximately the effects of random unresolved features. A leapfrog time stepping scheme is used, with a Robert filter. Flows simulated by the model agree well with those observed in the corresponding laboratory experiments

On the nature and variability of the east Greenland Spill Jet : a case study in Summer 2003

Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 41 (2011): 2307–2327, doi:10.1175/JPO-D-10-05004.1.Results from a high-resolution (~2 km) numerical simulation of the Irminger Basin during summer 2003 are presented. The focus is on the East Greenland Spill Jet, a recently discovered component of the circulation in the basin. The simulation compares well with observations of surface fields, the Denmark Strait overflow (DSO), and the hydrographic structure of typical sections in the basin. The model reveals new aspects of the circulation on scales of O(0.1–10) days and O(1–100) km. The model Spill Jet results from the cascade of dense waters over the East Greenland shelf. Spilling can occur in various locations southwest of the strait, and it is present throughout the simulation but exhibits large variations on periods of O(0.1–10) days. The Spill Jet sometimes cannot be distinguished in the velocity field from surface eddies or from the DSO. The vorticity structure of the jet confirms its unstable nature with peak relative and tilting vorticity terms reaching twice the planetary vorticity term. The average model Spill Jet transport is 4.9 ±1.7 Sv (1 Sv ≡ 106 m3 s−1) equatorward, about 2½ times larger than has been previously reported from a single ship transect in August 2001. Kinematic analysis of the model results suggests two different types of spilling events. In the first case (type I), a local perturbation results in dense waters descending over the shelf break into the Irminger Basin. In the second case (type II), surface cyclones associated with DSO deep domes initiate the spilling process. During summer 2003, more than half of the largest Spill Jet transport values are of type II.The research is supported by the National Science Foundation Grants OCE-0726393 and OCI-0904640 (MGM and TWNH) and OCE-0726640 (RSP).2012-06-0

Molecular Approaches to Identify Cryptic Species and Polymorphic Species within a Complex Community of Fig Wasps

Cryptic and polymorphic species can complicate traditional taxonomic research and both of these concerns are common in fig wasp communities. Species identification is very difficult, despite great effort and the ecological importance of fig wasps. Herein, we try to identify all chalcidoid wasp species hosted by one species of fig, using both morphological and molecular methods. We compare the efficiency of four different DNA regions and find that ITS2 is highly effective for species identification, while mitochondrial COI and Cytb regions appear less reliable, possibly due to the interference signals from either nuclear copies of mtDNA, i.e. NUMTs, or the effects of Wolbachia infections. The analyses suggest that combining multiple markers is the best choice for inferring species identifications as any one marker may be unsuitable in a given case

Comparative Phylogeography of a Coevolved Community: Concerted Population Expansions in Joshua Trees and Four Yucca Moths

Comparative phylogeographic studies have had mixed success in identifying common phylogeographic patterns among co-distributed organisms. Whereas some have found broadly similar patterns across a diverse array of taxa, others have found that the histories of different species are more idiosyncratic than congruent. The variation in the results of comparative phylogeographic studies could indicate that the extent to which sympatrically-distributed organisms share common biogeographic histories varies depending on the strength and specificity of ecological interactions between them. To test this hypothesis, we examined demographic and phylogeographic patterns in a highly specialized, coevolved community – Joshua trees (Yucca brevifolia) and their associated yucca moths. This tightly-integrated, mutually interdependent community is known to have experienced significant range changes at the end of the last glacial period, so there is a strong a priori expectation that these organisms will show common signatures of demographic and distributional changes over time. Using a database of >5000 GPS records for Joshua trees, and multi-locus DNA sequence data from the Joshua tree and four species of yucca moth, we combined paleaodistribution modeling with coalescent-based analyses of demographic and phylgeographic history. We extensively evaluated the power of our methods to infer past population size and distributional changes by evaluating the effect of different inference procedures on our results, comparing our palaeodistribution models to Pleistocene-aged packrat midden records, and simulating DNA sequence data under a variety of alternative demographic histories. Together the results indicate that these organisms have shared a common history of population expansion, and that these expansions were broadly coincident in time. However, contrary to our expectations, none of our analyses indicated significant range or population size reductions at the end of the last glacial period, and the inferred demographic changes substantially predate Holocene climate changes

Quantifying mesoscale variability in ocean transient tracer fields

A powerful way to test the realism of ocean general circulation models is to systematically compare observations of passive tracer concentration with model predictions. The general circulation models used in this way cannot resolve a full range of vigorous mesoscale activity (on length scales between 10–100 km). In the real ocean, however, this activity causes important variability in tracer fields. Thus, in order to rationally compare tracer observations with model predictions these unresolved fluctuations (the model variability error) must be estimated. We have analyzed this variability using an eddy‐resolving reduced‐gravity model in a simple midlatitude double‐gyre configuration. We find that the wave number spectrum of tracer variance is only weakly sensitive to the distribution of (large scale slowly varying) tracer sources and sinks. This suggests that a universal passive tracer spectrum may exist in the ocean. We estimate the spectral shape using high‐resolution measurements of potential temperature on an isopycnal in the upper northeast Atlantic Ocean, finding a slope near k −1.7 between 10 and 500 km. The typical magnitude of the variance is estimated by comparing tracer simulations using different resolutions. For CFC‐ and tritium‐type transient tracers the peak magnitude of the model variability saturation error may reach 0.20 for scales shorter than 100 km. This is of the same order as the time mean saturation itself and well over an order of magnitude greater than the instrumental uncertainty