Seasonal Changes in Chemical Composition of Eurasion Watermilfoil (Myriophyllum spicatum L.) and Water Temperature at Two Sites in northern California: implications for Herbivory

We compared seasonal changes in Eurasian watermilfoil (Myriophyllum spicatum L.) characteristics and water temperature for a shallow poind in Davis, CA, and the Truckee River, near Tahoe City, CA. Tissue C and N were 15% lower in plants from the Truckee River than in plants from the Davis pond. Seasonal fluctuations in tissue N were also different. Mean phenolic acid content of Truckee River palnts (162yM g-1) was less than those from the shallow pond (195 yM g-1). Phenolic acid content was positively related to tissue C for Truckee River and Davis pond plants and, tissue C:N ratio for Truckee River plants. Mean monthly water temperature (1990 to 1998) for the Truckee River site was less than 20 C. Water temperatures were warmer in August and September at this site. However, Eurasian watermilfoil collected during these months was characterized by lower levels of tissue N. During a 29-month period beginning January 1994, mean monthly water temperature for the Davis pond exceeded 20 C, only during July to September 1995. Tissue N was generally greater during summer for watermilfoil growing in the pond. These results imply that Eurasian watermilfoil biological control agents may have different developmental rates in these habitats, and thus different impacts on watermilfoil populations.

The Impact of Endothall on the Aquatic Plant Community of Kress Lake, Washington

CA dense mat-forming population of Eurasian watermilfoil ( Myriophyllum spicatum L . ) was interfering with fishing and recreation in a small western Washington lake. A low concentration (1.5 mg/L active ingredient) of the herbicide endothall formulated as Aquathol® K was used in 2000 to attempt to selectively control the Eurasian watermilfoil. Aquatic plant biomass and frequency data were collected before treatment, ten weeks after treatment and during the growing season for 3 additional years. Macrophyte data were analyzed to assess the herbicide’s impacts on Eurasian watermilfoil as well as the rest of the aquatic plant community. Results showed a significant decrease in Eurasian watermilfoil biomass and frequency 10 weeks after treatment. The Eurasian watermilfoil continued to be present, but at a significantly reduced level through the remainder of the study (3 years after treatment). Of the native plant species, large-leaf pondweed ( Potamogeton amplifolius Tucker . ) frequency and biomass was significantly reduced after treatment. Common elodea ( Elodea canadensis Rich.), muskgrass ( Chara sp. Vallaint.) and bladderwort ( Utricularia sp. L.) all increased significantly after treatment. (PDF has 6 pages.)

Effects of Triclopyr on Variable-Leaf Watermilfoil

The objective of the study described here was to determine the effect on variable-leaf watermilfoil of various combinations of triclopyr concentrations and exposure times using dosage rates that controlled Eurasian watermilfoil under laboratory and field conditions (Netherland and Getsinger 1992, Getsinger et al. 1997, Petty et al. 1998).

The weevil-watermilfoil interaction at different spatial scales: what we know and what we need to know

The North American weevil ( Euhrychiopsis lecontei (Dietz)) is being considered as a biological control agent for Eurasian watermilfoil ( Myriophyllum spicatum L.). This native insect damages watermilfoil plants and is frequently associated with declining watermilfoil populations

Relationship between water quality, watermilfoil frequency, and weevil distribution in the State of Washington

During the summer of 1997, we surveyed 50 waterbodies in Washington State to determine the distribution of the aquatic weevil Euhrychiopsis lecontei Dietz. We collected data on water quality and the frequency of occurrence of watermilfoil species within selected watermilfoil beds to compare the waterbodies and determine if they were related to the distribution E. lecontei . We found E. lecontei in 14 waterbodies, most of which were in eastern Washington. Only one lake with weevils was located in western Washington. Weevils were associated with both Eurasian ( Myriophyllum spicatum L.) and northern watermilfoil ( M. sibiricum K.). Waterbodies with E. lecontei had significantly higher ( P < 0.05) pH (8.7 ± 0.2) (mean ± 2SE), specific conductance (0.3 ± 0.08 mS cm -1 ) and total alkalinity (132.4 ± 30.8 mg CaCO 3 L -1 ). We also found that weevil presence was related to surface water temperature and waterbody location ( = 24.3, P ≤ 0.001) and of all the models tested, this model provided the best fit (Hosmer- Lemeshow goodness-of-fit = 4.0, P = 0.9). Our results suggest that in Washington State E. lecontei occurs primarily in eastern Washington in waterbodies with pH ≥ 8.2 and specific conductance ≥ 0.2 mS cm -1 . Furthermore, weevil distribution appears to be correlated with waterbody location (eastern versus western Washington) and surface water temperature.

Factors to consider when using native biological control organisms to manage exotic plants

