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.

Spatial distribution of bednet coverage under routine distribution through the public health sector in a rural district in Kenya.

Insecticide-treated nets (ITNs) are one of the most important and cost-effective tools for malaria control. Maximizing individual and community benefit from ITNs requires high population-based coverage. Several mechanisms are used to distribute ITNs, including health facility-based targeted distribution to high-risk groups; community-based mass distribution; social marketing with or without private sector subsidies; and integrating ITN delivery with other public health interventions. The objective of this analysis is to describe bednet coverage in a district in western Kenya where the primary mechanism for distribution is to pregnant women and infants who attend antenatal and immunization clinics. We use data from a population-based census to examine the extent of, and factors correlated with, ownership of bednets. We use both multivariable logistic regression and spatial techniques to explore the relationship between household bednet ownership and sociodemographic and geographic variables. We show that only 21% of households own any bednets, far lower than the national average, and that ownership is not significantly higher amongst pregnant women attending antenatal clinic. We also show that coverage is spatially heterogeneous with less than 2% of the population residing in zones with adequate coverage to experience indirect effects of ITN protection.

Distribution, ecology, life history, genetic variation, and risk of extinction of nonhuman primates from Costa Rica

We examined the association between geographic distribution, ecological traits, life history, genetic diversity, and risk of extinction in nonhuman primate species from Costa Rica. All of the current nonhuman primate species from Costa Rica are included in the study; spider monkeys (Ateles geoffroyi), howling monkeys (Alouatta palliata), capuchins (Cebus capucinus), and squirrel monkeys (Saimiri oerstedii). Geographic distribution was characterized accessing existing databases. Data on ecology and life history traits were obtained through a literature review. Genetic diversity was characterized using isozyme electrophoresis. Risk of extinction was assessed from the literature. We found that species differed in all these traits. Using these data, we conducted a Pearson correlation between risk of extinction and ecological and life history traits, and genetic variation, for widely distributed species. We found a negative association between risk of extinction and population birth and growth rates; indicating that slower reproducing species had a greater risk of extinction. We found a positive association between genetic variation and risk of extinction; i.e., species showing higher genetic variation had a greater risk of extinction. The relevance of these traits for conservation efforts is discussed.

The distribution of sales revenues from pharmaceutical innovation.

OBJECTIVE: This report updates our earlier work on the returns to pharmaceutical research and development (R&D) in the US (1980 to 1984), which showed that the returns distributions are highly skewed. It evaluates a more recent cohort of new drug introductions in the US (1988 to 1992) and examines how the returns distribution is emerging for drugs with life cycles concentrated in the 1990s versus the 1980s. DESIGN AND SETTING: Methods were described in detail in our earlier reports. The current sample included 110 new drug entities (including 28 orphan drugs), and sales data were obtained for the period 1988 to 1998, which represented between 7 and 11 years of sales for the drugs included. 20 years was chosen as the expected market life for this cohort, and a 2-step procedure was used to project future sales for the drugs--during the period until patent expiry and then beyond patent expiry until the 20-year time-horizon was completed. Thus, the values in the first half of the life cycle are essentially based on realised sales, while those in the second half are projected using information on patent expiry and other inputs. MAIN OUTCOME MEASURES AND RESULTS: Peak annual sales for the top decile of drugs introduced between 1988 and 1992 in the US amounted to almost $US1.1 billion compared with peak sales of less than $US175 million (1992 values) for the mean compound. In particular, the top decile accounted for 56% of overall sales revenue. Although the sales distributions were skewed in both our earlier and current analysis, the top decile in the later time-period exhibited more rapid rates of growth after launch, a peak that was more than 50% greater in real terms than for the 1980 to 1984 cohort, and a faster rate of expected decline in sales after patent expiry. One factor contributing to the distribution of sales revenues becoming more skewed over time is the orphan drug phenomenon (i.e. most of the orphan drugs are concentrated at the bottom of the distribution). CONCLUSION: The distribution of sales revenues for new drug compounds is highly skewed in nature. In this regard, the top decile of new drugs accounts for more than half of the total sales generated by the 1988 to 1992 cohort analysed. Furthermore, the distribution of sales revenues for this cohort is more skewed than that of the 1980 to 1984 cohort we analysed in previous research.

Convergent Surface Water Distributions in U.S. Cities

Earth's surface is rapidly urbanizing, resulting in dramatic changes in the abundance, distribution and character of surface water features in urban landscapes. However, the scope and consequences of surface water redistribution at broad spatial scales are not well understood. We hypothesized that urbanization would lead to convergent surface water abundance and distribution: in other words, cities will gain or lose water such that they become more similar to each other than are their surrounding natural landscapes. Using a database of more than 1 million water bodies and 1 million km of streams, we compared the surface water of 100 US cities with their surrounding undeveloped land. We evaluated differences in areal (A WB) and numeric densities (N WB) of water bodies (lakes, wetlands, and so on), the morphological characteristics of water bodies (size), and the density (D C) of surface flow channels (that is, streams and rivers). The variance of urban A WB, N WB, and D C across the 100 MSAs decreased, by 89, 25, and 71%, respectively, compared to undeveloped land. These data show that many cities are surface water poor relative to undeveloped land; however, in drier landscapes urbanization increases the occurrence of surface water. This convergence pattern strengthened with development intensity, such that high intensity urban development had an areal water body density 98% less than undeveloped lands. Urbanization appears to drive the convergence of hydrological features across the US, such that surface water distributions of cities are more similar to each other than to their surrounding landscapes. © 2014 The Author(s).

The distribution of early childhood memories.

The quantitative distribution of autobiographical memories for the first decade of life is described. The distribution, based on over 11,000 autobiographical memories from age 10 and younger from published studies, is nearly identical for males and females, for participants of different ages, and for different methods of collecting data, including using words to cue memories from anywhere in the lifespan or from just the childhood years, exhaustive listing of all early memories, and interviews.

The distribution of autobiographical memories across the lifespan.

Words were used to cue autobiographical memories from 20- and 70-year-old subjects. Both groups showed a decrease in memories from the childhood years and a power-function retention function for their most recent 10 years. Older subjects also had an increase in the number of memories from the ages 10 to 30. These results held for individual subjects as well as grouped data and held when either 124 or 921 memories were cued. Reaction times to produce memories were constant across decades except for childhood where they were longer.

Application of a rank-based genetic association test to age-at-onset data from the Collaborative Study on the Genetics of Alcoholism study.

Association studies of quantitative traits have often relied on methods in which a normal distribution of the trait is assumed. However, quantitative phenotypes from complex human diseases are often censored, highly skewed, or contaminated with outlying values. We recently developed a rank-based association method that takes into account censoring and makes no distributional assumptions about the trait. In this study, we applied our new method to age-at-onset data on ALDX1 and ALDX2. Both traits are highly skewed (skewness > 1.9) and often censored. We performed a whole genome association study of age at onset of the ALDX1 trait using Illumina single-nucleotide polymorphisms. Only slightly more than 5% of markers were significant. However, we identified two regions on chromosomes 14 and 15, which each have at least four significant markers clustering together. These two regions may harbor genes that regulate age at onset of ALDX1 and ALDX2. Future fine mapping of these two regions with densely spaced markers is warranted.