A rapidly melting ice sheet produces a distinctive geometry of sea level (SL) change. Thus, a network of SL observations may, in principle, be used to infer sources of meltwater flux. We outline a new method, based on a Kalman smoother, for using tide gauge observations to estimate the individual sources of global SL change. The Kalman smoother technique iteratively calculates the maximum likelihood estimate of Greenland and West Antarctic ice sheet melt rates at each time step, and it allows for data gaps while also permitting the estimation of non-linear trends. We have also implemented a fixed multi-model Kalman filter that allows us to rigorously account for additional contributions to SL changes, such as glacial isostatic adjustment and thermal expansion. We report on a series of detection experiments based on synthetic SL data that explore the feasibility of extracting source information from SL records before applying the new methodology to historical tide gauge records. In the historical tide gauge study we infer a global mean SL rise of ~1.5 ± 0.5 mm/yr up to 1970, followed by an acceleration to a rate of ~2.0 ± 0.5 mm/yr in 2008. In addition to its connection to SL, Greenland and its large ice sheet act as a barrier to storm systems traversing the North Atlantic. As a result of the interaction with Greenland, low-pressure systems located in the Irminger Sea, between Iceland and Greenland, often produce strong low-level winds. Through a combination of modeling and the analysis of rare in-situ observations, we explore the evolution of a lee cyclone that resulted in three high-speed-wind events in November 2004. Understanding Greenland’s role in these events is critical in our understanding of local weather in this region.

The Earth's magnetic field is generated in its fluid outer core through dynamo action. In this process, convection and differential rotation of an electrically conducting fluid maintain the magnetic field against its ohmic decay. Using numerical models, we can investigate planetary dynamo processes and the importance of various core properties on the dynamo. In this thesis, I use numerical dynamo models in Earth-like geometry in order to understand the influence of inner core electrical conductivity and the choice of thermal and velocity boundary conditions on the resulting magnetic field. I demonstrate how an electrically conducting inner core can reduce the frequency of reversals and produce axial-dipolar dominated fields in our models. I also demonstrate that a strong planetary magnetic field intensity does not imply that the dynamo operates in the strong field regime as is usually presumed. Through a scaling law analysis, I find that irrespective of the choice of thermal or velocity boundary conditions, the available power determines the magnetic and velocity field characteristics like the field strength, polarity and morphology. Also, whether a dynamo model is in a dipolar, transitional or multipolar regime is dependent on the force balance in the model. I demonstrate that the Lorentz force is balanced by the Coriolis force in the dipolar dynamo regime models resulting in magnetostrophically balanced dynamos whereas the Lorentz force is balanced by the Inertial force (and not the Coriolis force) in the multipolar dynamo regime models resulting in a non-magnetostrophically balanced dynamo. The generation mechanism differs between the regimes and depends on the velocity boundary conditions. The zonal flows of the stress-free models are stronger than in the no-slip models, and bistability is more prominent when stress-free boundary conditions are used. A single scaling law may be feasible for all the models, but there does appear to be some variation for models with different thermal and velocity boundary conditions. The results presented in this thesis are not only applicable to the geodynamo, but will also aid in understanding the dynamos of other planets and exoplanets.

The Earth's magnetic field is generated in its fluid outer core through dynamo action. In this process, convection and differential rotation of an electrically conducting fluid maintain the magnetic field against its ohmic decay. Using numerical models, we can investigate planetary dynamo processes and the importance of various core properties on the dynamo. In this thesis, I use numerical dynamo models in Earth-like geometry in order to understand the influence of inner core electrical conductivity and the choice of thermal and velocity boundary conditions on the resulting magnetic field. I demonstrate how an electrically conducting inner core can reduce the frequency of reversals and produce axial-dipolar dominated fields in our models. I also demonstrate that a strong planetary magnetic field intensity does not imply that the dynamo operates in the strong field regime as is usually presumed. Through a scaling law analysis, I find that irrespective of the choice of thermal or velocity boundary conditions, the available power determines the magnetic and velocity field characteristics like the field strength, polarity and morphology. Also, whether a dynamo model is in a dipolar, transitional or multipolar regime is dependent on the force balance in the model. I demonstrate that the Lorentz force is balanced by the Coriolis force in the dipolar dynamo regime models resulting in magnetostrophically balanced dynamos whereas the Lorentz force is balanced by the Inertial force (and not the Coriolis force) in the multipolar dynamo regime models resulting in a non-magnetostrophically balanced dynamo. The generation mechanism differs between the regimes and depends on the velocity boundary conditions. The zonal flows of the stress-free models are stronger than in the no-slip models, and bistability is more prominent when stress-free boundary conditions are used. A single scaling law may be feasible for all the models, but there does appear to be some variation for models with different thermal and velocity boundary conditions. The results presented in this thesis are not only applicable to the geodynamo, but will also aid in understanding the dynamos of other planets and exoplanets.