Biological control of exotic plant populations with native organisms appears to be increasing, even though its success to date has been limited. Although many researchers and managers feel that native organisms are easier to use and present less risk to the environment this may not be true. Developing a successful management program with a native insect is dependent on a number of critical factors that need to be considered. Information is needed on the feeding preference of the agent, agent effectiveness, environmental regulation of the agent, unique requirements of the agent, population maintenance of the agent, and time to desired impact. By understanding these factors, researchers and managers can develop a detailed protocol for using the native biological control agent for a specific target plant. . We found E. lecontei in 14 waterbodies, most of which were in eastern Washington. Only one lake with weevils was located in western Washington. Weevils were associated with both Eurasian ( Myriophyllum spicatum L.) and northern watermilfoil ( M. sibiricum K.). Waterbodies with E. lecontei had significantly higher ( P < 0.05) pH (8.7 ± 0.2) (mean ± 2SE), specific conductance (0.3 ± 0.08 mS cm -1 ) and total alkalinity (132.4 ± 30.8 mg CaCO 3 L -1 ). We also found that weevil presence was related to surface water temperature and waterbody location ( = 24.3, P ≤ 0.001) and of all the models tested, this model provided the best fit (Hosmer- Lemeshow goodness-of-fit = 4.0, P = 0.9). Our results suggest that in Washington State E. lecontei occurs primarily in eastern Washington in waterbodies with pH ≥ 8.2 and specific conductance ≥ 0.2 mS cm -1 . Furthermore, weevil distribution appears to be correlated with waterbody location (eastern versus western Washington) and surface water temperature.

Seismic structure above and below the core-mantle boundary

Seismic structure above and below the core-mantle boundary (CMB) has been studied through use of travel time and waveform analyses of several different seismic wave groups. Anomalous systematic trends in observables document mantle heterogeneity on both large and small scales. Analog and digital data has been utilized, and in many cases the analog data has been optically scanned and digitized prior to analysis.

Differential travel times of S - SKS are shown to be an excellent diagnostic of anomalous lower mantle shear velocity (V s) structure. Wavepath geometries beneath the central Pacific exhibit large S- SKS travel time residuals (up to 10 sec), and are consistent with a large scale 0(1000 km) slower than average V_s region (≥3%). S - SKS times for paths traversing this region exhibit smaller scale patterns and trends 0(100 km) indicating V_s perturbations on many scale lengths. These times are compared to predictions of three tomographically derived aspherical models: MDLSH of Tanimoto [1990], model SH12_WM13 of Suet al. [1992], and model SH.10c.17 of Masters et al. [1992]. Qualitative agreement between the tomographic model predictions and observations is encouraging, varying from fair to good. However, inconsistencies are present and suggest anomalies in the lower mantle of scale length smaller than the present 2000+ km scale resolution of tomographic models. 2-D wave propagation experiments show the importance of inhomogeneous raypaths when considering lateral heterogeneities in the lowermost mantle.

A dataset of waveforms and differential travel times of S, ScS, and the arrival from the D" layer, Scd, provides evidence for a laterally varying V_s velocity discontinuity at the base of the mantle. Two different localized D" regions beneath the central Pacific have been investigated. Predictions from a model having a V_s discontinuity 180 km above the CMB agree well with observations for an eastern mid-Pacific CMB region. This thickness differs from V_s discontinuity thicknesses found in other regions, such as a localized region beneath the western Pacific, which average near 280 km. The "sharpness" of the V_s jump at the top of D", i.e., the depth range over which the V_s increase occurs, is not resolved by our data, and our data can in fact may be modeled equally well by a lower mantle with the increase in V_s at the top of D" occurring over a 100 krn depth range. It is difficult at present to correlate D" thicknesses from this study to overall lower mantle heterogeneity, due to uncertainties in the 3-D models, as well as poor coverage in maps of D" discontinuity thicknesses.

P-wave velocity structure (V_p) at the base of the mantle is explored using the seismic phases SKS and SPdKS. SPdKS is formed when SKS waves at distances around 107° are incident upon the CMB with a slowness that allows for coupling with diffracted P-waves at the base of the mantle. The P-wave diffraction occurs at both the SKS entrance and exit locations of the outer core. SP_dKS arrives slightly later in time than SKS, having a wave path through the mantle and core very close to SKS. The difference time between SKS and SP_dKS strongly depends on V_p at the base of the mantle near SK Score entrance and exit points. Observations from deep focus Fiji-Tonga events recorded by North American stations, and South American events recorded by European and Eurasian stations exhibit anomalously large SP_dKS - SKS difference times. SKS and the later arriving SP_dKS phase are separated by several seconds more than predictions made by 1-D reference models, such as the global average PREM [Dziewonski and Anderson, 1981] model. Models having a pronounced low-velocity zone (5%) in V_p in the bottom 50-100 km of the mantle predict the size of the observed SP_dK S-SKS anomalies. Raypath perturbations from lower mantle V_s structure may also be contributing to the observed anomalies.

