Distribution of calcareous nannofossils in the Miocene interval of DSDP Hole 94-608 in the north Atlantic (Fig. 1)

DSDP North Atlantic Site 608 yielded an excellent Miocene pelagic section which affords a further opportunity for elucidating the chronology of the calcareous nannofossil succession in the framework of magnetostratigraphic control. Most of the conventional (zonal) markers have been documented for this site and some of the earlier results are confirmed and refined. In addition several unconventional and less known markers have been added. The first two are the highest (last) occurrence of Sphenolithus delphix and Sphenolithus capricornutus at 23.6 Ma, which is immediately above the Oligocene-Miocene boundary as identified by the last occurrence of Reticulofenestra bisecta at 23.7 Ma. The next unconventional datum is the highest (last) occurrence of Ilselithina fusa at 22.8 Ma, which is also the highest (last) occurrence of Helicosphaera recta. Calcidiscus tropicus' lowest (first) occurrence is at 19.5 Ma, which is also the lowest occurrence of Sphenolithus belemnos, and Calcidiscus leptoporus' lowest (first) occurrence coincides with that of Sphenolithus heteromorphus at 18.5 Ma. Sphenolithus dissimilis' highest (last) occurrence is at 18.2 Ma and the Calcidiscus premacintyrei lowest (first) and highest (last) occurrences are, respectively, at 17.7 and 11.7 Ma. Discoaster braarudii occurs from 11.6 to 11.3 Ma and its highest (last) occurrence corresponds to that of Cyclicargolithus floridanus. Minylitha convallis occurs from 9.0 to 6.9 Ma. Within the range of Minylitha, at 8.0 Ma, a major shift occurs in reticulofenestrid placoliths from dominantly large (Reticulofenestra pseudoumbilicus) and medium size (Reticulofenestra minutula) species below to significant numbers of very small species (Dictyococcites productus and Gephyrocapsa) above. This is interpreted to be a major, though perhaps seasonal, change of productivity of the North Atlantic at Site 608. A new genus and species Cryptococcolithus takayamae, is described and a variety, Reticulofenestra pseudoumbilicus var. amplus is identified.

Determining Spatial Distribution and Physical Properties of the Vashon Advance Outwash near Mountlake Terrace, Washington

This study is aimed at determining the spatial distribution, physical properties, and groundwater conditions of the Vashon advance outwash (Qva) in the Mountlake Terrace, WA area. The Qva is correlative with the Esperance Sand, as defined at its type section; however, local variations in the Qva are not well-characterized (Mullineaux, 1965). While the Qva is a dense glacial unit with low compressibility and high frictional shear strength (Gurtowski and Boirum, 1989), the strength of this unit can be reduced when it becomes saturated (Tubbs, 1974). This can lead to caving or flowing in excavations, and on a larger scale, can lead to slope failures and mass-wasting when intersected by steep slopes. By studying the Qva, we can better predict how it will behave under certain conditions, which will be beneficial to geologists, hydrogeologists, engineers, and environmental scientists during site assessments and early phases of project planning. In this study, I use data from 27 geotechnical borings from previous field investigations and C-Tech Corporation’s EnterVol software to create three-dimensional models of the subsurface geology in the study area. These models made it possible to visualize the spatial distribution of the Qva in relation to other geologic units. I also conducted a comparative study between data from the borings and generalized published data on the spatial distribution, relative density, soil classification, grain-size distribution, moisture content, groundwater conditions, and aquifer properties of the Qva. I found that the elevation of the top of the Qva ranges from 247 to 477 ft. I found that the Qva is thickest where the modern topography is high, and is thinnest where the topography is low. The thickness of the Qva ranges from absent to 242 ft. Along the northern, east-west trending transect, the Qva thins to the east as it rises above a ridge composed of Pre- Vashon glacial deposits. Along the southern, east-west trending transect, the Qva pinches out against a ridge composed of pre-Vashon interglacial deposits. Two plausible explanations for this ridge are paleotopography and active faulting associated with the Southern Whidbey Fault Zone. Further investigations should be done using geophysical methods and the modeling methods described in this study to determine the nature of this ridge. The relative density of the Qva in the study area ranges from loose to very dense, with the loose end of the spectrum probably relating to heave in saturated sands. I found subtle correlations between density and depth. Volumetric analysis of the soil groups listed in the boring logs indicate that the Qva in the study area is composed of approximately 9.5% gravel, 89.3% sand, and 1.2% silt and clay. The natural moisture content ranges from 3.0 to 35.4% in select samples from the Qva. The moisture content appears to increase with depth and fines content. The water table in the study area ranges in elevation from 231.9 to 458 ft, based on observations and measurements recorded in the boring logs. The results from rising-head and falling-head slug tests done at a single well in the study area indicate that the geometric mean of hydraulic conductivity is 15.93 ft/d (5.62 x 10-03 cm/s), the storativity is 3.28x10-03, and the estimated transmissivity is 738.58 ft2/d in the vicinity of this observation well. At this location, there was 1.73 ft of seasonal variation in groundwater elevation between August 2014 and March 2015.

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.