The combined finite-discrete element method (FEM/DEM) has been used to simulate processes of brittle fracturing and associated seismicity. With the newly extended FEM/DEM algorithm, two topics involving rock mechanics and geophysics are investigated. In the first topic, a velocity-weakening law is implemented to investigate the initiation of frictional slip, and an innovative method that incorporates surface roughness with varying friction coefficients is introduced to examine the influences of surface roughness. Simulated results revealed detailed responses of stresses to the propagation of the slip front. In the second topic, acoustic activities induced in confined compression tests are simulated and quantitatively studied using the internal monitoring algorithm in FEM/DEM. It is shown that with increasing confinement, AE events are spatially more concentrated and temporally more separated, accompanied by a decreasing b-value. Moreover, interesting correlation between orientations of cracks and the mechanical behavior of the rock was observed.

The combined finite-discrete element method (FEM/DEM) has been used to simulate processes of brittle fracturing and associated seismicity. With the newly extended FEM/DEM algorithm, two topics involving rock mechanics and geophysics are investigated. In the first topic, a velocity-weakening law is implemented to investigate the initiation of frictional slip, and an innovative method that incorporates surface roughness with varying friction coefficients is introduced to examine the influences of surface roughness. Simulated results revealed detailed responses of stresses to the propagation of the slip front. In the second topic, acoustic activities induced in confined compression tests are simulated and quantitatively studied using the internal monitoring algorithm in FEM/DEM. It is shown that with increasing confinement, AE events are spatially more concentrated and temporally more separated, accompanied by a decreasing b-value. Moreover, interesting correlation between orientations of cracks and the mechanical behavior of the rock was observed.

Gas hydrates are a solid, ice-like mixture of water and low molecular weight hydrocarbons. They are found under the permafrost and to a far greater extent under the ocean, usually at water depths greater than 300m. Hydrates are a potential energy resource, a possible factor in climate change, and a geohazard. For these reasons, it is critical that gas hydrate deposits are quantitatively assessed so that their concentrations, locations and distributions may be established. Due to their ice-like nature, hydrates are electrically insulating. Consequently, a method which remotely detects changes in seafloor electrical conductivity, such as marine controlled source electromagnetics (CSEM), is a useful geophysical tool for marine gas hydrate exploration. Hydrates are geometrically complex structures. Advanced electromagnetic modelling and imaging techniques are crucial for proper survey design and data interpretation. I develop a method to model thin resistive structures in conductive host media which may be useful in building approximate geological models of gas hydrate deposits using arrangements of multiple, bent sheets. I also investigate the possibility of interpreting diffusive electromagnetic data using seismic imaging techniques. To be processed in this way, such data must first be transformed into its non-diffusive, seismic-like counterpart. I examine such a transform from both an analytical and a numerical point of view, focusing on methods to overcome inherent numerical instabilities. This is the first step to applying seismic processing techniques to CSEM data to rapidly and efficiently image resistive gas hydrate structures. The University of Toronto marine electromagnetics group has deployed a permanent marine CSEM array offshore Vancouver Island, in the framework of the NEPTUNE Canada cabled observatory, for the purposes of monitoring gas hydrate deposits. In this thesis I also propose and examine a new CSEM survey technique for gas hydrate which would make use of the stationary seafloor transmitter already on the seafloor, along with a cabled receiver array, towed from a ship. I furthermore develop a modelling algorithm to examine the electromagnetic effects of conductive borehole casings which have been proposed to be placed in the vicinity of this permanent marine CSEM array, and make preliminary recommendations about their locations.

The goal of this research was to devise a new analytical technique, using Accelerator Mass Spectrometry (AMS), to measure Platinum Group Element (PGE) concentrations to the sup-ppb levels in natural, insulating, samples. The challenges were threefold. First, a method of sputtering an insulating sample to successfully produce a stable beam of anions needed to be devised. Second, a suitable standard of known PGE concentrations had to be found and third, spectral analysis of the beam had to verify any claims of PGE abundance. The first challenge was met by employing a modified high intensity negative ion source flooded with neutral caesium that successfully sputtered insulators to produce a beam of negative ions. The second challenge, that of finding a suitable standard, was fraught with difficulties, as no synthesized standards available were found to be appropriate for this work. As a result, direction is provided for future production of standards by ion implantation. The third challenge, successful spectral analysis, was accomplished using a newly designed gas ionization detector which allowed for resolution of the interfering molecular fragment from the PGE ions. Coupled with the use of the SRIM computer programme, positive identification of all peaks in the spectra of the analyzed samples was accomplished. The success of the first and third challenges lead to the qualitative analyses of geological samples for sub-ppb levels of PGE by AMS. Quantitative analyses await only for the appropriate standards and with those will come a whole new range of research possibilities for measuring sub-ppb levels of PGE in insulating samples by AMS.