Outer core structure is investigated using the family of SmKS (m=2,3,4) seismic waves. SmKS are waves that travel as S-waves in the mantle, P-waves in the core, and reflect (m-1) times on the underside of the CMB, and are well-suited for constraining outermost core V_p structure. This is due to closeness of the mantle paths and also the shallow depth range these waves travel in the outermost core. S3KS - S2KS and S4KS - S3KS differential travel times were measured using the cross-correlation method and compared to those from reflectivity synthetics created from core models of past studies. High quality recordings from a deep focus Java Sea event which sample the outer core beneath the northern Pacific, the Arctic, and northwestern North America (spanning 1/8th of the core's surface area), have SmKS wavepaths that traverse regions where lower mantle heterogeneity is pre- dieted small, and are well-modeled by the PREM core model, with possibly a small V_p decrease (1.5%) in the outermost 50 km of the core. Such a reduction implies chemical stratification in this 50 km zone, though this model feature is not uniquely resolved. Data having wave paths through areas of known D" heterogeneity (±2% and greater), such as the source-side of SmKS lower mantle paths from Fiji-Tonga to Eurasia and Africa, exhibit systematic SmKS differential time anomalies of up to several seconds. 2-D wave propagation experiments demonstrate how large scale lower mantle velocity perturbations can explain long wavelength behavior of such anomalous SmKS times. When improperly accounted for, lower mantle heterogeneity maps directly into core structure. Raypaths departing from homogeneity play an important role in producing SmKS anomalies. The existence of outermost core heterogeneity is difficult to resolve at present due to uncertainties in global lower mantle structure. Resolving a one-dimensional chemically stratified outermost core also remains difficult due to the same uncertainties. Restricting study to higher multiples of SmKS (m=2,3,4) can help reduce the affect of mantle heterogeneity due to the closeness of the mantle legs of the wavepaths. SmKS waves are ideal in providing additional information on the details of lower mantle heterogeneity.

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.

Geology of the Tick Canyon Area

The eastern part of the Ventura Basin contains great thicknesses of non-marine Tertiary sediments. The lower formations of the Tertiary strata outcrop in the Tick Canyon Area and are described in this report. Emphasis is placed on the description of the Vasquez formation which is the lowest Tertiary unit in the Tick Canyon Area and which contains the only Tertiary lavas found in the East Ventura Basin.

Viral metagenomic analysis of feces of wild small carnivores

Background: Recent studies have clearly demonstrated the enormous virus diversity that exists among wild animals. This exemplifies the required expansion of our knowledge of the virus diversity present in wildlife, as well as the potential transmission of these viruses to domestic animals or humans. Methods: In the present study we evaluated the viral diversity of fecal samples (n = 42) collected from 10 different species of wild small carnivores inhabiting the northern part of Spain using random PCR in combination with next-generation sequencing. Samples were collected from American mink (Neovison vison), European mink (Mustela lutreola), European polecat (Mustela putorius), European pine marten (Martes martes), stone marten (Martes foina), Eurasian otter (Lutra lutra) and Eurasian badger (Meles meles) of the family of Mustelidae; common genet (Genetta genetta) of the family of Viverridae; red fox (Vulpes vulpes) of the family of Canidae and European wild cat (Felis silvestris) of the family of Felidae. Results: A number of sequences of possible novel viruses or virus variants were detected, including a theilovirus, phleboviruses, an amdovirus, a kobuvirus and picobirnaviruses. Conclusions: Using random PCR in combination with next generation sequencing, sequences of various novel viruses or virus variants were detected in fecal samples collected from Spanish carnivores. Detected novel viruses highlight the viral diversity that is present in fecal material of wild carnivores.

Inferring Population Genetic Structure in Widely and Continuously Distributed Carnivores: The Stone Marten (Martes foina) as a Case Study

The stone marten is a widely distributed mustelid in the Palaearctic region that exhibits variable habitat preferences in different parts of its range. The species is a Holocene immigrant from southwest Asia which, according to fossil remains, followed the expansion of the Neolithic farming cultures into Europe and possibly colonized the Iberian Peninsula during the Early Neolithic (ca. 7,000 years BP). However, the population genetic structure and historical biogeography of this generalist carnivore remains essentially unknown. In this study we have combined mitochondrial DNA (mtDNA) sequencing (621 bp) and microsatellite genotyping (23 polymorphic markers) to infer the population genetic structure of the stone marten within the Iberian Peninsula. The mtDNA data revealed low haplotype and nucleotide diversities and a lack of phylogeographic structure, most likely due to a recent colonization of the Iberian Peninsula by a few mtDNA lineages during the Early Neolithic. The microsatellite data set was analysed with a) spatial and non-spatial Bayesian individual-based clustering (IBC) approaches (STRUCTURE, TESS, BAPS and GENELAND), and b) multivariate methods [discriminant analysis of principal components (DAPC) and spatial principal component analysis (sPCA)]. Additionally, because isolation by distance (IBD) is a common spatial genetic pattern in mobile and continuously distributed species and it may represent a challenge to the performance of the above methods, the microsatellite data set was tested for its presence. Overall, the genetic structure of the stone marten in the Iberian Peninsula was characterized by a NE-SW spatial pattern of IBD, and this may explain the observed disagreement between clustering solutions obtained by the different IBC methods. However, there was significant indication for contemporary genetic structuring, albeit weak, into at least three different subpopulations. The detected subdivision could be attributed to the influence of the rivers Ebro, Tagus and Guadiana, suggesting that main watercourses in the Iberian Peninsula may act as semi-permeable barriers to gene flow in stone martens. To our knowledge, this is the first phylogeographic and population genetic study of the species at a broad regional scale. We also wanted to make the case for the importance and benefits of using and comparing multiple different clustering and multivariate methods in spatial genetic analyses of mobile and continuously distributed species.

