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.

A Delay Vector Variance based marker for an output-only assessment of structural changes in tension leg platforms

Although aspects of power generation of many offshore renewable devices are well understood, their dynamic responses under high wind and wave conditions are still to be investigated to a great detail. Output only statistical markers are important for these offshore devices, since access to the device is limited and information about the exposure conditions and the true behaviour of the devices are generally partial, limited, and vague or even absent. The markers can summarise and characterise the behaviour of these devices from their dynamic response available as time series data. The behaviour may be linear or nonlinear and consequently a marker that can track the changes in structural situations can be quite important. These markers can then be helpful in assessing the current condition of the structure and can indicate possible intervention, monitoring or assessment. This paper considers a Delay Vector Variance based marker for changes in a tension leg platform tested in an ocean wave basin for structural changes brought about by single column dampers. The approach is based on dynamic outputs of the device alone and is based on the estimation of the nonlinearity of the output signal. The advantages of the selected marker and its response with changing structural properties are discussed. The marker is observed to be important for monitoring the as- deployed structural condition and is sensitive to changes in such conditions. Influence of exposure conditions of wave loading is also discussed in this study based only on experimental data.

AUTONOMOUS POWER DISTRIBUTION SYSTEMS

Using robotic systems for many missions that require power distribution can decrease the need for human intervention in such missions significantly. For accomplishing this capability a robotic system capable of autonomous navigation, power systems adaptation, and establishing physical connection needs to be developed. This thesis presents developed path planning and navigation algorithms for an autonomous ground power distribution system. In this work, a survey on existing path planning methods along with two developed algorithms by author is presented. One of these algorithms is a simple path planner suitable for implementation on lab-size platforms. A navigation hierarchy is developed for experimental validation of the path planner and proof of concept for autonomous ground power distribution system in lab environment. The second algorithm is a robust path planner developed for real-size implementation based on lessons learned from lab-size experiments. The simulation results illustrates that the algorithm is efficient and reliable in unknown environments. Future plans for developing intelligent power electronics and integrating them with robotic systems is presented. The ultimate goal is to create a power distribution system capable of regulating power flow at a desired voltage and frequency adaptable to load demands.

Advanced cell culture platforms: methods for drug testing with microfluidics and microstructured devices

Advanced cell cultures are developing rapidly in biomedical research. Nowadays, various approaches and technologies are being used, however, these culturing systems present limitations from increasing complexity, requiring high costs, and not easily customization. We present two versatile and cost-effective methods for developing culturing systems that integrate 3D cell culture and microfluidic platforms. Firstly, for drug screening applications, many high-quality cell spheres of homogeneous size and shape are required. Conventional approaches usually have a dearth of control over the size and geometry of cell spheres and require sample collection and manipulation. To overcome this difficulty, in this study, hundreds of spheroids of several cell lines were generated using multi-well plates that housed our microdevices. Tumor spheroids grow at a uniform rate (in scaffolded or scaffold-free environments) and can be harvested at will. Microscopy imaging are done in real time during or after the culture. After in situ immunostaining, fluorescence imaging can be conducted while keeping the spatial distribution of spheroids in the microwells. Drug effects were successfully observed through viability, growth, and morphologic investigations. Also, we fabricated a microfluidic device suitable for directed and selective cell culture treatments. The microfluidic device was used to reproduce and confirm in vitro investigations carried out using normal culture methods, using a microglia cell line. The device layout and the syringe pump system, entirely designed in our lab, successfully allowed culture growth and medium flow regulation. Solution flows can be finely controlled, allowing treatments and immunofluorescence in one single chamber selectively. To conclude, we propose the development of two culturing platforms (microstructured well devices and in-flow microfluidic chip), which are the result of separate scientific investigations but have the primary goal of performing treatments in a reproducible manner. Our devices shall improve future studies on drug exposure testing, representing adjustable and versatile cell culture systems.

Effect Of Platform Connection And Abutment Material On Stress Distribution In Single Anterior Implant-supported Restorations: A Nonlinear 3-dimensional Finite Element Analysis.

