SlideShare a Scribd company logo
Exploratory spatial analysis of  illegal oil discharges detected off Canada’s Pacific Coast Norma Serra Rosaline Canessa Patrick O’Hara Stefania Bertazzon  Marina Gavrilova   ICCSA Conference  2008
What are illegal vessel-source oil discharges? Photo courtesy of Gov. of Canada Photo courtesy of P. O’Hara Mystery Spill
The ecological impact in marine ecosystems of small-scale chronic-levels of oil pollution is shown to be greater than large accidental oil spills; especially for seabird populations.   ~ 80% global population of Cassin’s Auklets in BC ~   74% global population of Ancient Murrelets ~ 56% global population Rhinoceros Auklet   In British Columbia there is a rich diversity of seabirds at least 62 different species, including major proportions of world breeding populations: ~ 25% global population of Marbled Murrelets
National Aerial Surveillance Program  (NASP) Primary tool of Transport Canada for preventing and monitoring illegal oil pollution within Canada’s Economic Exclusive Zone, which generally extents 200 nautical miles (370 km) out from its coast . The Vancouver Sun. February 16 th , 2008
Research Goal Specify a multivariate regression model to predict  where illegal oil spills  are more likely to occur, to provide NASP with a  decision making tool  to maximize the effectiveness of their surveillance program. In order to identify a valid model, in necessary to explore the spatial properties of our data:  a) Spatial Heterogeneity  b) Spatial Dependence or Spatial    Autocorrelation (SAC)
Study Area Maximum bounding area covered by surveillance flights off the west coast of Vancouver Island, British Columbia. 1
Data sources From NASP : 53 detected oil spills in BC from 1997 to 2006 Surveillance flight path data in BC from 1997 to 2006. From Marine Communications Traffic Services - Vessel Traffic Operations Support System: Relative marine traffic densities in BC were estimated for different vessel type groups, from data collected in 2003.
Marine shipping traffic in  Canada’s West Coast Bulk carrier
Marine shipping traffic in  Canada’s West Coast Oil tanker
Marine shipping traffic in  Canada’s West Coast Tug vessel
Marine shipping traffic in  Canada’s West Coast Fishing vessels
Marine shipping traffic in  Canada’s West Coast Cruise ship
Methods Global measures of spatial autocorrelation   examine the nature and extent of the dependence within model variables and produce a single value for the entire data set. The two most common global measures: I ~ + 1    Positive SAC I ~ - 1    Negative SAC I =  0    No SAC Moran’s  I Geary's  c 0<C< 1    Positive SAC C >  1    Negative SAC C =  1    No SAC
Methods (Cont.) Local measures of spatial association – Local Moran’s  I identify the location and spatial scale of aggregations of unusual values, such as clusters of high values (hot spots) and low values (cold spots) and, because it generates an autocorrelation index for each data site, this can be mapped providing additional information about the pattern under study. E.g.  LISA Cluster Map  (GeoDA TM ) Outlier Outlier Cluster Cluster high values surrounded by low neighboring values High - Low low values surrounded by high neighboring values  Low – High low values surrounded by neighbors of similar low values  Low – Low high values surrounded by neighbors of similar high values  High – High
Results:   Spatial autocorrelation – Global indices 0.871 (0.01) 0.158 (0.01) Oil spills counts 0.455 (0.01) 0.518 (0.01) Cruise ship movement counts  0.371 (0.01) 0.553 (0.01) Tug vessel movement counts  0.395 (0.01) 0.553 (0.01) Carrier vessel movement counts  0.369 (0.01) 0.589 (0.01) Oil tankers movement counts  0.249 (0.01) 0.783 (0.01) Fishing vessel movement counts 0.172 (0.01) 0.805 (0.01) NASP flight counts Geary's  c  ( p -values) Moran's  I  ( p -values)
Results (Cont.): Spatial pattern description - LISA Cluster Maps   LISA cluster map of  surveillance flight counts.   Significance test using 999 permutations and applying a significance filter of  p  ≤ 0.01.
Results (Cont.): Spatial pattern description - LISA Cluster Maps   LISA cluster map of  detected oil spills counts.   Significance test using 999 permutations and applying a significance filter of  p  ≤ 0.01.
Results (Cont.): Spatial pattern description - LISA Cluster Maps   LISA cluster map of  carrier vessel movement counts.   Significance test using 999 permutations and applying a significance filter of  p  ≤ 0.01.
Results (Cont.): Spatial pattern description - LISA Cluster Maps   LISA cluster map of  oil tanker movement counts.   Significance test using 999 permutations and applying a significance filter of  p  ≤ 0.01.
Results (Cont.): Spatial pattern description - LISA Cluster Maps   LISA cluster map of  tug vessel movement counts.   Significance test using 999 permutations and applying a significance filter of  p  ≤ 0.01.
Results (Cont.): Spatial pattern description - LISA Cluster Maps   LISA cluster map of  fishing vessel movement counts.   Significance test using 999 permutations and applying a significance filter of  p  ≤ 0.01.
Results (Cont.): Spatial pattern description - LISA Cluster Maps   LISA cluster map of  cruise ships movement counts.   Significance test using 999 permutations and applying a significance filter of  p  ≤ 0.01.
Discussion Our exploratory analyses indicate that there is a positive spatial autocorrelation within datasets for all variables (dependent and independent).  LISA cluster maps allowed the visualization of significant aggregations of high values of oil spill counts, surveillance effort and shipping traffic.  To some extent LISA cluster maps show the effect of deterrence. For example, oil spill clusters were found in areas of Barkley Sound, whereas clusters of high counts of surveillance flights and marine vessels were mainly found in the entrance to the Strait of Juan de Fuca.
Conclusion (I) Based on our results:  Our variables exhibit spatial dependencies and non-stationarities,  therefore, spatially autoregressive (SAR) and  geographically weighted regression (GWR) models  are the most suitable to predict spatial patterns of illegal oil discharges based on marine commercial vessel type and movement patterns, while controlling for the spatial distribution of surveillance flights.
Conclusion (II) Based on our results:  We propose to calibrate a multivariate regression model, combining a spatial autoregressive and a local specification. Our model will estimate oil spill probabilities based only on areas poorly covered by the current NASP flight zone,  In order to provide effective and reliable recommendations for the planned surveillance program expansion.
Explore the degree of correlation between detected oil spills and  marine traffic densities for different vessel types (i.e. fishers, tugs, cargo vessels, tankers, cruise ships and ferries). Incorporate results from O’Hara  et al. * that defines the relationship between oil spill detection probabilities and surveillance effort.  Explore oil spill probabilities using a multivariate approach that will include variables such as season, distances between oil spills and shore and distances to nearest port or marina, and other characteristics that define the vessel type (i.e., flag state, inbound vs. outbound)  Future research * O’Hara, P.D., Serra-Sogas, N., Canessa, R., Keller, P., Pelot, P.: Estimating oil spill rates and deterrence based on aerial surveillance data  in Western Canadian marine waters.  Marine Pollution Bulletin, (submitted) (2008).
Thank you for your attention     Any questions? Acknowledgements:

More Related Content

PPTX
Australian seagrass habitats. Kathryn McMahon, ACEAS Grand 2014
PPT
Water Emergency Transportation Authority
PDF
HEC - RAS Final Project
PDF
Ship Emission Rating Index (SERI)
PPTX
Aswathy final(4)
PPTX
Rinkal.cpd.ppt
Australian seagrass habitats. Kathryn McMahon, ACEAS Grand 2014
Water Emergency Transportation Authority
HEC - RAS Final Project
Ship Emission Rating Index (SERI)
Aswathy final(4)
Rinkal.cpd.ppt

Viewers also liked (20)

PDF
Ijcatr04061005
PPTX
Data Philly Meetup - Big (Geo) Data
PDF
Web mapping with vector data. Is it the future ? 2012
PPTX
Spatial queries entity recognition and disambiguation
PPTX
Introduction to Oracle Spatial
PPT
3D Visibility with Vector GIS Data
PPTX
Spatial enhancement
PDF
Exploratory Spatial Analysis using GeoDa
KEY
Spatial Analytics, Where 2.0 2010
PDF
Components of Spatial Data Quality in GIS
PPT
Big Social Data: The Spatial Turn in Big Data (Video available soon on YouTube)
PDF
Spatial data analysis 1
PPT
Introduction To Gis With Employment Info
PDF
Dmitriy Kolesov - GIS as an environment for integration and analysis of spati...
PDF
QGIS Module 1
PPT
Vectors and Rasters
PPT
Spatial data mining
PPT
Vector analysis
PPT
4.2 spatial data mining
PPT
Raster Vs. Vector Presentation1
Ijcatr04061005
Data Philly Meetup - Big (Geo) Data
Web mapping with vector data. Is it the future ? 2012
Spatial queries entity recognition and disambiguation
Introduction to Oracle Spatial
3D Visibility with Vector GIS Data
Spatial enhancement
Exploratory Spatial Analysis using GeoDa
Spatial Analytics, Where 2.0 2010
Components of Spatial Data Quality in GIS
Big Social Data: The Spatial Turn in Big Data (Video available soon on YouTube)
Spatial data analysis 1
Introduction To Gis With Employment Info
Dmitriy Kolesov - GIS as an environment for integration and analysis of spati...
QGIS Module 1
Vectors and Rasters
Spatial data mining
Vector analysis
4.2 spatial data mining
Raster Vs. Vector Presentation1
Ad