Evolutionary trends of the mitochondrial lineage differentiation in species of genera Martes and Mustela

We compared partial sequences (402 bp) of the mitochondrial cytochrome b gene in 68 individuals of martens (Martes), weasels (Mustela) and their relatives from the Northern Hemisphere to identify the modes of geographic differentiation in each species. We then compared complete sequences (1140 bp) of the gene in 17 species of the family Mustelidae to know the spatial and temporal modes of speciation, constructing linearized trees with transversional substitutions for deeper lineage divergences and with transversions and transitions for younger lineages. Our data suggested that these lineages of Martes and Mustela differentiated in a stepwise fashion with five radiation stages from the generic divergences (stage I) to the intraspecific divergences (stage V), during the last 10 or 20 million years as the fossil evidence suggests. In the lineage of Martes, the first offshoots are of Martes flavigula, M. pennanti, and Gulo gulo (stage II), the second is M. foina (stage III), and the third are M. americana, M. martes, M. melampus, and M. zibellina (stage IV). The divergence of the lineages of Mustela is likely to have taken place concurrently with the radiations of the Martes. These divergence processes are attributable in part to the geographic allocation along the two continents, North America and Eurasia, as well as among peripheral insular domains, such as Taiwan and the Japanese Islands. In addition, the Eurasian continent itself was shown to have been involved in the species diversification in the martens and weasels.

Molecular phylogenetics and biogeography of Lepus in Eastern Asia based on mitochondrial DNA sequences

In spite of several classification attempts among taxa of the genus Lepus, phylogenetic relationships still remain poorly understood. Here, we present molecular genetic evidence that may resolve some of the current incongruities in the phylogeny of the leporids. The complete mitochondrial cytb, 12S genes, and parts of ND4 and control region fragments were sequenced to examine phylogenetic relationships among Chinese hare taxa and other leporids throughout the World using maximum parsimony, maximum likelihood, and Bayesian phylogenetic reconstruction approaches. Using reconstructed phylogenies, we observed that the Chinese hare is not a single monophyletic group as originally thought. Instead, the data infers that the genus Lepus is monophyletic with three unique species groups: North American, Eurasian, and African. Ancestral area analysis indicated that ancestral Lepus arose in North America and then dispersed into Eurasia via the Bering Land Bridge eventually extending to Africa. Brooks Parsimony analysis showed that dispersal events followed by subsequent speciation have occurred in other geographic areas as well and resulted in the rapid radiation and speciation of Lepus. A Bayesian relaxed molecular clock approach based on the continuous autocorrelation of evolutionary rates along branches estimated the divergence time between the three major groups within Lepus. The genus appears to have arisen approximately 10.76 MYA (+/- 0.86 MYA), with most speciation events occurring during the Pliocene epoch (5.65 +/- 1.15 MYA similar to 1.12 +/- 10.47 MYA). (c) 2005 Elsevier Inc. All rights reserved.

Phylogeny of mitochondrial DNA macrohaplogroup N in India, based on complete sequencing: Implications for the peopling of South Asia

To resolve the phylogeny of the autochthonous mitochondrial DNA (mtDNA) haplogroups of India and determine the relationship between the Indian and western Eurasian mtDNA pools more precisely, a diverse subset of 75 macrohaplogroup N lineages was chosen fo