Professor Corey J. A. Bradshaw
Global Ecology | Partuyarta Ngadluku Wardli Kuu, Flinders University
e-mail

Team:

• Associate Professor Larissa Schneider, Australian National University
• Adjunct Professor Patrice de Caritat, Curtin University
• Professor Simon Haberle, Australian National University
• James Taylor, Australian National University
• Dr Olha Furman, Department of Climate Change, Energy, the Environment and Water, Commonwealth of Australia

• identify external and indirect determinants of mercury (Hg)
• understand environmental conditions that influence mercury retention and mobility
• predict continental distribution of soil mercury

All scripts and data support the following paper:

Schneider, L, P de Caritat, JR Taylor, OS Furman, SG Haberle, CJA Bradshaw. 2026. Predicting continental-scale soil mercury concentrations in Australia to refine global frameworks. Environmental Science and Technology doi:10.1021/acs.est.5c11189

• HgGH.R: all required R code combined

• geochem.csv: geochemical data
• field.csv: sample point characteristics
• hgTSID.csv: re-analysed [Hg] estimates (ng/g)
• gs.csv: grain-size category percentages

Most of these files are too large to store in this repository directly, so in most cases the links refer to the original repository URLs where you can download the datasets.

• aus.shp: Australia boundary shapefile (zipped)

• NLUM_v7_250_ALUMV8_2020_21_alb.tif: land use of Australia 2010–11 to 2020–21 (download geotif raster from original site)

• wwf_terr_ecos.shp: WWF ecoregions shapefile (download zipped file from original site)
• OzWALD.GPP.AnnualMeans.nc: vegetation carbon uptake (gross primary production) NetCDF (download from original site)
• OzWALD.LAI.AnnualMeans.nc: leaf area index NetCDF (download from original site)

• GeologicUnitPolygons1M.shp: 1:1,000,000 geological unit polygon shapefile (download zipped file from original site)
• lithreclass.csv: reclassified lithology groups text file
• radmap_v4_2019_filtered_ML_KThU_RGB_24bit.tif: potassium:thorium:uranium geotif raster (download data package from original site)
• radmap_v4_2019_filtered_ML_ppmTh_32bitfloat_grid.tif: thorium ppm geotif raster (download data package from original site)
• radmap_v4_2019_filtered_ML_ppmU_32bitfloat_grid.tif: uranium ppm geotif raster (download data package from original site)
• radmap_v4_2019_filtered_ML_pctk_32bitfloat_grid.tif: % potassium geotif raster (download data package from original site)
• NTO_000_005_EV_N_P_AU_NAT_C_20231101.tif: soil nitrogen (0-5 cm) geotif raster (download from original site)
• PTO_000_005_EV_N_P_AU_NAT_C_20231101.tif: soil phosphorus (0-5 cm) geotif raster (download from original site)
• pHc_000_005_EV_N_P_AU_NAT_C_20140801.tif: soil pH (0-5 cm) geotif raster (download from original site)
• CLY_000_005_EV_N_P_AU_TRN_N_20210902.tif: % soil clay content geotif raster (download from original site)
• SLT_000_005_EV_N_P_AU_TRN_N_20210902.tif: % soil silt content geotif raster (download from original site)
• ferric2rsmp.rds: principal component 2 of blue, red, NIR, SWIR1 from the enhanced barest Earth (proxy for soil categories on continuous scale) (based on Loughin 1991); layer pre-prepared for appropriate projection and clipped area. Due to Github file-size constraints, this is the resampled .rds file (ferric2.rsmp) indicated in the script. We have broken the file into similar-sized chunks, and then compressed them using the fast, lossless compression algorithm zstd. First, decompress each .zst chunk using the following command in Terminal (or equivalent): zstd -d 'ferric2rsmp_chunk_**.zst', and then combine chunks a to h using the following Terminal (or equivalent) command: cat ferric2rsmp_chunk_* > ferric2rsmp.rds (individual .zst files are in .../data/spatial/barest earth/ferricPC2/ in this repository). Once chunks are recombined, import the .rds file in R with this code: ferric2.rsmp <- readRDS(ferric2rsmp.rds)
• ferric4rsmp.rds: principal component 4 of blue, red, NIR, SWIR1 from the enhanced barest Earth (proxy for soil categories on continuous scale) (based on Loughin 1991); layer pre-prepared for appropriate projection and clipped area. Due to Github file-size constraints, this is the resampled .rds file (ferric4.rsmp) indicated in the script. We have broken the file into similar-sized chunks, and then compressed them using the fast, lossless compression algorithm zstd. First, decompress each .zst chunk using the following command in Terminal (or equivalent): zstd -d 'ferric4rsmp_chunk_**.zst', and then combine chunks a to h using the following Terminal (or equivalent) command: cat ferric4rsmp_chunk_* > ferric4rsmp.rds (individual .zst files are in .../data/spatial/barest earth/ferricPC4/ in this repository). Once chunks are recombined, import the .rds file in R with this code: ferric4.rsmp <- readRDS(ferric4rsmp.rds)
• Aluminium_oxide_prediction_median.tif: aluminium oxide (Al₂O₃)
• Feox9473.tif: iron oxide (Fe₂O₃); pre-prepared for correct projection. Due to Github file-size constraints, we have broken the file into similar-sized chunks, and then compressed them using the fast, lossless compression algorithm zstd. First, decompress each .zst chunk using the following command in Terminal (or equivalent): zstd -d 'Feox9473_chunk_**.zst', and then combine chunks aa to jg using the following Terminal (or equivalent) command: cat Feox9473_chunk_* > Feox9473.tif (individual .zst files are in .../data/spatial/oxides/Fe/part1 ... part19 in this repository).
• Pox9473.tif: phosphorus oxide (P₂O₅); pre-prepared for correct projection. Due to Github file-size constraints, we have broken the file into similar-sized chunks, and then compressed them using the fast, lossless compression algorithm zstd. First, decompress each .zst chunk using the following command in Terminal (or equivalent): zstd -d 'Pox9473_chunk_**.zst', and then combine chunks aa to jg using the following Terminal (or equivalent) command: cat Pox9473_chunk_* > Pox9473.tif (individual .zst files are in .../data/spatial/oxides/P/part1 ... part19 in this repository).

