Geostatistics analyzes and predicts the values associated with spatial or spatial-temporal phenomena. It incorporates the spatial (and in some cases temporal) coordinates of the data within the analyses. It is a practical means of describing spatial patterns and interpolating values for locations where samples were not taken (and measures the uncertainty of those values, which is critical to informed decision making). This archive contains results of geostatistical analysis of COVID-19 case counts for all available US counties. Test results were obtained with ArcGIS Pro (ESRI). Sources are state health departments, which are scraped and aggregated by the Johns Hopkins Coronavirus Resource Center and then pre-processed by MappingSupport.com.

This update of the Zenodo dataset (version 6) consists of three compressed archives containing geostatistical analyses of SARS-CoV-2 testing data. This dataset utilizes many of the geostatistical techniques used in previous versions of this Zenodo archive, but has been significantly expanded to include analyses of up-to-date U.S. COVID-19 case data (from March 24th to September 8th, 2020):

Archive #1: “1.Geostat. Space-Time analysis of SARS-CoV-2 in the US (Mar24-Sept6).zip” – results of a geostatistical analysis of COVID-19 cases incorporating spatially-weighted hotspots that are conserved over one-week timespans. Results are reported starting from when U.S. COVID-19 case data first became available (March 24th, 2020) for 25 consecutive 1-week intervals (March 24th through to September 6th, 2020). Hotspots, where found, are reported in each individual state, rather than the entire continental United States.

Archive #2: "2.Geostat. Spatial analysis of SARS-CoV-2 in the US (Mar24-Sept8).zip" – the results from geostatistical spatial analyses only of corrected COVID-19 case data for the continental United States, spanning the period from March 24th through September 8th, 2020. The geostatistical techniques utilized in this archive includes ‘Hot Spot’ analysis and ‘Cluster and Outlier’ analysis.

Archive #3: "3.Kriging and Densification of SARS-CoV-2 in LA and MA.zip" – this dataset provides preliminary kriging and densification analysis of COVID-19 case data for certain dates within the U.S. states of Louisiana and Massachusetts.