Beluga whale distribution, migration, and behavior in a changing Pacific Arctic

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

Implications of seasonal climate forecasts on world wheat trade: a stochastic, dynamic analysis

Improvements in seasonal climate forecasts have potential economic implications for international agriculture. A stochastic, dynamic simulation model of the international wheat economy is developed to estimate the potential effects of seasonal climate forecasts for various countries' wheat production, exports and world trade. Previous studies have generally ignored the stochastic and dynamic aspects of the effects associated with the use of climate forecasts. This study shows the importance of these aspects. In particular with free trade, the use of seasonal forecasts results in increased producer surplus across all exporting countries. In fact, producers appear to capture a large share of the economic surplus created by using the forecasts. Further, the stochastic dimensions suggest that while the expected long-run benefits of seasonal forecasts are positive, considerable year-to-year variation in the distribution of benefits between producers and consumers should be expected. The possibility exists for an economic measure to increase or decrease over a 20-year horizon, depending on the particular sequence of years.

Three putative types of El Nino revealed by spatial variability in impact on Australian wheat yield

The El Nino-Southern Oscillation (ENSO) phenomenon significantly impacts rainfall and ensuing crop yields in many parts of the world. In Australia, El Nino events are often associated with severe drought conditions. However, El Nino events differ spatially and temporally in their manifestations and impacts, reducing the relevance of ENSO-based seasonal forecasts. In this analysis, three putative types of El Nino are identified among the 24 occurrences since the beginning of the twentieth century. The three types are based on coherent spatial patterns (footprints) found in the El Nino impact on Australian wheat yield. This bioindicator reveals aligned spatial patterns in rainfall anomalies, indicating linkage to atmospheric drivers. Analysis of the associated ocean-atmosphere dynamics identifies three types of El Nino differing in the timing of onset and location of major ocean temperature and atmospheric pressure anomalies. Potential causal mechanisms associated with these differences in anomaly patterns need to be investigated further using the increasing capabilities of general circulation models. Any improved predictability would be extremely valuable in forecasting effects of individual El Nino events on agricultural systems.

Potential to increase yield in lucerne (Medicago sativa subsp sativa) through introgression of Medicago sativa subsp falcata into Australian adapted material

The effect of interspecific heterosis in crosses between Medicago sativa subsp. sativa and M. sativa subsp. falcata was assessed. Three sativa and 3 falcata plants were crossed in a diallel design. Progeny dry matter yield and natural plant height were assessed in a replicated field experiment at Gatton, Queensland. Yield data were analysed using the method of residual maximum likelihood (REML) and Griffing's model 1. There were significant differences between the reciprocal, general combining ability (GCA), and specific combining ability (SCA) effects. As expected, S-1 populations were lower yielding than their respective intraspecific cross and falcata x falcata crosses were significantly lower yielding than sativa x sativa crosses. Some of the interspecific crosses indicated substantial SCA effects, yielding at least as well as the best sativa x sativa crosses. We have demonstrated the potential usefulness of unselected M. sativa subsp. falcata as a heterotic group in the improvement of yield in northern Australian adapted lucerne material, and discuss how it could be incorporated into future breeding to overcome the yield stagnation currently being experienced in Australian programs.

DNA markers linked to yield, yield components, and morphological traits in autotetraploid lucerne (Medicago sativa L.)