Similar to Exploratory Spatial Analysis Norma (20)

PPTX
Embedding Risk in Everything we do
PPTX
DA 2 NDMM.pptx
PPT
igarss staples july 2011.ppt
PPTX
Extra_Li_XF_2011_IGARSS_OilSpill.pptx
DOCX
Meeting the Threat of MaritimeImprovised Explosive Devices.docx
PDF
Hydrographic Surveying : Safe Navigation
PDF
FVCOM's MAIN CHARACTERISTICS AND FEATURES
PPTX
IBA Presentation Combo Updated 3-28
PDF
PPTX
Introduction of Motion Control and requirements of Navigations systems techno...
PPTX
Introduction of Motion Control and requirements of Navigations systems techno...
PPTX
Group Final Presentation.ppt (1)
PPT
Exceeding Expectations: Port Readiness Gary smith
PDF
LindseyHamilton_NocturnalSurfScoterDistributionPoster_2015
PPTX
WAVE Project Overview
PDF
D3602026041
PDF
Reserve Estimation of Initial Oil and Gas by using Volumetric Method in Mann ...
PDF
Optimal routing development based on real voyage data presented by_sewonkim
PPTX
WAVE Project Presentation at the University of Gibraltar
PDF
G04702074087
Embedding Risk in Everything we do
DA 2 NDMM.pptx
igarss staples july 2011.ppt
Extra_Li_XF_2011_IGARSS_OilSpill.pptx
Meeting the Threat of MaritimeImprovised Explosive Devices.docx
Hydrographic Surveying : Safe Navigation
FVCOM's MAIN CHARACTERISTICS AND FEATURES
IBA Presentation Combo Updated 3-28
Introduction of Motion Control and requirements of Navigations systems techno...
Introduction of Motion Control and requirements of Navigations systems techno...
Group Final Presentation.ppt (1)
Exceeding Expectations: Port Readiness Gary smith
LindseyHamilton_NocturnalSurfScoterDistributionPoster_2015
WAVE Project Overview
D3602026041
Reserve Estimation of Initial Oil and Gas by using Volumetric Method in Mann ...
Optimal routing development based on real voyage data presented by_sewonkim
WAVE Project Presentation at the University of Gibraltar
G04702074087
Ad

More from Beniamino Murgante (20)