These archives consist of map files (as both static images and as animations) and data files (including text files which contain the underlying data of said map files [where applicable]) which were generated when performing the following Geostatistical analyses: Hot Spot analysis (Getis-Ord Gi*) [‘Archive #1’: consecutive weeklong Space-Time Hot Spot analysis; ‘Archive #2’: daily Hot Spot Analysis], Cluster and Outlier analysis (Anselin Local Moran's I) [‘Archive #2’], Spatial Autocorrelation (Global Moran's I) [‘Archive #2’], and point-to-point comparisons with Kriging and Densification analysis [‘Archive #3’].

The Word document provided ("Description-of-Archive.Updated-Geostatistical-Analysis-of-SARS-CoV-2 (version 6).docx") details the contents of each file and folder within these three archives and gives general interpretations of these results.

The top of the Floridan Aquifer System was created by using FGS well core and cuttings data. The surface was interpolated by using the kriging method of the ArcGIS Geostatistical Analyst extension.

The top of IAS grid was created by using Florida Geological Survey well core and cuttings data. The surface was interpolated by using the kriging method from the ArcGIS Geostatistical Analyst package.

This raster in this tile service portrays the combined kriging prediction standard errors and depth to bedrock measurement uncertainties at 100-meter resolution. The kriging prediction standard errors for each grid cell are produced as output from the empirical Bayesian kriging process. The measurement uncertainties are assigned to the data by the authors based on who acquired the depth to bedrock data, the specific source of the data, or the method of data acquisition (drill hole or geophysical method).Map service is also available.See full metadata page.

This submission contains a shapefile of heat flow contour lines around the FORGE site located in Milford, Utah. The model was interpolated from data points in the Milford_wells shapefile. This heat flow model was interpolated from 66 data points using the kriging method in Geostatistical Analyst tool of ArcGIS. The resulting model was smoothed 100%. The well dataset contains 59 wells from various sources, with lat/long coordinates, temperature, quality, basement depth, and heat flow. This data was used to make models of the specific characteristics.

This is the Zenodo archive for the manuscript "Likely community transmission of COVID-19 infections between neighboring, persistent hotspots in Ontario, Canada" (Mucaki EJ, Shirley BC and Rogan PK. F1000Research 2021, 10:1312, DOI: 10.12688/f1000research.75891.1). This study aimed to produce community-level geo-spatial mapping of patterns and clusters of symptoms, and of confirmed COVID-19 cases, in near real-time in order to support decision-making. This was accomplished by area-to-area geostatistical analysis, space-time integration, and spatial interpolation of COVID-19 positive individuals. This archive will contain data and image files from this study, which were too numerous to be included in the manuscript for this study. It also provides all program files pertaining to the Geostatistical Epidemiology Toolbox (Geostatistical analysis software package to be used in ArcGIS), as well as all other scripts described in this manuscript and other software developed (cluster, outlier, streak identification and pairing)..

We also provide a guide which provides a general description of the contents of the four sections in this archive (Documentation_for_Sections_of_Zenodo_Archive.docx). If you have any intent to utilize the data provided in Section 3, we greatly advise you to review this document as it describes the output of all geostatistical analyses performed in this study in detail.

Data Files:

Section 1. "Section_1.Tables_S1_S7.Figures_S1_S11.zip"

This section contains all additional tables and figures described in the manuscript "Likely community transmission of COVID-19 infections between neighboring, persistent hotspots in Ontario, Canada". Additional tables S1 to S7 are presented in an Excel document. These 7 tables provide summary statistics of various geostatistical tests described in the study (“Section 1 – Tables S1-S4”) and lists all identified single and paired high-case cluster streaks (“Section 1 – Tables S5-S7”). This section also contains 11 additional figures referred to in the manuscript (“Section 1 – Figures S1-S11”) both individually and within a Word document which describes them.

Section 2. "Section_2.Localized_Hotspot_Lists.zip"

All localized hotspots (identified through kriging analysis) were catalogued for each municipality evaluated (Hamilton, Kitchener/Waterloo, London, Ottawa, Toronto, Windsor/Essex). These files indicate the FSA in which the hotspot was identified, the date in which it was identified (utilizing 3-day case data at the postal code level), the amount of cases which occurred within the FSA within these 3 dates, the range of cases interpolated by kriging analysis (between 5-10, 10-15, 15-20, 20-25, 25-30, 30-35, 35-40, 40-50, >50), and whether or not the FSA was deemed a hotspot by Gi* relative to the rest of Ontario on any of the three dates evaluated. Please see Section 4 for map images of these localized hotspots.

Section 3. "Section_3.All-Data_Files.Kriging_GiStar_Local_and_GlobalMorans.2020_2021"

Section 3 – All output files from the geostatistical tests performed in this study are provided in this section. This includes the output from Ontario-wide FSA-level Gi* and Cluster and Outlier analyses, and PC-level Cluster and Outlier, Spatial Autocorrelation, and kriging analysis of 6 municipal regions. It also includes kriging analysis of 7 other municipal regions adjacent to Toronto (Ajax, Brampton, Markham, Mississauga, Pickering, Richmond Hill and Vaughan). This section also provides data files from our analyses of stratified case data (by age, gender, and at-risk condition). All coordinates presented in these data files are given in “PCS_Lambert_Conformal_Conic” format. Case values between 1-5 were masked (appear as “NA”).

Section 4. "Section_4.All_Map_Images_of_Geostat_Analyses.zip"

Sets of image files which map the results of our geostatistical analyses onto a map of Ontario or within the municipalities evaluated (Hamilton, Kitchener/Waterloo, London, Ottawa, Toronto, Windsor/Essex) are provided. This includes: Kriging analysis (PC-level), Local Moran's I cluster and outlier analysis (FSA and PC-level), normal and space-time Gi* analysis, and all images for all analyses performed on stratified data (by age, gender and at-risk condition). Kriging contour maps are also included for 7 other municipal regions adjacent to Toronto (Ajax, Brampton, Markham, Mississauga, Pickering, Richmond Hill and Vaughan).

Software:

This Zenodo archive also provides all program files pertaining to the Geostatistical Epidemiology Toolbox (Geostatistical analysis software package to be used in ArcGIS), as well as all other scripts described in this manuscript. This geostatistical toolbox was developed by CytoGnomix Inc., London ON, Canada and is distributed freely under the terms of the GNU General Public License v3.0. It can be easily modified to accommodate other Canadian provinces and, with some additional effort, other countries.

This distribution of the Geostatistical Epidemiology Toolbox does not include postal code (PC) boundary files (which are required for some of the tools included in the toolbox). The PC boundary shapefiles used to test the toolbox were obtained from DMTI (https://www.dmtispatial.com/canmap/) through the Scholar's Geoportal at the University of Western Ontario (http://geo2.scholarsportal.info/). The distribution of these files (through sharing, sale, donation, transfer, or exchange) is strictly prohibited. However, any equivalent PC boundary shape file should suffice, provided it contains polygon boundaries representing postal code regions (see guide for more details).

Software File 1. "Software.GeostatisticalEpidemiologyToolbox.zip"

The Geostatistical Epidemiology Toolbox is a set of custom Python-based geoprocessing tools which function as any built-in tool in the ArcGIS system. This toolbox implements data preprocessing, geostatistical analysis and post-processing software developed to evaluate the distribution and progression of COVID-19 cases in Canada. The purpose of developing this toolbox is to allow external users without programming knowledge to utilize the software scripts which generated our analyses and was intended to be used to evaluate Canadian datasets. While the toolbox was developed for evaluating the distribution of COVID-19, it could be utilized for other purposes.

The toolbox was developed to evaluate statistically significant distributions of COVID-19 case data at Canadian Forward Sortation Area (FSA) and Postal Code-level in the province of Ontario utilizing geostatistical tools available through the ArcGIS system. These tools include: 1) Standard Gi* analysis (finds areas where cases are significantly spatially clustered), 2) spacetime based Gi* analysis (finds areas where cases are both spatially and temporally clustered), 3) cluster and outlier analysis (determines if high case regions are an regional outlier or part of a case cluster), 4) spatial autocorrelation (determines the cases in a region are clustered overall) and, 5) Empirical Bayesian Kriging analysis (creates contour maps which define the interpolation of COVID-19 cases in measured and unmeasured areas). Post-processing tools are included that import these all of the preceding results into the ArcGIS system and automatically generate PNG images.

