Spatial distribution of key zooplankton species using continuous plankton recorder (CPR) data from the North Atlantic (2000-2009)

The continuous plankton recorder (CPR) survey is an upper layer plankton monitoring program that has regularly collected samples, at monthly intervals, in the North Atlantic and adjacent seas since 1946. Water from approximately 6 m depth enters the CPR through a small aperture at the front of the sampler and travels down a tunnel where it passes through a silk filtering mesh of 270 µm before exiting at the back of the CPR. The plankton filtered on the silk is analyzed in sections corresponding to 10 nautical miles (approx. 3 m**3 of seawater filtered) and the plankton microscopically identified (Richardson et al., 2006 and reference therein). In the present study we used the CPR data to investigate the current basin scale distribution of C. finmarchicus (C5-C6), C. helgolandicus (C5-C6), C. hyperboreus (C5-C6), Pseudocalanus spp. (C6), Oithona spp. (C1-C6), total Euphausiida, total Thecosomata and the presence/absence of Cnidaria and the Phytoplankton Colour Index (PCI). The PCI, which is a visual assessment of the greenness of the silk, is used as an indicator of the distribution of total phytoplankton biomass across the Atlantic basin (Batten et al., 2003). Monthly data collected between 2000 and 2009 were gridded using the inverse-distance interpolation method, in which the interpolated values were the nodes of a 2 degree by 2 degree grid. The resulting twelve monthly matrices were then averaged within the year and in the case of the zooplankton the data were log-transformed (i.e. log10 (x+1).

Spectral analysis of seabed elevation from the Danish Wadden Sea

Bedforms both reflect and influence shallow water hydrodynamics and sediment dynamics. A correct characterization of their spatial distribution and dimensions is required for the understanding, assessment and prediction of numerous coastal processes. A method to parameterize geometrical characteristics using two-dimensional (2D) spectral analysis is presented and tested on seabed elevation data from the Knudedyb tidal inlet in the Danish Wadden Sea, where large compound bedforms are found. The bathymetric data were divided into 20x20 m areas on which a 2D spectral analysis was applied. The most energetic peak of the 2D spectrum was found and its energy, frequency and direction were calculated. A power-law was fitted to the average of slices taken through the 2D spectrum; its slope and y-intercept were calculated. Using these results the test area was morphologically classified into 4 distinct morphological regions. The most energetic peak and the slope and intercept of the power-law showed high values above the crest of the primary bedforms and scour holes, low values in areas without bedforms, and intermediate values in areas with secondary bedforms. The secondary bedform dimensions and orientations were calculated. An area of 700x700 m was used to determine the characteristics of the primary bedforms. However, they were less distinctively characterized compared to the secondary bedforms due to relatively large variations in their orientations and wavelengths. The method is thus appropriate for morphological classification of the seabed and for bedform characterization, being most efficient in areas characterized by bedforms with regular dimensions and directions.

Analysis of the Hydrogeological Conditions and the Effects of Development on the Hydraulic Head at Redmond Ridge East

Redmond Ridge East (RRE) is a large-scale master plan community in East King County, WA. In this report, I evaluate the spatial variability of the Quaternary Advance Outwash (Qva) at RRE and the time-series data for 16 water wells with the intent to better understand groundwater below the RRE area. I investigate changes between pre- and post-development conditions through the determination of temporal changes in annual water level, annual water level fluctuations, hydraulic head response to precipitation, and ambient drainage of the aquifer. I also perform a basic analysis of the annual aquifer recharge and a determination for the storage through the implementation of the water table fluctuation (WTF) method. Associated Earth Sciences (AESI) was tasked with monitoring the geological and environmental impacts during the development of RRE and collected the data I use in this report. AESI involvement in monitoring began in 1998 and extends to the present. Sixteen wells were identified in the RRE area with adequate temporal data to conduct the analysis. A comparison of the well logs and aquifer testing data allowed local variations in the Qva to be mapped. The WTF was used to determine a range of reasonable specific yield values for locations where the Qva was unconfined. Yearly average of the seasonal water level high and lows, and the fluctuations were quantified. Temporal relationships were established through linear regression. The average water level was found to be increasing in some locations, and the corresponding fluctuations were found to decrease. However, no clear change between pre- and post-development was observed. The response of hydraulic head to precipitation was investigated through an analysis of hydrographs for ten wells. Periods of consistent response and the corresponding precipitation during each period were delineated. A linear relationship between precipitation and water level change was determined. The threshold precipitation under which there is a positive response in the hydraulic head was established. No observable changes were apparent between pre- and post-development conditions. The ambient drainage for the Qva was calculated using recessional periods on the hydrograph. The transmissivity of Qva varies with thickness of the overlying lodgment till and thickness of the Qva, itself. Water level fluctuations observed in the Qva are consistent with regional observations. Localized areas in the Qva display the large 10 foot fluctuations and these anomalies are likely due to a combination of the local variability in the storativity as well as the concentration and channeling of water due to geographical variations in the Qva and the overlying topography. All trends seen in the RRE area remained relatively constant through time. There was no evidence showing an effect of development on the hydraulic head at RRE. This implies that the style and distribution of infiltration has not changed as a result of development, and that any measures in place are properly mitigating the effects of development on the RRE region.

Exploring the spatial variation of tides for seafloor mapping data corrections

Senior thesis written for Oceanography 445

Essays on the Analysis of Dynamic Games

Thesis (Ph.D.)--University of Washington, 2016-06

A vulnerability analysis of the temperate forests of south central Chile A vulnerability analysis of the temperate forests of south central Chile. Biological Conservation, 122 2005: 9-21.

Areas of the landscape that are priorities for conservation should be those that are both vulnerable to threatening processes and that if lost or degraded, will result in conservation targets being compromised. While much attention is directed towards understanding the patterns of biodiversity, much less is given to determining the areas of the landscape most vulnerable to threats. We assessed the relative vulnerability of remaining areas of native forest to conversion to plantations in the ecologically significant temperate rainforest region of south central Chile. The area of the study region is 4.2 million ha and the extent of plantations is approximately 200000 ha. First, the spatial distribution of native forest conversion to plantations was determined. The variables related to the spatial distribution of this threatening process were identified through the development of a classification tree and the generation of a multivariate. spatially explicit, statistical model. The model of native forest conversion explained 43% of the deviance and the discrimination ability of the model was high. Predictions were made of where native forest conversion is likely to occur in the future. Due to patterns of climate, topography, soils and proximity to infrastructure and towns, remaining forest areas differ in their relative risk of being converted to plantations. Another factor that may increase the vulnerability of remaining native forest in a subset of the study region is the proposed construction of a highway. We found that 90% of the area of existing plantations within this region is within 2.5 km of roads. When the predictions of native forest conversion were recalculated accounting for the construction of this highway, it was found that: approximately 27000 ha of native forest had an increased probability of conversion. The areas of native forest identified to be vulnerable to conversion are outside of the existing reserve network. (C) 2004 Elsevier Ltd. All tights reserved.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Identification, detection, and spatial resolution of Clostridium populations responsible for cellulose degradation in a methanogenic landfill leachate bioreactor

An anaerobic landfill leachate bioreactor was operated with crystalline cellulose and sterile landfill leacbate until a steady state was reached. Cellulose hydrolysis, acidogenesis, and methanogenesis were measured. Microorganisms attached to the cellulose surfaces were hypothesized to be the cellulose hydrolyzers. 16S rRNA gene clone libraries were prepared from this attached fraction and also from the mixed fraction (biomass associated with cellulose particles and in the planktonic phase). Both clone libraries were dominated by Firmicutes phylum sequences (100% of the attached library and 90% of the mixed library), and the majority fell into one of five lineages of the clostridia. Clone group 1 (most closely related to Clostridium stercorarium), clone group 2 (most closely related to Clostridium thermocellum), and clone group 5 (most closely related to Bacteroides cellulosolvens) comprised sequences in Clostridium group III. Clone group 3 sequences were in Clostridium group XIVa (most closely related to Clostridium sp. strain XB90). Clone group 4 sequences were affiliated with a deeply branching clostridial lineage peripherally associated with Clostridium group VI. This monophyletic group comprises a new Clostridium cluster, designated cluster VIa. Specific fluorescence in situ hybridization (FISH) probes for the five groups were designed and synthesized, and it was demonstrated in FISH experiments that bacteria targeted by the probes for clone groups 1, 2, 4, and 5 were very abundant on the surfaces of the cellulose particles and likely the key cellulolytic microorganisms in the landfill bioreactor. The FISH probe for clone group 3 targeted cells in the planktonic phase, and these organisms were hypothesized to be glucose fermenters.

A broad-scale analysis of links between coastal fisheries production and mangrove extent: A case-study for northeastern Australia

The paradigm that mangroves are critical for sustaining production in coastal fisheries is widely accepted, but empirical evidence has been tenuous. This study showed that links between mangrove extent and coastal fisheries production could be detected for some species at a broad regional scale (1000s of kilometres) on the east coast of Queensland, Australia. The relationships between catch-per-unit-effort for different commercially caught species in four fisheries (trawl, line, net and pot fisheries) and mangrove characteristics, estimated from Landsat images were examined using multiple regression analyses. The species were categorised into three groups based on information on their life history characteristics, namely mangrove-related species (banana prawns Penaeus merguiensis, mud crabs Scylla serrata and barramundi Lates calcarifer), estuarine species (tiger prawns Penaeus esculentus and Penaeus semisulcatus, blue swimmer crabs Portunus pelagicus and blue threadfin Eleutheronema tetradactylum) and offshore species (coral trout Plectropomus spp.). For the mangrove-related species, mangrove characteristics such as area and perimeter accounted for most of the variation in the model; for the non-mangrove estuarine species, latitude was the dominant parameter but some mangrove characteristics (e.g. mangrove perimeter) also made significant contributions to the models. In contrast, for the offshore species, latitude was the dominant variable, with no contribution from mangrove characteristics. This study also identified that finer scale spatial data for the fisheries, to enable catch information to be attributed to a particular catchment, would help to improve our understanding of relationships between mangroves and fisheries production. (C) 2005 Elsevier B.V. All rights reserved.

Demonstration of the spatial separation of the entangled quantum sidebands of an optical field

Quantum optics experiments on bright beams are based on the spectral analysis of field fluctuations and typically probe correlations between radio-frequency sideband modes. However, the extra degree of freedom represented by this dual-mode picture is generally ignored. We demonstrate the experimental operation of a device which can be used to separate the quantum sidebands of an optical field. We use this device to explicitly demonstrate the quantum entanglement between the sidebands of a squeezed beam.

Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model

This paper describes a biventricular model, which couples the electrical and mechanical properties of the heart, and computer simulations of ventricular wall motion and deformation by means of a biventricular model. In the constructed electromechanical model, the mechanical analysis was based on composite material theory and the finite-element method; the propagation of electrical excitation was simulated using an electrical heart model, and the resulting active forces were used to calculate ventricular wall motion. Regional deformation and Lagrangian strain tensors were calculated during the systole phase. Displacements, minimum principal strains and torsion angle were used to describe the motion of the two ventricles. The simulations showed that during the period of systole, (1) the right ventricular free wall moves towards the septum, and at the same time, the base and middle of the free wall move towards the apex, which reduces the volume of the right ventricle; the minimum principle strain (E3) is largest at the apex, then at the middle of the free wall and its direction is in the approximate direction of the epicardial muscle fibres; (2) the base and middle of the left ventricular free wall move towards the apex and the apex remains almost static; the torsion angle is largest at the apex; the minimum principle strain E3 is largest at the apex and its direction on the surface of the middle wall of the left ventricle is roughly in the fibre orientation. These results are in good accordance with results obtained from MR tagging images reported in the literature. This study suggests that such an electromechanical biventricular model has the potential to be used to assess the mechanical function of the two ventricles, and also could improve the accuracy ECG simulation when it is used in heart torso model-based body surface potential simulation studies.

Analysis of transient eddy currents in MRI using a cylindrical FDTD method

Most magnetic resonance imaging (MRI) spatial encoding techniques employ low-frequency pulsed magnetic field gradients that undesirably induce multiexponentially decaying eddy currents in nearby conducting structures of the MRI system. The eddy currents degrade the switching performance of the gradient system, distort the MRI image, and introduce thermal loads in the cryostat vessel and superconducting MRI components. Heating of superconducting magnets due to induced eddy currents is particularly problematic as it offsets the superconducting operating point, which can cause a system quench. A numerical characterization of transient eddy current effects is vital for their compensation/control and further advancement of the MRI technology as a whole. However, transient eddy current calculations are particularly computationally intensive. In large-scale problems, such as gradient switching in MRI, conventional finite-element method (FEM)-based routines impose very large computational loads during generation/solving of the system equations. Therefore, other computational alternatives need to be explored. This paper outlines a three-dimensional finite-difference time-domain (FDTD) method in cylindrical coordinates for the modeling of low-frequency transient eddy currents in MRI, as an extension to the recently proposed time-harmonic scheme. The weakly coupled Maxwell's equations are adapted to the low-frequency regime by downscaling the speed of light constant, which permits the use of larger FDTD time steps while maintaining the validity of the Courant-Friedrich-Levy stability condition. The principal hypothesis of this work is that the modified FDTD routine can be employed to analyze pulsed-gradient-induced, transient eddy currents in superconducting MRI system models. The hypothesis is supported through a verification of the numerical scheme on a canonical problem and by analyzing undesired temporal eddy current effects such as the B-0-shift caused by actively shielded symmetric/asymmetric transverse x-gradient head and unshielded z-gradient whole-body coils operating in proximity to a superconducting MRI magnet.

Coupled thermal, structural and vibrational analysis of a hypersonic engine for flight test

This paper describes a relatively simple and quick method for implementing aerodynamic heating models into a finite element code for non-linear transient thermal-structural and thermal-structural-vibrational analyses of a Mach 10 generic HyShot scramjet engine. The thermal-structural-vibrational response of the engine was studied for the descent trajectory from 60 to 26 km. Aerodynamic heating fluxes, as a function of spatial position and time for varying trajectory points, were implemented in the transient heat analysis. Additionally, the combined effect of varying dynamic pressure and thermal loads with altitude was considered. This aero-thermal-structural analysis capability was used to assess the temperature distribution, engine geometry distortion and yielding of the structural material due to aerodynamic heating during the descent trajectory, and for optimising the wall thickness, nose radius of leading edge, etc. of the engine intake. A structural vibration analysis was also performed following the aero-thermal-structural analysis to determine the changes in natural frequencies of the structural vibration modes that occur at the various temperatures associated with the descent trajectory. This analysis provides a unique and relatively simple design strategy for predicting and mitigating the thermal-structural-vibrational response of hypersonic engines. (C) 2006 Elsevier SAS. All rights reserved.

Definition and spatial annotation of the dynamic secretome during early kidney development

The term secretome has been defined as a set of secreted proteins (Grimmond et al. [2003] Genome Res 13:1350-1359). The term secreted protein encompasses all proteins exported from the cell including growth factors, extracellular proteinases, morphogens, and extracellular matrix molecules. Defining the genes encoding secreted proteins that change in expression during organogenesis, the dynamic secretome, is likely to point to key drivers of morphogenesis. Such secreted proteins are involved in the reciprocal interactions between the ureteric bud (UB) and the metanephric mesenchyme (AM) that occur during organogenesis of the metanephros. Some key metanephric secreted proteins have been identified, but many remain to be determined. In this study, microarray expression profiling of E10.5, E11.5, and E13.5 kidney and consensus bioinformatic analysis were used to define a dynamic secretome of early metanephric development. In situ hybridisation was used to confirm microarray results and clarify spatial expression patterns for these genes. Forty-one secreted factors were dynamically expressed between the E10.5 and E13.5 timeframe profiled, and 25 of these factors had not previously been implicated in kidney development. A text-based anatomical ontology was used to spatially annotate the expression pattern of these genes in cultured metanephric explants.