We have mapped and identifed DNA markers linked to morphology, yield, and yield components of lucerne, using a backcross population derived from winter-active parents. The high-yielding and recurrent parent, D, produced individual markers that accounted for up to 18% of total yield over 6 harvests, at Gatton, south-eastern Queensland. The same marker, AC/TT8, was consistently identified at each individual harvest, and in individual harvests accounted for up to 26% of the phenotypic variation for yield. This marker was located in linkage group 2 of the D map, and several other markers positively associated with yield were consistently identified in this linkage group. Similarly, markers negatively associated with yield were consistently identified in the W116 map, W116 being the low-yielding parent. Highly significant positive correlations were observed between total yield and yield for harvests 1-6, and between total yield and stem length, tiller number, leaf yield/plant, leaf yield/5 stems, stem yield/plant, and stem yield/5 stems. Highly significant QTL were located for all these characters as well as for leaf shape and pubescence.

Seasonal growth of foliose lichens in successive years in South Gwynedd, Wales

The aim of this study was to test the hypothesis that differences in the pattern of seasonal growth in foliose lichens from year to year were determined by yearly differences in the distribution of rainfall, shortwave radiation and temperature. Hence, the radial growth of Parmelia conspersa (Ehrh. Ex Ach.) Ach. , P. glabratula ssp. fuliginosa (Fr. ex Duby) Laund. and Physcia orbicularis (Neck) Poetsch. was studied on slate fragments over 34 successive months in an area of South Gwynedd, Wales. U.K. Similarities and differences were observed in the pattern of seasonal growth in the three species. Periods of maximum growth of a species occurred in different seasons in successive years. Correlation and multiple regression analysis suggested that total rainfall per month was the most important climatic variable positively correlated with monthly growth. Significant positive correlations were found in some growth periods with number of raindays per month, average wind speed and maximum and minimum temperature. Total number of sunshine hours per month and the frequency of ground frosts were negatively correlated with monthly growth in some growth periods. For each species, monthly radial growth was correlated with different climatic variables in each growth period. Hence, the results support the hypothesis in that periods of maximum growth can occur in any season in South Gwynedd and depend on (1) the distribution of periods of high total rainfall and (2) whether or not these periods coincide with periods of maximum sunlight.

Principal component and second generation wavelet analysis of Treasury yield curve evolution

Prices of U.S. Treasury securities vary over time and across maturities. When the market in Treasurys is sufficiently complete and frictionless, these prices may be modeled by a function time and maturity. A cross-section of this function for time held fixed is called the yield curve; the aggregate of these sections is the evolution of the yield curve. This dissertation studies aspects of this evolution. ^ There are two complementary approaches to the study of yield curve evolution here. The first is principal components analysis; the second is wavelet analysis. In both approaches both the time and maturity variables are discretized. In principal components analysis the vectors of yield curve shifts are viewed as observations of a multivariate normal distribution. The resulting covariance matrix is diagonalized; the resulting eigenvalues and eigenvectors (the principal components) are used to draw inferences about the yield curve evolution. ^ In wavelet analysis, the vectors of shifts are resolved into hierarchies of localized fundamental shifts (wavelets) that leave specified global properties invariant (average change and duration change). The hierarchies relate to the degree of localization with movements restricted to a single maturity at the base and general movements at the apex. Second generation wavelet techniques allow better adaptation of the model to economic observables. Statistically, the wavelet approach is inherently nonparametric while the wavelets themselves are better adapted to describing a complete market. ^ Principal components analysis provides information on the dimension of the yield curve process. While there is no clear demarkation between operative factors and noise, the top six principal components pick up 99% of total interest rate variation 95% of the time. An economically justified basis of this process is hard to find; for example a simple linear model will not suffice for the first principal component and the shape of this component is nonstationary. ^ Wavelet analysis works more directly with yield curve observations than principal components analysis. In fact the complete process from bond data to multiresolution is presented, including the dedicated Perl programs and the details of the portfolio metrics and specially adapted wavelet construction. The result is more robust statistics which provide balance to the more fragile principal components analysis. ^

Plasma and cold sprayed aluminum carbon nanotube composites: Quantification of nanotube distribution and multi-scale mechanical properties

