An evaluation of home range models for marine fish tracking and fine scale habitat use and movement patterns of age 1 Greenland cod (Gadus macrocephalus ogac)

My thesis examines fine-scale habitat use and movement patterns of age 1 Greenland cod (Gadus macrocephalus ogac) tracked using acoustic telemetry. Recent advances in tracking technologies such as GPS and acoustic telemetry have led to increasingly large and detailed datasets that present new opportunities for researchers to address fine-scale ecological questions regarding animal movement and spatial distribution. There is a growing demand for home range models that will not only work with massive quantities of autocorrelated data, but that can also exploit the added detail inherent in these high-resolution datasets. Most published home range studies use radio-telemetry or satellite data from terrestrial mammals or avian species, and most studies that evaluate the relative performance of home range models use simulated data. In Chapter 2, I used actual field-collected data from age-1 Greenland cod tracked with acoustic telemetry to evaluate the accuracy and precision of six home range models: minimum convex polygons, kernel densities with plug-in bandwidth selection and the reference bandwidth, adaptive local convex hulls, Brownian bridges, and dynamic Brownian bridges. I then applied the most appropriate model to two years (2010-2012) of tracking data collected from 82 tagged Greenland cod tracked in Newman Sound, Newfoundland, Canada, to determine diel and seasonal differences in habitat use and movement patterns (Chapter 3). Little is known of juvenile cod ecology, so resolving these relationships will provide valuable insight into activity patterns, habitat use, and predator-prey dynamics, while filling a knowledge gap regarding the use of space by age 1 Greenland cod in a coastal nursery habitat. By doing so, my thesis demonstrates an appropriate technique for modelling the spatial use of fish from acoustic telemetry data that can be applied to high-resolution, high-frequency tracking datasets collected from mobile organisms in any environment.

Iron speciation and magnetic susceptibility in bottom sediments of the Black Sea, Core SMG91-903

12 cores of Late Pleistocene - Holocene deposits were studied. They were collected by gravity cores on the continental slope and in the deep-water part of the Black Sea within the Adler-Tuapse polygon. In four of them in New Euxinian deposits at the base of a packet of hydrotroilite laminae paleomagnetic anomalies likely resulting from the Gothenburg magnetic excursion occur. Comparison with results of similar studies in the western Black Sea, where the Gothenburg magnetic excursion was previously found, let to validate stratigraphic synchronism of the hydrotroilite horizon in the eastern and western parts of the Black Sea and to confirm the authors' views about peculiarities of paleogeographical development of the Black Sea basin in the Late Pleistocene - Holocene.

Composition of bottom sediments from the Clarion Fracture Zone

New data on lithology and stratigraphy of Cenozoic sediments from the Clarion Transform Fault Zone (Pacific Ocean) have been obtained on the base of polygon studies. It has been established that on different blocks (uplifted and subsided) of the Clarion tectonic structure deposits of different age (Eocene to Quaternary) occur. Unconsolidated sediments have been deposited under pelagic conditions since Eocene (probably, since Early Cretaceous) until now. Their mineral composition and content of different ore components are given.

Optimized square Fresnel zone plates for microoptics applications

Polygonal Fresnel zone plates with a low number of sides have deserved attention in micro and nanoptics, because they can be straightforwardly integrated in photonic devices, and, at the same time, they represent a balance between the high-focusing performance of a circular zone plate and the easiness of fabrication at micro and nano-scales of polygons. Among them, the most representative family are Square Fresnel Zone Plates (SFZP). In this work, we propose two different customized designs of SFZP for optical wavelengths. Both designs are based on the optimization of a SFZP to perform as close as possible as a usual Fresnel Zone Plate. In the first case, the criterion followed to compute it is the minimization of the difference between the area covered by the angular sector of the zone of the corresponding circular plate and the one covered by the polygon traced on the former. Such a requirement leads to a customized polygon-like Fresnel zone. The simplest one is a square zone with a pattern of phases repeating each five zones. On the other hand, an alternative SFZP can be designed guided by the same criterion but with a new restriction. In this case, the distance between the borders of different zones remains unaltered. A comparison between the two lenses is carried out. The irradiance at focus is computed for both and suitable merit figures are defined to account for the difference between them.