This archive also includes a guide ("UserManual_GeostatisticalEpidemiologyToolbox_CytoGnomix.pdf") which describes in detail how to set up the toolbox, how to format input case data, and how to use each tool (describing both the relevant input parameters and the structure of the resultant output files).

Software File 2: “Software.Additional_Programs_for_Cluster_Outlier_Streak_Idendification_and_Pairing.zip"

In the manuscript associated with this archive, Perl scripts were utilized to evaluate postal code-level Cluster and Outlier analysis to identify significantly, highly clustered postal codes over consecutive periods (i.e., high-case cluster “streaks”). The identified streaks are then paired to those in close proximity, based on the neighbors of each postal code from PC centroid data ("paired streaks"). Multinomial logistic regression models were then derived in the R programming language to measure the correlation between the number of cases reported in each paired streak, the interval of time separating each streak, and the physical distance between the two postal codes. Here, we provide the 3 Perl scripts and the R markdown file which perform these tasks:

“Ontario_City_Closest_Postal_Code_Identification.pl”

Using an input file with postal code coordinates (by centroid), this program identifies the nearest neighbors to all postal codes for a given municipal region (the name of this region is entered on the command line). Postal code centroids were calculated in ArcGIS using the “Calculate Geometry” function against DMTI postal code boundary files (not provided). Input from other sources could be used, however, as long as the input includes a list of coordinates with a unique label associated with a particular municipality.

The output of this program (for the same municipal region being evaluated) is required for the following two Perl

This is the data repository for the PLOS ONE Manuscript: "Meeting radiation dosimetry capacity requirements of population-scale exposures by geostatistical sampling". This repository contains the following data:

1. "State-and-Subdivision-Boundary-Files.Edited-for-ArcMap-10.4.KML-Format.zip":

This file contains modified U.S. state and sub-division boundary files [in KML format], which can be imported into ArcMap using its KMLtoLayer function. These files have been modified to prevent sub-division naming issues that we encountered when importing boundary data into ArcMap: A) State sub-divisions with identical names are considered a single sub-division by ArcMap (corrected by adding a letter after each sub-division of the same name, i.e. CenterA, CenterB, etc), and; B) ArcMap would only identify the sub-division by its first word if sub-division name contained spaces (corrected by converting all spaces into dashes).

2. "HPAC-Plumes.Processed.zip" and "HPAC-Plumes.Unprocessed.zip":