Carbon nanotubes (CNT) could serve as potential reinforcement for metal matrix composites for improved mechanical properties. However dispersion of carbon nanotubes (CNT) in the matrix has been a longstanding problem, since they tend to form clusters to minimize their surface area. The aim of this study was to use plasma and cold spraying techniques to synthesize CNT reinforced aluminum composite with improved dispersion and to quantify the degree of CNT dispersion as it influences the mechanical properties. Novel method of spray drying was used to disperse CNTs in Al-12 wt.% Si prealloyed powder, which was used as feedstock for plasma and cold spraying. A new method for quantification of CNT distribution was developed. Two parameters for CNT dispersion quantification, namely Dispersion parameter (DP) and Clustering Parameter (CP) have been proposed based on the image analysis and distance between the centers of CNTs. Nanomechanical properties were correlated with the dispersion of CNTs in the microstructure. Coating microstructure evolution has been discussed in terms of splat formation, deformation and damage of CNTs and CNT/matrix interface. Effect of Si and CNT content on the reaction at CNT/matrix interface was thermodynamically and kinetically studied. A pseudo phase diagram was computed which predicts the interfacial carbide for reaction between CNT and Al-Si alloy at processing temperature. Kinetic aspects showed that Al4C3 forms with Al-12 wt.% Si alloy while SiC forms with Al-23wt.% Si alloy. Mechanical properties at nano, micro and macro-scale were evaluated using nanoindentation and nanoscratch, microindentation and bulk tensile testing respectively. Nano and micro-scale mechanical properties (elastic modulus, hardness and yield strength) displayed improvement whereas macro-scale mechanical properties were poor. The inversion of the mechanical properties at different scale length was attributed to the porosity, CNT clustering, CNT-splat adhesion and Al 4C3 formation at the CNT/matrix interface. The Dispersion parameter (DP) was more sensitive than Clustering parameter (CP) in measuring degree of CNT distribution in the matrix.

Physical factors influencing the distribution of a top predator in a subtropical oligotrophic estuary

We used longline fishing to determine the effects of distance from the ocean, season, and short-term variation in abiotic conditions on the abundance of juvenile bull sharks (Carcharhinus leucas) in an estuary of the Florida Everglades, U.S.A. Logistic regression revealed that young-of-the-year sharks were concentrated at a protected site 20 km upstream and were present in greater abundance when dissolved oxygen (DO) levels were high. For older juvenile sharks (age 1+), DO levels had the greatest influence on catch probabilities followed by distance from the ocean; they were most likely to be caught at sites with .3.5 mg L21 DO and on the main branch of the river 20 km upstream. Salinity had a relatively small effect on catch rates and there were no seasonal shifts in shark distribution. Our results highlight the importance of considering DO as a possible driver of top predator distributions in estuaries, even in the absence of hypoxia. In Everglades estuaries hydrological drivers that affect DO levels (e.g., groundwater discharge, modification of primary productivity through nutrient fluxes) will be important in determining shark distributions, and the effects of planned ecosystem restoration efforts on bull sharks will not simply be mediated by changing salinity regimes and the location of the oligohaline zone. More generally, variation in DO levels could structure the nature and spatiotemporal pattern of top predator effects in the coastal Everglades, and other tropical and subtropical estuaries, because of interspecific variation in reliance on DO within the top predator guild.

Distribution of Diatoms and Development of Diatom-Based Models for Inferring Salinity and Nutrient Concentrations in Florida Bay and Adjacent Coastal Wetlands of South Florida (USA)