Although various abutment connections and materials have recently been introduced, insufficient data exist regarding the effect of stress distribution on their mechanical performance. The purpose of this study was to investigate the effect of different abutment materials and platform connections on stress distribution in single anterior implant-supported restorations with the finite element method. Nine experimental groups were modeled from the combination of 3 platform connections (external hexagon, internal hexagon, and Morse tapered) and 3 abutment materials (titanium, zirconia, and hybrid) as follows: external hexagon-titanium, external hexagon-zirconia, external hexagon-hybrid, internal hexagon-titanium, internal hexagon-zirconia, internal hexagon-hybrid, Morse tapered-titanium, Morse tapered-zirconia, and Morse tapered-hybrid. Finite element models consisted of a 4×13-mm implant, anatomic abutment, and lithium disilicate central incisor crown cemented over the abutment. The 49 N occlusal loading was applied in 6 steps to simulate the incisal guidance. Equivalent von Mises stress (σvM) was used for both the qualitative and quantitative evaluation of the implant and abutment in all the groups and the maximum (σmax) and minimum (σmin) principal stresses for the numerical comparison of the zirconia parts. The highest abutment σvM occurred in the Morse-tapered groups and the lowest in the external hexagon-hybrid, internal hexagon-titanium, and internal hexagon-hybrid groups. The σmax and σmin values were lower in the hybrid groups than in the zirconia groups. The stress distribution concentrated in the abutment-implant interface in all the groups, regardless of the platform connection or abutment material. The platform connection influenced the stress on abutments more than the abutment material. The stress values for implants were similar among different platform connections, but greater stress concentrations were observed in internal connections.

Censored Linear Regression Models For Irregularly Observed Longitudinal Data Using The Multivariate-t Distribution.

In acquired immunodeficiency syndrome (AIDS) studies it is quite common to observe viral load measurements collected irregularly over time. Moreover, these measurements can be subjected to some upper and/or lower detection limits depending on the quantification assays. A complication arises when these continuous repeated measures have a heavy-tailed behavior. For such data structures, we propose a robust structure for a censored linear model based on the multivariate Student's t-distribution. To compensate for the autocorrelation existing among irregularly observed measures, a damped exponential correlation structure is employed. An efficient expectation maximization type algorithm is developed for computing the maximum likelihood estimates, obtaining as a by-product the standard errors of the fixed effects and the log-likelihood function. The proposed algorithm uses closed-form expressions at the E-step that rely on formulas for the mean and variance of a truncated multivariate Student's t-distribution. The methodology is illustrated through an application to an Human Immunodeficiency Virus-AIDS (HIV-AIDS) study and several simulation studies.

Photoelastic Analysis Of Fixed Partial Prosthesis Crown Height And Implant Length On Distribution Of Stress In Two Dental Implant Systems.

The aim of this study was to evaluate by photoelastic analysis stress distribution on short and long implants of two dental implant systems with 2-unit implant-supported fixed partial prostheses of 8 mm and 13 mm heights. Sixteen photoelastic models were divided into 4 groups: I: long implant (5 × 11 mm) (Neodent), II: long implant (5 × 11 mm) (Bicon), III: short implant (5 × 6 mm) (Neodent), and IV: short implants (5 × 6 mm) (Bicon). The models were positioned in a circular polariscope associated with a cell load and static axial (0.5 Kgf) and nonaxial load (15°, 0.5 Kgf) were applied to each group for both prosthetic crown heights. Three-way ANOVA was used to compare the factors implant length, crown height, and implant system (α = 0.05). The results showed that implant length was a statistically significant factor for both axial and nonaxial loading. The 13 mm prosthetic crown did not result in statistically significant differences in stress distribution between the implant systems and implant lengths studied, regardless of load type (P > 0.05). It can be concluded that short implants showed higher stress levels than long implants. Implant system and length was not relevant factors when prosthetic crown height were increased.

Geographical Distribution Of The Association Between El Niño South Oscillation And Dengue Fever In The Americas: A Continental Analysis Using Geographical Information System-based Techniques.

El Niño South Oscillation (ENSO) is one climatic phenomenon related to the inter-annual variability of global meteorological patterns influencing sea surface temperature and rainfall variability. It influences human health indirectly through extreme temperature and moisture conditions that may accelerate the spread of some vector-borne viral diseases, like dengue fever (DF). This work examines the spatial distribution of association between ENSO and DF in the countries of the Americas during 1995-2004, which includes the 1997-1998 El Niño, one of the most important climatic events of 20(th) century. Data regarding the South Oscillation index (SOI), indicating El Niño-La Niña activity, were obtained from Australian Bureau of Meteorology. The annual DF incidence (AIy) by country was computed using Pan-American Health Association data. SOI and AIy values were standardised as deviations from the mean and plotted in bars-line graphics. The regression coefficient values between SOI and AIy (rSOI,AI) were calculated and spatially interpolated by an inverse distance weighted algorithm. The results indicate that among the five years registering high number of cases (1998, 2002, 2001, 2003 and 1997), four had El Niño activity. In the southern hemisphere, the annual spatial weighted mean centre of epidemics moved southward, from 6° 31' S in 1995 to 21° 12' S in 1999 and the rSOI,AI values were negative in Cuba, Belize, Guyana and Costa Rica, indicating a synchrony between higher DF incidence rates and a higher El Niño activity. The rSOI,AI map allows visualisation of a graded surface with higher values of ENSO-DF associations for Mexico, Central America, northern Caribbean islands and the extreme north-northwest of South America.

