Tag: CORS solution

  • NGS will soon compute third multi-year CORS solution

    NGS will soon compute third multi-year CORS solution

    On Aug. 5, the National Geodetic Survey (NGS) stated it will be updating the NOAA CORS to be aligned with the latest International Terrestrial Reference frame, ITRF2020 (see below). As stated in the announcement, NGS will soon compute a third multi-year continuously operating reference station (CORS) solution, MYCS3.

    The last multi-year CORS solution, MYCS2, was performed by NGS in 2019. I discussed the MYCS2 in my February 2019 and April 2019 columns. This new multi-year CORS solution will be important to the 2022 modernized National Spatial Reference System (NSRS), because NGS will establish a strict mathematical relationship between the 2022 NSRS frames and the ITRF2020 frame. This will allow direct access to the NSRS (NOAA Technical Report NOS NGS 67).

    NGS Aligns National System to Global Reference Frame

    August 5, 2022

    The International Global Navigation Satellite System (GNSS) Service, which provides GNSS data products globally, recently released a new GNSS-only version of the International Terrestrial Reference Frame. This provides GNSS users access to the reference frame through coordinate functions for a global set of reference stations. In response, NGS will soon compute the multi-year Continuously Operating Reference Station (CORS) Solution 3, which will modernize the National Spatial Reference System. Aligning the National Spatial Reference System with the updated global reference frame will allow greater access for the global community of scientists, educators, and commercial users of location science.

    For more information, contact: Phillip McFarland

    As in the past, the multi-year CORS solution will mean that the NOAA CORS coordinates will be updated to be consistent with the latest International Terrestrial Reference Frame of 2020 (ITRF2020). The International GNSS Service provides information about its GNSS products and services. Readers can find information on the latest International Terrestrial Reference Frame 2020 here. This column will provide basic information on the ITRF2020. Please note: NGS stated that it will soon start computing the third multi-year CORS solution, but — as of October — all NOAA CORS coordinates are still based on MYCS2 and provide coordinates in ITRF2014 epoch 2010.00 and NAD 83 (2011, MA11, PA11) epoch 2010.00. As in the past, NGS will provide advance notice before publishing the results of its third multi-year CORS solution.

    A document on the ITRF website stated the ITRF2020 is expected to be an improved solution compared to the previous solution, ITRF2014. It listed several innovations introduced in the ITRF2020 processing.

    Description from ITRF2020 Document

    ITRF2020 is the new realization of the International Terrestrial Reference System. Following the procedure already used for previous ITRF solutions, the ITRF2020 uses as input data time series of station positions and Earth Orientation Parameters (EOPs) provided by the Technique Centers of the four space geodetic techniques (VLBI, SLR, GNSS and DORIS), as well as local ties at colocation sites. Based on completely reprocessed solutions of the four techniques, the ITRF2020 is expected to be an improved solution compared to ITF2014. A number of innovations were introduced in the ITRF2020 processing, including:

    • The time series of the four techniques were stacked all together, adding local ties and equating station velocities and seasonal signals at colocation sites;
    • Annual and semi-annual terms were estimated for stations of the 4 techniques with sufficient time spans;
    • Post-Seismic Deformation (PSD) models for stations subject to major earthquakes were determined by fitting GNSS/IGS data. The PSD models were then applied to the 3 other technique time series at earthquake colocation sites.

    The box below provides a good summary of the International Reference Frame and why it’s important to the scientific community as well as the surveying and mapping community. Readers can download the article from the June 2022 International GNSS Service Issue 4 newsletter. Users also can sign up to receive notices and newsletters from the International GNSS Service.

    ITRF2020: A new release of the International Terrestrial Reference Frame By Zuheir Altamimi

    What is the current rate of sea level rise in different regions of the globe? How does our Earth deform under the effect of plate tectonics, seismic phenomena, or the melting of ice caps? How the Earth’s center of mass is varying? How to determine the position of a point on the surface of a constantly deforming Earth and compare it to positions estimated decades apart? The answers to these fundamental questions for understanding the dynamics of our planet require the availability of a global, long-term stable terrestrial reference frame, but preferably a standard reference so to ensure interoperability and consistency of various measurements collected by sensors on the ground, or via artificial satellites. The International Terrestrial Reference Frame (ITRF) is the standard reference recommended by a number of international scientific organizations, including the International Union of Geodesy and Geophysics (IUGG) and the International Association of Geodesy (IAG) for earth science, satellite navigation and operational geodesy applications. The ITRF is an international effort that is built on the investments of space and mapping agencies, universities and research groups in operating geodetic observatories, archiving and analyzing the collected geodetic observations to derive not only the ITRF, but also critical geodetic products for science and society.

    The ITRF integrates and unifies technique-specific reference frames provided by the four IAG’s international services of space geodetic technique (DORIS/IDS, GNSS/IGS, SLR/ILRS, VLBI/ IVS). It is supplied to the users in the form of temporal coordinates of more than 1500 stations, Earth Orientation Parameters, as well as parametric functions describing nonlinear station motions: seasonal signals due to mainly loading effects and post-seismic deformations for sites subject to major earthquakes. It is necessary to regularly update the ITRF (approximately every 5 years) in order to benefit from continuous observations so to improve its accuracy, considering station position temporal variations due to geophysical phenomena.

    The ITRF is maintained by a research group at IGN-France and IPGP (Institut de Physique de Globe de Paris), and whose new release called ITRF2020 was published on April 15 and accessible here: https://itrf.ign.fr/en/solutions/ITRF2020. The ITRF2020 brings significant improvements compared to previous achievements: it confirms the estimate of the position of the center of mass of the Earth as it was determined in 2016, but also provides its seasonal variations; it improves the accuracy of the scale of the frame at the millimeter level, which represents a gain in precision of a factor of 8 on the measurement of the size of the Earth (compared to that determined in 2016); it provides a precise quantification of co- and post-seismic displacements caused by devastating earthquakes, such as that of Sumatra in 2004, Chile in 2010 and Japan in 2011. The IAG Services rely on the ITRF to align their geodetic products to it, and therefore disseminate it widely among the various users. In particular, using the IGS products, such as the orbits, allows a universal access in space and time to the ITRF.

    As stated in the article by Zuheir Altamimi, ITRF2020 involves IAG’s international services of four space geodetic techniques: DORIS/IDS, GNSS/IGS, SLR/ILRS, VLBI/ IVS. Computing an International Terrestrial Frame is very complex and requires analyses of difference types of geodetic and geophysical data. It is beyond the scope of this column, but online is more detailed technical information.

    For this column, I downloaded the station lists from the four space geodetic techniques and provided a few plots that depict the location and velocities of these sites. The box below depicts the location of the space geodetic techniques around the world. As indicated in the plot, some locations have more than one technique collocated at the same site.

    Plot of the Four Different Space Geodetic Techniques

    Image: David Zilkoski
    Image: Dave Zilkoski

    The following plots depict the locations using each space geodetic techniques: GNSS sites, DORIS sites, SLR sites and VLBI sites.

    Plot of GNSS Sites

    Image: David Zilkoski
    Image: Dave Zilkoski

     

    Plot of DORIS Sites

    Image: David Zilkoski
    Image: Dave Zilkoski

    Plot of SLR Sites

    Image: Dave Zilkoski
    Image: Dave Zilkoski

    Plot of VLBI Sites

    Image: Dave Zilkoski
    Image: Dave Zilkoski

    The box below shows the location of the techniques in the conterminous United States.

    Plot of the Four Different Space Geodetic Techniques in the CONUS

    Image: Dave Zilkoski
    Image: Dave Zilkoski

    The plot below depicts the sites in the state of Alaska.

    Plot of the Four Different Space Geodetic Techniques in the Alaska

    Image: Dave Zilkoski
    Image: Dave Zilkoski

    The images below depict each of the four space geodetic techniques in the conterminous United States.

    Plots of the Space Geodetic Techniques by Technique in the CONUS

    Image: Dave Zilkoski
    Plot of GNSS Sites in CONUS Image: Dave Zilkoski
    Image: Dave Zilkoski
    Plot of DORIS Sites in CONUS (Image: Dave Zilkoski)
    Image: Dave Zilkoski
    Plot of SLR Sites in CONUS (Image: Dave Zilkoski)
    Image: Dave Zilkoski
    Plot of VLBI Sites in CONUS (Image: Dave Zilkoski)

    Altamimi’s article on the ITRF2020 stated it is “necessary to regularly update the ITRF (approximately every 5 years) to account for station position temporal variations due to geophysical phenomena.” My February 2022 column discussed the tectonic plates and why is it necessary to account for movement in a geodetic reference frame. As I stated then, coordinates basically change because the Earth’s surface is moving due to the movement of major tectonic plates. See the box titled “What is Tectonic Shift?” for information about why it is called plate movement or tectonic shift. The world’s geodesists understand this and are attempting to manage the changing coordinates by providing a time-dependent component of the international terrestrial reference frame.

    Image: National Ocean Service Website
    Image: National Ocean Service website
    Image: National Ocean Service Website
    Image: National Ocean Service website

    The box below depicts the horizontal velocity based on the ITRF2020 velocities (downloaded on 08/12/2022).

    Plot of the Horizontal Velocity Vectors based on the ITRF2020 Velocities

    (Image: Dave Zilkoski)
    Image: Dave Zilkoski

    The box below depicts the horizontal velocities in the North America. These vectors look very similar to the velocities reported in my February 2022 column.

    Plot of the Horizontal Velocity Vectors in North America based on the ITRF2020 Velocities

    Image: Dave Zilkoski
    Image: Dave Zilkoski

    For a comparison to North America vectors, the box below depicts the velocity vectors in Europe.

    Plot of the Horizontal Velocity Vectors in Europe based on the ITRF2020 Velocities

    Image: Dave Zilkoski
    Image: Dave Zilkoski

    They are similar in magnitude, but not in direction. Once again, looking at the map of tectonic plates, North America is located mostly on the North American plate and Europe is on the Eurasian plate.

    Australia is on the Indo-Australian plate and has some fairly large horizontal velocities vectors. See the box below.

    Plot of the Horizontal Velocity Vectors in Australia based on the ITRF2020 Velocities

    Image: Dave Zilkoski
    Image: Dave Zilkoski

    So, what’s the difference between ITRF2014 and the new ITRF2020? The box below provides the 14 transformation parameters from ITRF2020 to ITRF2014. These transformation parameters have been estimated using 131 stations located at 105 sites. See the box “Plot of the Stations used in the Transformation Parameters from ITRF2020 to ITRF2014” for the location of these stations. Notice that the translation values in X,Y,Z are very small (<1.5 mm) between the two reference frames.

    Transformation Parameters from ITRF2020 to ITRF2014

    (https://itrf.ign.fr/en/solutions/ITRF2020)
    (https://itrf.ign.fr/en/solutions/ITRF2020)

    Transformation parameters at epoch 2015.0 and their rates from ITRF2020 to ITRF2014 (ITRF2014 minus ITRF2020)

    (https://itrf.ign.fr/en/solutions/ITRF2020)
    (https://itrf.ign.fr/docs/solutions/itrf2020/Transfo-ITRF2020_TRFs.txt)

    X,Y,Z are the coordinates in ITRF2020, and XS,YS,ZS are the coordinates in ITRF2014.

    Plot of the Stations used in the Transformation Parameters from ITRF2020 to ITRF2014

    Image: Dave Zilkoski
    Image: Dave Zilkoski

    The transformation parameters from ITRF2020 and past ITRFs are provided in the table below. As indicated in the table, most of the changes in X,Y and Z are very small since ITRF2005.

    Transformation Parameters from ITRF2020 to Past ITRFs

    (https://itrf.ign.fr/docs/solutions/itrf2020/Transfo-ITRF2020_TRFs.txt )
    (https://itrf.ign.fr/docs/solutions/itrf2020/Transfo-ITRF2020_TRFs.txt)

    As previously stated, the third multi-year CORS solution will be important to the new 2022 modernized National Spatial Reference System (NSRS) because NGS will establish a strict mathematical relationship between the 2022 NSRS frames and the ITRF2020 frame. This will allow direct access to the NSRS, according to NOAA Technical Report NOS NGS 67. Again, there will not be any changes to NGS’s NOAA CORS coordinates due to ITRF2020 until NGS completes its third multi-year CORS solution.

    Users can receive emails about the latest NGS News by signing up for NGS’s newsletters. These notices will highlight the release of new products, updates to existing services, progress reports for major projects, information about upcoming NGS-sponsored events, and job opportunities at NGS.

  • NGS releases beta coordinates and multi-year CORS solution

    NGS releases beta coordinates and multi-year CORS solution

    My last column discussed the preliminary results of NGS’ second Multi-Year CORS Beta Solution of the National CORS. Since my last column, NGS announced the release of the beta version of the hybrid geoid model GEOID18 and, on Feb. 15, NGS officially released the Beta CORS ITRF2014 coordinates and velocities.

    This column provides the official links to NGS website that provide the beta coordinates and information about the latest multi-year CORS solution. Below is the NGS announcement of the beta release of the updated coordinates.

    Excerpt from Feb. 15, 2019, "NOTICE: New BETA Coordinates Available for CORS and OPUS". (Screenshot: NGS)
    Excerpt from Feb. 15, 2019, “NOTICE: New BETA Coordinates Available for CORS and OPUS”. (Screenshot: NGS)

    NGS also provides a notice of the new beta coordinates on the National Geodetic Survey homepage, with a link to the Beta CORS ITRF14 coordinates (see the highlighted section below).

    National Geodetic Survey homepage. (Screenshot: NGS)
    National Geodetic Survey homepage. (Screenshot: NGS)

    Clicking on the hyperlink labeled BETA CORS ITRF14 Coordinates directs you to the Multi-Years CORS Solution informational homepage.

    Information on Multi-Year CORS Solution 2. (Screenshot: NGS)
    Information on Multi-Year CORS Solution 2. (Screenshot: NGS)

    By clicking on the CORS Home button, the user is directed to the Beta CORS page.

    Beta CORS release page. (Screenshot: NGS)
    Beta CORS release page. (Screenshot: NGS)

    This page clearly states that the ITRF2014 reference frame for CORS is available as a beta product. It also implies that these coordinates are being used in other beta products such as OPUS. I’ll address this later in this column.

    Users can obtain information about the MYCS and other related products and services such as Beta OPUS by clicking on links provided on the Beta CORS homepage.

    Accessing information about ITRF2014 frame in NGS beta products. (Screenshot: NGA)
    Accessing information about ITRF2014 frame in NGS beta products. (Screenshot: NGA)

    It should be noted that these values are considered “beta” and are available to users for testing and feedback. NGS provides a statement about its beta release products. Basically, it states that users should only use beta products to test their workflows and never for official or production work.

    The NGS beta release statement. (Screenshot: NGS)
    The NGS beta release statement. (Screenshot: NGS)

    To facilitate testing of the beta CORS coordinates and velocities, NGS provides links to other beta products that will use the MYCS 2 coordinates and velocities.

    By clicking on the link labeled BETA OPUS on the beta CORS homepage, the user is directed to the BETA OPUS webpage. This page clearly states that the beta OPUS routine uses the new ITRF2014 reference frame for CORS.

    The NGS beta OPUS webpage. (Screenshot: NGS)
    The NGS beta OPUS webpage. (Screenshot: NGS)

    NGS also provides a link to Beta OPUS Projects that use the MCYS2 coordinates and velocities.

    Beta OPUS Projects webpage. (Screenshot: NGS)
    Beta OPUS Projects webpage. (Screenshot: NGS)

    Once again, the Beta OPUS Projects website clearly states that the beta version is using the CORS coordinates and velocities from the MYCS2. It also states that, at this time, NGS will not accept ITRF2014 submissions for publication. As previously stated, NGS’ beta products are for users to test their workflows and should never be used for official or production work.

    The Beta CORS webpage provides a lot of valuable information on the processing and establishment of the multi-years CORS solution. I’ve highlighted several of the sections below.

    First, by clicking on the link MYCS2 Processing, the user is directed to the section that describes the data used and the processing strategy.

    Excerpt from Beta CORS Webpage – MYCS2 Processing. (Screenshot: NGS)
    Excerpt from Beta CORS Webpage – MYCS2 Processing. (Screenshot: NGS)

    The following are highlights from the section:

    • The processing included data spanning 1996 to 2016 and involved around 3050 CORS, IGS and other (e.g., NGA) stations.
    • The corresponding input and output data occupied about 25 TB on the NGS computers.
    • The residual time series in the early 1990s showed exceptionally noisy behavior at times, which were deleted in the alignment/velocity computation stage.
    • The processing was performed in 3 steps:

    1. The global processing step solves for orbits, Earth Orientation Parameters (EOPs), hourly tropospheric delay parameters and weekly global (IGS) station positions in an IGS-NNR frame.

    2. The CORS processing step ties the remaining CORS to global, backbone, sites holding fixed estimated orbits, troposphere, EOPs and IGS station coordinates. This leads to estimated CORS coordinates in a no net rotation (NNR) frame.

    3. The last step is the alignment of the estimated coordinates with ITRF2014 and velocity estimation. This process was done in 15 iterations to achieve rigorous quality control and discontinuity detection.

    Linear velocities for all stations are estimated in the NGS realization of ITRF2014. NGS explains how this was implemented in the section titled “The velocity field relative to ITRF2014” (see box titled “Section Describing the Velocity Field Relative to ITRF2014”). The website provides figures that depict the horizontal and vertical velocities used in the processing.

    The following are a few highlights from the section:

    • Unless an earthquake or a post seismic adjustment occurred, the velocities of a station in between discontinuities are constrained to have the same value.
    • Stations that experience earthquakes, post seismic adjustment and in a few cases, non-uniform vertical motion, are allowed to have different velocities in between events as dictated by the data.
    • The webpage provides figures that depict the estimated horizontal and vertical CORS velocities.
    Section describing the velocity field relative to ITRF2014. (Screenshot: NGS)
    Section describing the velocity field relative to ITRF2014. (Screenshot: NGS)

    What users usually want to know is how much the coordinates have changed and what it means to their surveying activities. The section titled “Main Changes Compared to Previous Reference Frames” provides information and plots that depict the changes of coordinates.

    Section on changes in coordinates. (Screenshot: NGS)
    Section on changes in coordinates. (Screenshot: NGS)

    This section provides NAD83 (MYCS2) coordinate values minus NAD83 (MYCS1) coordinate values.

    The following are a few highlights from the section:

    • The ITRF2014 coordinates of all computed CORS coordinates from MYCS2 processing are converted to NAD83 (2011) using HTDP.
    • The resulting NAD83 (2011) coordinates are then compared to those obtained from MYCS1 at all common sites.
    • The coordinate differences are compared at epoch 2010.0 (MYCS2 – MYCS1).
    • The differences are less than 5 mm in most areas with some exceptions.
    • The largest differences are seen in southern Alaska.
    • Other visible changes are seen in areas of significant and real subsidence and in places where the time series are too short, such as in Iowa where almost all time series are three years long.
    • Vertical coordinates (ellipsoidal heights) are compared using the same criteria.
    • The stations with the HTDP estimated velocities from MYCS1 (no vertical velocities) show the largest differences. In addition, non-secular subsidence areas also show larger differences.

    By clicking on the plots, the user is directed to a larger figure that is easier to interpret. (See boxes titled “NAD83 (MYCS2) – NAD83 (MYCS1) Horizontal Position Differences” and “NAD83 (MYCS2) – NAD83 (MYCS1) Vertical Position Differences.”)

    NAD83 (MYCS2) - NAD83 (MYCS1) Horizontal Position Differences. (Screenshot: NGS)
    NAD83 (MYCS2) – NAD83 (MYCS1) Horizontal Position Differences. (Screenshot: NGS)
    NAD83 (MYCS2) - NAD83 (MYCS1) Vertical Position Differences. (Screenshot: NGS)
    NAD83 (MYCS2) – NAD83 (MYCS1) Vertical Position Differences. (Screenshot: NGS)

    NGS has done a tremendous job of explaining the MYCS2 process and results. As the results indicate, most differences between the MYCS1 and MYCS2 are small. Saying that, I would encourage all users to look at the NGS Beta webpages and obtain an understanding of the MYCS2 process and results. Users should also use the beta products and compare their results to the current production products to evaluate the CORS beta coordinates and velocities in their region of interest.

    Notice announcing beta version of Geoid18 on NGS homepage. (Screenshot: NGS)
    Notice announcing beta version of Geoid18 on NGS homepage. (Screenshot: NGS)

    It should also be noted that in late February, NGS released a beta version of the latest hybrid geoid model, Geoid18. This model can be accessed here; the site provides an opportunity for users to compute a beta Geoid18 value for a particular station.

    Excerpt from beta Geoid18 website. (Screenshot: NGS)
    Excerpt from beta Geoid18 website. (Screenshot: NGS)

    I would encourage all users to obtain an understanding of the new hybrid model. Once again, it should be noted that this model is a beta model for users to test their workflows and should never be used for official or production work.

    My next column will discuss the beta hybrid Geoid18 model, and the differences between the beta model and the official hybrid geoid model, Geoid12B.