In the late Pleistocene, most of British Columbia and northern Washington was covered by the Cordilleran ice sheet. The weight of the ice sheet caused up to several hundred metres of depression of the Earth’s crust. This caused relative sea level to be higher in southwestern British Columbia despite lower global eustatic sea level. After deglaciation, postglacial rebound of the crust caused sea level to quickly drop to below present levels. The rate of sea-level fall is used here to determine the rheology of the mantle in southwestern British Columbia. The first section of this study deals with determination of the postglacial sea-level history in the Victoria area. Constraints on sea-level position come from isolation basin cores collected in 2000 and 2001, as well as from previously published data from the past 45 years. The position of sea-level is well constrained at elevations greater than -4 m, and there are only loose constraints below that. The highstand position in the Victoria area is between 75-80 m. Sea level fell rapidly from the highstand position to below 0 m between 14.3 and 13.2 thousand calendar years before present (cal kyr BP). The magnitude of the lowstand position was between -11 and -40 m. Though there are few constraints on the lowstand position, analysis of the crustal response favours larger lowstand. Well constrained sea-level histories from Victoria, central Strait of Georgia and northern Strait of Georgia are used to model the rheology of the mantle in southwestern British Columbia. A new ice sheet model for the southwestern Cordillera was developed as older models systematically underpredicted the magnitude of sea level in late glacial times. Radiocarbon dates are compiled to provide constraints on ice sheet advance and retreat. The Cordillera ice sheet reached maximum extent between 17 and 15.4 cal kyr BP. After 15.4 cal kyr, the ice sheet retreated, and by 13.7 cal kyr BP Puget Sound, Juan de Fuca Strait and Strait of Georgia were ice free. By 10.7 cal kyr BP, ice was restricted to mountain glaciers at levels similar to present. With the new ice model, and using an Earth model with a 60 km lithosphere, asthenosphere with variable viscosity and thickness, and transitional and lower mantle viscosity based on the VM2 Earth model, predicted sea level matches the observed sea level constraints in southwestern British Columbia. Nearly identical predicted sea-level curves are found using asthenosphere thicknesses between 140-380 km with viscosity values between 3x10^18 and 4x10^19 Pa s. Predicted sea level is almost completely insensitive to the mantle below the asthenosphere. Modeled present day postglacial uplift rates are less than 0.5 mm yr^-1. Despite the tight fit of the predicted sea level to observed late-glacial sea level observations, the modelling was not able to fit the early Holocene rise of sea level to present levels in the central and northern Strait of Georgia.

The study focuses mainly on two gas hydrate-related targets, located on the Northern Cascadia Margin, offshore Vancouver Island: (1) a recently identified 70-80-m high carbonate mound, Cucumber Ridge, located ~3.5-km west of Ocean Drilling Program (ODP) Site 889 and Integrated Ocean Drilling Program (IODP) Site U1327, and (2) a large cold vent, Bullseye vent, which is up to ~500 m in diameter and was drilled by IODP at Site U1328. The objective of this thesis is to analyze seismic data that provide indicators of locally focused fluid flow and characteristics of the gas hydrate occurrence associated with these two features. A grid of closely-spaced single channel seismic (SCS) data was collected at Cucumber Ridge in July/August 2001, and deep-towed multichannel seismic (MCS) lines were collected using Deep-towed Acoustics and Geophysics System (DTAGS) at the Bullseye vent area and at Cucumber Ridge in October 2002. The high-resolution SCS data, with a frequency bandpass of 40-150 Hz, recorded coherent reflectivity down to about 400 m beneath the seafloor, and provide excellent images of the subseafloor structure of Cucumber Ridge and of the gas hydrate bottom-simulating reflector (BSR) beneath it. Cucumber Ridge is interpreted to have developed as a structural topographic high in the hanging wall of a large reverse fault formed at the base of the current seaward slope. The fault zone provides pathways for fluids including gas to migrate to the seafloor where diagenetic carbonate forms and cements the near-surface sediments. Over the seismic grid, heat flow was derived from the depth of the BSR. A simple 2-D analytical correction for theoretical heat flow variations due to topography is applied to the data. Across the mound, most of the variability in heat flow is explained by topographic effects, including a local 6 mW/m2 negative anomaly over the central mound and a large 20 mW/m2 positive anomaly over the mound steep side slope. However, just south of the mound, there is a 6-7 mW/m2 positive anomaly in a 2-km-long band that has predominantly flat seafloor. Most of this anomaly is probably unrelated to topographic effects, but rather likely due to warm upward fluid flow along faults or fracture zones. Towed ~300 m above seafloor, the high frequency (220-1k Hz) DTAGS signal can provide high vertical resolution images with increased lateral resolution. The major problems of DTAGS are significant nonlinear variations of the source depths and receivers locations. New routines were developed for optimal DTAGS data processing, mainly including (1) cable geometry estimation by node depths, direct arrivals and seasurface reflections using a Genetic Algorithm inversion method, (2) acoustic image stitching based on accurate relative-source positioning by crosscorrelation of redundant data between two adjacent shots, and (3) velocity inversion of wide-angle traveltimes using a nonlinear global grid search method. The final processed DTAGS images resolve multiple seismic blanking zones and fine details of subseafloor features in the slope sediments. At Bullseye vent, where a 35-m-thick near-surface massive hydrate layer was drilled at U1328, the DTAGS data resolved the upper part of layer as a dipping diffraction zone, likely corresponding to a fracture zone. The inverted velocity structure in upper 100 m sediments successfully revealed a 17-m-thick layer of high velocity (~1650 m/s) just below seafloor, probably related to carbonate presence. A local high velocity zone, with a positive velocity anomaly of ~40-80 m/s in the upper 50 m beneath seafloor, was observed over the ~100-m wide region between U1328 and the deepest part of a seafloor depression; the high velocity zone is consistent with the dipping diffraction zone in the DTAGS image and with the massive hydrate drilled at U1328.