Distribution Of Lutzomyia Longipalpis Chemotype Populations In São Paulo State, Brazil.

American visceral leishmaniasis (AVL) is an emerging disease in the state of São Paulo, Brazil. Its geographical expansion and the increase in the number of human cases has been linked to dispersion of Lutzomyia longipalpis into urban areas. To produce more accurate risk maps we investigated the geographic distribution and routes of expansion of the disease as well as chemotype populations of the vector. A database, containing the annual records of municipalities which had notified human and canine AVL cases as well as the presence of the vector, was compiled. The chemotypes of L. longipalpis populations from municipalities in different regions of São Paulo State were determined by Coupled Gas Chromatography - Mass Spectrometry. From 1997 to June 2014, L. longipalpis has been reported in 166 municipalities, 148 of them in the Western region. A total of 106 municipalities were identified with transmission and 99 were located in the Western region, where all 2,204 autochthonous human cases occurred. Both the vector and the occurrence of human cases have expanded in a South-easterly direction, from the Western to central region, and from there, a further expansion to the North and the South. The (S)-9-methylgermacrene-B population of L. longipalpis is widely distributed in the Western region and the cembrene-1 population is restricted to the Eastern region. The maps in the present study show that there are two distinct epidemiological patterns of AVL in São Paulo State and that the expansion of human and canine AVL cases through the Western region has followed the same dispersion route of only one of the two species of the L. longipalpis complex, (S)-9-methylgermacrene-B. Entomological vigilance based on the routes of dispersion and identification of the chemotype population could be used to identify at-risk areas and consequently define the priorities for control measures.

Microbial distribution in the root canal system after periapical lesion induction using different methods

The aim of this study was to evaluate the microbial distribution in the root canal system after periapical lesion induction in dogs' teeth using different methods. Fifty-two root canals were assigned to 4 groups (n=13). Groups I and II: root canals were exposed to the oral cavity for 180 days; groups III and IV: root canals were exposed for 7 days and then the coronal openings were sealed for 53 days. The root apices of groups I and III were perforated, while those of groups II and IV remained intact. After the experimental periods, the animals were euthanized and the anatomic pieces containing the roots were processed and stained with the Brown & Brenn method to assess the presence and distribution of microorganisms. The incidence of microorganisms at different sites of the roots and periapical lesions was analyzed statistically by the chi-square test at 5% significance level. All groups presented microorganisms in the entire root canal system. A larger number of microorganisms was observed on the root canal walls, apical delta and dentinal tubules (p<0.05), followed by cementum and cemental resorption areas. In spite of the different periods of exposure to the oral environment, the methods used for induction of periapical periodontitis yielded similar distribution of microorganisms in the root canal system.

Distribution of biotypes and leukotoxic activity of Aggregatibacter actinomycetemcomitans isolated from Brazilian patients with chronic periodontitis

Aggregatibacter actinomycetemcomitans is an important etiologic agent of the periodontitis and is associated with extra-oral infections. In this study, the detection of the ltxA gene as well as the ltx promoter region from leukotoxic A. actinomycetemcomitans isolated from 50 Brazilian patients with periodontitis and 50 healthy subjects was performed. The leukotoxic activity on HL-60 cells was also evaluated. Leukotoxic activity was determined using a trypan blue exclusion method. The 530 bp deletion in the promoter region was evaluated by PCR using a PRO primer pair. A. actinomycetemcomitans was detected by culture and directly from crude subgingival biofilm by PCR using specific primers. By culture, A. actinomycetemcomitans was detected in nine (18%) of the periodontal patients and one (2%) healthy subject. However, by PCR, this organism was detected in 44% of the periodontal patients and in 16% of the healthy subjects. It was verified a great discrepancy between PCR detection of the ltx operon promoter directly from crude subgingival biofilm and from bacterial DNA. Only one periodontal sample harbored highly leukotoxic A. actinomycetemcomitans. Moreover, biotype II was the most prevalent and no correlation between biotypes and leukotoxic activity was observed. The diversity of leukotoxin expression by A. actinomycetemcomitans suggests a role of this toxin in the pathogenesis of periodontal disease and other infectious diseases.

