Tag: orthometric height

  • Estimating heights with subsidence changes using NGS data and tools

    Estimating heights with subsidence changes using NGS data and tools

    This column details the potential effects of crustal movement on published heights in various regions of the United States.

    In my last column (in the April 2021 Survey Scene), I mentioned that the National Geodetic Survey (NGS) announced that it is suppressing height information in Southeast Texas.

    The April column also highlighted one of NGS’ four use cases – “Use Case 1: Flood Mapping.” The case study discusses the Elevation Certificate (CE) Example, Flood Insurance Rate Map (FIRM) and Flood Insurance Study (FIS).

    The column highlighted the potential effects of subsidence on published heights in the Houston region, which implied that most of the published heights that are based on older surveys in the region are not current or accurate.

    This column will provide more details of the suppression of heights in the Southeast Texas region, and potential effects of crustal movement on published heights in other regions of the United States.

    NGS announcement to suppress height information for Southeast Texas. (Image: NGS)
    NGS announcement that it suppressed height information for Southeast Texas. (Image: NGS)

    According to NGS’ announcement, only 28 marks will have publicly available orthometric heights on NGS datasheets in Southeast Texas.

    The “Link to Map: SE TX Valid Ortho. Heights” button provides the benchmarks available to users (see the box titled “Link to Map SE TX Valid Ortho Heights”). The website provides links to the published stations.

    Link to Map SE TX Valid Ortho Heights. (Image: NGS website)
    Link to Map SE TX Valid Ortho Heights. (Image: NGS website)

    Clicking on an icon provides the PID and name of the station with a link to a datasheet. Click  “Get Datasheet” for a datasheet of the station. Below is an excerpt from the datasheet of Station P 1200.

    Excerpt from Datasheet of Station P 1200.(Image: NGS Website)
    Excerpt from Datasheet of Station P 1200.(Image: NGS Website)

    Let’s address why NGS is suppressing the stations in Southeast Texas. My last column provided plots depicting the amount of movement in the Harris-Galveston, Texas, region. See the box titled “Estimate of Amount of Subsidence in 5 Years in Harris-Galveston, Texas, Region – Units Feet.”

    As indicated in the plot, some of the marks are estimated to have moved almost ½ foot (approximately 0.15 meters) in 5 years. In addition, some of the relative height differences approach 1/3 of a foot (approximately 0.1 meter) between neighboring stations. See the highlighted stations in the box titled “Estimate of Amount of Subsidence in 5 Years in Harris-Galveston, Texas, Region – Units Feet.”

    Estimate of Amount of Subsidence in 5 Years in the Harris-Galveston, Texas, Region – Units Feet. (Image: David Zilkoski)
    Estimate of Amount of Subsidence in 5 Years in the Harris-Galveston, Texas, Region – Units Feet. (Image: David Zilkoski)

    The last major releveling incorporated into NGS’ Database in the Harris-Galveston, Texas, region was performed more than 30 years ago in the 1986/1987 timeframe. Therefore, some of the published stations in the region could have subsided more than three feet (or about a meter).

    As stated in NGS’ Blueprint 3, “Most leveling data in NGS archives comes from the mid-20th century, in support of the NAVD 88 project.” Of course, most regions of the United States are not subsiding at the same rates as in the Houston-Galveston, Texas, region.

    In a previous newsletter, I discussed NGS’ second Multi-Year CORS Solution of the National CORS (MYCS2). I downloaded the coordinates and velocities from NGS’ website and created a plot of the vertical velocities. For those who prefer to use feet as opposed to meters, I provided velocities with units in feet/year and mm/year.

    See the boxes titled “Estimate of Velocity Rates Based on MYCS2 – CONUS (feet/year),” “Estimate of Velocity Rates Based on MYCS2 – Alaska (feet/year),” “Estimate of Velocity Rates Based on MYCS2 – CONUS (mm/year)” and “Estimate of Velocity Rates Based on MYCS2 – Alaska (mm/year).”

    It should be noted that the intent of these four plots is to provide a wide-ranging view of the values and some of the variation in rates across the United States.

    Estimate of Velocity Rates Based on MYCS2 – CONUS (feet/year). (Image: David Zilkoski)
    Estimate of Velocity Rates Based on MYCS2 – CONUS (feet/year). (Image: David Zilkoski)
    Estimate of Velocity Rates Based on MYCS2 – CONUS (feet/year). (Image: David Zilkoski)
    Estimate of Velocity Rates Based on MYCS2 – CONUS (feet/year). (Image: David Zilkoski)
    Estimate of Velocity Rates Based on MYCS2 – CONUS (mm/year). (Image: David Zilkoski)
    Estimate of Velocity Rates Based on MYCS2 – CONUS (mm/year). (Image: David Zilkoski)
    Estimate of Velocity Rates Based on MYCS2 – Alaska (mm/year). (Image: David Zilkoski)
    Estimate of Velocity Rates Based on MYCS2 – Alaska (mm/year). (Image: David Zilkoski)

    The rates appear to be small in most regions of the United States. As an example, the rates are all less than -0.0062 feet/year (-0.0019 meters/year) in the Lake Norman region in North Carolina (see the box titled “Potential Subsidence Rates in the Lake Norman Region in North Carolina). It would take many years for the crustal movement to make a difference to some projects in this region.

    Potential Subsidence Rates in the Lake Norman Region in North Carolina. (Image: David Zilkoski)
    Potential Subsidence Rates in the Lake Norman Region in North Carolina. (Image: David Zilkoski)

    That said, let’s look at another region of the country. For example, in the vicinity of Maryville, Missouri, the rate of subsidence is around -0.0187 feet/year (-0.0057 meters/year). See the box titled “Potential Subsidence Rates in the Maryville, Missouri, Region.” These subsidence rates don’t appear to be large values but if you take into account the last time the height of a mark was established by leveling data it could result in a large difference from the true orthometric height.

    Potential Subsidence Rates in the Maryville, Missouri, Region. (Image: David Zilkoski)
    Potential Subsidence Rates in the Maryville, Missouri, Region. (Image: David Zilkoski)

    According to NGS’ database, it appears that many of the marks in the Maryville, Missouri, region were last leveled in 1935. I used NGS’ Passive Mark Lookup tool and Leveling Project Page tool to identify the marks and associated leveling lines in the area of the CORS stations in the Maryville, Missouri, region.

    I described the Passive Mark Lookup webtool in a previous column. As previously mentioned, these subsidence rates all seem very small, but if you take into account the last time the height of mark was established by leveling data, the subsidence value can be very large.

    See the box titled “Potential Subsidence in 86 Years in the Maryville, Missouri, Region.” The box indicates that, if you account for the last 86 years (2021 – 1935), the potential subsidence exceeds 1½ feet (-1.6082 feet, -0.4902 meters).

    Potential Subsidence in 86 Years in the Maryville, Missouri, Region. (Image: David Zilkoski)
    Potential Subsidence in 86 Years in the Maryville, Missouri, Region. (Image: David Zilkoski)

    Continuing across the country to Colorado, the box titled “Potential Subsidence Rates in the Grand Junction Region, Colorado,” provides the estimate of subsidence rates in Mesa County, Colorado. As the plot indicates, the rates vary between -0.0046 feet/year (-1.4 mm/year) and -0.0128 feet/year (-3.9 mm/year). Once again, these rates all seem relatively small but many of the marks near CORS MC06 were last leveled in 1985. This means the potential change in height could be as large as 0.2592 feet (0.0792 meters).

    Potential Subsidence Rates in the Grand Junction. Colorado, Region. (Image: David Zilkoski)
    Potential Subsidence Rates in the Grand Junction Region, Colorado. (Image: David Zilkoski)

    Obviously, this is only an estimate of the subsidence in the region and the actual amount of subsidence is unknown since the last time the mark was leveled. These estimates are based on the MYCS2, which uses current data to estimate the velocity. The processing included data spanning 1996 to 2016 (week 0834 to 1933), 1099 weeks or about 21 years in total.

    The point of this column is not to provide the exact change in height of a mark, but to highlight that the publicly available orthometric height on a NGS datasheet may not be up to date based on crustal movement. The new modernized National Spatial Reference System will enable users to determine an accurate, current height on a mark and be able to efficiently and effectively monitor changes in a mark’s height.

    As stated in NGS’ announcement to suppress the heights in Southeast Texas, the agency has developed tools to assist users in submitted data to NGS. See the box titled “Excerpt from NGS Announcement to Suppresses Height Information for Southeast Texas.”

    Excerpt from NGS Announcement to Suppresses Height Information for Southeast Texas. (Image: NGS website)

    This assistance is for every user, not just for individuals performing surveys in Southeast Texas. NGS has Regional Geodetic Advisors throughout the United States.

    NGS Regional Geodetic Advisors. (Image: NGS Website)
    NGS Regional Geodetic Advisors. (Image: NGS Website)

    The Regional Geodetic Advisors provide guidance and assistance to constituents within their region. They are subject-matter experts in geodesy and regional geodetic issues. These individuals can assist users that are planning GNSS campaigns to re-densify the network.

    NGS also provides a website detailing how users can help densify the network to prepare for the new, modernized North American-Pacific Geopotential Datum of 2022 (NAPGD2022). See the box titled “NGS GPS on Bench Marks Webpage.”

    As mentioned in previous newsletters, a benefit of the new modernized National Spatial Reference System (NSRS) will facilitate the establishment of consistent, accurate NAPGD2022 GNSS-derived orthometric heights.

    NGS GPS on Bench Marks webpage. (Image: NGS Website)
    NGS GPS on Bench Marks webpage. (Image: NGS Website)

    This column provided details on the suppression of heights in the Southeast Texas region, and potential effects of crustal movement on published heights in other regions of the United States. NGS suppressed the heights in the Southeast Texas region because of the large amount of crustal movement since the last time the heights of the marks were established.

    As indicated by NGS’ MYCS2 velocities, every mark could be affected by crustal movement. In my opinion, the question a user should be asking is “How much has the height of the mark changed since it was last determined? Not, “Has the height of the mark changed?”

  • NGS releases beta tool for obtaining geodetic information

    NGS releases beta tool for obtaining geodetic information

    NGS has developed a new beta tool for obtaining geodetic information about a passive mark in their database. This column will highlight some features (available as of Oct. 5, 2020) that may be of interest to GNSS users. It provides all of the information about a station in a more user-friendly format. The box titled “Passive Mark Lookup Tool” is an example of the webtool. The tool provides a lot of information so I have separated the output of the tool into several boxes titled “Passive Mark Lookup Tool — A through D.”

    I will highlight several attributes that I believe will be very useful to users, especially users of leveling-derived and GNSS-derived orthometric heights. I’ve highlighted several attributes in the box titled “Passive Mark Lookup Tool — A” that are important to users such as published coordinates, their datum and source, Geoid18 value, GNSS Useable, and the date of last recovery. All of these values are available on a NGS datasheet but, in my opinion, this provides the information in a more user-friendly format.

    Passive Mark Lookup Tool — A

    (https://beta.ngs.noaa.gov/datasheets/passive-marks/index.html)

    Image: National Geodetic Survey
    Image: National Geodetic Survey

    One calculation that the user can easily compute for marks that have been leveled to and occupied by GNSS equipment, is the difference between the published leveling-derived orthometric height and the computed GNSS-derived orthometric height. This may indicate that the mark has moved since the last time it was leveled to or that its height coordinate has been readjusted since the creation of the published geoid model.

    The table below provides the calculation using the data from the box titled “Passive Mark Lookup Tool — A.” The calculation [HGNSS = hGNSS — NGeoid18; Difference = HGNSS — HNAVD 88] has been described in several of my previous columns (this one, for example).

    Data: National Geodetic Survey
    Data: National Geodetic Survey

    In this example, the difference between the GNSS-derived orthometric height and the Published NAVD 88 height is 6.1 cm. NGS is looking for comments on this beta webtool so if users would like this computation added to the tool, they should send a comment to NGS using the link provided on the site (This is a beta product. NGS is interested in your feedback concerning its function and usability as well as how users would like to interact with NGS datasheet information in the future. Email us at [email protected].) So, the user should ask the question, did the station move since the last time it was leveled?

    Another attribute that would be nice to be part of this tool is which station was used to create the hybrid geoid model. As of Oct. 5, 2020, users have to go to the Geoid18 webpage to get the information. The Excel file and shapefiles provide whether the station was used to create the Geoid18 model. In the case of this example, KK1531, CHAMBERS, the mark was not used in the creation of Geoid18 so NGS felt that the station may have moved and/or the GPS on Bench Mark residual was large relative to its neighbors. See NGS’s technical report on Geoid18 for more information on the creation of Geoid18. The GPS on Bench Mark residual analysis was described in several of my previous columns (see “The differences between Geoid18 values and NAD 83, NAVD 88 values” and “NGS 2018 GPS on BMs program in support of NAPGD2022 — Part 6” for examples).

    The webtool provides a map depicting the location of the station, photos (if available), and previously published, superceded values of the mark. See the box titled “Passive Mark Lookup Tool — B.”

    Passive Mark Lookup Tool — B

    Image: National Geodetic Survey
    Image: National Geodetic Survey

    In the example of Chambers, KK1531, no photos were available. It would be helpful if a user would provide photos to NGS when visiting this station. (Note: NGS has a webtool for users to submit recovery information about a mark as well as to provide current photos of the station.) The new Passive Mark webtool also provides information about the survey projects that the mark has been involved in such as leveling and GNSS projects.

    In this example, mark CHAMBERS was leveled to in a 1984 first-order, class 2 leveling project (Leveling Line number L24838/6) and, in 1995, the mark was part of a GNSS project (GNSS Project GPS1010). It also provides all the descriptive text and recovery information (See boxes titled “Passive Mark Lookup Tool – C” and “Passive Mark Lookup Tool – D”).

    Passive Mark Lookup Tool — C

    Data: National Geodetic Survey
    Data: National Geodetic Survey

    Passive Mark Lookup Tool — D

    Image: National Geodetic Survey
    Image: National Geodetic Survey

    I want to highlight a few other attributes of this webtool. The station, PID AA3862, has an interesting attribute that users should take note of; that is, the NAD 83 (2011) position source is NO CHECK. See box titled “Passive Mark Page for PID AA3862.”

    This means that the mark’s NAD 83 (2011) coordinates were determined without redundant observations. This is not a good survey practice but there are times that a project may contain check observations for some purpose or, more likely, the mark did contain other GNSS vector but they were rejected in the final adjustment. Either way, a good survey practice would be for users to verify the coordinates of these marks before using them.

    Passive Mark Page for PID AA3862

    Data: National Geodetic Survey
    Data: National Geodetic Survey

    As previously mentioned, the tool provides the location of the station on a map and photos if they are available. This is a really nice feature for anyone searching for the mark. The map can be enlarged as well reduced by clicking on the box. See boxes titled “Passive Mark Page for PID AA3862” and “Photos of Mark PID AA3862.” The box titled “Photos of Mark PID AA3862” provides all three photos of mark PID AA3862.

    Photos of Mark PID AA3862

    Photo: National Geodetic Survey
    Photo: National Geodetic Survey
    Photo: National Geodetic Survey
    Photo: National Geodetic Survey

    Photo: National Geodetic Survey
    Photo: National Geodetic Survey

    It should be noted, according to the Geoid18 GPS on BMs dataset that users can download, this station, AA3862, was not used in the creation of Geoid18. The table below provides the difference between the GNSS-derived orthometric height and the published NAVD 88 height.

    In this example, the difference between the GNSS-derived orthometric height and the published NAVD 88 height is 9.9 cm. Also, the webtool provides the network accuracy values for the station. In this example, the horizontal network accuracy is 20.65 cm and the vertical network accuracy value is 14.50 cm (see highlighted values in box titled “Passive Mark Page for PID AA3862”). These are very large network accuracy values. This should be a flag to anyone that is using this station as control.

    Data: National Geodetic Survey
    Data: National Geodetic Survey

    As I previously mentioned, as a beta site, users should verify all information from the site. NGS is requesting feedback on this tool so they can improve it and make it an operational webtool. I encourage everyone to access the tool and check out a few of their favorite marks, and then send an email to NGS informing them of what you like, what you would like to change, and what you would like to see added to the tool.

    NGS is releasing this tool as a beta product to get feedback from users. As NGS states in the heading of the tool, they are interested in your feedback concerning its function and usability as well as how users would like to interact with NGS datasheet information in the future. Email NGS at [email protected].

    One last item that may be of interest to GNSS users is that NGS, working with the University Corporation for Atmospheric Research (UCAR), developed another online GNSS lesson (see box titled “New GNSS Lesson by NGS and UCAR”). These lessons are free but users must sign up to access the website and lesson.

    New GNSS Lesson by NGS and UCAR

    Image: National Geodetic Survey
    Image: National Geodetic Survey