Optimal phase distributions for polygonal Fresnel lenses

Polygonal Fresnel zone plates can be configured in a variety of forms depending on the number of sides of the polygon and the number of phase steps used. This contribution deals with some specific polygonal designs that tessellate the plane: triangles, squares, and hexagons. The phase distribution is chosen as a continuous one to form a polygonal kinoform. The selected designs have been simulated and its behaviour compared. Although their performance is worse than the circular Fresnel plate, they may present some other advantages as the tessellation capability, and the possibility to fabricate them as extruded profiles.

On Provable Algorithms for Determination of Continuous Protein Interdomain Motions from Residual Dipolar Couplings

Dynamics of biomolecules over various spatial and time scales are essential for biological functions such as molecular recognition, catalysis and signaling. However, reconstruction of biomolecular dynamics from experimental observables requires the determination of a conformational probability distribution. Unfortunately, these distributions cannot be fully constrained by the limited information from experiments, making the problem an ill-posed one in the terminology of Hadamard. The ill-posed nature of the problem comes from the fact that it has no unique solution. Multiple or even an infinite number of solutions may exist. To avoid the ill-posed nature, the problem needs to be regularized by making assumptions, which inevitably introduce biases into the result.

Here, I present two continuous probability density function approaches to solve an important inverse problem called the RDC trigonometric moment problem. By focusing on interdomain orientations we reduced the problem to determination of a distribution on the 3D rotational space from residual dipolar couplings (RDCs). We derived an analytical equation that relates alignment tensors of adjacent domains, which serves as the foundation of the two methods. In the first approach, the ill-posed nature of the problem was avoided by introducing a continuous distribution model, which enjoys a smoothness assumption. To find the optimal solution for the distribution, we also designed an efficient branch-and-bound algorithm that exploits the mathematical structure of the analytical solutions. The algorithm is guaranteed to find the distribution that best satisfies the analytical relationship. We observed good performance of the method when tested under various levels of experimental noise and when applied to two protein systems. The second approach avoids the use of any model by employing maximum entropy principles. This 'model-free' approach delivers the least biased result which presents our state of knowledge. In this approach, the solution is an exponential function of Lagrange multipliers. To determine the multipliers, a convex objective function is constructed. Consequently, the maximum entropy solution can be found easily by gradient descent methods. Both algorithms can be applied to biomolecular RDC data in general, including data from RNA and DNA molecules.

Advanced Applications of 3D Dosimetry and 3D Printing in Radiation Therapy

As complex radiotherapy techniques become more readily-practiced, comprehensive 3D dosimetry is a growing necessity for advanced quality assurance. However, clinical implementation has been impeded by a wide variety of factors, including the expense of dedicated optical dosimeter readout tools, high operational costs, and the overall difficulty of use. To address these issues, a novel dry-tank optical CT scanner was designed for PRESAGE 3D dosimeter readout, relying on 3D printed components and omitting costly parts from preceding optical scanners. This work details the design, prototyping, and basic commissioning of the Duke Integrated-lens Optical Scanner (DIOS).

The convex scanning geometry was designed in ScanSim, an in-house Monte Carlo optical ray-tracing simulation. ScanSim parameters were used to build a 3D rendering of a convex ‘solid tank’ for optical-CT, which is capable of collimating a point light source into telecentric geometry without significant quantities of refractive-index matched fluid. The model was 3D printed, processed, and converted into a negative mold via rubber casting to produce a transparent polyurethane scanning tank. The DIOS was assembled with the solid tank, a 3W red LED light source, a computer-controlled rotation stage, and a 12-bit CCD camera. Initial optical phantom studies show negligible spatial inaccuracies in 2D projection images and 3D tomographic reconstructions. A PRESAGE 3D dose measurement for a 4-field box treatment plan from Eclipse shows 95% of voxels passing gamma analysis at 3%/3mm criteria. Gamma analysis between tomographic images of the same dosimeter in the DIOS and DLOS systems show 93.1% agreement at 5%/1mm criteria. From this initial study, the DIOS has demonstrated promise as an economically-viable optical-CT scanner. However, further improvements will be necessary to fully develop this system into an accurate and reliable tool for advanced QA.