Inequalities in the distribution of dental caries among 12-year-old Brazilian schoolchildren

We assessed the inequality in the distribution of dental caries and the association between indicators of socioeconomic status and caries experience in a representative sample of schoolchildren. This study followed a cross-sectional design, with a sample of 792 schoolchildren aged 12 years, representative of this age group in Santa Maria, RS, Brazil. Guardians answered questions on socioeconomic status and a dental examination provided information on the dental caries experience (DMF-T). Inequality in dental caries distribution was measured by the Gini coefficient and the Significant Caries Index (SiC). The assessment of association used Poisson regression models. Socioeconomic factors were associated with prevalence of dental caries for the whole sample and also for individuals with a high-caries level. Children from low-income households had the highest prevalence of dental caries. The Gini coefficient was 0.7 and the SiC Index 2.5. The percentage of caries prevalence was 39.3% (95% CI: 35.8%-42.8%) and the mean for DMF-T was 0.9 (± SD 1.5). Inequalities in the distribution of dental caries were observed and socioeconomic factors were found to be strong predictors of the prevalence of oral disease in children of this age group.

A new species of the genus Pimelerodius Vanin from the Amazon Region, with notes on the geografic distribution of Pimelerodius motacilla (Boheman) (Coleoptera, Curculionidae, Otidocephalini)

Pimelerodius punctiventris sp. nov. (type locality Brazil, Amazonas, Itacoatiara) is described and illustrated. The new taxon is compared with similar species, being distinguished from the other 12 known species of the genus by the presence of punctures in ventrite I. The available published key for identification of species of Pimelerodius is adapted to include the new species. A modification of the generic description of the aedeagus of Pimelerodius is provided, a necessity due to the differences observed in the aedeagus of the new species. The occurrence of P. motacilla (Boheman, 1843) in the Amazon Region, recorded in sympatry with P. punctiventris in Itacoatiara, AM, is discussed and confirmed, based on the study of 41 available specimens.

Glycosaminoglycan distribution in the rat uterine cervix during the estrous cycle

OBJECTIVE: To analyze the amount of glycosaminoglycans in the uterine cervix during each phase of the rat estrous cycle. DESIGN: Based on vaginal smears, forty female, regularly cycling rats were divided into four groups (n = 10 for each group): GI - proestrous, GII - estrous, GIII - metaestrous and GIV - diestrous. Animals were sacrificed at each phase of the cycle, and the cervix was immediately removed and submitted to biochemical extraction and determination of sulfated glycosaminoglycans and hyaluronic acid. The results were analyzed by ANOVA followed by the Bonferroni post-hoc test. RESULTS: The uterine cervix had the highest amount of total sulfated glycosaminoglycans and dermatan sulfate during the estrous phase (8.90 ± 0.55 mg/g of cetonic extract, p<0.001; and 8.86 ± 0.57 mg/g of cetonic extract, p<0.001). In addition, there was more heparan sulfate at the cervix during the proestrous phase (0.185 ± 0.03 mg/g of cetonic extract) than during any other phase (p<0.001). There were no significant changes in the concentration of hyaluronic acid in the uterine cervix during the estrous cycle. CONCLUSION: Our data suggest that the amount of total sulfated glycosaminoglycans may be influenced by hormonal fluctuations related to the estrous cycle, with dermatan sulfate and heparan sulfate being the glycosaminoglycans most sensitive to hormonal change.

Stress distribution in the cervical region of an upper central incisor in a 3D finite element model

The aim of this study was to evaluate the stress distribution in the cervical region of a sound upper central incisor in two clinical situations, standard and maximum masticatory forces, by means of a 3D model with the highest possible level of fidelity to the anatomic dimensions. Two models with 331,887 linear tetrahedral elements that represent a sound upper central incisor with periodontal ligament, cortical and trabecular bones were loaded at 45º in relation to the tooth's long axis. All structures were considered to be homogeneous and isotropic, with the exception of the enamel (anisotropic). A standard masticatory force (100 N) was simulated on one of the models, while on the other one a maximum masticatory force was simulated (235.9 N). The software used were: PATRAN for pre- and post-processing and Nastran for processing. In the cementoenamel junction area, tensile forces reached 14.7 MPa in the 100 N model, and 40.2 MPa in the 235.9 N model, exceeding the enamel's tensile strength (16.7 MPa). The fact that the stress concentration in the amelodentinal junction exceeded the enamel's tensile strength under simulated conditions of maximum masticatory force suggests the possibility of the occurrence of non-carious cervical lesions such as abfractions.