In this thesis, acoustic inversion methods are employed to estimate array element locations and the geoacoustic properties of the seabed using measured acoustic data consisting of noise from a surface ship in the Gulf of Mexico. The array element localization utilizes relative travel-time information obtained by cross-correlating the recorded time series of ship noise received at spatially separated hydrophones. The relative travel-time data are used in an inversion, based on the regularized least-squares method and the acoustic ray tracing equations, to obtain improved estimates of the receiver and source positions and their uncertainties. Optimization and Bayesian matched-field inversion methods are employed to estimate seabed geoacoustic properties and their uncertainties in the vicinity of a bottom-moored vertical line array using the recorded surface ship noise. This study is used to test the feasibility of matched-field methods to detect temporal changes in the geoacoustic properties of the seabed near a known gas hydrate mound in the Gulf of Mexico. Finally, a synthetic study is performed that demonstrates how ignoring environmental range dependence of seabed sound speed and water depth in matched-field inversion can lead to biases in the estimated geoacoustic parameters. The study considers the distributions of optimal parameter estimates obtained from a large number of range-independent inversions of synthetic data generated for random range-dependent environments. Range-independent Bayesian inversions are also performed on selected data sets and the marginal parameter distributions are examined. Both hard- and soft-bottom environments are examined at a number of scales of variability in sound speed and water depth.

A 616-year Douglas-fir (Pseudotsuga menziesii) chronology was developed to examine the history of drought and moisture variability in the Sooke Watershed, near Victoria, British Columbia. Ring-width chronologies were compared to historical precipitation, air temperature and drought variables (Palmer Drought Severity Index (PDSI) and Standardized Precipitation Index (SPI)) to determine the climate/radial-growth response to moisture stress on the sampled stands. Correlations between the ring-width chronologies and climate variables revealed that May to July precipitation, May-June SPI and July PDSI were significant limiting factors to radial-width growth. A transfer function was established for each of these variables to create a proxy climate reconstruction of drought in the watershed. The summer precipitation model provided the most accurate representation of past moisture variability (R2 = 0.20) and reveals substantial variation in precipitation over the past six centuries. Evidence from the periodicity of the tree-ring record to suggest that some modes of atmospheric circulation are influencing precipitation supply to the watershed.