• OzWALD.annual.Pg.AnnualSums.nc: annual precipitation NetCDF (download from original site)
• PrescottIndex_01_3s_lzw.tif: Prescott index geotif raster (download from original site)
• OzWALD.Ssoil.AnnualMeans.nc: soil water availability NetCDF (download from original site)

• XXXHgPredSpatRFlog10.nc: This is the NetCDF file for the top (XXX = TOS) and bottom (XXX = BOS) outlet sediment sample output maps of Australia-wide Hg concentration (log₁₀ values) predicted from the random forest model (resolution = 0.005° × 0.005° latitude/longitude ≈ 0.55 km × 0.55 km ≈ 0.306 km²). Due to Github file-size constraints, we have broken the file into similar-sized chunks, and then compressed them using the fast, lossless compression algorithm zstd. First, decompress each .zst chunk using the following command in Terminal (or equivalent): zstd -d 'XXXHgPredSpatRFlog10.nc_chunk_a*.zst', and then combine chunks a to n using the following Terminal (or equivalent) command: cat XXXHgPredSpatRFlog10.nc_chunk_* > XXXHgPredSpatRFlog10.nc. You can import the NetCDF file in R using the ncdf4 package and its function nc_open. This produces a ncdf4 object that can be converted to a SpatRaster object using the rast function in package terra.
• XXXHgPredSpatRFlin.nc: This is the NetCDF file for the top (XXX = TOS) and bottom (XXX = BOS) outlet sediment sample output maps of Australia-wide Hg concentration (back-transformed to linear values) predicted from the random forest model (resolution = 0.005° × 0.005° latitude/longitude ≈ 0.55 km × 0.55 km ≈ 0.306 km²). Due to Github file-size constraints, we have broken the file into similar-sized chunks, and then compressed them using the fast, lossless compression algorithm zstd. First, decompress each .zst chunk using the following command in Terminal (or equivalent): zstd -d 'XXXHgPredSpatRFlin.nc_chunk_a*.zst', and then combine chunks a to n using the following Terminal (or equivalent) command: cat XXXHgPredSpatRFlin.nc_chunk_* > XXXHgPredSpatRFlin.nc. You can import the NetCDF file in R using the ncdf4 package and its function nc_open. This produces a ncdf4 object that can be converted to a SpatRaster object using the rast function in package terra.

• ade4, adegraphics, adespatial, boot, dismo, dplyr, fields, gbm, geodata, geosphere, ggplot2, kableExtra, ncdf4, patchwork, performance, pdp, randomForestExplainer, ranger, rnaturalearthdata, sf, sjPlot, sp, spatialRF, spatstat, spdep, terra, tidyverse, usdm