Pre-clinical animal studies are used as a conventional means of translational research, as a midpoint between in-vitro cell studies and clinical implementation. However, modern small animal radiotherapy platforms are primitive in comparison with conventional linear accelerators. This work also investigates a series of 3D printed tools to expand the treatment capabilities of the X-RAD 225Cx orthovoltage irradiator, and applies them to a feasibility study of hippocampal avoidance in rodent whole-brain radiotherapy.

As an alternative material to lead, a novel 3D-printable tungsten-composite ABS plastic, GMASS, was tested to create precisely-shaped blocks. Film studies show virtually all primary radiation at 225 kVp can be attenuated by GMASS blocks of 0.5cm thickness. A state-of-the-art software, BlockGen, was used to create custom hippocampus-shaped blocks from medical image data, for any possible axial treatment field arrangement. A custom 3D printed bite block was developed to immobilize and position a supine rat for optimal hippocampal conformity. An immobilized rat CT with digitally-inserted blocks was imported into the SmART-Plan Monte-Carlo simulation software to determine the optimal beam arrangement. Protocols with 4 and 7 equally-spaced fields were considered as viable treatment options, featuring improved hippocampal conformity and whole-brain coverage when compared to prior lateral-opposed protocols. Custom rodent-morphic PRESAGE dosimeters were developed to accurately reflect these treatment scenarios, and a 3D dosimetry study was performed to confirm the SmART-Plan simulations. Measured doses indicate significant hippocampal sparing and moderate whole-brain coverage.

On surfaces with prescribed shape operator

The problem of immersing a simply connected surface with a prescribed shape operator is discussed. I show that, aside from some special degenerate cases, such as when the shape operator can be realized by a surface with one family of principal curves being geodesic, the space of such realizations is a convex set in an affine space of dimension at most 3. The cases where this maximum dimension of realizability is achieved are analyzed and it is found that there are two such families of shape operators, one depending essentially on three arbitrary functions of one variable and another depending essentially on two arbitrary functions of one variable. The space of realizations is discussed in each case, along with some of their remarkable geometric properties. Several explicit examples are constructed.

Optimal Capital Requirements in Financial Networks with Fire Sales

I explore and analyze a problem of finding the socially optimal capital requirements for financial institutions considering two distinct channels of contagion: direct exposures among the institutions, as represented by a network and fire sales externalities, which reflect the negative price impact of massive liquidation of assets.These two channels amplify shocks from individual financial institutions to the financial system as a whole and thus increase the risk of joint defaults amongst the interconnected financial institutions; this is often referred to as systemic risk. In the model, there is a trade-off between reducing systemic risk and raising the capital requirements of the financial institutions. The policymaker considers this trade-off and determines the optimal capital requirements for individual financial institutions. I provide a method for finding and analyzing the optimal capital requirements that can be applied to arbitrary network structures and arbitrary distributions of investment returns.

In particular, I first consider a network model consisting only of direct exposures and show that the optimal capital requirements can be found by solving a stochastic linear programming problem. I then extend the analysis to financial networks with default costs and show the optimal capital requirements can be found by solving a stochastic mixed integer programming problem. The computational complexity of this problem poses a challenge, and I develop an iterative algorithm that can be efficiently executed. I show that the iterative algorithm leads to solutions that are nearly optimal by comparing it with lower bounds based on a dual approach. I also show that the iterative algorithm converges to the optimal solution.

Finally, I incorporate fire sales externalities into the model. In particular, I am able to extend the analysis of systemic risk and the optimal capital requirements with a single illiquid asset to a model with multiple illiquid assets. The model with multiple illiquid assets incorporates liquidation rules used by the banks. I provide an optimization formulation whose solution provides the equilibrium payments for a given liquidation rule.

I further show that the socially optimal capital problem using the ``socially optimal liquidation" and prioritized liquidation rules can be formulated as a convex and convex mixed integer problem, respectively. Finally, I illustrate the results of the methodology on numerical examples and

discuss some implications for capital regulation policy and stress testing.

Le processus de thermo-érosion du pergélisol dans la zone de pergélisol continu