These files contain HPAC plume coordinate (WGS1984) and dose (in cGy) values for all scenarios discussed in the manuscript. We provide "processed" and "unprocessed" HPAC plume data files. The "unprocessed" HPAC plume data is provided in its original XML format, which cannot be imported into ArcMap directly. The "processed" HPAC plumes are provided in tab-delimited X,Y,Z format (Latitude, Longitude, and Dose). We have also added a "0 cGy" contour in the "processed" plumes (surrounding the HPAC plume), as the presence of unirradiated data points adjacent to the plume was found to be crucial for accurate kriging, since these points served as boundaries for kriging.

3. "Final-Derived-Plumes.Data-Points.zip":

This file contains geostatistically-derived plume coordinate (WGS1984) and dose (in cGy) values for all scenarios discussed in the manuscript. Data is in comma-delimited format (Latitude, Longitude, and Dose). Data points consist of a set of initial coordinates generated at random locations within each Census sub-division using the ArcMap tool, ‘CreateRandomPoints_management’, and subsequent points generated by densification (the geostatistical procedure that targets and localizes an additional small cohort of irradiated individuals to mitigate uncertainty in environmental measurements). These data points were assigned radiation level values corresponding to the adjacent outer HPAC contour by a script comparing each sample with its location within the HPAC plume of the same scenario.

4. “Intermediate-Derived-Plumes.Data-Points.zip”

This archive contains coordinate data (WGS1984) and dose values (in cGy) for all intermediary steps of plume development (using our geostatistical method) for all scenarios. Like (3), the data is comma-delimited (Latitude, Longitude, and Dose), and were assigned radiation level values by a script comparing sampling locations with the location of the HPAC plume of the same scenario. Scenario replicate folders contains text files for each iteration step of the plume derivation process, including a file containing just the initial random sampling (“Iteration-1”), a file containing initial sampling and sampling locations selected by the first densification step (“Iteration-2”), a file containing initial sampling and sampling locations selected by the first and second densification steps (“Iteration-3”), and so on.

This archive also contains a Table (“Progression-of-New-Densification-Selected-Sampling-Locations-For-All-Scenarios.xslx”) which provides a categorical breakdown of how many unique densification-selected sampling locations occur within the irradiated region (i.e. overlap the HPAC plume) for each iteration of all scenario replicates. The fraction of irradiated to unirradiated sampling locations varies among each scenario and individual replicates for the same scenario. Our analysis shows that these results depend on the population densities and exact topography of the HPAC plume which is different among each scenario.

5. "Geostatistical-Sampling-Project.All-Scripts.zip"