This thesis reports eight measured sections and >400 new detrital zircon U-Pb SHRIMP-II ages from the Mesoproterozoic (~1.4 Ga) upper Purcell Supergroup of southeastern British Columbia, Canada. The goal of my study is to constrain the depositional, tectonic and paleogeographic setting of the upper Purcell Supergroup. Stratigraphic sections across the Purcell Anticlinorium, constructed from measured sections, reveal three syn-depositional growth faults: (1) paleo-Hall Lake, (2) paleo-Larchwood Lake, and (3) paleo-Moyie. Stratigraphic sections were combined into a fence diagram, revealing a large north-northeast trending graben bound to the east by the paleo-Larchwood Lake fault and to the west by the paleo-Hall Lake fault. Five samples were collected for detrital zircon analysis along the eastern extent of exposed Purcell strata; one sample was collected from the western limit of strata. All samples are characterized by subordinate numbers of detrital zircons that yield Paleoproterozoic and Archean ages. Detrital zircon ages from the Sheppard Formation are dominated by 1500, 1700, 1750, and 1850 Ma grains. The overlying Gateway Formation is dominated by 1400-1450, 1700, 1850, and 1900 Ma zircon grains. The overlying Phillips, Roosville (east), and Mount Nelson formations are dominated by detrital zircon ages between 1375-1450 Ma and 1650-1800 Ma. Detrital zircon ages from the Roosville Formation (west) are dominated by 1500-1625 Ma grains. Based on the margin perpendicular orientation of the long axis of syn-depositional grabens relative to Laurentia, and on the presence of syn-depositional aged zircons through the entire sedimentary succession, we interpret the upper Purcell Supergroup to have been deposited in a transpressional pull-apart basin setting, adjacent to a convergent/translational plate margin bound to the west by terranes now located in northeastern Australia.

Planktonic ecosystems provide a key mechanism for the transfer of CO2 from the atmosphere to the deep ocean via the so-called "biological pump". Mathematical models of these ecosystems have been used to predict CO2 uptake in surface waters, and more recently have been embedded in global climate models. While the equilibrium properties of these models are well studied, less attention has been paid to their response to external perturbations, despite the fact that as a result of the variability of environmental forcing such ecosystems are rarely, if ever, in equilibrium. Human induced perturbations to these ecosystems, namely the addition of limiting nutrients (e.g. iron) to areas where nitrate is plentiful to accelerate the biological pump, have been proposed as a solution to reduce atmospheric CO2. Linear theory is used to determine the structure of "unit-norm" perturbations (size in mmol N m^-3) to state variables of an ecosystem model in steady state, describing Ocean Station P (50N 145W) in summer, that optimize either instantaneous export flux of organic matter at fixed times or integrated export as the ecosystem relaxes towards equilibrium. For all perturbations, the flux to higher trophic levels is the primary contributor to export flux, the contribution of aggregation is negligible, and (sinking) detritus increases significantly in the transient dynamics. Two perturbations considered optimize instantaneous export flux; both perturbations synchronize P1 and Z1 relative to their predator prey cycle, resulting in a maximum instantaneous export flux of 4.4 mmol N m^-2 d^-1, and also increased integrated export above that at steady state (6 g C m^-2 over 150 days). An increase in larger phytoplankton (P2), representing diatoms, results in the highest integrated export (7 g C m^-2). The perturbations in which P2 persist the longest give the highest integrated export, and these perturbations are primarily increases in P2. The additional integrated export in response to a proportional increase to steady state concentrations of both large and small phytoplankton is positive, but much lower than the optimal perturbations. However, the additional integrated export in response to an increase in only P1 is negligible. The linear and nonlinear ecosystem and export responses to two perturbations are compared; for perturbations of magnitude 0.5 mmol N m^-3, the linearization of the ecosystem dynamics, rather than of the export flux, is the primary cause for differences between the fully linear and fully nonlinear cases.

Current meter records from a long term mooring site on the continental slope off the west coast of Vancouver Island, British Columbia, Canada are used to investigate the scales of variability of the subsurface California Undercurrent and its relation to possible driving mechanisms. Observed along the west coast of North America from Baja California to Vancouver Island, the California Undercurrent is part of the California Current System, a typical basin-scale eastern boundary circulation system. Of the four instruments at nominal depths of 35, 100, 175, and 400 m, the upper two show seasonally reversing flow, while the 175 m instrument registers a year-round poleward flow. The deepest current meter, located approximately 100 m above the bottom, reflects the influence of a nearby submarine canyon. The flow at 100 and 175 m depths, as well as the water properties sampled in the region with CTD casts, are characteristic of the temporal and spatial variability of the California Undercurrent over the continental slope off central and southern Vancouver Island. The correlation of the 175 m flow with local atmospheric forcing (wind stress) in the low-frequency band (periods of months) is higher than with ocean-wide climatic indices, suggesting that regional processes play a key role in the forcing of the subsurface flow.

