View on GitHub

# Data Products and Interpolation

## Oak Ridges Moraine Groundwater Program

Spatial interpolation of point data

# Interpolation of Hourly scalars to Sub-daily Basins

Bolow is an itemized description of the interpolation process, with open data access. All interpolated (i.e., “vector”) data are automatically updated and maintained using the ORMGP-FEWS system system.

workflow: preprocessMSCtoBasinsHourly.xml

## Meteorological Service of Canada (MSC) scalars

1. Scrape MSC for recent data, executed from FEWS.
2. Import scraped MSC hourly scalars into FEWS.
3. Export hourly MSC NetCDF file (*.nc) from FEWS, from 1989-10-01

_exportMSChourlyNetcdf.nc $T_a, p_a, r, \text{vis}, u, u_\alpha$

4. Interpolate to 10km sub-watersheds, using a python script executed from config file pyMSChourliesToBasin.xml that executes: ncMSCtoHourlyBasinNetCDF.py.
• Air temperature $(T_a)$, relative humidity $(r)$ and wind speeds $(u)$ are interpolated using a radial basis function (RBF) (with a cubic kernel and a smoothing factor $\lambda=1/1000$ to prevent singular matrices).
• Air/barometric pressure $(p_a)$ are first corrected for elevation then interpolated using the same cubic RBF.
• Wind directions $(u_\alpha)$ are split into their x-y components, each interpolated separately using a cubic RBF before returned to an angle. A sample result is shown below:
5. Interpolated $(T_a, r, u)$ are applied to compute potential evaporation flux [m/s] (Novák, 2012):
$E_a=7.46\times 10^{-9} \cdot (a + ub) d_a$ $d_a=(1-r) \cdot e_s(T_a)$
1. Save to NetCDF (.nc) for import back to FEWS.

_exportMSChourlyNetcdf_interp.nc $T_a, p_a, r, u, E_a$ hourly basins

## Hourly Basin to 6-hourly Basin interpolation

1. Hourly aggregation to 6-hourly time intervals (00:00 06:00 12:00 18:00 UTC) is performed in FEWS using the:

These data have a set expiry.

2. Export 6-hourly, basin-interpolated $P_\text{HRDPA}, T_a, p_a, r, u$ to NetCDF.

yyyyMMddHHmm-6hourlyBasin.nc

to be altered to $P_\text{HRDPA}, P_R, P_S, T_a, p_a, r, u, E_a$

# 6-hourly Precipitation to Basins

The 6-hourly CaPA-RDPA precipitation $(P)$ field is a gridded raster that is routinely scraped of open web resources and proportioned to the sub-watersheds using our ORMGP-FEWS system, using the Interpolation: SpatialAverage transformation.

# References

Novák, V., 2012. Evapotranspiration in the Soil-Plant-Atmosphere System. Springer Science+Business Media Dordrecht. 253pp.