PDF
Analyzing and assessing ecological transition in building sustainable cities
PDF
Smart Cities: New Science for the Cities
PDF
The evolution of spatial analysis and modeling in decision processes
PDF
Smart City or Urban Science?
PDF
Involving citizens in smart energy approaches: the experience of an energy pa...
PPTX
Programmazione per la governance territoriale in tema di tutela della biodive...
PDF
Involving Citizens in a Participation Process for Increasing Walkability
PDF
Presentation of ICCSA 2019 at the University of Saint petersburg
PPTX
RISCHIO TERRITORIALE NEL GOVERNO DEL TERRITORIO: Ricerca e formazione nelle s...
PDF
Presentation of ICCSA 2017 at the University of trieste
PPT
GEOGRAPHIC INFORMATION – NEED TO KNOW (GI-N2K) Towards a more demand-driven g...
PPTX
Focussing Energy Consumers’ Behaviour Change towards Energy Efficiency and Lo...
PPTX
Socio-Economic Planning profiles: Sciences VS Daily activities in public sector 
PPT
GEOGRAPHIC INFORMATION – NEED TO KNOW (GI-N2K) Towards a more demand-driven g...
PPT
Garden in motion. An experience of citizens involvement in public space regen...
PPT
Planning and Smartness: the true challenge
PDF
GeoSDI: una piattaforma social di dati geografici basata sui principi di INSP...
PPT
Murgante smart energy
PPT
Informazione Geografica, Città, Smartness
PPT
Tecnologie, Territorio, Smartness
Analyzing and assessing ecological transition in building sustainable cities
Smart Cities: New Science for the Cities
The evolution of spatial analysis and modeling in decision processes
Smart City or Urban Science?
Involving citizens in smart energy approaches: the experience of an energy pa...
Programmazione per la governance territoriale in tema di tutela della biodive...
Involving Citizens in a Participation Process for Increasing Walkability
Presentation of ICCSA 2019 at the University of Saint petersburg
RISCHIO TERRITORIALE NEL GOVERNO DEL TERRITORIO: Ricerca e formazione nelle s...
Presentation of ICCSA 2017 at the University of trieste
GEOGRAPHIC INFORMATION – NEED TO KNOW (GI-N2K) Towards a more demand-driven g...
Focussing Energy Consumers’ Behaviour Change towards Energy Efficiency and Lo...
Socio-Economic Planning profiles: Sciences VS Daily activities in public sector 
GEOGRAPHIC INFORMATION – NEED TO KNOW (GI-N2K) Towards a more demand-driven g...
Garden in motion. An experience of citizens involvement in public space regen...
Planning and Smartness: the true challenge
GeoSDI: una piattaforma social di dati geografici basata sui principi di INSP...
Murgante smart energy
Informazione Geografica, Città, Smartness
Tecnologie, Territorio, Smartness

Recently uploaded (20)

PDF
Web App vs Mobile App What Should You Build First.pdf
PDF
Unlocking AI with Model Context Protocol (MCP)
PDF
August Patch Tuesday
PDF
Transform Your ITIL® 4 & ITSM Strategy with AI in 2025.pdf
PDF
From MVP to Full-Scale Product A Startup’s Software Journey.pdf
PPTX
Chapter 5: Probability Theory and Statistics
PDF
A comparative analysis of optical character recognition models for extracting...
PDF
Video forgery: An extensive analysis of inter-and intra-frame manipulation al...
PPTX
TLE Review Electricity (Electricity).pptx
PPTX
A Presentation on Touch Screen Technology
PPTX
Digital-Transformation-Roadmap-for-Companies.pptx
PPTX
A Presentation on Artificial Intelligence
PDF
gpt5_lecture_notes_comprehensive_20250812015547.pdf
PDF
Zenith AI: Advanced Artificial Intelligence
PDF
Agricultural_Statistics_at_a_Glance_2022_0.pdf
PPTX
Tartificialntelligence_presentation.pptx
PPTX
SOPHOS-XG Firewall Administrator PPT.pptx
PDF
Approach and Philosophy of On baking technology
PDF
Microsoft Solutions Partner Drive Digital Transformation with D365.pdf
PDF
Mushroom cultivation and it's methods.pdf
Web App vs Mobile App What Should You Build First.pdf
Unlocking AI with Model Context Protocol (MCP)
August Patch Tuesday
Transform Your ITIL® 4 & ITSM Strategy with AI in 2025.pdf
From MVP to Full-Scale Product A Startup’s Software Journey.pdf
Chapter 5: Probability Theory and Statistics
A comparative analysis of optical character recognition models for extracting...
Video forgery: An extensive analysis of inter-and intra-frame manipulation al...
TLE Review Electricity (Electricity).pptx
A Presentation on Touch Screen Technology
Digital-Transformation-Roadmap-for-Companies.pptx
A Presentation on Artificial Intelligence
gpt5_lecture_notes_comprehensive_20250812015547.pdf
Zenith AI: Advanced Artificial Intelligence
Agricultural_Statistics_at_a_Glance_2022_0.pdf
Tartificialntelligence_presentation.pptx
SOPHOS-XG Firewall Administrator PPT.pptx
Approach and Philosophy of On baking technology
Microsoft Solutions Partner Drive Digital Transformation with D365.pdf
Mushroom cultivation and it's methods.pdf