Gas hydrates are a solid, ice-like mixture of water and low molecular weight hydrocarbons. They are found under the permafrost and to a far greater extent under the ocean, usually at water depths greater than 300m. Hydrates are a potential energy resource, a possible factor in climate change, and a geohazard. For these reasons, it is critical that gas hydrate deposits are quantitatively assessed so that their concentrations, locations and distributions may be established. Due to their ice-like nature, hydrates are electrically insulating. Consequently, a method which remotely detects changes in seafloor electrical conductivity, such as marine controlled source electromagnetics (CSEM), is a useful geophysical tool for marine gas hydrate exploration. Hydrates are geometrically complex structures. Advanced electromagnetic modelling and imaging techniques are crucial for proper survey design and data interpretation. I develop a method to model thin resistive structures in conductive host media which may be useful in building approximate geological models of gas hydrate deposits using arrangements of multiple, bent sheets. I also investigate the possibility of interpreting diffusive electromagnetic data using seismic imaging techniques. To be processed in this way, such data must first be transformed into its non-diffusive, seismic-like counterpart. I examine such a transform from both an analytical and a numerical point of view, focusing on methods to overcome inherent numerical instabilities. This is the first step to applying seismic processing techniques to CSEM data to rapidly and efficiently image resistive gas hydrate structures. The University of Toronto marine electromagnetics group has deployed a permanent marine CSEM array offshore Vancouver Island, in the framework of the NEPTUNE Canada cabled observatory, for the purposes of monitoring gas hydrate deposits. In this thesis I also propose and examine a new CSEM survey technique for gas hydrate which would make use of the stationary seafloor transmitter already on the seafloor, along with a cabled receiver array, towed from a ship. I furthermore develop a modelling algorithm to examine the electromagnetic effects of conductive borehole casings which have been proposed to be placed in the vicinity of this permanent marine CSEM array, and make preliminary recommendations about their locations.

A rapidly melting ice sheet produces a distinctive geometry of sea level (SL) change. Thus, a network of SL observations may, in principle, be used to infer sources of meltwater flux. We outline a new method, based on a Kalman smoother, for using tide gauge observations to estimate the individual sources of global SL change. The Kalman smoother technique iteratively calculates the maximum likelihood estimate of Greenland and West Antarctic ice sheet melt rates at each time step, and it allows for data gaps while also permitting the estimation of non-linear trends. We have also implemented a fixed multi-model Kalman filter that allows us to rigorously account for additional contributions to SL changes, such as glacial isostatic adjustment and thermal expansion. We report on a series of detection experiments based on synthetic SL data that explore the feasibility of extracting source information from SL records before applying the new methodology to historical tide gauge records. In the historical tide gauge study we infer a global mean SL rise of ~1.5 ± 0.5 mm/yr up to 1970, followed by an acceleration to a rate of ~2.0 ± 0.5 mm/yr in 2008. In addition to its connection to SL, Greenland and its large ice sheet act as a barrier to storm systems traversing the North Atlantic. As a result of the interaction with Greenland, low-pressure systems located in the Irminger Sea, between Iceland and Greenland, often produce strong low-level winds. Through a combination of modeling and the analysis of rare in-situ observations, we explore the evolution of a lee cyclone that resulted in three high-speed-wind events in November 2004. Understanding Greenland’s role in these events is critical in our understanding of local weather in this region.

Master of Science / Department of Geology / Abdelmoneam Raef / Seismic reservoir characterization and prospect evaluation based 3D seismic attributes analysis in Kansas has been successful in contributing to the tasks of building static and dynamic reservoir models and in identifying commercial hydrocarbon prospects. In some areas, reservoir heterogeneities introduce challenges, resulting in some wells with poor economics. Analysis of seismic attributes gives insight into hydrocarbon presence, fluid movement (in time lapse mode), porosity, and other factors used in evaluating reservoir potential. This study evaluates a producing lease using seismic attributes analysis of an area covered by a 2010 3D seismic survey in the Morrison Northeast field and Morrison field of Clark County, KS. The target horizon is the Viola Limestone, which continues to produce from seven of twelve wells completed within the survey area. In order to understand reservoir heterogeneities, hydrocarbon entrapment settings and the implications for future development plans, a seismic attributes extraction and analysis, guided with geophysical well-logs, was conducted with emphasis on instantaneous attributes and amplitude anomalies. Investigations into tuning effects were conducted in light of amplitude anomalies to gain insight into what seismic results led to the completion of the twelve wells in the area drilled based on the seismic survey results. Further analysis was conducted to determine if the unsuccessful wells completed could have been avoided. Finally the study attempts to present a set of 3D seismic attributes associated with the successful wells, which will assist in placing new wells in other locations within the two fields, as well as promote a consistent understanding of entrapment controls in this field.