Understanding the role of forests in the cycling of methane (CH4) and carbon dioxide (CO2) is of importance to the elucidation of global greenhouse gas budgets. Previous studies have shown aerated forests soils to be net sinks of atmospheric CH4 and sources of carbon dioxide. While much research has focused on the role of forest soils as CO2 sources and CH4 sinks, few studies have utilized 13C-isotope studies to clarify the nature of subsurface CO2 production and CH4 consumption. The present study, carried out in 3 temperate forest environments on Vancouver Island during 2006 and 2007, and a boreal forest in northern Quebec in 2005, is intended to address this paucity of information. The isotope and concentration data corroborates previous studies in suggesting that both temperate and boreal forest environments act as net CH4 sinks and CO2 sources. No clear evidence of methanogenesis is apparent in either Vancouver Island or northern Quebec, where the isotopic composition of subsurface CH4 is influenced by diffusive and biological fractionation. Near-surface photosynthetic uptake may have a strong influence on the isotopic composition of soil CO2 and the resultant fluxes, acting to reduce apparent fluxes due to CO2 consumption. Intra-site variability of CH4 and CO2 fluxes indicates that the use of two static chambers in a single site, while sufficient for the confirmation of gas uptake or emission, may be less adequate in the determination of actual rates of efflux/influx. Future studies should address this by either sampling a larger area, installing a greater number of chambers, or by utilizing entirely different methods, such as the use of eddy covariance techniques.

To understand the mechanisms of energy exchange between the tropics and the midlatitudes, it is necessary to develop simplified climate models. Motivated by linear wave theory, one such model is derived below. It captures the nonlinear interaction between barotropic and first baroclinic modes. In particular, it allows for the study of the barotropic response to a baroclinic forcing. Numerical methods for handling this nonlinear system are carefully developed and validated. The response generated by a physically realistic Kelvin wave forcing is studied and is found to consist mainly of one eastward propagating wave (phase-locked to the forcing) and two westward propagating (Rossby) waves. The Rossby waves are shown to be highly constrained by the initial parameters of the forcing and an explanation of this result is proposed.

Generalized linear stability theory is used to calculate the optimal initial conditions that result in transient amplification of the thermohaline circulation (THC) in a zonally-averaged single basin ocean model. The eigenmodes of the tangent linear model verify that the system is asymptotically stable but the nonnormality of the system permits the growth of perturbations for a finite period through the interference of nonorthogonal eigenmodes. It is found that the maximum amplification of the THC anomalies occurs after 6 years with both the thermally driven and salinity driven components playing major roles in the amplification process. The transient amplification of THC anomalies is due to the constructive and destructive interference of a large number of eigenmodes and the evolution over time is determined by how the interference pattern evolves. It is found that five of the most highly nonnormal eigenmodes are critical to the initial cancellation of the salinity and temperature contributions to the THC while 11 oscillating modes with decay timescales ranging from 2 to 6 years are the major contributors at the time of maximum amplification. This analysis demonstrates that the different dynamics of salinity and temperature anomalies allows the dramatic growth of perturbations to the THC on relatively short (interannual to decadal) timescales. In addition the ideas of generalized stability theory are used to calculate the stochastic optimals which are the spatial patterns of stochastic forcing that are most efficient at generating variance growth in the THC. It is found that the optimal stochastic forcing occurs at high latitudes and induces low-frequency THC variability by exciting the salinity-dominated modes of the THC. The first stochastic optimal is found to have its largest projection on the same five highly nonnormal eigenmodes found to be critical to the structure of the optimal initial conditions. The model’s response to stochastic forcing is not controlled by the least damped eigenmodes of the tangent linear model but rather by the linear interference of these highly nonnormal eigenmodes. The process of pseudoresonance suggests that the nonnormal eigenmodes are excited and sustained by stochastically induced perturbations which in turn lead to maximum THC variance. Finally, it was shown that the addition of wind stress did not have a large impact on the nonnormal dynamics of the linearised system. Adding wind allowed the value of the vertical diffusivity to be reduced to achieve the same maximum linearised THC amplitude as was used in the case with no wind stress.