Exploratory Spatial Analysis Norma

  • 1. Exploratory spatial analysis of illegal oil discharges detected off Canada’s Pacific Coast Norma Serra Rosaline Canessa Patrick O’Hara Stefania Bertazzon Marina Gavrilova ICCSA Conference 2008
  • 2. What are illegal vessel-source oil discharges? Photo courtesy of Gov. of Canada Photo courtesy of P. O’Hara Mystery Spill
  • 3. The ecological impact in marine ecosystems of small-scale chronic-levels of oil pollution is shown to be greater than large accidental oil spills; especially for seabird populations. ~ 80% global population of Cassin’s Auklets in BC ~ 74% global population of Ancient Murrelets ~ 56% global population Rhinoceros Auklet In British Columbia there is a rich diversity of seabirds at least 62 different species, including major proportions of world breeding populations: ~ 25% global population of Marbled Murrelets
  • 4. National Aerial Surveillance Program (NASP) Primary tool of Transport Canada for preventing and monitoring illegal oil pollution within Canada’s Economic Exclusive Zone, which generally extents 200 nautical miles (370 km) out from its coast . The Vancouver Sun. February 16 th , 2008
  • 5. Research Goal Specify a multivariate regression model to predict where illegal oil spills are more likely to occur, to provide NASP with a decision making tool to maximize the effectiveness of their surveillance program. In order to identify a valid model, in necessary to explore the spatial properties of our data: a) Spatial Heterogeneity b) Spatial Dependence or Spatial Autocorrelation (SAC)
  • 6. Study Area Maximum bounding area covered by surveillance flights off the west coast of Vancouver Island, British Columbia. 1
  • 7. Data sources From NASP : 53 detected oil spills in BC from 1997 to 2006 Surveillance flight path data in BC from 1997 to 2006. From Marine Communications Traffic Services - Vessel Traffic Operations Support System: Relative marine traffic densities in BC were estimated for different vessel type groups, from data collected in 2003.
  • 8. Marine shipping traffic in Canada’s West Coast Bulk carrier
  • 9. Marine shipping traffic in Canada’s West Coast Oil tanker
  • 10. Marine shipping traffic in Canada’s West Coast Tug vessel
  • 11. Marine shipping traffic in Canada’s West Coast Fishing vessels
  • 12. Marine shipping traffic in Canada’s West Coast Cruise ship
  • 13. Methods Global measures of spatial autocorrelation examine the nature and extent of the dependence within model variables and produce a single value for the entire data set. The two most common global measures: I ~ + 1  Positive SAC I ~ - 1  Negative SAC I = 0  No SAC Moran’s I Geary's c 0<C< 1  Positive SAC C > 1  Negative SAC C = 1  No SAC
  • 14. Methods (Cont.) Local measures of spatial association – Local Moran’s I identify the location and spatial scale of aggregations of unusual values, such as clusters of high values (hot spots) and low values (cold spots) and, because it generates an autocorrelation index for each data site, this can be mapped providing additional information about the pattern under study. E.g. LISA Cluster Map (GeoDA TM ) Outlier Outlier Cluster Cluster high values surrounded by low neighboring values High - Low low values surrounded by high neighboring values Low – High low values surrounded by neighbors of similar low values Low – Low high values surrounded by neighbors of similar high values High – High
  • 15. Results: Spatial autocorrelation – Global indices 0.871 (0.01) 0.158 (0.01) Oil spills counts 0.455 (0.01) 0.518 (0.01) Cruise ship movement counts 0.371 (0.01) 0.553 (0.01) Tug vessel movement counts 0.395 (0.01) 0.553 (0.01) Carrier vessel movement counts 0.369 (0.01) 0.589 (0.01) Oil tankers movement counts 0.249 (0.01) 0.783 (0.01) Fishing vessel movement counts 0.172 (0.01) 0.805 (0.01) NASP flight counts Geary's c ( p -values) Moran's I ( p -values)
  • 16. Results (Cont.): Spatial pattern description - LISA Cluster Maps LISA cluster map of surveillance flight counts. Significance test using 999 permutations and applying a significance filter of p ≤ 0.01.
  • 17. Results (Cont.): Spatial pattern description - LISA Cluster Maps LISA cluster map of detected oil spills counts. Significance test using 999 permutations and applying a significance filter of p ≤ 0.01.
  • 18. Results (Cont.): Spatial pattern description - LISA Cluster Maps LISA cluster map of carrier vessel movement counts. Significance test using 999 permutations and applying a significance filter of p ≤ 0.01.
  • 19. Results (Cont.): Spatial pattern description - LISA Cluster Maps LISA cluster map of oil tanker movement counts. Significance test using 999 permutations and applying a significance filter of p ≤ 0.01.
  • 20. Results (Cont.): Spatial pattern description - LISA Cluster Maps LISA cluster map of tug vessel movement counts. Significance test using 999 permutations and applying a significance filter of p ≤ 0.01.
  • 21. Results (Cont.): Spatial pattern description - LISA Cluster Maps LISA cluster map of fishing vessel movement counts. Significance test using 999 permutations and applying a significance filter of p ≤ 0.01.
  • 22. Results (Cont.): Spatial pattern description - LISA Cluster Maps LISA cluster map of cruise ships movement counts. Significance test using 999 permutations and applying a significance filter of p ≤ 0.01.
  • 23. Discussion Our exploratory analyses indicate that there is a positive spatial autocorrelation within datasets for all variables (dependent and independent). LISA cluster maps allowed the visualization of significant aggregations of high values of oil spill counts, surveillance effort and shipping traffic. To some extent LISA cluster maps show the effect of deterrence. For example, oil spill clusters were found in areas of Barkley Sound, whereas clusters of high counts of surveillance flights and marine vessels were mainly found in the entrance to the Strait of Juan de Fuca.
  • 24. Conclusion (I) Based on our results: Our variables exhibit spatial dependencies and non-stationarities, therefore, spatially autoregressive (SAR) and geographically weighted regression (GWR) models are the most suitable to predict spatial patterns of illegal oil discharges based on marine commercial vessel type and movement patterns, while controlling for the spatial distribution of surveillance flights.
  • 25. Conclusion (II) Based on our results: We propose to calibrate a multivariate regression model, combining a spatial autoregressive and a local specification. Our model will estimate oil spill probabilities based only on areas poorly covered by the current NASP flight zone, In order to provide effective and reliable recommendations for the planned surveillance program expansion.
  • 26. Explore the degree of correlation between detected oil spills and marine traffic densities for different vessel types (i.e. fishers, tugs, cargo vessels, tankers, cruise ships and ferries). Incorporate results from O’Hara et al. * that defines the relationship between oil spill detection probabilities and surveillance effort. Explore oil spill probabilities using a multivariate approach that will include variables such as season, distances between oil spills and shore and distances to nearest port or marina, and other characteristics that define the vessel type (i.e., flag state, inbound vs. outbound) Future research * O’Hara, P.D., Serra-Sogas, N., Canessa, R., Keller, P., Pelot, P.: Estimating oil spill rates and deterrence based on aerial surveillance data in Western Canadian marine waters. Marine Pollution Bulletin, (submitted) (2008).
  • 27. Thank you for your attention Any questions? Acknowledgements: