Polar stereographic projection with latitude of standard parallel

Q&A
2

Not sure of this but your second listing uses the WKT2 whilst the first uses WKT1, and it says

PARAMETER["latitude_of_origin",60],

So maybe (likely) it was applied (that’s a work of GDAL/PROJ that are called by GMT) and is the WKT1 report that is less explicit.

And BTW, EPSG:4313 is for Belgium

julia> println(epsg2wkt(4313))
GEOGCS["BD72",
    DATUM["Reseau_National_Belge_1972",
        SPHEROID["International 1924",6378388,297,
            AUTHORITY["EPSG","7022"]],
        AUTHORITY["EPSG","6313"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST],
    AUTHORITY["EPSG","4313"]]