The composition and distribution of diatom algae inhabiting estuaries and coasts of the subtropical Americas are poorly documented, especially relative to the central role diatoms play in coastal food webs and to their potential utility as sentinels of environmental change in these threatened ecosystems. Here, we document the distribution of diatoms among the diverse habitat types and long environmental gradients represented by the shallow topographic relief of the South Florida, USA, coastline. A total of 592 species were encountered from 38 freshwater, mangrove, and marine locations in the Everglades wetland and Florida Bay during two seasonal collections, with the highest diversity occurring at sites of high salinity and low water column organic carbon concentration (WTOC). Freshwater, mangrove, and estuarine assemblages were compositionally distinct, but seasonal differences were only detected in mangrove and estuarine sites where solute concentration differed greatly between wet and dry seasons. Epiphytic, planktonic, and sediment assemblages were compositionally similar, implying a high degree of mixing along the shallow, tidal, and storm-prone coast. The relationships between diatom taxa and salinity, water total phosphorus (WTP), water total nitrogen (WTN), and WTOC concentrations were determined and incorporated into weighted averaging partial least squares regression models. Salinity was the most influential variable, resulting in a highly predictive model (r apparent 2  = 0.97, r jackknife 2  = 0.95) that can be used in the future to infer changes in coastal freshwater delivery or sea-level rise in South Florida and compositionally similar environments. Models predicting WTN (r apparent 2  = 0.75, r jackknife 2  = 0.46), WTP (r apparent 2  = 0.75, r jackknife 2  = 0.49), and WTOC (r apparent 2  = 0.79, r jackknife 2  = 0.57) were also strong, suggesting that diatoms can provide reliable inferences of changes in solute delivery to the coastal ecosystem.

Plasma And Cold Sprayed Aluminum Carbon Nanotube Composites: Quantification Of Nanotube Distribution And Multi-Scale Mechanical Properties

Carbon nanotubes (CNT) could serve as potential reinforcement for metal matrix composites for improved mechanical properties. However dispersion of carbon nanotubes (CNT) in the matrix has been a longstanding problem, since they tend to form clusters to minimize their surface area. The aim of this study was to use plasma and cold spraying techniques to synthesize CNT reinforced aluminum composite with improved dispersion and to quantify the degree of CNT dispersion as it influences the mechanical properties. Novel method of spray drying was used to disperse CNTs in Al-12 wt.% Si pre-alloyed powder, which was used as feedstock for plasma and cold spraying. A new method for quantification of CNT distribution was developed. Two parameters for CNT dispersion quantification, namely Dispersion parameter (DP) and Clustering Parameter (CP) have been proposed based on the image analysis and distance between the centers of CNTs. Nanomechanical properties were correlated with the dispersion of CNTs in the microstructure. Coating microstructure evolution has been discussed in terms of splat formation, deformation and damage of CNTs and CNT/matrix interface. Effect of Si and CNT content on the reaction at CNT/matrix interface was thermodynamically and kinetically studied. A pseudo phase diagram was computed which predicts the interfacial carbide for reaction between CNT and Al-Si alloy at processing temperature. Kinetic aspects showed that Al4C3 forms with Al-12 wt.% Si alloy while SiC forms with Al-23wt.% Si alloy. Mechanical properties at nano, micro and macro-scale were evaluated using nanoindentation and nanoscratch, microindentation and bulk tensile testing respectively. Nano and micro-scale mechanical properties (elastic modulus, hardness and yield strength) displayed improvement whereas macro-scale mechanical properties were poor. The inversion of the mechanical properties at different scale length was attributed to the porosity, CNT clustering, CNT-splat adhesion and Al4C3 formation at the CNT/matrix interface. The Dispersion parameter (DP) was more sensitive than Clustering parameter (CP) in measuring degree of CNT distribution in the matrix.

Data to Decision in a Dynamic Ocean: Robust Species Distribution Models and Spatial Decision Frameworks

Human use of the oceans is increasingly in conflict with conservation of endangered species. Methods for managing the spatial and temporal placement of industries such as military, fishing, transportation and offshore energy, have historically been post hoc; i.e. the time and place of human activity is often already determined before assessment of environmental impacts. In this dissertation, I build robust species distribution models in two case study areas, US Atlantic (Best et al. 2012) and British Columbia (Best et al. 2015), predicting presence and abundance respectively, from scientific surveys. These models are then applied to novel decision frameworks for preemptively suggesting optimal placement of human activities in space and time to minimize ecological impacts: siting for offshore wind energy development, and routing ships to minimize risk of striking whales. Both decision frameworks relate the tradeoff between conservation risk and industry profit with synchronized variable and map views as online spatial decision support systems.

For siting offshore wind energy development (OWED) in the U.S. Atlantic (chapter 4), bird density maps are combined across species with weights of OWED sensitivity to collision and displacement and 10 km2 sites are compared against OWED profitability based on average annual wind speed at 90m hub heights and distance to transmission grid. A spatial decision support system enables toggling between the map and tradeoff plot views by site. A selected site can be inspected for sensitivity to a cetaceans throughout the year, so as to capture months of the year which minimize episodic impacts of pre-operational activities such as seismic airgun surveying and pile driving.

