EPSG:10016 DEPRECATED

COORDINATEOPERATION["NAD83 to NAVD88 height (4)",
    VERSION["NGS-US Conus NE"],
    SOURCECRS[
        GEOGCRS["NAD83",
            DATUM["North American Datum 1983",
                ELLIPSOID["GRS 1980",6378137,298.257222101,
                    LENGTHUNIT["metre",1]]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["degree",0.0174532925199433]],
            CS[ellipsoidal,2],
                AXIS["geodetic latitude (Lat)",north,
                    ORDER[1],
                    ANGLEUNIT["degree",0.0174532925199433]],
                AXIS["geodetic longitude (Lon)",east,
                    ORDER[2],
                    ANGLEUNIT["degree",0.0174532925199433]],
            ID["EPSG",4269]]],
    TARGETCRS[
        VERTCRS["NAVD88 height",
            VDATUM["North American Vertical Datum 1988"],
            CS[vertical,1],
                AXIS["gravity-related height (H)",up,
                    LENGTHUNIT["metre",1]],
            ID["EPSG",5703]]],
    METHOD["Geographic3D to GravityRelatedHeight (gtx)",
        ID["EPSG",9665]],
    PARAMETERFILE["Geoid (height correction) model file","g2003u04.bin"],
    OPERATIONACCURACY[0.05],
    USAGE[
        SCOPE["Derivation of approximate gravity-related heights from GNSS observations."],
        AREA["United States (USA) - CONUS onshore north of 41°N and east of 78°W - Connecticut, Maine, Massachusetts, New Hampshire, Rhode Island and Vermont; New Jersey north of 41°N; New York and Pennsylvania north of 41°N and east of 78°W."],
        BBOX[41,-78,47.47,-66.91]],
    ID["EPSG",10016],
    REMARK["Uses Geoid03 hybrid model. See information source for further information. Note: Source CRS is 2D, used in this application of the method as a pseudo-3D CRS."]]

proj4.defs("EPSG:10016","+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +inv +proj=vgridshift +grids=g2003u04.bin +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 to NAVD88 height (4)",
    "source_crs": {
        "type": "GeographicCRS",
        "name": "NAD83",
        "datum": {
            "type": "GeodeticReferenceFrame",
            "name": "North American Datum 1983",
            "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"
                }
            ]
        },
        "id": {
            "authority": "EPSG",
            "code": 4269
        },
        "remarks": "Longitude is POSITIVE EAST. The adjustment included connections to Greenland and Mexico but the system was not adopted there. For applications with an accuracy of better than 1m replaced by NAD83(HARN) in the US and PRVI and by NAD83(CSRS) in Canada."
    },
    "target_crs": {
        "type": "VerticalCRS",
        "name": "NAVD88 height",
        "datum": {
            "type": "VerticalReferenceFrame",
            "name": "North American Vertical Datum 1988"
        },
        "coordinate_system": {
            "subtype": "vertical",
            "axis": [
                {
                    "name": "Gravity-related height",
                    "abbreviation": "H",
                    "direction": "up",
                    "unit": "metre"
                }
            ]
        },
        "id": {
            "authority": "EPSG",
            "code": 5703
        }
    },
    "method": {
        "name": "Geographic3D to GravityRelatedHeight (gtx)",
        "id": {
            "authority": "EPSG",
            "code": 9665
        }
    },
    "parameters": [
        {
            "name": "Geoid (height correction) model file",
            "value": "g2003u04.bin",
            "id": {
                "authority": "EPSG",
                "code": 8666
            }
        }
    ],
    "accuracy": "0.05",
    "scope": "Derivation of approximate gravity-related heights from GNSS observations.",
    "area": "United States (USA) - CONUS onshore north of 41°N and east of 78°W - Connecticut, Maine, Massachusetts, New Hampshire, Rhode Island and Vermont; New Jersey north of 41°N; New York and Pennsylvania north of 41°N and east of 78°W.",
    "bbox": {
        "south_latitude": 41,
        "west_longitude": -78,
        "north_latitude": 47.47,
        "east_longitude": -66.91
    },
    "id": {
        "authority": "EPSG",
        "code": 10016
    },
    "remarks": "Uses Geoid03 hybrid model. See information source for further information. Note: Source CRS is 2D, used in this application of the method as a pseudo-3D CRS."
}

<?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=g2003u04.bin +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=g2003u04.bin +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 ( 10016, 'EPSG', 10016, '+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +inv +proj=vgridshift +grids=g2003u04.bin +multiplier=1 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1', 'None');