Les polygones à coin de glace sont très répandus dans la zone du pergélisol continu. Lorsque le ruissellement d’eau de fonte nivale s’infiltre de façon concentrée dans une cavité, il peut initier le processus de thermo-érosion du pergélisol (notamment des coins de glace) pouvant mener à la formation de ravins. Dans la vallée de Qalikturvik sur l’Ile Bylot (NU, Canada), le développement de ravins de thermo-érosion dans un milieu de polygones à coins de glace entraîne comme impact : i. la réorganisation des réseaux de drainage impliquant un assèchement des milieux humides en marge des chenaux d’érosion, ii. des variations dans le régime thermique et de l’humidité de proche-surface et iii. la prise en charge et le déplacement des sédiments vers l’extérieur du bassin-versant. L’objectif de cette thèse vise à approfondir les connaissances géomorphologiques propres au ravinement par thermo-érosion, d’examiner, caractériser et quantifier les impacts du ravinement (tel que sus-mentionné en i. ii. iii.) et le rôle de celui-ci dans une optique d’évolution du paysage périglaciaire à l’échelle temporelle de l’année à la décennie. Les ravins sont dynamiques : un ravin en particulier déclenché en 1999 et étudié depuis s’érodait à une vitesse de 38 à 50 m/a durant sa première décennie d’existence, pour atteindre une longueur totale de ~750 m et une surface érodée de ~25 000 m² en 2009. Des puits sont localisés près des zones de ravinement actives ; des levées alluviale, mares et polygones effondrés dans les zones stabilisées post-perturbation. Sur la terrasse de polygones recouvrant le plancher de la vallée au site à l’étude, 35 ravins furent identifiés et 1401 polygones furent perturbés avec 200 000 m³ de sols transportés. Une amélioration du drainage, une dégradation de la capacité de rétention de l’humidité, une transition d’un écoulement de ruissellement vers un écoulement canalisé caractérise les aires ravinées et leurs environs. Les polygones intacts sont homogènes d’un à l’autre et dans leurs centres ; les polygones perturbés ont une réponse hétérogène (flore, humidité et régime thermique). Les milieux érodés hétérogènes succèdent aux milieux homogènes et deviennent le nouvel état d’équilibre pour plusieurs décennies.

Radiocarbon dating, carbon storage and soil properties of samples from Tulemalu Lake, central Canadian Arctic

We investigated total storage and landscape partitioning of soil organic carbon (SOC) in continuous permafrost terrain, central Canadian Arctic. The study is based on soil chemical analyses of pedons sampled to 1 m depth at 35 individual sites along three transects. Radiocarbon dating of cryoturbated soil pockets, basal peat and fossil wood shows that cryoturbation processes have been occurring since the Middle Holocene and that peat deposits started to accumulate in a forest-tundra environment where spruce was present (~6000 cal yrs BP). Detailed partitioning of SOC into surface organic horizons, cryoturbated soil pockets and non-cryoturbated mineral soil horizons is calculated (with storage in active layer and permafrost calculated separately) and explored using principal component analysis. The detailed partitioning and mean storage of SOC in the landscape are estimated from transect vegetation inventories and a land cover classification based on a Landsat satellite image. Mean SOC storage in the 0-100 cm depth interval is 33.8 kg C/m**2, of which 11.8 kg C/m**2 is in permafrost. Fifty-six per cent of the total SOC mass is stored in peatlands (mainly bogs), but cryoturbated soil pockets in Turbic Cryosols also contribute significantly (17%). Elemental C/N ratios indicate that this cryoturbated soil organic matter (SOM) decomposes more slowly than SOM in surface O-horizons.

Description of manganese deposits from the Papago mine, Saipan, Mariana Islands

Saipan, situated about 15° N. and 146° E., is one of the larger and more southerly of the Mariana Islands. The 15 small islands of this chain are strung along an eastwardly convex ridge for more than 400 miles north to south, midway between Honshu and New Guinea and about 1,200 miles east of the Philippines. Paralleling this ridge 60 to 100 miles further east is a deep submarine trench, beyond which lies the Pacific Basin proper. To the west is the Philippine Sea, generally deeper than 2,000 fathoms. The trench coincides with a zone of negative gravity anomalies, earthquake foci occur at increasing depths westward from it, and silica- and alumina-rich volcanic rocks characterize the emergent island chain itself. The contrast between these features and those of the Pacific Basin proper to the east is held to favor the conclusion that the Mariana island arc and trench define the structural and petrographic front of Asia.