Routing ships to avoid whale strikes (chapter 5) can be similarly viewed as a tradeoff, but is a different problem spatially. A cumulative cost surface is generated from density surface maps and conservation status of cetaceans, before applying as a resistance surface to calculate least-cost routes between start and end locations, i.e. ports and entrance locations to study areas. Varying a multiplier to the cost surface enables calculation of multiple routes with different costs to conservation of cetaceans versus cost to transportation industry, measured as distance. Similar to the siting chapter, a spatial decisions support system enables toggling between the map and tradeoff plot view of proposed routes. The user can also input arbitrary start and end locations to calculate the tradeoff on the fly.

Essential to the input of these decision frameworks are distributions of the species. The two preceding chapters comprise species distribution models from two case study areas, U.S. Atlantic (chapter 2) and British Columbia (chapter 3), predicting presence and density, respectively. Although density is preferred to estimate potential biological removal, per Marine Mammal Protection Act requirements in the U.S., all the necessary parameters, especially distance and angle of observation, are less readily available across publicly mined datasets.

In the case of predicting cetacean presence in the U.S. Atlantic (chapter 2), I extracted datasets from the online OBIS-SEAMAP geo-database, and integrated scientific surveys conducted by ship (n=36) and aircraft (n=16), weighting a Generalized Additive Model by minutes surveyed within space-time grid cells to harmonize effort between the two survey platforms. For each of 16 cetacean species guilds, I predicted the probability of occurrence from static environmental variables (water depth, distance to shore, distance to continental shelf break) and time-varying conditions (monthly sea-surface temperature). To generate maps of presence vs. absence, Receiver Operator Characteristic (ROC) curves were used to define the optimal threshold that minimizes false positive and false negative error rates. I integrated model outputs, including tables (species in guilds, input surveys) and plots (fit of environmental variables, ROC curve), into an online spatial decision support system, allowing for easy navigation of models by taxon, region, season, and data provider.

For predicting cetacean density within the inner waters of British Columbia (chapter 3), I calculated density from systematic, line-transect marine mammal surveys over multiple years and seasons (summer 2004, 2005, 2008, and spring/autumn 2007) conducted by Raincoast Conservation Foundation. Abundance estimates were calculated using two different methods: Conventional Distance Sampling (CDS) and Density Surface Modelling (DSM). CDS generates a single density estimate for each stratum, whereas DSM explicitly models spatial variation and offers potential for greater precision by incorporating environmental predictors. Although DSM yields a more relevant product for the purposes of marine spatial planning, CDS has proven to be useful in cases where there are fewer observations available for seasonal and inter-annual comparison, particularly for the scarcely observed elephant seal. Abundance estimates are provided on a stratum-specific basis. Steller sea lions and harbour seals are further differentiated by ‘hauled out’ and ‘in water’. This analysis updates previous estimates (Williams & Thomas 2007) by including additional years of effort, providing greater spatial precision with the DSM method over CDS, novel reporting for spring and autumn seasons (rather than summer alone), and providing new abundance estimates for Steller sea lion and northern elephant seal. In addition to providing a baseline of marine mammal abundance and distribution, against which future changes can be compared, this information offers the opportunity to assess the risks posed to marine mammals by existing and emerging threats, such as fisheries bycatch, ship strikes, and increased oil spill and ocean noise issues associated with increases of container ship and oil tanker traffic in British Columbia’s continental shelf waters.

Starting with marine animal observations at specific coordinates and times, I combine these data with environmental data, often satellite derived, to produce seascape predictions generalizable in space and time. These habitat-based models enable prediction of encounter rates and, in the case of density surface models, abundance that can then be applied to management scenarios. Specific human activities, OWED and shipping, are then compared within a tradeoff decision support framework, enabling interchangeable map and tradeoff plot views. These products make complex processes transparent for gaming conservation, industry and stakeholders towards optimal marine spatial management, fundamental to the tenets of marine spatial planning, ecosystem-based management and dynamic ocean management.