EPSG:9986

COORDINATEOPERATION["NAD83(CSRS)v6 to CGVD28 height (1)",
    VERSION["NRC-Can CGG2000 2010"],
    SOURCECRS[
        GEOGCRS["NAD83(CSRS)v6",
            DATUM["North American Datum of 1983 (CSRS) version 6",
                ELLIPSOID["GRS 1980",6378137,298.257222101,
                    LENGTHUNIT["metre",1]],
                ANCHOREPOCH[2010]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["degree",0.0174532925199433]],
            CS[ellipsoidal,3],
                AXIS["geodetic latitude (Lat)",north,
                    ORDER[1],
                    ANGLEUNIT["degree",0.0174532925199433]],
                AXIS["geodetic longitude (Lon)",east,
                    ORDER[2],
                    ANGLEUNIT["degree",0.0174532925199433]],
                AXIS["ellipsoidal height (h)",up,
                    ORDER[3],
                    LENGTHUNIT["metre",1]],
            ID["EPSG",8251]]],
    TARGETCRS[
        VERTCRS["CGVD28 height",
            VDATUM["Canadian Geodetic Vertical Datum of 1928"],
            CS[vertical,1],
                AXIS["gravity-related height (H)",up,
                    LENGTHUNIT["metre",1]],
            ID["EPSG",5713]]],
    METHOD["Geographic3D to GravityRelatedHeight (NRCan byn)",
        ID["EPSG",1060]],
    PARAMETERFILE["Geoid (height correction) model file","HT2_2010v70.byn"],
    OPERATIONACCURACY[0.05],
    USAGE[
        SCOPE["Derivation of gravity-related heights from GNSS observations."],
        AREA["Canada - onshore - Alberta; British Columbia; Manitoba south of 59°N; New Brunswick; Newfoundland, Labrador between 52°50’N and 54°30’N; Northwest Territories south west of a line between 60°40’N, 110°W and the coast at 132°W; Nova Scotia; Ontario south of 52°20'N; Prince Edward Island; Quebec - mainland south of 55°N and west of 64°W, north of 55°N between 70°W and 66°W, coastal area between 66°W and 64°W, and Anticosta island; Saskatchewan south of 58°30’N; Yukon."],
        BBOX[41.67,-141.01,69.81,-52.54]],
    ID["EPSG",9986],
    REMARK["Hybrid geoid model valid at epoch 2010.0. Derived at epoch 1997.0 through NAD83(CSRS)v3 and modified to include correction derived from the Canada velocity grid v7 for propagation of height between 1997 and 2010. Also used as CT to CGVD28(v2.0) height."]]

proj4.defs("EPSG:9986","+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +inv +proj=vgridshift +grids=ca_nrc_HT2_2010v70.tif +multiplier=1 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1");

{
    "$schema": "https://proj.org/schemas/v0.7/projjson.schema.json",
    "type": "Transformation",
    "name": "NAD83(CSRS)v6 to CGVD28 height (1)",
    "source_crs": {
        "type": "GeographicCRS",
        "name": "NAD83(CSRS)v6",
        "datum": {
            "type": "GeodeticReferenceFrame",
            "name": "North American Datum of 1983 (CSRS) version 6",
            "anchor_epoch": 2010,
            "ellipsoid": {
                "name": "GRS 1980",
                "semi_major_axis": 6378137,
                "inverse_flattening": 298.257222101
            }
        },
        "coordinate_system": {
            "subtype": "ellipsoidal",
            "axis": [
                {
                    "name": "Geodetic latitude",
                    "abbreviation": "Lat",
                    "direction": "north",
                    "unit": "degree"
                },
                {
                    "name": "Geodetic longitude",
                    "abbreviation": "Lon",
                    "direction": "east",
                    "unit": "degree"
                },
                {
                    "name": "Ellipsoidal height",
                    "abbreviation": "h",
                    "direction": "up",
                    "unit": "metre"
                }
            ]
        },
        "id": {
            "authority": "EPSG",
            "code": 8251
        }
    },
    "target_crs": {
        "type": "VerticalCRS",
        "name": "CGVD28 height",
        "datum": {
            "type": "VerticalReferenceFrame",
            "name": "Canadian Geodetic Vertical Datum of 1928"
        },
        "coordinate_system": {
            "subtype": "vertical",
            "axis": [
                {
                    "name": "Gravity-related height",
                    "abbreviation": "H",
                    "direction": "up",
                    "unit": "metre"
                }
            ]
        },
        "id": {
            "authority": "EPSG",
            "code": 5713
        }
    },
    "method": {
        "name": "Geographic3D to GravityRelatedHeight (NRCan byn)",
        "id": {
            "authority": "EPSG",
            "code": 1060
        }
    },
    "parameters": [
        {
            "name": "Geoid (height correction) model file",
            "value": "HT2_2010v70.byn",
            "id": {
                "authority": "EPSG",
                "code": 8666
            }
        }
    ],
    "accuracy": "0.05",
    "scope": "Derivation of gravity-related heights from GNSS observations.",
    "area": "Canada - onshore - Alberta; British Columbia; Manitoba south of 59°N; New Brunswick; Newfoundland, Labrador between 52°50’N and 54°30’N; Northwest Territories south west of a line between 60°40’N, 110°W and the coast at 132°W; Nova Scotia; Ontario south of 52°20'N; Prince Edward Island; Quebec - mainland south of 55°N and west of 64°W, north of 55°N between 70°W and 66°W, coastal area between 66°W and 64°W, and Anticosta island; Saskatchewan south of 58°30’N; Yukon.",
    "bbox": {
        "south_latitude": 41.67,
        "west_longitude": -141.01,
        "north_latitude": 69.81,
        "east_longitude": -52.54
    },
    "id": {
        "authority": "EPSG",
        "code": 9986
    },
    "remarks": "Hybrid geoid model valid at epoch 2010.0. Derived at epoch 1997.0 through NAD83(CSRS)v3 and modified to include correction derived from the Canada velocity grid v7 for propagation of height between 1997 and 2010. Also used as CT to CGVD28(v2.0) height."
}

<?xml version="1.0" encoding="utf-8"?>
<Map srs="+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +inv +proj=vgridshift +grids=ca_nrc_HT2_2010v70.tif +multiplier=1 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1">
	<Layer srs="+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +inv +proj=vgridshift +grids=ca_nrc_HT2_2010v70.tif +multiplier=1 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1">
	</Layer>
</Map>

INSERT into spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext) values ( 9986, 'EPSG', 9986, '+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +inv +proj=vgridshift +grids=ca_nrc_HT2_2010v70.tif +multiplier=1 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1', 'None');