This archive contains all programs required for this project. This includes Python scripts meant to be run within the ArcMap software environment (for random point generation and data extraction), and Perl scripts used to process HPAC and U.S. State and Sub-division boundary files, and to assign radiation values to sample locations based on a modified HPAC plume. A java program, “CompareReplicates.jar”, compares the overlapping areas between a pair of polygons that overlap one other using the ArcMap software environment, and requires access to the ArcGIS Runtime SDK (https://developers.arcgis.com/arcgis-runtime/).

This data set archives all inputs, outputs and scripts needed to reproduce the findings of W.H. Farmer and G.F. Koltun in the 2017 Journal of Hydrology Regional Studies article entitled “Geospatial Tools Effectively Estimated Nonexceedance Probabilities of Daily Streamflow at Ungauged and Intermittently Gauged Locations in Ohio”. Input data includes observed streamflow values, in cubic feet per second, for 152 streamgages in and around Ohio from 01 January 2009 through 31 August 2015. Data from the Ohio Environmental Protection Agency on where and when water quality samples were taken are also provided. Geospatial locations are provided for all streamgages and sampling sites considered. ESRI ArcGIS shapefiles are available for all maps produced in the original publication. Comma-separated-value files contain the output data required to reproduce every figure in the report. This archive also includes an R script capable of reading the input files and producing output files and figures. See the README.txt file for a full description of model application.

Michigan DNR Fisheries Division hydroacoustic survey data (2003-2011) were collected annually at survey stations in Lake St. Claire to document the abundance and distribution of submerged aquatic vegetation. The hydrouacoustic surveys at each station consisted of 11 parallel east-west boat transects spaced 10 m apart within a square hectare plot. Not all survey stations were sampled every year; to allow for sufficient data for mapping, sonar records across multiple years were compiled and averaged on a per station basis. Plant height, percent of lake bottom covered by aquatic plants and percent of water column containing plant biomass interpolated layers were generated from this data using ordinary kriging and the ArcGIS Geostatistical Analyst. The MDNR Fisheries Research Report 2099, documenting hydroacoustic survey methods and results from the years 2003-2007. Kriging parameters are as follows: SAV_mean_per_cover_2008_2011; no transformations, no trend removal; Gaussian semivariogram model, nugget = 999, major range = 9063, partial sill = 759, number of lags = 12, lag size = 1133; smooth search neighborhood, smoothing factor = 0.1

This dataset contains base-flow recession time constant (tau) contours that are interpreted from tau values calculated at streamgages in the Niobrara National Scenic River study area. The contours were created by interpolating the calculated tau values using geostatistical kriging methods. Kriging is a geostatistical method that can be used to determine optimal weights for measurements at sampled locations (streamgages) for the estimation of values at unsampled locations (ungaged sites). The kriged tau map could be used (1) as the basis for identifying areas with different hydrologic responsiveness, and (2) in the development of regional low-flow regression equations. The Geostatistical Analyst tools in ArcGIS Pro version 2.5.2 (Environmental Systems Research Institute, 2012) were used to create the kriged tau map and perform cross validation to determine the root mean square error (RMSE) of the tau map. The results of the kriging were then exported from ArcGIS to contours.

The Intermediate Aquifer System (IAS) consists of highly-variable, fine-grained, siliciclastic deposits and low permeability limestone and dolostone. It refers to a statewide hydrostratigraphic unit which provides variable confinement to the Floridan Aquifer system. The IAS grid was created by using Florida Geological Survey well core and cuttings data. The surface was interpolated by using the kriging method from the ArcGIS Geostatistical Analyst package.

Annual and five-year (5YA) average wet deposition maps for the nitrate ion are available. The file formats include geodatabase files (*.gdb) compatible with geospatial software (e.g. ESRI ArcGIS) and KMZ files compatible with virtual globe software (e.g. Google Earth™). Maps can also be viewed online via Open Maps and the ArcGIS online viewer. Annual deposition from each site was screened for completeness using the following criteria: (1) precipitation amounts were recorded for >90% of the year and >60% of each quarter, and (2) nitrate concentrations were reported for >70% of the precipitation measured over the year and for >60% of each quarter. Five-year average wet deposition values are averaged annual deposition values with a completeness criterion >60% for the five-year period. Units for wet deposition fluxes are in kg of NO3 per hectare per year (kg ha-1 y-1). Sources of measurement data and spatial interpolation method are described here: https://doi.org/10.18164/e8896575-1fb8-4e53-8acd-8579c3c055c2. Recommended citation: Environment and Climate Change Canada, [year published]. NO3 Wet Deposition Maps. Air Quality Research Division, Environment and Climate Change Canada, Toronto, Ontario, Canada. [URL/DOI], accessed [date]. Recommended acknowledgement: The author(s) acknowledge Environment and Climate Change Canada for the provision of Canada-U.S. wet deposition kriging maps accessed from the Government of Canada Open Government Portal at open.canada.ca, and the data providers referenced therein.

This 1m Digital Surface Model (DSM) shaded relief is derived from first-stop Light Detection and Ranging (LiDAR) point cloud data from September 2005 for the Green Lakes Valley, near Boulder Colorado. The DSM was created from LiDAR point cloud tiles subsampled to 1-meter postings, acquired by the National Center for Airborne Laser Mapping (NCALM) project. This data was collected in collaboration between the University of Colorado, Institute of Arctic and Alpine Research (INSTAAR) and NCALM, which is funded by the National Science Foundation (NSF). The DSM shaded relief has the functionality of a map layer for use in Geographic Information Systems (GIS) or remote sensing software. Total area imaged is 35 km^2. The LiDAR point cloud data was acquired with an Optech 1233 Airborne Laser Terrain Mapper (ALTM) and mounted in a twin engine Piper Chieftain (N931SA) with Inertial Measurement Unit (IMU) at a flying height of 600 m. Data from two GPS (Global Positioning System) ground stations were used for aircraft trajectory determination. The continuous DSM surface was created by mosaicing and then kriging 1 km2 LiDAR point cloud LAS-formated tiles using Golden Software's Surfer 8 Kriging algorithm. Horizontal accuracy and vertical accuracy is unknown. cm RMSE at 1 sigma. The layer is available in GEOTIF format approx. 265 MB of data. It has a UTM zone 13 projection, with a NAD83 horizonal datum and a NAVD88 vertical datum computed using NGS GEOID03 model, with FGDC-compliant metadata. This shaded relief model was also generated. A similar layer, the Digital Terrain Model (DTM), is a ground-surface elevation dataset better suited for derived layers such as slope angle, aspect, and contours. A processing report and readme file are included with this data release. The DSM dataset is available through an unrestricted public license. The LiDAR DEMs will be of interest to land managers, scientists, and others for study of topography, ecosystems, and environmental change. NOTE: This EML metadata file does not contain important geospatial data processing information. Before using any NWT LTER geospatial data read the arcgis metadata XML file in either ISO or FGDC compliant format, using ArcGIS software (ArcCatalog > description), or by viewing the .xml file provided with the geospatial dataset.

Annual and five-year (5YA) average wet deposition maps for the non-sea-salt sulfate ion are available. The file formats include geodatabase files (*.gdb) compatible with geospatial software (e.g. ESRI ArcGIS) and KMZ files compatible with virtual globe software (e.g. Google Earth™). Maps can also be viewed online via Open Maps and the ArcGIS online viewer. Annual deposition from each site was screened for completeness using the following criteria: (1) precipitation amounts were recorded for >90% of the year and >60% of each quarter, and (2) sulfate concentrations were reported for >70% of the precipitation measured over the year and for >60% of each quarter. Five-year average wet deposition values are averaged annual deposition values with a completeness criterion >60% for the five-year period. Units for wet deposition fluxes are in kg of xSO4 per hectare per year (kg ha-1 y-1). Sources of measurement data and spatial interpolation method are described here: https://doi.org/10.18164/e8896575-1fb8-4e53-8acd-8579c3c055c2. Recommended citation: Environment and Climate Change Canada, [year published]. xSO4 Wet Deposition Maps. Air Quality Research Division, Environment and Climate Change Canada, Toronto, Ontario, Canada. [URL/DOI], accessed [date]. Recommended acknowledgement: The author(s) acknowledge Environment and Climate Change Canada for the provision of Canada-U.S. wet deposition kriging maps accessed from the Government of Canada Open Government Portal at open.canada.ca, and the data providers referenced therein.

Processed results from of surface x-ray florescence (XRF) analysis of the sediment grab samples recovered as part of the Long Island Sound mapping project Phase II. Sediment grab samples have been collected in the summers of 2017 and 2018 using a modified van Veen grab sampler. A subsample of the top two centimeter was taken for further lab analysis. Dried and homogenized splits of the samples were analyzed for chemical composition using an Innov-X Alpha series 4000 XRF (Innov-X Systems, Woburn, MA). The results of the measurements are presented as ppm. The XRF analytical protocol included the following elements: P, S, Cl, K, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn As, Se, Br, Rb, Sr, Zr, Mo, Ag, Cd, Sn, Sb, I, Ba, Hg, Pb, Bi, Th, and U. However, only Cl, K, Ca, Ti, Cr, Mn, Fe, Co, Cu, Zn, As, Br, Rb, Sr, Zr and Pb were consistently present at levels above background detection in surficial sediments collected in the LIS Phase II area. ArcGIS Pro empirical kriging has been used to interpolate values for selected elements onto a 50 m raster. The interpolated raster has been clipped to fit the extent of the phase 2 survey area. The final raster data are in GeoTiff format with UTM 18 N projection.Time period of content: 2017-08-01 to 2022-11-16Attribute accuracy: The attribute accuracy has not been determined. This raster dataset shown mainly the major trends and patterns of the value distribution in the Phase 2 study area.Completeness: The dataset is complete.Positional accuracy: The raster resolution is 50 m.Attributes: XRF Zn raster: Interpolated Zinc (Zn) distribution (in ppm)

This 1m Digital Terrain Model (DTM) is derived from bare-ground Light Detection and Ranging (LiDAR) point cloud data from September 2005 for the Green Lakes Valley, near Boulder Colorado. This dataset is better suited for derived layers such as slope angle, aspect, and contours. The DTM was created from LiDAR point cloud tiles subsampled to 1-meter postings, acquired by the National Center for Airborne Laser Mapping (NCALM) project. This data was collected in collaboration between the University of Colorado, Institute of Arctic and Alpine Research (INSTAAR) and NCALM, which is funded by the National Science Foundation (NSF). The DTM has the functionality of a map layer for use in Geographic Information Systems (GIS) or remote sensing software. Total area imaged is 35 km^2. The LiDAR point cloud data was acquired with an Optech 1233 Airborne Laser Terrain Mapper (ALTM) and mounted in a twin engine Piper Chieftain (N931SA) with Inertial Measurement Unit (IMU) at a flying height of 600 m. Data from two GPS (Global Positioning System) ground stations were used for aircraft trajectory determination. The continuous DTM surface was created by mosaicing and then kriging 1 km2 LiDAR point cloud LAS-formated tiles using Golden Software's Surfer 8 Kriging algorithm. Horizontal accuracy and vertical accuracy is unknown. The layer is available in GEOTIF format approx. 265 MB of data. It has a UTM zone 13 projection, with a NAD83 horizonal datum and a NAVD88 vertical datum computed using NGS GEOID03 model, with FGDC-compliant metadata. A shaded relief model was also generated. A similar layer, the Digital Surface Model (DSM), is a first-stop elevation layer. A processing report and readme file are included with this data release. The DTM is available through an unrestricted public license. The LiDAR DEMs will be of interest to land managers, scientists, and others for study of topography, ecosystems, and environmental change. NOTE: This EML metadata file does not contain important geospatial data processing information. Before using any NWT LTER geospatial data read the arcgis metadata XML file in either ISO or FGDC compliant format, using ArcGIS software (ArcCatalog > description), or by viewing the .xml file provided with the geospatial dataset.

This 1m Digital Surface Model (DSM) is derived from first-stop Light Detection and Ranging (LiDAR) point cloud data from September 2005 for the Green Lakes Valley, near Boulder Colorado. The DSM was created from LiDAR point cloud tiles subsampled to 1-meter postings, acquired by the National Center for Airborne Laser Mapping (NCALM) project. This data was collected in collaboration between the University of Colorado, Institute of Arctic and Alpine Research (INSTAAR) and NCALM, which is funded by the National Science Foundation (NSF). The DSM has the functionality of a map layer for use in Geographic Information Systems (GIS) or remote sensing software. Total area imaged is 35 km^2. The LiDAR point cloud data was acquired with an Optech 1233 Airborne Laser Terrain Mapper (ALTM) and mounted in a twin engine Piper Chieftain (N931SA) with Inertial Measurement Unit (IMU) at a flying height of 600 m. Data from two GPS (Global Positioning System) ground stations were used for aircraft trajectory determination. The continuous DSM surface was created by mosaicing and then kriging 1 km2 LiDAR point cloud LAS-formated tiles using Golden Software's Surfer 8 Kriging algorithm. Horizontal accuracy and vertical accuracy is unknown. cm RMSE at 1 sigma. The layer is available in GEOTIF format approx. 265 MB of data. It has a UTM zone 13 projection, with a NAD83 horizonal datum and a NAVD88 vertical datum computed using NGS GEOID03 model, with FGDC-compliant metadata. A shaded relief model was also generated. A similar layer, the Digital Terrain Model (DTM), is a ground-surface elevation dataset better suited for derived layers such as slope angle, aspect, and contours. A processing report and readme file are included with this data release. The DSM is available through an unrestricted public license. The LiDAR DEMs will be of interest to land managers, scientists, and others for study of topography, ecosystems, and environmental change. NOTE: This EML metadata file does not contain important geospatial data processing information. Before using any NWT LTER geospatial data read the arcgis metadata XML file in either ISO or FGDC compliant format, using ArcGIS software (ArcCatalog > description), or by viewing the .xml file provided with the geospatial dataset.

Geoscience Australia has been deriving raster sediment datasets for the continental Australian Exclusive Economic Zone (AEEZ) using existing marine samples collected by Geoscience Australia and external organisations. Since seabed sediment data are collected at sparsely and unevenly distributed locations, spatial interpolation methods become essential tools for generating spatially continuous information. Previous studies have examined a number of factors that affect the performance of spatial interpolation methods. These factors include sample density, data variation, sampling design, spatial distribution of samples, data quality, correlation of primary and secondary variables, and interaction among some of these factors. Apart from these factors, a spatial reference system used to define sample locations is potentially another factor and is worth investigating. In this study, we aim to examine the degree to which spatial reference systems can affect the predictive accuracy of spatial interpolation methods in predicting marine environmental variables in the continental AEEZ. Firstly, we reviewed spatial reference systems including geographic coordinate systems and projected coordinate systems/map projections, with particular attention paid to map projection classification, distortion and selection schemes; secondly, we selected eight systems that are suitable for the spatial prediction of marine environmental data in the continental AEEZ. These systems include two geographic coordinate systems (WGS84 and GDA94) and six map projections (Lambert Equal-area Azimuthal, Equidistant Azimuthal, Stereographic Conformal Azimuthal, Albers Equal-Area Conic, Equidistant Conic and Lambert Conformal Conic); thirdly, we applied two most commonly used spatial interpolation methods, i.e. inverse distance squared (IDS) and ordinary kriging (OK) to a marine dataset projected using the eight systems. The accuracy of the methods was assessed using leave-one-out cross validation in terms of their predictive errors and, visualization of prediction maps. The difference in the predictive errors between WGS84 and the map projections were compared using paired Mann-Whitney test for both IDW and OK. The data manipulation and modelling work were implemented in ArcGIS and R. The result from this study confirms that the little shift caused by the tectonic movement between WGS84 and GDA94 does not affect the accuracy of the spatial interpolation methods examined (IDS and OK). With respect to whether the unit difference in geographical coordinates or distortions introduced by map projections has more effect on the performance of the spatial interpolation methods, the result shows that the accuracies of the spatial interpolation methods in predicting seabed sediment data in the SW region of AEEZ are similar and the differences are considered negligible, both in terms of predictive errors and prediction map visualisations. Among the six map projections, the slightly better prediction performance from Lambert Equal-Area Azimuthal and Equidistant Azimuthal projections for both IDS and OK indicates that Equal-Area and Equidistant projections with Azimuthal surfaces are more suitable than other projections for spatial predictions of seabed sediment data in the SW region of AEEZ. The outcomes of this study have significant implications for spatial predictions in environmental science. Future spatial prediction work using a data density greater than that in this study may use data based on WGS84 directly and may not have to project the data using certain spatial reference systems. The findings are applicable to spatial predictions of both marine and terrestrial environmental variables.

You can also purchase hard copies of Geoscience Australia data and other products at http://www.ga.gov.au/products-services/how-to-order-products/sales-centre.html

This study subdivides the Weddell Sea, Antarctica, into seafloor regions using multivariate statistical methods. These regions are categories used for comparing, contrasting and quantifying biogeochemical processes and biodiversity between ocean regions geographically but also regions under development within the scope of global change. The division obtained is characterized by the dominating components and interpreted in terms of ruling environmental conditions. The analysis uses 28 environmental variables for the sea surface, 25 variables for the seabed and 9 variables for the analysis between surface and bottom variables. The data were taken during the years 1983-2013. Some data were interpolated. The statistical errors of several interpolation methods (e.g. IDW, Indicator, Ordinary and Co-Kriging) with changing settings have been compared for the identification of the most reasonable method. The multivariate mathematical procedures used are regionalized classification via k means cluster analysis, canonical-correlation analysis and multidimensional scaling. Canonical-correlation analysis identifies the influencing factors in the different parts of the cove. Several methods for the identification of the optimum number of clusters have been tested. For the seabed 8 and 12 clusters were identified as reasonable numbers for clustering the Weddell Sea. For the sea surface the numbers 8 and 13 and for the top/bottom analysis 8 and 3 were identified, respectively. Additionally, the results of 20 clusters are presented for the three alternatives offering the first small scale environmental regionalization of the Weddell Sea. Especially the results of 12 clusters identify marine-influenced regions which can be clearly separated from those determined by the geological catchment area and the ones dominated by river discharge. Supplement to: Jerosch, Kerstin; Pehlke, Hendrik; Weber, Lukas; Teschke, Katharina; Heidemann, Teresa; Scharf, Frauke Katharina (in prep.): Comparing the surface and the bottom of the Southern Ocean using multivariate cluster analysis: regional effects of environmental parameters.

Climate variable. Map with the average annual mean temperature for the period 1991-2019 obtained by interpolation (kriging with external drift) from data observed at manual weather stations.