St Michael Bay

**2023 June 6, I am still working to update this use case with the new code.”


at01 is the only tide gauge in this region of Alaska.


Station Name: at01

Location: St Michael, Alaska

Archive: UNAVCO

Ellipsoidal Coordinates:

  • Latitude: 63.4840

  • Longitude: -162.0064

  • Height: 21.7 m

Station Page at UNAVCO

Station Page at Nevada Geodetic Laboratory

Google Map Link

Data Summary

at01 observes all constellation signals and high-rate data are available. There is good visibility over a large azimuthal region. The monument is about a meter taller than normal geodetic monuments to improve visibility of the sea surface. From the geoid app you can see the general location of the antenna with respect to the coast. You will also note that it is at ~12 meters above sea level.

Use the webapp to set a possible mask. In this example we used the option for elevation angles between 5 and 12 degrees.

Take a Quick Look at the Data

First examine the spectral characteristics of the SNR data. (For details on quickLook output.). Make a SNR file:

rinex2snr at01 2020 109 -archive unavco

This will generate a SNR file at a sampling rate of 15 seconds. One second data are available if you would like to use them (-rate high).

If you use the default settings with quickLook you will mistakenly think it is not a useful site. Nearly every single retrieval is set as bad (i.e. it is gray rather than blue):

quickLook at01 2020 109


Remember that the site is 12 meters above sea level and the default restricts the reflector height region to < 8 meters. Try again using a reflector height region that includes the water surface (h1 = 8 and h2 = 15) and better elevation angles 5 and 13 degrees:

quickLook at01 2020 109 -e1 5 -e2 13 -h1 8 -h2 15


Now you see good retrievals at azimuths sweeping from true north to about 220 degrees.
You will also see strong retrievals in the Lomb Scargle periodograms:


This site has modern GPS signals, Galileo signal, and Glonass signals. Here are some sample results for L2C:

quickLook at01 2020 109 -e1 5 -e2 13 -h1 8 -h2 15 -fr 20



quickLook at01 2020 109 -e1 5 -e2 13 -h1 8 -h2 15 -fr 101


and Galileo:

quickLook at01 2020 109 -e1 5 -e2 13 -h1 8 -h2 15 -fr 205


Analyze the Data

Next we analyze data for two months in the fall of 2020. First make the SNR files:

rinex2snr at01 2020 230 -archive unavco -doy_end 290

Now set up the analysis instructions (assume database receiver coordinates are correct) using our new utility (extending h2 a bit):

gnssir_input at01 -h1 8 -h2 17 -e1 5 -e2 13 -ampl 4 -frlist 1 20 5 101 102 201 205 206 207 -azlist2 20 220

Next estimate reflector height (RH) for the two month time period

gnssir at01 2020 230 -doy_end 290

We have written some code to help you look at these subdaily files - it is a work in progress, but you can certainly give it a try. (Note: These figures were generated with an older version of the code and should be updated.) Concatenates the RH files for this period:

subdaily at01 2020 -doy1 230 -doy2 290

You can see that there are a very large number of RH retrievals per day:


This preliminary version of the code removes outliers and makes an effort to compute the RH dot correction if rhdot is set to true. It uses a cubic spline to fit the RH data which allows a first order estimate for the surface rate of change. That, along with geometrical information as to the elevation angle rate of change, is used to make the RH dot correction (for more information).


This term is very important for sites with large tidal ranges, but is of less importance at sites like at01. Nevertheless, you can see here that it does help a bit:

RMS no RHdot correction (m)  0.084
RMS w/ RHdot correction (m)  0.071

Please note: the corrected RH value is not written to column 3. Please look into the file.

Finally, we attempt to remove the inter-frequency biases and fit a new spline, yielding a RMS agreement of about 5 cm with respect to the spline:


Again: the corrected RH values are written to new columns in the output files, not to column 3. Please look into the file.

I would like to include Simon Williams’ RH retrieval/tidal estimation code in this package. Simon has been kind enough to make the Matlab code open source. If someone is willing to convert it to python, that would be fabulous.