This dataset is part of the 3D hydrostratigraphic model of the Calgary-Lethbridge Corridor (CLC), southwestern Alberta. It represents the structure top of hydrostratigraphic unit (HSU) S1, which contains coarse-grained fluvial materialsl lying directly on bedrock and constrained to paleovalleys and some upland surfaces in the Neogene-Quaternary succession of the CLC.
This dataset supplements Alberta Geological Survey (AGS) Report 91, which includes a full description of steps taken to produce the structure grid.
Numerous water wells from the AWWID were used in the study to broadly define hydrostratigraphic units in the Neogene-Quaternary succession. The AWWID has abundant, relatively low quality data records, which include location inaccuracies resulting in up to 10 wells with duplicate locations. In cases of collocated data, where three or more wells with the same location information were encountered in the database, the deepest and shallowest well logs were selected to include wells that were completed in bedrock and to also characterize near surface sediments.The remaining duplicate wells were omitted.
The original input data were classified using unpublished hydrostratigraphic assignments to create a laterally extensive and mappable geobody. Classification of lithological data into hydrostratigraphic units may incorporate assumptions and human error when simplifying complex geological data, as unit geometry, extent, and description are a product of the author’s conceptual understanding. The HSU S1 grid is for use in regional-scale groundwater studies, as these mappable hydrostratigraphic units may be more complex or absent at the scale of a local investigation.
This grid represents all areas within the geographic boundary of the HSU S1.
All borehole data from the AWWID were sampled to the 25 m Provincial Digital Elevation Model (DEM) to provide a borehole top elevation (m asl). This process was necessary because the borehole top elevation of AWWID water wells, in many cases, was not available/recorded or verified, and if reported was estimated from a topographic map, field investigation or handheld Global Positional System (GPS) device. Elevation values of moderate- and high-quality data from AGS boreholes and field sites, geophysical logs, and previously published geological and hydrogeological maps were also checked with the DEM to verify elevation accuracy and modified to match the DEM, if necessary.
The vertical accuracy of the 25 m DEM is stated as 5 m @ 90% in the CLC derived from 1:60 000 aerial photographs using photogrammetric methods.
Methodology to create this grid dataset is fully documented in AGS Report 91.
Modelling Program used to create grid: ArcMap (grid interpolation) and RockWorks16 (deterministic alterations)
Number of data points available: 3595
Data type: A variety of low- to high-quality borehole data were used including lithologs submitted by water well contractors and compiled by the AWWID, downhole geophysical logs collected from water wells, and AGS borehole logs and field data.
Cell Size: 400 m
The Geostatistical Analyst ArcMap extension was used to interpolate the elevation value for HSU S1. A regional trend surface was modelled from the HSU S1 data points and removed prior to interpolation. An ordinary kriging algorithm was then used in the ArcMap Geostatistical Analyst extension to interpolate the residual elevation values and create a grid of the structure top of HSU S1.
The following Method Report and Prediction Errors provide information about the parameters used to create the grid, and any errors associated with the gridding process.
Method Report:
Method: Kriging
Type: Ordinary
Output type: Prediction
Variogram: Semivariogram
Number of lags: 12
Lag size: 569.03
Nugget: 3.28
Measurement error percentage: 100
Model type: spherical
Range: 6334.23
Anisotropy: Yes
Minor range: 3821.68
Direction: 102.31
Partial sill: 377.47
Prediction Errors:
Total points: 14793 (points include 0 m thickness for HSU S1, inserted at each geographic coordinate (X,Y locations) where the unit does not exist)
Mean: -0.07
Root mean square error: 9.53
Mean standardized: -0.00
Root mean square error, standardized: 1.05
Average standard error: 8.75
Error exclusion: 0
Number of points excluded: 0
Percent error exclusion: 0
Deterministic Grid Alteration Information:
All HSU gridded data in the CLC model (HSU S2, C, and S1) as well as the CLC bedrock topography were exported as ESRI ascii grids from ArcMap and imported into RockWorks16 for integration into the full hydrostratigraphic model.
Modelling constraints applied: The grid of the top of HSU S1 was constrained to the land surface (i.e., the 25 m Provincial Digital Elevation Model) and to the underlying CLC bedrock topography where the grids intersected.
Root mean square error (after deterministic grid alterations): 3.94 m
The RMSE compares how similar interpolated grid values are to the original input data. This is a global estimate of error for the HSU S1 grid and provides neither indication of how error is distributed spatially in the model domain nor the accuracy of the input data.