ConceptioArchiveUSGS Publications
USGS Publicationspublic full text

The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Jacob T. Murchek, Benjamin J. Drenth, James J. Reitman, Eric D. Anderson, Benjamin P. Magnin, James M. DeGraff · U.S. Geological Survey
USGS Publications · Papers · License: Public Domain
Open Source ↗Direct PDF ↓
datadigitalelevationfilereportgravitylidarmodelsusgsopen
usgs, geological survey, united states, science, government publication, usgs open-file report, 2025, 1019, feasibility, using, lidar, derived, digital, elevation, models, gravity, data, reduction

Mineral Resources Program

The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Open-File Report 2025–1019 Version 1.1. July 2025

U.S. Department of the Interior U.S. Geological Survey

The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction By Jacob T. Murchek, Benjamin J. Drenth, James J. Reitman, Eric D. Anderson, Benjamin P. Magnin, and James M. DeGraff

Mineral Resources Program

Open-File Report 2025–1019 Version 1.1. July 2025

U.S. Department of the Interior U.S. Geological Survey

U.S. Geological Survey, Reston, Virginia: 2025

First release: 2025 Revised: July 2025 (ver. 1.1)

For more information on the USGS—the Federal source for science about the Earth, its natural and living resources, natural hazards, and the environment—visit https://www.usgs.gov or call 1–888–392–8545. For an overview of USGS information products, including maps, imagery, and publications, visit https://store.usgs.gov/ or contact the store at 1–888–275–8747. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Although this information product, for the most part, is in the public domain, it also may contain copyrighted materials as noted in the text. Permission to reproduce copyrighted items must be secured from the copyright owner. Suggested citation: Murchek, J.T., Drenth, B.J., Reitman, J.J., Anderson, E.D., Magnin, B.P., and DeGraff, J.M., 2025, The feasibility of using lidar-derived digital elevation models for gravity data reduction (ver. 1.1, July 2025): U.S. Geological Survey Open-File Report 2025–1019, 33 p., https://doi.org/10.3133/ofr20251019. ISSN 2331-1258 (online)

iii

Acknowledgments The authors would like to thank Dr. Snehamoy Chatterjee of the Michigan Technological University for his helpful discussion of statistical data analysis; Joshua Nimetz of the U.S. Geological Survey (USGS) for his helpful discussion on lidar data acquisition, processing, and formats; and Anji Shah and Phil Brown of the USGS for their thorough reviews of our manuscript. Val Chandler of the Minnesota Geological Survey was an early proponent of using lidar-derived digital elevation models for gravity work and challenged us to undertake this study.

v

Contents Acknowledgments����������������������������������������������������������������������������������������������������������������������������������������iii Abstract�����������������������������������������������������������������������������������������������������������������������������������������������������������1 Introduction����������������������������������������������������������������������������������������������������������������������������������������������������1 Gravity Data Acquisition and Reduction����������������������������������������������������������������������������������������������������1 Lidar Acquisition and Processing���������������������������������������������������������������������������������������������������������������2 Study Design���������������������������������������������������������������������������������������������������������������������������������������������������3 Results�������������������������������������������������������������������������������������������������������������������������������������������������������������9 Discussion�����������������������������������������������������������������������������������������������������������������������������������������������������12 Recommended Field Practices���������������������������������������������������������������������������������������������������������12 Conclusion����������������������������������������������������������������������������������������������������������������������������������������������������12 Tables 3–7.......................................................................................................................................................13 References Cited�����������������������������������������������������������������������������������������������������������������������������������������32

Figures 1.

2.

3.

4.

5.

A schematic of light detection and ranging data acquisition showing an aircraft scanning the Earth’s surface with lidar and revealing lidar topography information while simultaneously receiving global navigation satellite system and inertial measurement unit data�������������������������������������������������������������������������������������������2 Three maps (A–C) showing the 237 gravity stations in the Upper Peninsula of Michigan for which location information was collected as part of regional gravity surveys������������������������������������������������������������������������������������������������������������������������������4 Two maps (A, B) showing the 198 gravity stations in the western United States (Wyoming and Colorado) for which location information was collected as part of regional gravity surveys��������������������������������������������������������������������������������������������7 A flowchart showing how global navigational satellite system positional data and light detection and ranging data were acquired, processed, and prepared and then how the resulting data for both sets were selected, extracted, and interpreted�������������������������������������������������������������������������������������������������������������������������������������9 Two bar graphs (A, B) showing elevation distributions�������������������������������������������������������10

Tables 1.

2. 3. 4. 5.

The minimum number of checkpoints recommended per square kilometer of a light detection and ranging survey from the American Society for Photogrammetry and Remote Sensing (2023)��������������������������������������������������������������������������3 The requirements for different quality-level light detection and ranging data; quality level 2 data are used in this report�������������������������������������������������������������������������������3 Location information for the 94 gravity stations located at west-central Upper Peninsula of Michigan from Drenth and Others (2024)��������������������������������������������������������13 Location information for the 87 gravity stations located on the central Upper Peninsula of Michigan from Drenth and others (2024)���������������������������������������������������������17 Location information for the 139 gravity stations located on the Medicine Bow Mountains of Wyoming from Brown and others (2025)��������������������������������������������������������21

vi

6.

Location information for the 59 gravity stations located on the Wet Mountains of Colorado from Magnin and Anderson (2024)���������������������������������������������������������������������26 7. Location information for the 56 gravity stations located on the Keweenaw Peninsula of Michigan from Murchek and others (2025)�����������������������������������������������������29 8. A statistical breakdown of each study area in this project that shows the average elevation difference observed and the number of points with a difference greater than 1 meter�����������������������������������������������������������������������������������������������11

Conversion Factors International System of Units to U.S. customary units

Multiply

By

To obtain Length

centimeter (cm)

0.3937

inch (in.)

meter (m)

3.281

foot (ft)

kilometer (km)

0.6214

mile (mi)

meter (m)

1.094

yard (yd)

square meter (m2)

10.76

square foot (ft2)

square kilometer (km2)

0.3861

square mile (mi2)

Area

Datums Vertical coordinate information is referenced to the North American Vertical Datum of 1988 (NAVD 88). Horizontal coordinate information is referenced to the North American Datum of 1983 (NAD 83). Elevation, as used in this report, refers to distance above the vertical datum.

vii

Abbreviations 3DEP

3D Elevation Program

ASPRS

American Society for Photogrammetry and Remote Sensing

cm

centimeter

CO–WM

Wet Mountains of Colorado

DEM

digital elevation model

dGNSS

differential global navigational satellite system

GLONASS Globalnaya Navigazionnaya Sputnikovaya Sistema [Global Navigation Satellite System] GNSS

global navigation satellite system

GPS

Global Positioning System

IDW

inverse distance weighted

IMU

inertial measurement unit

km

kilometer

km2

square kilometer

lidar

light detection and ranging

m

meter

mGal

milligal

MI–HR

west-central Upper Peninsula of Michigan

MI–KP0

Keweenaw Peninsula of Michigan

MI–TM

central Upper Peninsula of Michigan

NAD 83

North American Datum of 1983

NAVD 88

North American Vertical Datum of 1988

NGS

National Geodetic Survey

NOAA

National Oceanic and Atmospheric Administration

NPD

nominal pulse density

NPS

nominal pulse spacing

QL

quality level

RMSEz

root mean square error in the vertical (z) direction

TIN

triangulated irregular network

UP

upper peninsula

USGS

U.S. Geological Survey

WY–MB

Medicine Bow Mountains of Wyoming

The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction By Jacob T. Murchek,1,2 Benjamin J. Drenth,1 James J. Reitman,1 Eric D. Anderson,1 Benjamin P. Magnin,1 and James M. DeGraff2

Abstract Gravity data require submeter elevation accuracy for data processing, and differential global navigation satellite system (dGNSS) equipment is commonly used to acquire three-dimensional positional data to achieve such accuracy. However, lidar (light detection and ranging) data are commonly used to develop digital elevation models (DEMs) of Earth’s surface. Therefore, using elevations from lidar-derived DEMs for gravity-data acquisition and reduction may improve field efficiency and reduce cost. This study examines the feasibility of using DEMs for gravity-data reduction by comparing dGNSS elevation data from 435 gravity stations in Michigan, Wyoming, and Colorado with their respective DEM elevations. The results show that the average difference between DEM and dGNSS elevations is 13 centimeters (cm) and that 93 percent of those differences are less than 50 cm, even in areas with steep terrain. Because an elevation discrepancy of 50 cm corresponds to an error of roughly 0.1 milligals (mGal) in the simple Bouguer gravity anomaly, the results suggest that lidar-derived DEMs are a viable source for acquiring the elevation data needed to process gravity data, thus improving both the cost and efficiency of data collection for regional surveys where an accuracy of less than 1.0 mGal is desired.

Introduction The gravity geophysical method is a useful, often essential, tool for mapping subsurface geology. Gravity data reduction (processing) relies heavily on accurate and precise elevation measurements, which are often achieved using high-precision (submeter) differential global navigation satellite system (dGNSS) equipment that can determine gravity-station elevations within less than 1 meter. This sort of “surveying grade” dGNSS equipment typically costs tens of thousands of dollars and requires the daily setup of a local base station that must be secured. If the requirement for expensive dGNSS equipment and the corresponding need for a local 1U.S. Geological Survey. 2Michigan Technological University.

base station can be eliminated, gravity data acquisition would become less expensive and more efficient. This report examines the feasibility of using elevations taken from lidar-derived DEMs for use in the reduction of gravity data as a possible replacement for measurements that require the use of expensive dGNSS equipment. The assessment involves comparing gravity-station elevations acquired using high-precision dGNSS equipment with elevations sampled from lidarderived DEMs at the same locations. Nondifferentially corrected horizontal coordinates provided by the dGNSS measurements were used to simulate the effect of data acquisition without using a GNSS base station.

Gravity Data Acquisition and Reduction Gravity data are acquired using gravimeters, which are usually spring-type balances that measure relative gravity by the change of strain (length) on an internal spring (for example, Hinze and others [2012]). A change in gravity causes the displacement of a test mass within these devices that can be nullified by adjusting the spring length to compensate for the test mass displacement. The spring length needed to nullify the balance is then used to calculate the gravity at a specific location. Gravity measurements are taken relative to a base station, where the absolute value of the gravity field intensity is typically known to within 0.1 milligal (mGal). Measurements of variations in the Earth’s gravity field (anomalies) are produced by lateral variations of density within the subsurface (as in Hinze and others [2012]). The measured gravity field is affected by multiple factors, including time-varying factors (such as tides and instrument drift), latitude, elevation, the terrain surrounding the measurement locations (stations), and the density of the Earth. To isolate anomalies produced by local density variations relatable to geology, a series of corresponding corrections must be applied to the observed gravity (for example, Longman [1959] and Blakely [1995]). Accurate elevation measurements are crucial for calculating the free-air and Bouguer gravity anomalies. So, the quality of a gravity survey is heavily dependent on precise elevation control, as an elevation error of 50 centimeters (cm) results in a 0.1 mGal error in the simple Bouguer anomaly.

2   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Lidar Acquisition and Processing Lidar is a remote-sensing method that uses light as a pulsed laser to measure distances to the Earth from an airborne platform (for example, National Oceanic and Atmospheric Administration [NOAA], [2023]). Data are collected from an aircraft with three main components: a laser scanner unit, a global navigation satellite system (GNSS) unit, and an inertial measurement unit (IMU) (Habib and others, 2005; Hollaus and others, 2005; Reutebuch and others, 2005; Webster and Dias, 2006; Pfeifer and Briese, 2007; Liu, 2008). The laser scanner emits pulses around 1,000 nanometers at a near-infrared wavelength (for example, Elaksher [2016]) and contains a receiver that detects the time it takes for the pulsed laser to reach the Earth and return (fig. 1). These reflections are recorded as individual points that define a point cloud, which can be processed to represent locations on the surface of

the Earth (Sugarbaker and others, 2014). The recorded reflections occur from vegetation, the ground surface, and even human-made objects (see Barber and Shortrudge [2004] and Stoker and others [2006]). Therefore, a critical step in lidar processing is ensuring that unwanted artifacts (nonground) are extracted from the data before DEM construction (Liu, 2008). The GNSS unit records the aircraft’s trajectory, and the IMU measures the aircraft’s altitude; both directly influence the accuracy of the lidar points (Webster and Dias, 2006). Analysts supplement and validate the lidar data with groundcontrol checkpoints having known horizontal and vertical positions to ensure horizontal and vertical accuracy. Table 1 gives the minimum number of checkpoints recommended by the American Society for Photogrammetry and Remote Sensing (ASPRS) based on the area of the DEM being constructed; however, these checkpoints may be altered based on the desired quality level (QL) of the DEM (ASPRS, 2023).

Z GNSS

Y

X

Pitch

IMU Roll

Laser scanner

Yaw

θ = full scan angle Single laser shot GNSS base station

Figure 1. A schematic of light detection and ranging (lidar) data acquisition, modified from National Oceanic and Atmospheric Administration (2012), showing an aircraft scanning the Earth’s surface with lidar and revealing lidar topography information while simultaneously receiving global navigation satellite system (GNSS) and inertial measurement unit data. The full scan angle (θ) from the aircraft and a single laser shot within that scan are shown. A GNSS base station is shown on Earth’s surface.

Lidar Acquisition and Processing   3 The quality level (QL) of a lidar survey is determined by the nominal pulse spacing (NPS) and the vertical positional accuracy, as specified in table 2 (U. S. Geological Survey [USGS], undated.). QL2 lidar-derived DEMs were used for this analysis. Table 1. The minimum number of checkpoints recommended per square kilometer (km2) of a light detection and ranging survey from the American Society for Photogrammetry and Remote Sensing (2023). [≤, less than or equal to]

Area (km2)

Number of checkpoints

≤ 500

30

501–750

35

751–1,000

40

1,001–1,250

45

1,251–1,500

50

1,501–1,750

55

1,751–2,000

60

2,001–2,250

65

2,251–2,500

70

Study Design

Table 2. The requirements for different quality level (QL) light detection and ranging (lidar) data; quality level 2 (QL2) data are used in this report. [This table is modified from the U.S. Geological Survey (USGS) 3D Elevation Program (3DEP) (USGS, undated). cm, centimeter; DEM, digital elevation model; m, meter; NPD, nominal pulse density; NPS, nominal pulse spacing; pts per m2, points per square meter; RMSEz, root mean square error in the vertical (z) direction, QL; quality level; ≤, less than or equal to; ≥, greater than or equal to]

Quality Data level1 source

algorithm to create a surface model from the point-cloud data. The height of each point is compared with its neighbor within the triangle, and points having a significantly higher height than their neighbors are removed (Peucker and others, 1976). After extracting nonground points from the lidar data, a DEM is interpolated based on the remaining ground points. Interpolation is predicting values at an unsampled location using the measured values nearby (for example, Burrough and McDonnell [1998]). DEM interpolation typically uses inverse distance weighted (IDW), spline-based, or geostatistical methods (such as kriging) (Liu, 2008). Once the interpolation method is chosen, a DEM grid is constructed based on the QL of the data.

Vertical accuracy RMSEz (cm)

Nominal pulse spacing (m)

Nominal pulse density (pts per m2)

DEM cell size (m)

QL0

Lidar

5

≤ 0.35

≥8

0.5

QL1

Lidar

10

≤ 0.35

≥8

0.5

QL2

Lidar

10

≤ 0.71

≥2

1

QL3

Lidar

20

≤ 1.41

≥ 0.5

2

1Quality levels are explained in U.S. Geological Survey (undated).

Several methods can be used for extracting nonground points from the raw lidar point-cloud, such as triangulated irregular network (TIN) filtering, slope-based filtering, mathematical morphological filtering, interpolation-based filtering, and machine-learning-based filtering (Cai and others, 2020). The USGS 3D Elevation Program (3DEP) typically requires TIN filtering to assess both the vegetated and nonvegetated vertical accuracy of lidar data (USGS, 2024). TIN filtering uses a triangulation

To study the feasibility of using lidar-derived elevations for gravity-data reduction, we acquired dGNSS elevation data for 435 gravity stations in the Upper Peninsula of Michigan (237 stations, fig. 2A–C), the Wet Mountains of Colorado (59 stations, fig. 3A), and the Medicine Bow Mountains of Wyoming (139 stations, fig. 3B). For each station, Leica Viva GS16 GNSS equipment or a Trimble Geo7x handheld GNSS receiver was used to determine differentially corrected locations (latitude and longitude) and elevations, typically accurate to within 10 cm. To compute differential corrections for the Leica Viva GS16, a GNSS base station located within 25 kilometers (km) of the stations was used. For the Trimble Geo7x, the NOAA Continuously Operating Reference Stations located within 65 km of the gravity stations were used to correct the data differentially. Additionally, GNSS-derived locations were recorded without applying differential corrections to simulate data acquisition conditions without using a GNSS base station. For the lidar comparisons, 1-meter DEM tiles (QL2) were obtained using the USGS LidarExplorer (USGS, 2022; USGS, 2023), the Michigan Technological University Geospatial Research Facility DEM downloader (Sanborn Map Company, Inc., 2020), or the Colorado Hazard Mapping & Risk MAP Portal (Merrick & Company, 2016; Quantum Spatial, Inc., 2020) and were “mosaiced” together in the ArcGIS Pro software package to encompass the entirety of each gravity survey.3 The nondifferentially corrected gravity station coordinates (North American Datum of 1983) were added to ArcGIS Pro, and the DEM was sampled at each station for comparison with the differentially corrected elevation provided by the dGNSS equipment. The flowchart shown in figure 4 describes each data acquisition, processing, and interpretation step. However, it should be noted that the lidar data were acquired and processed before this study. 3Datasets are “mosaiced” or merged using the “Mosaic” tool in the ArcGIS Pro software package.

4   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction A

Base map from Esri and its licensors, copyright 2022

Figure 2. Three maps (A–C) showing the 237 gravity stations in the Upper Peninsula of Michigan for which location information was collected as part of regional gravity surveys. The stations overlay the mosaiced 1-meter (m) digital elevation model (DEM) rasters obtained from the U.S. Geological Survey (USGS) 3D Elevation Program (3DEP) light detection and ranging (lidar) data in the National Map Downloader (USGS, 2022) or the Michigan Technological University Geospatial Research Facility DEM tool (Sanborn Map Company, Inc., 2020). (A) Map showing the 56 gravity stations in the Keweenaw Peninsula of Michigan (MI–KP0) where location data were collected with a Global Navigational Satellite System (GNSS) receiver. Elevation data are from 0 to 450 m. (B) Map showing the 94 gravity stations in the west-central Upper Peninsula of Michigan (MI–HR) where location information was collected with GNSS equipment. Elevation data are from 400 to 600 m. (C) Map showing the 87 gravity stations in the central Upper Peninsula of Michigan (MI–TM) where location information was collected with GNSS equipment. Elevation data are from 150 to 300 m. The prefixes discussed here—MI–KP0, MI–HR, and MI–TM—represent the location aspect of the gravity station numbers in tables 3, 4, and 7 (at the end of this report). ft, foot.

Lidar Acquisition and Processing   5 B

Figure 2.—Continued

6   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction C

Base map from Esri and its licensors, copyright 2022

Figure 2.—Continued

0

10

20 KILOMETERS

Lidar Acquisition and Processing   7 A

Base map from Esri and its licensors, copyright 2022 0

10

20 KILOMETERS

Figure 3. Two maps (A, B) showing the 198 gravity stations in the western United States (Wyoming and Colorado) for which location information was collected as part of regional gravity surveys. The stations overlay the mosaiced 1-meter (m) digital elevation model (DEM) rasters. Light detection and ranging (lidar) data were obtained from two sources. U.S. Geological Survey (USGS) 3D Elevation Program lidar data were retrieved from the USGS National Map Downloader (USGS, 2024). Lidar data from Merrick & Company (2016) and Quantum Spatial, Inc. (2020) were retrieved from the Colorado Hazard Mapping & Risk MAP Portal (Merrick & Company, 2016). (A) Map showing the 59 gravity stations in the Wet Mountains of Colorado (CO–WM) where location information was collected with Global Navigational Satellite System (GNSS) equipment. Elevation data are from 2,000 to 4,000 m. (B) Map showing the 139 gravity stations in the Medicine Bow Mountains of Wyoming (WY–MB) where location information was collected with GNSS equipment. Elevation data are from 2,000 to 4,000 m. The prefixes discussed here—CO–WM and WY– MB—represent the location aspect of the gravity station numbers in tables 5 and 6 (at the end of this report). ft, foot.

8   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction B

Base map from Esri and its licensors, copyright 2022

Figure 3.—Continued

0

10

20 KILOMETERS

Lidar Acquisition and Processing   9

GNSS Posi�onal Data Data Acquisi�on Eleva�on data collected at gravity sta�on using differen�al GNSS equipment

Data Processing GNSS data post-processed to determine both nondifferen�ally and differen�ally corrected loca�ons

Data Prepara�on Tables of data prepared and loaded as point data (NAD83)

Lidar Data Acquisi�on Airborne lidar point cloud data collected prior to study Data Interpreta�on Data were exported for sta�s�cal analysis

Data Extrac�on The DEMs were sampled at each gravity sta�on loca�on and a column with the DEM eleva�on for each sta�on was created

Data Selec�on Data points selected for each gravity survey. Data outside DEM coverage filtered out

Data Processing Ar�facts removed and bare Earth surface iden�fied. 1meter DEM created

Data Prepara�on 1-meter DEM �les for each study area were downloaded and mosaiced

Figure 4. A flowchart showing how global navigational satellite system (GNSS) positional data and light detection and ranging (lidar) data were acquired, processed, and prepared and then how the resulting data for both sets were selected, extracted, and interpreted. DEM, digital elevation model; NAD 83, North American Datum of 1983.

Results The elevation differences between the dGNSS measurements and DEMs were consistently (93 percent) less than 50 cm, corresponding to a 0.1 mGal error in the simple Bouguer gravity anomaly. Horizontal and vertical positions for all 435 gravity stations are shown in tables 3–7 (at the end of the report) with their respective lidar-derived elevation and the differences between the GNSS and lidar-derived elevations. The distribution of differences between DEM and GNSS elevations is shown in figure 5A, and the distribution of their absolute differences is shown in figure 5B. The mean difference between the GNSS and DEMderived elevations is 13 cm with a standard deviation of 46 cm. Approximately 93 percent (N = 406) of the data

points fall within 1 standard deviation of the mean (minimum –33 cm; maximum 59 cm). The lower limit for the 95 percent confidence interval is 9 cm, and the upper limit is 17 cm. A maximum positive difference of 3.9 meters (m) and a maximum negative difference of –2.1 m occur in the central Upper Peninsula. However, the mean elevation difference in the Upper Peninsula is 8 cm, a measure less than the average in Colorado and Wyoming (20 cm), which contain steep terrain. Overall, 97 percent (N = 423) of the 435 stations exhibit an elevation difference of less than 1 m. The statistical breakdown for each study area is shown in table 8.

10   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction A

200

GNSS Elevations - DEM Elevations Mean = 0.13 Std Dev = 0.46 96% Confidence interval for mean [0.08, 0.17]

46%

Number of measurements

−σ

150 29%

100

15%

50

3%

2% ~1%

~1%

0 −2.5

−2.0

−1.5

−1.0

−0.5

0

µ

0.5

1.0

1.5

2.0

2.5

3.0

3.5

4.0

Elevation differences, in meters

Figure 5. Two bar graphs (A, B) showing elevation distributions. (A) Bar graph that shows the distribution (mean = 0.13, standard deviation = 0.46) of elevation differences (global navigational satellite system minus light detection and ranging [GNSS – lidar]) and indicates that 93 percent of the data fall within 1 standard deviation (σ) of the mean (μ) (13 centimeters [cm]). (B) Bar graph that shows the distribution (mean = 0.22, standard deviation = 0.42) of the absolute value of the elevation differences (|GNSS – lidar|) and a mean difference of 22 cm in the absolute differences observed within the data. %, percent; m, meter.

Lidar Acquisition and Processing   11 GNSS Elevations - DEM Elevations

B

Mean = 0.22 Std Dev = 0.42 95% Confidence interval for mean [0.18, 0.26]

75%

300 +σ

200

100

29%

3% 0

µ

~1%

0.5

1.0

1.5

2.0

2.5

3.0

3.5

4.0

Number of measurements

−σ

0

Absolute elevation differences, in meters

Figure 5.—Continued

Table 8. A statistical breakdown of each study area in this project that shows the average elevation difference observed and the number of points with a difference greater than 1 meter (m). [UP, upper peninsula ; > greater than]

Mean elevation difference (m)

Standard deviation

Minimum (m)

Maximum (m)

No. of points >1 m elevation difference

Medicine Bow Mountains, Wyoming

0.25

0.37

–0.39

3.1

4

Wet Mountains, Colorado

0.11

0.53

–0.93

3.7

1

Keweenaw Peninsula, Michigan

0.26

0.47

–0.07

2.7

2

West-central UP, Michigan

–0.07

0.17

–0.66

0.54

0

Central UP, Michigan

0.11

0.62

–2.1

3.9

5

Total UP, Michigan (less-steep terrain)

0.08

0.44

–2.1

3.9

7

Total Wyoming and Colorado (steep terrain)

0.20

0.42

–0.93

3.7

5

Location

Regional totals

12   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Discussion Of the 435 gravity stations, 406 (approximately 93 percent) lidar-derived elevations fall within ±50 cm of the dGNSS elevations, meaning that their use would produce an error of less than 0.1 mGal in the simple Bouguer anomaly. Regional gravity surveys commonly focus on mapping anomalies greater than 1 mGal, so an error of 0.1 mGal is not a cause for concern in most cases. Moreover, most absolute base (reference) stations used for gravity-data reduction have an uncertainty of ±0.1 mGal (for example, Morelli and others [1972]). Therefore, this study shows that lidar-derived elevations are sufficiently accurate for regional gravity-data reduction, and, in most cases, their use produces errors comparable to the inherent absolute accuracy of standard regional gravity surveys. For detailed gravity surveys focused on smaller anomalies—those with an accuracy of less than 0.1 mGal—dGNSS equipment may be necessary. Only 12 gravity stations (less than 3 percent) had elevation differences greater than 1 m between dGNSS and lidar-derived elevations. Seven were in the Upper Peninsula of Michigan, four were in the Medicine Bow Mountains of Wyoming, and one was in the Wet Mountains of Colorado. To identify the sources of these discrepancies, we examined the elevation differences between the lidar-derived elevations and dGNSS elevations using the differentially corrected horizontal coordinates in the Medicine Bow Mountains. The average elevation difference between the measurements was roughly 1-cm less than when the nondifferentially corrected coordinates were used, and the same four stations showed an elevation difference of greater than 1 m. This observation indicates that the four largest elevation discrepancies are not likely related to the lack of differential processing for the horizontal coordinates. Instead, these discrepancies may result from dense vegetation causing a multipath in the GNSS signal, a locally steep elevation gradient, or an unknown error in the dGNSS-derived elevations.

Recommended Field Practices The use of several recommendations may optimize positional data acquisition for gravity data processing and yield acceptable uncertainty in simple Bouguer gravity values. GNSS equipment (for example, a handheld, non-dGNSS unit) with the ability to detect both Global Positioning System (GPS) and Globalnaya Navigatsionnaya Sputnikovaya Sistema (GLONASS)4 satellites is recommended for obtaining the best possible horizontal position data, especially at higher latitudes where GPS readings are less accurate. Horizontal positional errors may result in the gravity station being mislocated by more than several meters, thus producing elevation inaccuracies. For this 4The “Global Navigation Satellite System” is a Russian satellite-based navigation system that provides position information.

reason, placing gravity stations in areas more than 5 m away from steep topographic slopes (such as a locally flat surface) is recommended. The reduction or elimination of multipaths is also necessary for optimizing positional accuracy during data acquisition. A multipath can result from the primary GNSS signal reflecting off buildings, vegetation, or mountains, or from atmospheric scattering. A reduction of this effect can be achieved by placing the GNSS antenna in an elevated position, such as on the roof of a vehicle or a range pole and away from areas with dense vegetation. Placing stations on small platforms above the ground’s surface can also cause inaccuracies in the lidar-derived elevation and should be avoided.

Conclusion This study evaluated whether lidar-derived DEMs could serve as elevation control for gravity surveys to reduce cost and improve field efficiency. Elevation data from 435 gravity stations in the Upper Peninsula of Michigan, the Medicine Bow Mountains of Wyoming, and the Wet Mountains of Colorado were compared with elevations from a 1-m horizontal-resolution DEM for each station. An average difference of 13 cm between dGNSS and DEM elevations was observed, with a standard deviation of 46 cm. Over 93 percent of the DEM data were within 50 cm of GNSS data, corresponding to a simple Bouguer gravity error of approximately 0.1 mGal. These results indicate that lidar-derived DEMs provide acceptable elevation control for gravity data reduction, particularly for regional gravity surveys where anomalies of interest are usually greater than 1 mGal in amplitude. In detailed surveys requiring greater accuracy and precision, other means of determining elevations may be required. Twelve stations exhibited elevation differences between dGNSS and DEM data greater than 1 m. The source of this relatively rare discrepancy is unknown but may be related to dense vegetation cover, locally steep elevation gradients, or an unknown error in the dGNSS elevations. This study suggests that using QL2 lidar-derived DEMs for gravity-data reduction is acceptable for most regional surveys.

Tables 3–7 Table 3. Location information for the 94 gravity stations located at west-central Upper Peninsula of Michigan (MI–HR) from Drenth and Others (2024). [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

Gravity station

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

46.42494

–88.16299

46.42493

–88.16300

523.87

523.74

0.13

0.13

46.43683

–88.16355

46.43681

–88.16355

538.11

538.10

0.01

0.01

MI–HR211

46.44739

–88.16216

46.44738

–88.16216

531.44

531.10

0.35

0.35

MI–HR212

46.46448

–88.17829

46.46448

–88.17830

544.78

544.73

0.04

0.04

MI–HR213

46.43243

–88.18091

46.43242

–88.18092

525.93

525.81

0.12

0.12

MI–HR214

46.43118

–88.19405

46.43117

–88.19406

514.27

514.27

0.00

0.00

MI–HR215

46.42438

–88.21213

46.42439

–88.21215

515.77

515.53

0.24

0.24

MI–HR216

46.42267

–88.22513

46.42266

–88.22513

536.88

536.80

0.09

0.09

MI–HR217

46.41882

–88.24124

46.41881

–88.24127

534.57

534.72

–0.15

0.15

MI–HR218

46.43243

–88.22989

46.43243

–88.22991

543.09

543.01

0.08

0.08

MI–HR219

46.44121

–88.23626

46.44121

–88.23627

540.62

540.19

0.43

0.43

MI–HR220

46.45163

–88.23310

46.45164

–88.23311

534.12

534.74

–0.62

0.62

MI–HR221

46.45329

–88.21923

46.45330

–88.21922

525.97

525.43

0.54

0.54

MI–HR222

46.43917

–88.21494

46.43918

–88.21494

528.55

528.60

–0.05

0.05

MI–HR223

46.44226

–88.18349

46.44227

–88.18349

529.38

529.38

0.00

0.00

MI–HR239

46.44133

–88.08456

46.44131

–88.08455

510.37

510.46

–0.10

0.10

MI–HR240

46.44881

–88.07866

46.44879

–88.07865

511.92

511.90

0.02

0.02

MI–HR241

46.46879

–88.08052

46.46877

–88.08051

478.68

478.89

–0.21

0.21

MI–HR242

46.45943

–88.08535

46.45941

–88.08534

472.61

472.67

–0.06

0.06

MI–HR243

46.45848

–88.09835

46.45848

–88.09833

479.39

479.48

–0.09

0.09

MI–HR244

46.46113

–88.10863

46.46113

–88.10862

498.09

498.04

0.06

0.06

MI–HR245

46.46627

–88.11848

46.46626

–88.11848

521.55

521.60

–0.05

0.05

MI–HR246

46.47889

–88.12047

46.47888

–88.12046

525.69

525.75

–0.06

0.06

MI–HR247

46.45086

–88.10303

46.45085

–88.10303

492.13

492.27

–0.14

0.14

MI–HR248

46.44016

–88.11112

46.44014

–88.11112

486.16

486.13

0.02

0.02

Tables 3–7  13

MI–HR209 MI–HR210

[“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–HR249

46.43465

–88.12686

46.43464

–88.12685

508.44

508.33

0.11

0.11

MI–HR250

46.44574

–88.12259

46.44573

–88.12259

508.40

508.40

0.00

0.00

MI–HR251

46.45507

–88.12062

46.45507

–88.12063

512.48

512.41

0.08

0.08

MI–HR252

46.45011

–88.06530

46.45011

–88.06529

471.30

471.16

0.14

0.14

MI–HR253

46.44357

–88.05018

46.44357

–88.05017

465.11

465.32

–0.21

0.21

MI–HR254

46.42801

–88.05931

46.42801

–88.05930

497.72

497.72

0.00

0.00

MI–HR255

46.43435

–88.03451

46.43434

–88.03450

457.50

457.70

–0.20

0.20

MI–HR256

46.42357

–88.02001

46.42358

–88.02000

457.56

457.47

0.09

0.09

MI–HR257

46.41524

–88.00702

46.41523

–88.00701

461.07

461.20

–0.13

0.13

MI–HR258

46.43268

–88.09584

46.43266

–88.09583

519.07

519.20

–0.13

0.13

MI–HR259

46.42517

–88.10726

46.42517

–88.10726

509.49

509.43

0.06

0.06

MI–HR260

46.41580

–88.11487

46.41580

–88.11487

511.05

511.11

–0.07

0.07

MI–HR261

46.40543

–88.11949

46.40543

–88.11949

503.56

503.58

–0.02

0.02

MI–HR262

46.40569

–88.11962

46.40569

–88.11961

503.74

503.75

0.00

0.00

MI–HR263

46.39486

–88.11466

46.39486

–88.11466

490.43

490.56

–0.13

0.13

MI–HR265

46.42661

–88.09738

46.42660

–88.09738

524.72

524.96

–0.24

0.24

MI–HR266

46.41892

–88.08438

46.41892

–88.08436

508.71

508.78

–0.07

0.07

MI–HR267

46.43559

–88.05164

46.43558

–88.05165

499.46

499.57

–0.11

0.11

MI–HR268

46.41399

–88.02718

46.41398

–88.02718

475.55

475.52

0.04

0.04

MI–HR269

46.41310

–88.01249

46.41309

–88.01248

472.11

472.25

–0.13

0.13

MI–HR270

46.41278

–87.99475

46.41277

–87.99474

457.07

457.03

0.04

0.04

MI–HR271

46.40503

–87.98644

46.40502

–87.98643

456.36

456.40

–0.04

0.04

MI–HR272

46.39359

–87.99155

46.39358

–87.99154

462.14

462.16

–0.01

0.01

MI–HR273

46.39923

–88.01748

46.39922

–88.01748

457.80

457.90

–0.10

0.10

MI–HR274

46.38052

–88.02036

46.38051

–88.02037

461.31

461.31

0.00

0.00

MI–HR275

46.38799

–88.00558

46.38799

–88.00558

456.66

456.70

–0.04

0.04

MI–HR276

46.39172

–88.09548

46.39171

–88.09546

505.30

505.39

–0.09

0.09

MI–HR277

46.39956

–88.08134

46.39956

–88.08134

514.50

514.63

–0.12

0.12

Gravity station

14   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Table 3. Location information for the 94 gravity stations located at west-central Upper Peninsula of Michigan (MI–HR) from Drenth and Others (2024).—Continued

Table 3. Location information for the 94 gravity stations located at west-central Upper Peninsula of Michigan (MI–HR) from Drenth and Others (2024).—Continued [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–HR278

46.39919

–88.06436

46.39920

–88.06434

498.56

498.47

0.09

0.09

MI–HR279

46.40759

–88.05725

46.40758

–88.05723

488.88

488.88

0.01

0.01

MI–HR280

46.41484

–88.06662

46.41484

–88.06661

493.27

493.31

–0.04

0.04

MI–HR281

46.40704

–88.04444

46.40703

–88.04443

477.41

477.68

–0.27

0.27

MI–HR282

46.38617

–88.04290

46.38617

–88.04290

474.39

474.54

–0.15

0.15

MI–HR283

46.39784

–88.04632

46.39783

–88.04632

476.50

476.69

–0.18

0.18

MI–HR284

46.38457

–88.08025

46.38456

–88.08025

507.73

507.89

–0.16

0.16

MI–HR285

46.37563

–88.07300

46.37561

–88.07301

501.87

502.09

–0.22

0.22

MI–HR286

46.37335

–88.05530

46.37333

–88.05530

482.75

482.83

–0.08

0.08

MI–HR287

46.37545

–88.03577

46.37543

–88.03578

471.54

471.62

–0.08

0.08

MI–HR288

46.37665

–88.08606

46.37664

–88.08606

507.56

507.64

–0.08

0.08

MI–HR289

46.37673

–88.09903

46.37672

–88.09903

495.34

495.37

–0.03

0.03

MI–HR290

46.37274

–88.10606

46.37273

–88.10606

488.82

488.87

–0.05

0.05

MI–HR291

46.37188

–88.11852

46.37187

–88.11853

479.97

480.05

–0.08

0.08

MI–HR301

46.41593

–87.98031

46.41594

–87.98032

482.42

482.41

0.01

0.01

MI–HR302

46.41762

–87.96272

46.41761

–87.96274

474.48

474.78

–0.30

0.30

MI–HR303

46.42835

–87.95810

46.42835

–87.95812

466.70

466.94

–0.24

0.24

MI–HR304

46.42658

–87.93938

46.42657

–87.93940

465.76

465.84

–0.08

0.08

MI–HR305

46.42620

–87.92179

46.42620

–87.92181

462.74

462.75

0.00

0.00

MI–HR306

46.42973

–87.90791

46.42970

–87.90793

461.32

461.42

–0.11

0.11

MI–HR307

46.44103

–87.91421

46.44099

–87.91423

477.03

476.99

0.04

0.04

MI–HR308

46.42258

–87.89031

46.42254

–87.89032

456.78

456.92

–0.14

0.14

MI–HR309

46.42578

–87.87663

46.42575

–87.87665

459.60

459.74

–0.15

0.15

MI–HR310

46.41063

–87.89066

46.41060

–87.89066

457.74

457.91

–0.17

0.17

MI–HR311

46.41775

–87.92782

46.41772

–87.92783

460.52

460.70

–0.18

0.18

MI–HR312

46.40566

–87.92574

46.40565

–87.92576

464.99

465.32

–0.33

0.33

MI–HR313

46.40483

–87.96428

46.40484

–87.96430

465.96

466.39

–0.43

0.43

MI–HR314

46.39421

–87.95238

46.39421

–87.95241

467.12

467.40

–0.27

0.27

Tables 3–7  15

UC LAT (NAD 83)

Gravity station

[“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–HR315

46.38814

–87.93001

46.38815

–87.93003

465.65

465.94

–0.29

0.29

MI–HR316

46.38068

–87.91778

46.38070

–87.91779

459.97

460.14

–0.16

0.16

MI–HR317

46.36793

–87.90407

46.36795

–87.90407

454.28

454.45

–0.17

0.17

MI–HR318

46.35638

–87.89669

46.35640

–87.89667

456.48

457.15

–0.66

0.66

MI–HR319

46.34743

–87.87883

46.34744

–87.87881

451.22

451.53

–0.31

0.31

MI–HR320

46.33948

–87.89772

46.33949

–87.89771

454.00

454.09

–0.09

0.09

MI–HR321

46.34497

–87.91396

46.34498

–87.91394

458.56

458.65

–0.09

0.09

MI–HR322

46.34791

–87.93149

46.34791

–87.93146

452.12

452.42

–0.30

0.30

MI–HR323

46.36046

–87.93067

46.36046

–87.93065

455.94

456.20

–0.26

0.26

MI–HR324

46.37647

–87.93349

46.37646

–87.93346

460.19

460.06

0.14

0.14

MI–HR325

46.40201

–87.99925

46.40202

–87.99923

481.08

481.11

–0.03

0.03

MI–HR326

46.37824

–87.98663

46.37825

–87.98662

451.77

451.76

0.01

0.01

MI–HR327

46.36282

–87.97087

46.36284

–87.97084

447.92

448.19

–0.27

0.27

Gravity station

16   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Table 3. Location information for the 94 gravity stations located at west-central Upper Peninsula of Michigan (MI–HR) from Drenth and Others (2024).—Continued

Table 4. Location information for the 87 gravity stations located on the central Upper Peninsula of Michigan (MI–TM) from Drenth and others (2024). [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–TM001

46.05761

–86.76523

46.05760

–86.76523

228.46

228.55

–0.10

0.10

MI–TM002

46.05005

–86.77984

46.05004

–86.77984

237.26

237.29

–0.04

0.04

MI–TM003

46.04471

–86.79793

46.04470

–86.79792

235.38

235.29

0.09

0.09

MI–TM004

46.03561

–86.82491

46.03561

–86.82492

231.49

231.54

–0.05

0.05

MI–TM005

46.05917

–86.82250

46.05919

–86.82249

231.44

233.53

–2.09

2.09

MI–TM006

46.10298

–86.76241

46.10298

–86.76243

240.67

240.74

–0.07

0.07

MI–TM007

46.12079

–86.75999

46.12079

–86.76001

250.59

250.52

0.07

0.07

MI–TM008

46.13699

–86.75775

46.13700

–86.75778

243.03

243.05

–0.01

0.01

MI–TM009

46.14718

–86.78123

46.14719

–86.78124

240.63

240.60

0.02

0.02

MI–TM010

46.14721

–86.80917

46.14722

–86.80918

242.80

242.83

–0.02

0.02

MI–TM011

46.15373

–86.83689

46.15373

–86.83691

237.25

237.29

–0.04

0.04

MI–TM012

46.13682

–86.83782

46.13683

–86.83782

248.12

247.11

1.01

1.01

MI–TM013

46.11941

–86.83631

46.11941

–86.83633

243.01

243.11

–0.10

0.10

MI–TM014

46.10464

–86.83439

46.10466

–86.83442

228.64

228.71

–0.07

0.07

MI–TM015

46.09690

–86.86844

46.09690

–86.86845

214.65

214.68

–0.04

0.04

MI–TM016

46.08374

–86.84393

46.08374

–86.84394

231.24

231.25

–0.02

0.02

MI–TM017

46.06476

–86.85307

46.06475

–86.85306

232.13

231.87

0.26

0.26

MI–TM018

46.04313

–86.85153

46.04314

–86.85154

230.66

230.67

–0.01

0.01

MI–TM019

46.02522

–86.86199

46.02522

–86.86200

224.94

225.03

–0.09

0.09

MI–TM020

46.02190

–86.83746

46.02190

–86.83748

228.33

228.41

–0.08

0.08

MI–TM021

46.09679

–86.77856

46.09679

–86.77857

243.10

242.94

0.16

0.16

MI–TM022

46.07565

–86.77473

46.07564

–86.77473

237.34

237.23

0.11

0.11

MI–TM023

46.04142

–86.75639

46.04141

–86.75639

215.18

215.18

0.00

0.00

MI–TM024

46.02985

–86.73892

46.02984

–86.73893

231.27

231.18

0.09

0.09

MI–TM025

46.03102

–86.71163

46.03102

–86.71163

229.00

228.94

0.06

0.06

MI–TM026

46.04420

–86.71386

46.04421

–86.71387

230.28

230.34

–0.07

0.07

MI–TM027

46.02166

–86.69368

46.02168

–86.69369

230.49

230.47

0.02

0.02

MI–TM028

46.01128

–86.67806

46.01129

–86.67808

230.95

230.94

0.01

0.01

Tables 3–7  17

UC LAT (NAD 83)

Gravity station

[“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–TM029

46.01631

–86.65150

46.01633

–86.65150

229.18

229.18

0.00

0.00

MI–TM030

46.02011

–86.62713

46.02013

–86.62713

238.58

234.70

3.89

3.89

MI–TM031

46.03034

–86.60562

46.03034

–86.60562

232.86

232.85

0.01

0.01

MI–TM032

46.04037

–86.58267

46.04037

–86.58268

233.50

233.46

0.04

0.04

MI–TM033

46.04585

–86.56003

46.04586

–86.56004

233.23

233.06

0.17

0.17

MI–TM034

46.04844

–86.52518

46.04844

–86.52518

231.52

231.51

0.01

0.01

MI–TM035

46.06206

–86.54437

46.06206

–86.54437

234.53

234.51

0.02

0.02

MI–TM036

46.07869

–86.54662

46.07871

–86.54663

237.41

237.31

0.10

0.10

MI–TM037

46.09413

–86.54977

46.09414

–86.54978

239.60

239.56

0.04

0.04

MI–TM038

46.11218

–86.55485

46.11219

–86.55486

244.43

244.38

0.05

0.05

MI–TM039

46.13041

–86.55865

46.13043

–86.55866

247.42

247.39

0.02

0.02

MI–TM040

46.14352

–86.55224

46.14354

–86.55226

255.16

255.02

0.13

0.13

MI–TM041

46.15762

–86.55333

46.15763

–86.55333

253.88

253.87

0.01

0.01

MI–TM042

46.15343

–86.58112

46.15343

–86.58114

242.02

242.12

–0.10

0.10

MI–TM043

46.15213

–86.60946

46.15214

–86.60948

240.03

240.05

–0.02

0.02

MI–TM044

46.07875

–86.74031

46.07874

–86.74032

237.59

236.63

0.95

0.95

MI–TM045

46.07712

–86.69008

46.07711

–86.69009

278.00

278.33

–0.32

0.32

MI–TM046

46.09977

–86.67117

46.09976

–86.67117

249.36

249.29

0.07

0.07

MI–TM047

46.11141

–86.67072

46.11143

–86.67073

239.36

239.24

0.11

0.11

MI–TM048

46.06748

–86.71794

46.06747

–86.71797

231.27

231.28

0.00

0.00

MI–TM049

46.07899

–86.72260

46.07900

–86.72261

235.34

235.34

0.00

0.00

MI–TM050

46.06587

–86.74520

46.06587

–86.74521

223.32

223.30

0.02

0.02

MI–TM051

46.02779

–86.66780

46.02781

–86.66781

242.08

242.04

0.04

0.04

MI–TM052

46.04770

–86.65981

46.04770

–86.65981

249.52

249.51

0.01

0.01

MI–TM053

46.04221

–86.63906

46.04224

–86.63907

240.68

240.27

0.41

0.41

MI–TM054

46.06501

–86.64277

46.06502

–86.64278

245.66

245.55

0.11

0.11

MI–TM055

46.08670

–86.64203

46.08671

–86.64204

248.09

248.09

0.00

0.00

MI–TM056

46.10423

–86.63742

46.10425

–86.63741

253.69

252.71

0.98

0.98

Gravity station

18   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Table 4. Location information for the 87 gravity stations located on the central Upper Peninsula of Michigan (MI–TM) from Drenth and others (2024).—Continued

Table 4. Location information for the 87 gravity stations located on the central Upper Peninsula of Michigan (MI–TM) from Drenth and others (2024).—Continued [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–TM057

46.16162

–86.75468

46.16161

–86.75469

245.85

245.96

–0.12

0.12

MI–TM058

46.17326

–86.72937

46.17326

–86.72939

241.60

241.64

–0.04

0.04

MI–TM059

46.15155

–86.72839

46.15156

–86.72840

241.30

241.38

–0.09

0.09

MI–TM060

46.12253

–86.69312

46.12254

–86.69314

258.60

258.57

0.04

0.04

MI–TM061

46.13349

–86.68783

46.13351

–86.68784

239.11

239.13

–0.02

0.02

MI–TM062

46.14831

–86.70171

46.14833

–86.70173

237.16

237.11

0.05

0.05

MI–TM063

46.16396

–86.69830

46.16397

–86.69831

238.12

238.16

–0.04

0.04

MI–TM064

46.15843

–86.67672

46.15844

–86.67673

255.73

255.73

–0.01

0.01

MI–TM065

46.16060

–86.65028

46.16063

–86.65029

245.13

245.13

0.01

0.01

MI–TM066

46.14994

–86.62531

46.14995

–86.62532

259.00

259.06

–0.06

0.06

MI–TM067

46.12999

–86.60514

46.13001

–86.60516

259.66

258.78

0.88

0.88

MI–TM068

46.13003

–86.57811

46.13004

–86.57811

251.35

251.37

–0.02

0.02

MI–TM069

46.11426

–86.58261

46.11428

–86.58263

248.29

248.27

0.02

0.02

MI–TM070

46.11179

–86.60224

46.11181

–86.60225

252.11

252.11

0.00

0.00

MI–TM071

46.10062

–86.58400

46.10063

–86.58401

245.12

245.04

0.09

0.09

MI–TM072

46.08423

–86.59097

46.08424

–86.59097

240.72

240.70

0.02

0.02

MI–TM073

46.07043

–86.57232

46.07044

–86.57232

236.64

236.62

0.02

0.02

MI–TM074

46.05627

–86.59073

46.05627

–86.59074

235.12

235.84

–0.71

0.71

MI–TM075

46.11936

–86.63511

46.11935

–86.63514

243.36

243.68

–0.32

0.32

MI–TM076

46.14429

–86.65346

46.14430

–86.65349

243.00

242.59

0.41

0.41

MI–TM077

46.14126

–86.63864

46.14127

–86.63868

252.68

251.23

1.45

1.45

MI–TM078

46.12989

–86.80151

46.12989

–86.80151

245.00

245.02

–0.02

0.02

MI–TM079

46.16822

–86.83967

46.16821

–86.83969

239.68

239.76

–0.08

0.08

MI–TM080

46.18263

–86.83808

46.18263

–86.83810

251.87

248.97

2.90

2.90

MI–TM081

46.19654

–86.83722

46.19654

–86.83725

247.93

247.95

–0.03

0.03

MI–TM082

46.16161

–86.79492

46.16161

–86.79493

246.17

246.19

–0.02

0.02

MI–TM083

46.17095

–86.77776

46.17096

–86.77777

245.68

245.66

0.02

0.02

MI–TM084

46.18688

–86.78758

46.18689

–86.78760

245.01

245.08

–0.07

0.07

Tables 3–7  19

UC LAT (NAD 83)

Gravity station

[“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–TM085

46.17720

–86.75986

46.17722

–86.75988

241.53

241.41

0.12

0.12

MI–TM086

46.19373

–86.77394

46.19374

–86.77395

245.13

245.18

–0.05

0.05

MI–TM087

46.20733

–86.79234

46.20736

–86.79234

248.92

248.74

0.19

0.19

Gravity station

20   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Table 4. Location information for the 87 gravity stations located on the central Upper Peninsula of Michigan (MI–TM) from Drenth and others (2024).—Continued

Table 5. Location information for the 139 gravity stations located on the Medicine Bow Mountains of Wyoming (WY–MB) from Brown and others (2025). [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

WY–MB501

41.19046

–106.11896

41.19047

–106.11896

2460.88

2460.79

0.09

0.09

WY–MB502

41.19209

–106.09834

41.19209

–106.09833

2429.51

2429.49

0.02

0.02

WY–MB503

41.18895

–106.08313

41.18895

–106.08313

2437.69

2437.67

0.02

0.02

WY–MB504

41.17132

–106.06800

41.17133

–106.06801

2500.66

2500.68

–0.02

0.02

WY–MB505

41.16023

–106.05934

41.16023

–106.05934

2472.42

2472.48

–0.06

0.06

WY–MB506

41.15066

–106.05078

41.15065

–106.05078

2442.59

2442.57

0.02

0.02

WY–MB507

41.12864

–106.00738

41.12862

–106.00738

2275.38

2275.50

–0.12

0.12

WY–MB508

41.14603

–105.97609

41.14602

–105.97608

2262.97

2263.02

–0.05

0.05

WY–MB509

41.15448

–105.96188

41.15447

–105.96188

2258.65

2258.61

0.03

0.03

WY–MB510

41.17454

–105.96686

41.17453

–105.96686

2247.40

2247.43

–0.03

0.03

WY–MB511

41.18156

–105.98720

41.18155

–105.98720

2282.29

2282.41

–0.12

0.12

WY–MB512

41.20351

–106.08888

41.20350

–106.08889

2410.12

2410.13

–0.02

0.02

WY–MB513

41.21824

–106.08941

41.21823

–106.08941

2391.92

2391.99

–0.07

0.07

WY–MB514

41.15853

–106.14494

41.15853

–106.14494

2795.83

2795.68

0.15

0.15

WY–MB515

41.14441

–106.16480

41.14441

–106.16481

2779.63

2779.50

0.13

0.13

WY–MB516

41.13696

–106.17559

41.13695

–106.17559

2758.67

2758.61

0.06

0.06

WY–MB517

41.14246

–106.32488

41.14245

–106.32489

2753.32

2753.25

0.07

0.07

WY–MB518

41.13765

–106.28225

41.13764

–106.28225

2858.08

2857.98

0.09

0.09

WY–MB519

41.12641

–106.30779

41.12640

–106.30779

2750.82

2750.57

0.25

0.25

WY–MB520

41.23346

–106.08523

41.23346

–106.08524

2390.43

2390.38

0.06

0.06

WY–MB521

41.25531

–106.07788

41.25530

–106.07789

2398.72

2398.82

–0.10

0.10

WY–MB522

41.26883

–106.07623

41.26882

–106.07624

2377.67

2377.53

0.14

0.14

WY–MB523

41.37340

–106.14062

41.37340

–106.14062

2992.11

2991.83

0.28

0.28

WY–MB524

41.38964

–106.13893

41.38964

–106.13892

3011.35

3011.26

0.09

0.09

WY–MB525

41.39842

–106.15353

41.39843

–106.15353

3011.69

3011.50

0.19

0.19

WY–MB526

41.40553

–106.13729

41.40554

–106.13730

2993.47

2993.19

0.28

0.28

WY–MB527

41.42013

–106.13554

41.42014

–106.13554

2948.85

2948.73

0.13

0.13

WY–MB528

41.43273

–106.14276

41.43274

–106.14277

2951.91

2951.76

0.16

0.16

Tables 3–7  21

UC LAT (NAD 83)

Gravity station

[“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

WY–MB529

41.44584

–106.12669

41.44585

–106.12670

2859.09

2858.94

0.15

0.15

WY–MB530

41.45834

–106.13831

41.45835

–106.13831

2819.46

2819.35

0.11

0.11

WY–MB531

41.46409

–106.12843

41.46408

–106.12843

2728.25

2726.51

1.74

1.74

WY–MB532

41.47037

–106.13498

41.47037

–106.13499

2712.18

2711.95

0.22

0.22

WY–MB533

41.44986

–106.14798

41.44985

–106.14799

3003.16

3002.90

0.26

0.26

WY–MB534

41.46325

–106.16424

41.46323

–106.16423

3002.84

3000.83

2.01

2.01

WY–MB535

41.35344

–106.14914

41.35343

–106.14915

3004.88

3004.53

0.36

0.36

WY–MB536

41.34814

–106.13838

41.34813

–106.13839

2988.23

2988.12

0.11

0.11

WY–MB537

41.34938

–106.12493

41.34936

–106.12494

2934.99

2935.39

–0.39

0.39

WY–MB538

41.36627

–106.15873

41.36628

–106.15873

2911.80

2911.48

0.33

0.33

WY–MB539

41.35579

–106.16766

41.35579

–106.16765

2783.06

2782.97

0.09

0.09

WY–MB540

41.34026

–106.16665

41.34026

–106.16665

2715.32

2715.31

0.01

0.01

WY–MB541

41.32412

–106.15787

41.32412

–106.15788

2591.96

2591.90

0.06

0.06

WY–MB542

41.28275

–106.10119

41.28275

–106.10119

2398.06

2398.07

–0.01

0.01

WY–MB543

41.22071

–106.07606

41.22071

–106.07605

2424.21

2424.16

0.04

0.04

WY–MB544

41.29736

–106.09520

41.29734

–106.09519

2421.13

2421.07

0.06

0.06

WY–MB545

41.31386

–106.13360

41.31385

–106.13360

2503.76

2503.76

0.00

0.00

WY–MB546

41.30799

–106.19738

41.30798

–106.19738

2918.25

2918.09

0.16

0.16

WY–MB547

41.30010

–106.21714

41.30009

–106.21715

3031.58

3031.27

0.31

0.31

WY–MB548

41.30087

–106.23800

41.30087

–106.23801

3051.26

3050.86

0.40

0.40

WY–MB549

41.31434

–106.25226

41.31435

–106.25228

3175.54

3175.29

0.26

0.26

WY–MB550

41.29642

–106.25626

41.29643

–106.25626

3037.98

3037.74

0.24

0.24

WY–MB551

41.29071

–106.27609

41.29071

–106.27609

3074.98

3074.79

0.18

0.18

WY–MB552

41.28869

–106.23471

41.28869

–106.23470

2960.95

2960.86

0.09

0.09

WY–MB553

41.36680

–106.17763

41.36680

–106.17764

2863.04

2862.91

0.13

0.13

WY–MB554

41.37821

–106.18809

41.37821

–106.18811

2961.79

2961.66

0.13

0.13

WY–MB555

41.39324

–106.19073

41.39322

–106.19074

3036.08

3035.80

0.28

0.28

WY–MB556

41.40661

–106.19289

41.40661

–106.19292

3071.54

3071.25

0.29

0.29

Gravity station

22   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Table 5. Location information for the 139 gravity stations located on the Medicine Bow Mountains of Wyoming (WY–MB) from Brown and others (2025).—Continued

Table 5. Location information for the 139 gravity stations located on the Medicine Bow Mountains of Wyoming (WY–MB) from Brown and others (2025).—Continued [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

WY–MB557

41.41760

–106.20245

41.41760

–106.20246

3053.63

3053.41

0.22

0.22

WY–MB558

41.43099

–106.19507

41.43097

–106.19507

2995.17

2994.94

0.23

0.23

WY–MB559

41.44403

–106.20143

41.44402

–106.20143

2984.57

2984.45

0.12

0.12

WY–MB560

41.45465

–106.21267

41.45464

–106.21267

2990.60

2990.29

0.31

0.31

WY–MB561

41.47833

–106.20236

41.47833

–106.20237

2948.91

2948.56

0.35

0.35

WY–MB562

41.46460

–106.21658

41.46458

–106.21659

2981.29

2980.69

0.60

0.60

WY–MB563

41.45965

–106.23325

41.45964

–106.23326

2973.63

2972.78

0.84

0.84

WY–MB564

41.47559

–106.23098

41.47558

–106.23098

3071.38

3071.14

0.24

0.24

WY–MB565

41.48384

–106.21930

41.48382

–106.21931

3006.00

3005.99

0.01

0.01

WY–MB566

41.46829

–106.24426

41.46830

–106.24427

3078.87

3078.57

0.30

0.30

WY–MB567

41.46409

–106.25900

41.46410

–106.25901

3077.83

3077.51

0.32

0.32

WY–MB568

41.45625

–106.27278

41.45626

–106.27279

3073.55

3073.28

0.28

0.28

WY–MB569

41.47083

–106.28379

41.47083

–106.28379

3055.43

3055.10

0.33

0.33

WY–MB570

41.48219

–106.29583

41.48219

–106.29583

3035.92

3035.62

0.30

0.30

WY–MB571

41.48175

–106.26286

41.48175

–106.26286

3049.84

3049.53

0.31

0.31

WY–MB572

41.48593

–106.27830

41.48593

–106.27829

2946.91

2943.84

3.07

3.07

WY–MB573

41.49861

–106.28389

41.49859

–106.28388

2940.61

2940.15

0.46

0.46

WY–MB574

41.35981

–106.18889

41.35980

–106.18890

2905.54

2905.28

0.26

0.26

WY–MB575

41.37285

–106.20701

41.37284

–106.20703

3062.39

3062.32

0.06

0.06

WY–MB576

41.34232

–106.18478

41.34233

–106.18479

2782.28

2782.24

0.04

0.04

WY–MB577

41.32670

–106.18456

41.32671

–106.18457

2715.24

2715.07

0.17

0.17

WY–MB578

41.34118

–106.20559

41.34119

–106.20560

2889.53

2889.33

0.19

0.19

WY–MB579

41.34941

–106.21585

41.34942

–106.21586

2987.38

2987.24

0.14

0.14

WY–MB580

41.37390

–106.24701

41.37391

–106.24702

3225.84

3225.55

0.28

0.28

WY–MB581

41.35780

–106.23310

41.35781

–106.23311

3095.83

3095.61

0.21

0.21

WY–MB582

41.35124

–106.25958

41.35124

–106.25958

3194.78

3194.64

0.15

0.15

WY–MB583

41.35149

–106.28097

41.35150

–106.28097

3232.63

3232.42

0.21

0.21

WY–MB584

41.36052

–106.26446

41.36051

–106.26446

3218.29

3218.16

0.14

0.14

Tables 3–7  23

UC LAT (NAD 83)

Gravity station

[“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

WY–MB585

41.35868

–106.29448

41.35868

–106.29448

3284.27

3283.91

0.36

0.36

WY–MB586

41.34155

–106.30579

41.34154

–106.30579

3297.49

3297.19

0.30

0.30

WY–MB587

41.33727

–106.32016

41.33726

–106.32014

3217.13

3216.98

0.15

0.15

WY–MB588

41.31288

–106.35900

41.31287

–106.35901

3193.29

3193.05

0.24

0.24

WY–MB589

41.32837

–106.36387

41.32837

–106.36388

3101.15

3100.74

0.40

0.40

WY–MB590

41.31746

–106.34631

41.31745

–106.34631

3237.34

3237.02

0.32

0.32

WY–MB591

41.34466

–106.37755

41.34466

–106.37756

2967.09

2966.92

0.17

0.17

WY–MB592

41.36871

–106.37738

41.36872

–106.37738

3165.50

3165.29

0.21

0.21

WY–MB593

41.35236

–106.38354

41.35238

–106.38354

3051.39

3051.08

0.31

0.31

WY–MB594

41.35572

–106.39594

41.35571

–106.39594

3062.38

3062.04

0.35

0.35

WY–MB595

41.34407

–106.39819

41.34405

–106.39819

2891.16

2890.84

0.32

0.32

WY–MB596

41.33876

–106.41694

41.33875

–106.41695

2813.67

2813.49

0.18

0.18

WY–MB597

41.33429

–106.43329

41.33428

–106.43327

2760.35

2759.99

0.37

0.37

WY–MB598

41.33547

–106.44952

41.33545

–106.44952

2710.88

2710.85

0.03

0.03

WY–MB599

41.33790

–106.46655

41.33789

–106.46655

2655.43

2655.29

0.14

0.14

WY–MB600

41.33987

–106.48643

41.33986

–106.48643

2594.88

2593.65

1.23

1.23

WY–MB601

41.34448

–106.50441

41.34447

–106.50442

2532.86

2532.66

0.20

0.20

WY–MB602

41.34786

–106.52291

41.34786

–106.52291

2473.84

2473.60

0.24

0.24

WY–MB603

41.33354

–106.51849

41.33353

–106.51850

2501.09

2500.91

0.18

0.18

WY–MB604

41.32386

–106.50193

41.32386

–106.50194

2542.62

2542.48

0.13

0.13

WY–MB605

41.32934

–106.47678

41.32934

–106.47679

2651.04

2651.03

0.01

0.01

WY–MB606

41.32362

–106.45199

41.32361

–106.45200

2795.31

2795.19

0.11

0.11

WY–MB607

41.32444

–106.42427

41.32445

–106.42427

2924.09

2923.92

0.17

0.17

WY–MB608

41.34529

–106.36044

41.34529

–106.36044

3139.11

3138.94

0.17

0.17

WY–MB609

41.32228

–106.38720

41.32226

–106.38721

3002.52

3002.33

0.20

0.20

WY–MB610

41.32165

–106.40788

41.32164

–106.40789

2963.94

2963.83

0.12

0.12

WY–MB611

41.30854

–106.38885

41.30853

–106.38886

3091.15

3090.82

0.33

0.33

WY–MB612

41.29757

–106.38299

41.29756

–106.38301

3145.04

3144.82

0.23

0.23

Gravity station

24   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Table 5. Location information for the 139 gravity stations located on the Medicine Bow Mountains of Wyoming (WY–MB) from Brown and others (2025).—Continued

Table 5. Location information for the 139 gravity stations located on the Medicine Bow Mountains of Wyoming (WY–MB) from Brown and others (2025).—Continued [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

WY–MB613

41.28690

–106.37650

41.28688

–106.37651

3101.85

3101.70

0.15

0.15

WY–MB614

41.28625

–106.36023

41.28622

–106.36024

3000.34

3000.32

0.02

0.02

WY–MB615

41.29944

–106.35908

41.29941

–106.35908

3041.05

3040.98

0.07

0.07

WY–MB616

41.31065

–106.43329

41.31064

–106.43329

2844.89

2844.62

0.27

0.27

WY–MB617

41.30018

–106.44522

41.30019

–106.44522

2758.25

2757.99

0.26

0.26

WY–MB618

41.28856

–106.45303

41.28856

–106.45302

2683.76

2683.56

0.20

0.20

WY–MB619

41.27598

–106.45759

41.27598

–106.45759

2639.70

2639.35

0.36

0.36

WY–MB620

41.26670

–106.46745

41.26671

–106.46746

2645.14

2645.02

0.12

0.12

WY–MB621

41.26679

–106.44476

41.26681

–106.44476

2736.19

2735.88

0.31

0.31

WY–MB622

41.27006

–106.43131

41.27006

–106.43133

2882.87

2882.68

0.19

0.19

WY–MB623

41.27354

–106.41756

41.27355

–106.41759

2952.08

2951.80

0.28

0.28

WY–MB624

41.27963

–106.40188

41.27964

–106.40189

3063.82

3063.56

0.26

0.26

WY–MB625

41.28087

–106.39031

41.28089

–106.39033

2996.80

2996.54

0.27

0.27

WY–MB626

41.26106

–106.45549

41.26107

–106.45551

2661.72

2661.57

0.14

0.14

WY–MB627

41.24914

–106.46006

41.24914

–106.46008

2630.29

2630.16

0.13

0.13

WY–MB628

41.25510

–106.44214

41.25509

–106.44216

2607.63

2607.44

0.19

0.19

WY–MB629

41.25491

–106.42496

41.25491

–106.42498

2589.58

2589.53

0.05

0.05

WY–MB630

41.26188

–106.40894

41.26187

–106.40896

2607.66

2607.64

0.02

0.02

WY–MB631

41.26583

–106.39063

41.26583

–106.39064

2630.06

2630.01

0.05

0.05

WY–MB632

41.27255

–106.37663

41.27255

–106.37665

2654.22

2654.14

0.08

0.08

WY–MB633

41.27468

–106.35713

41.27467

–106.35714

2709.69

2709.54

0.15

0.15

WY–MB634

41.28110

–106.34158

41.28111

–106.34158

2786.58

2786.38

0.20

0.20

WY–MB635

41.28853

–106.32510

41.28853

–106.32510

2869.34

2869.14

0.20

0.20

WY–MB636

41.29911

–106.31249

41.29912

–106.31250

2944.56

2944.28

0.28

0.28

WY–MB637

41.28596

–106.47379

41.28599

–106.47380

2667.80

2667.65

0.15

0.15

WY–MB638

41.29813

–106.46541

41.29815

–106.46541

2713.53

2712.81

0.72

0.72

WY–MB639

41.30993

–106.45597

41.30994

–106.45597

2781.09

2780.78

0.31

0.31

Tables 3–7  25

UC LAT (NAD 83)

Gravity station

[“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

CO–WM300

38.31006

–105.35527

38.31004

–105.35526

2101.07

2100.79

0.28

0.28

CO–WM301

38.29837

–105.35721

38.29836

–105.35720

2121.65

2121.24

0.41

0.41

CO–WM302

38.29417

–105.35912

38.29414

–105.35911

2129.39

2129.27

0.12

0.12

CO–WM303

38.28391

–105.35495

38.28388

–105.35494

2148.09

2147.80

0.29

0.29

CO–WM304

38.27815

–105.35453

38.27814

–105.35454

2161.73

2161.46

0.27

0.27

CO–WM305

38.26982

–105.35141

38.26981

–105.35141

2179.96

2179.67

0.29

0.29

CO–WM306

38.26719

–105.35480

38.26718

–105.35478

2190.84

2190.52

0.31

0.31

CO–WM307

38.26307

–105.35934

38.26307

–105.35934

2201.92

2201.58

0.34

0.34

CO–WM308

38.26102

–105.36307

38.26101

–105.36307

2205.89

2205.52

0.37

0.37

CO–WM309

38.25873

–105.37197

38.25871

–105.37198

2216.51

2216.19

0.32

0.32

CO–WM310

38.25520

–105.37646

38.25519

–105.37647

2225.66

2225.29

0.37

0.37

CO–WM311

38.24934

–105.38073

38.24933

–105.38074

2234.73

2235.09

–0.36

0.36

CO–WM312

38.24533

–105.38533

38.24530

–105.38534

2243.33

2243.24

0.09

0.09

CO–WM313

38.23901

–105.38846

38.23899

–105.38848

2251.57

2251.65

–0.08

0.08

CO–WM314

38.23475

–105.38846

38.23472

–105.38848

2255.81

2255.86

–0.05

0.05

CO–WM315

38.23029

–105.39197

38.23027

–105.39198

2268.27

2268.04

0.23

0.23

CO–WM316

38.22717

–105.39525

38.22715

–105.39526

2274.55

2274.18

0.37

0.37

CO–WM317

38.22334

–105.40223

38.22333

–105.40227

2283.91

2284.06

–0.15

0.15

CO–WM318

38.23332

–105.40288

38.23332

–105.40289

2408.71

2408.94

–0.23

0.23

CO–WM319

38.24616

–105.39248

38.24616

–105.39249

2265.00

2264.96

0.04

0.04

CO–WM320

38.25602

–105.38334

38.25603

–105.38334

2280.00

2279.92

0.08

0.08

CO–WM321

38.25907

–105.39394

38.25907

–105.39394

2395.34

2395.25

0.09

0.09

CO–WM322

38.26950

–105.38806

38.26949

–105.38806

2398.60

2398.58

0.02

0.02

CO–WM323

38.27286

–105.39011

38.27285

–105.39010

2366.05

2366.03

0.02

0.02

CO–WM324

38.28735

–105.40667

38.28734

–105.40668

2274.35

2274.34

0.01

0.01

CO–WM325

38.29275

–105.39243

38.29273

–105.39243

2380.41

2380.41

0.00

0.00

CO–WM326

38.28718

–105.38645

38.28716

–105.38644

2326.71

2326.70

0.01

0.01

CO–WM327

38.28938

–105.37284

38.28937

–105.37284

2224.21

2224.36

–0.15

0.15

Gravity station

26   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Table 6. Location information for the 59 gravity stations located on the Wet Mountains of Colorado (CO–WM) from Magnin and Anderson (2024).

Table 6. Location information for the 59 gravity stations located on the Wet Mountains of Colorado (CO–WM) from Magnin and Anderson (2024).—Continued [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

CO–WM328

38.26465

–105.37548

38.26464

–105.37548

2345.54

2345.57

–0.03

0.03

CO–WM329

38.25471

–105.36188

38.25469

–105.36189

2257.82

2257.71

0.11

0.11

CO–WM330

38.24246

–105.37071

38.24243

–105.37072

2302.59

2302.50

0.09

0.09

CO–WM331

38.24149

–105.37685

38.24148

–105.37686

2271.77

2271.80

–0.03

0.03

CO–WM332

38.11421

–105.18721

38.11420

–105.18722

2861.52

2861.53

–0.01

0.01

CO–WM333

38.11834

–105.19500

38.11833

–105.19501

2844.62

2844.62

0.00

0.00

CO–WM334

38.11539

–105.20536

38.11538

–105.20537

2851.18

2851.05

0.13

0.13

CO–WM335

38.09331

–105.21065

38.09329

–105.21067

3263.29

3262.77

0.52

0.52

CO–WM336

38.09336

–105.17546

38.09335

–105.17546

2882.04

2881.89

0.15

0.15

CO–WM337

38.09126

–105.17045

38.09125

–105.17046

2955.79

2956.01

–0.22

0.22

CO–WM338

38.08223

–105.15881

38.08223

–105.15881

3160.31

3160.21

0.10

0.10

CO–WM339

38.07603

–105.16797

38.07602

–105.16798

3161.92

3161.82

0.10

0.10

CO–WM340

38.07120

–105.17146

38.07119

–105.17146

3105.18

3105.85

–0.67

0.67

CO–WM341

38.06252

–105.18589

38.06251

–105.18590

3159.78

3159.74

0.04

0.04

CO–WM342

38.07947

–105.18030

38.07944

–105.18029

2925.35

2921.68

3.67

3.67

CO–WM343

38.09203

–105.14929

38.09200

–105.14930

2941.98

2941.82

0.16

0.16

CO–WM344

38.08606

–105.13990

38.08604

–105.13992

2934.02

2934.11

–0.09

0.09

CO–WM345

38.06968

–105.22355

38.06966

–105.22357

3171.49

3171.41

0.08

0.08

CO–WM346

38.07352

–105.22973

38.07349

–105.22974

3189.36

3189.64

–0.28

0.28

CO–WM347

38.07574

–105.24332

38.07573

–105.24332

3115.80

3115.97

–0.17

0.17

CO–WM348

38.06952

–105.23591

38.06952

–105.23593

3068.66

3068.71

–0.05

0.05

CO–WM349

38.06894

–105.22953

38.06894

–105.22955

3119.76

3119.79

–0.03

0.03

CO–WM350

38.07174

–105.23154

38.07174

–105.23155

3142.02

3141.75

0.27

0.27

CO–WM351

38.08607

–105.24776

38.08608

–105.24777

3170.15

3169.97

0.18

0.18

CO–WM352

38.08224

–105.24561

38.08225

–105.24560

3140.92

3140.93

–0.01

0.01

CO–WM353

38.08720

–105.22464

38.08718

–105.22464

3220.52

3220.55

–0.03

0.03

CO–WM354

38.09102

–105.22718

38.09100

–105.22718

3260.35

3260.33

0.02

0.02

CO–WM355

38.09067

–105.23085

38.09065

–105.23085

3287.47

3288.41

–0.94

0.94

Tables 3–7  27

UC LAT (NAD 83)

Gravity station

[“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

CO–WM356

38.08886

–105.23514

38.08884

–105.23515

3280.50

3280.75

–0.25

0.25

CO–WM357

38.08606

–105.23944

38.08604

–105.23944

3224.57

3224.24

0.33

0.33

CO–WM358

38.09156

–105.23900

38.09156

–105.23900

3313.85

3313.63

0.22

0.22

Gravity station

28   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Table 6. Location information for the 59 gravity stations located on the Wet Mountains of Colorado (CO–WM) from Magnin and Anderson (2024).—Continued

Table 7. Location information for the 56 gravity stations located on the Keweenaw Peninsula of Michigan (MI–KP0) from Murchek and others (2025). [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–KP0–1

47.19247

–88.48937

47.19248

–88.48937

357.97

357.83

0.14

0.14

MI–KP0–2

47.2386

–88.61182

47.2386

–88.61181

186.12

185.82

0.31

0.31

MI–KP0–3

47.23998

–88.59468

47.23998

–88.59468

189.02

188.48

0.53

0.53

MI–KP0–4

47.22843

–88.59448

47.22843

–88.59447

198.65

198.31

0.34

0.34

MI–KP0–5

47.23928

–88.59002

47.23927

–88.59001

194.57

194.42

0.15

0.15

MI–KP0–6

47.22849

–88.58479

47.22849

–88.58478

207.38

206.85

0.52

0.52

MI–KP0–7

47.22858

–88.57424

47.22857

–88.57423

221.14

220.83

0.31

0.31

MI–KP0–8

47.23942

–88.57388

47.23942

–88.57387

193.96

193.64

0.32

0.32

MI–KP0–9

47.21424

–88.56344

47.21424

–88.56344

283.33

283.08

0.24

0.24

MI–KP0–10

47.18516

–88.56329

47.18514

–88.56329

318.8

318.72

0.07

0.07

MI–KP0–11

47.21431

–88.55265

47.21431

–88.55265

300.8

300.71

0.09

0.09

MI–KP0–12

47.2138

–88.54236

47.21379

–88.54235

323.6

323.45

0.15

0.15

MI–KP0–13

47.18488

–88.53965

47.18487

–88.53965

333.12

333.07

0.05

0.05

MI–KP0–14

47.21352

–88.53145

47.21351

–88.53145

327.83

327.78

0.05

0.05

MI–KP0–15

47.2067

–88.52091

47.2067

–88.52091

351.73

351.52

0.21

0.21

MI–KP0–16

47.21363

–88.52081

47.21362

–88.52080

351.01

350.86

0.14

0.14

MI–KP0–17

47.17455

–88.51855

47.17455

–88.51855

333.56

333.58

–0.02

0.02

MI–KP0–18

47.20738

–88.50895

47.20738

–88.50894

360.93

360.62

0.3

0.3

MI–KP0–19

47.19935

–88.49976

47.19935

–88.49976

359.46

359.34

0.12

0.12

MI–KP0–20

47.19036

–88.49965

47.19035

–88.49964

353.43

353.33

0.11

0.11

MI–KP0–21

47.19392

–88.49585

47.19391

–88.49584

356.93

356.81

0.13

0.13

MI–KP0–22

47.18917

–88.49367

47.18916

–88.49366

357.74

357.75

–0.01

0.01

MI–KP0–23

47.13179

–88.48612

47.13179

–88.48611

190.49

190.48

0.01

0.01

MI–KP0–24

47.18894

–88.48568

47.18893

–88.48569

365.34

365.31

0.03

0.03

MI–KP0–25

47.19946

–88.48134

47.19945

–88.48133

359.63

359.64

–0.01

0.01

MI–KP0–26

47.18842

–88.47950

47.18843

–88.47950

354.53

354.47

0.06

0.06

MI–KP0–27

47.14506

–88.47867

47.14505

–88.47866

269.47

269.54

–0.07

0.07

MI–KP0–28

47.1853

–88.47486

47.1853

–88.47485

354.27

354.18

0.09

0.09

Tables 3–7  29

UC LAT (NAD 83)

Gravity station

[“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

Gravity station

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–KP0–29

47.18223

–88.47155

47.18223

–88.47155

337.38

337.18

0.2

0.2

MI–KP0–30

47.20707

–88.47135

47.20707

–88.47134

365.23

365.08

0.16

0.16

MI–KP0–31

47.17804

–88.46743

47.17804

–88.46742

324.81

324.66

0.15

0.15

MI–KP0–32

47.17419

–88.46209

47.17417

–88.46208

314.46

313.5

0.96

0.96

MI–KP0–33

47.14882

–88.45806

47.14882

–88.45805

190.87

190.83

0.04

0.04

MI–KP0–34

47.1709

–88.45737

47.1709

–88.45735

306.2

305.98

0.23

0.23

MI–KP0–35

47.1674

–88.45202

47.16739

–88.45201

273.23

273.09

0.14

0.14

MI–KP0–36

47.20788

–88.44713

47.20787

–88.44713

356.72

356.54

0.18

0.18

MI–KP0–37

47.18497

–88.44653

47.18497

–88.44652

311.7

311.74

–0.04

0.04

MI–KP0–38

47.19592

–88.44639

47.19592

–88.44638

337

337.07

–0.07

0.07

MI–KP0–39

47.16467

–88.44616

47.16467

–88.44616

226.5

226.36

0.14

0.14

MI–KP0–40

47.17307

–88.44471

47.17307

–88.44472

267.49

267.5

–0.02

0.02

MI–KP0–41

47.16355

–88.43915

47.16355

–88.43915

187.27

187.07

0.2

0.2

MI–KP0–42

47.16107

–88.43507

47.16106

–88.43507

184.11

184.01

0.09

0.09

MI–KP0–43

47.15894

–88.43253

47.15893

–88.43253

184.07

183.39

0.68

0.68

MI–KP0–44

47.15447

–88.42535

47.15447

–88.42535

198.16

197.89

0.27

0.27

MI–KP0–45

47.14855

–88.42529

47.14855

–88.42530

216.99

214.63

2.36

2.36

MI–KP0–46

47.14989

–88.41428

47.14988

–88.41427

185.32

185.14

0.18

0.18

MI–KP0–47

47.14911

–88.39405

47.14911

–88.39405

196.15

195.98

0.17

0.17

MI–KP0–48

47.17075

–88.38358

47.17075

–88.38358

226.35

226.3

0.05

0.05

MI–KP0–49

47.14186

–88.38341

47.14186

–88.38341

218.89

218.71

0.19

0.19

MI–KP0–50

47.12744

–88.36215

47.12744

–88.36215

221.22

221

0.22

0.22

MI–KP0–51

47.11808

–88.34597

47.11808

–88.34597

202.97

202.92

0.05

0.05

MI–KP0–52

47.1147

–88.34358

47.11471

–88.34357

195.46

194.97

0.49

0.49

MI–KP0–53

47.17086

–88.34093

47.17086

–88.34092

226.63

226.62

0.01

0.01

30   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

Table 7. Location information for the 56 gravity stations located on the Keweenaw Peninsula of Michigan (MI–KP0) from Murchek and others (2025).—Continued

Table 7. Location information for the 56 gravity stations located on the Keweenaw Peninsula of Michigan (MI–KP0) from Murchek and others (2025).—Continued [“Gravity station” refers to a unique identifier that is based on the station number used in the data source for this table. “UC LAT” and “UC LON” refer to the horizontal coordinates before differential correction. “DC LAT” and “DC LON” refer to the differentially corrected horizontal coordinates. “DGNSS ELEV” refers to the elevation data from the differential global navigation satellite system (dGNSS) equipment from the differentially corrected coordinates. “LIDAR ELEV” refers to the light detection and ranging (lidar) derived elevation from the 1-meter (m) digital elevation model (DEM). “ELEV DIFF” refers to the elevation difference (dGNSS – lidar), and “ABS ELEV DIFF” refers to the absolute value of this difference. NAD 83, North American Datum of 1983; NAVD 88, North American Vertical Datum of 1988]

Gravity station

UC LAT (NAD 83)

UC LON (NAD 83)

DC LAT (NAD 83)

DC LON (NAD 83)

DGNSS ELEV (m) (NAVD 88)

LIDAR ELEV (m) (NAVD 88)

ELEV DIFF (m)

ABS ELEV DIFF (m)

MI–KP0–54

47.13504

–88.29974

47.13504

–88.29974

192.78

192.76

0.02

0.02

MI–KP0–55

47.12553

–88.29006

47.12553

–88.29006

189.78

187.11

2.67

2.67

1MI–KP0–M65

47.2436

–88.44747

47.2436

–88.44747

380.18

379.69

0.49

0.49

1The data for this National Geodetic Survey (NGS) station are from NGS (2025). This gravity station identifier was created for use in this report.

Tables 3–7  31

32   The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction

References Cited American Society for Photogrammetry and Remote Sensing [ASPRS], 2023, ASPRS positional accuracy standards for digital geospatial data (2d ed., ver. 1): ASPRS web page, accessed August 20, 2023, at http​s://public​ documents.​asprs.org/​Po​sitionalAc​curacyStd-​Ed2-​V1. Barber, C.P., and Shortrudge, A.M., 2004, Light detection and ranging (lidar)-derived elevation data for surface hydrology applications: East Lansing, Mich., Institute of Water Resources, Michigan State University, 17 p. Blakely, R.J., 1995, Gravity anomalies, section 7.3 in chap. 7 of Potential theory in gravity and magnetic applications: New York, Cambridge University Press, p. 136–146. Brown, P.J., Reitman, J.J., Drenth, B.J., and Lynds, R.M., 2025, Principal facts of regional gravity data in the Medicine Bow Mountains, Wyoming, 2022–2024: U.S. Geological Survey data release, accessed April 15, 2025, at https://doi.org/​ 10.5066/​P1XAVOXP. Burrough, P.A., and McDonnell, R.A., 1998, Principles of geographical information systems: New York, Oxford University Press, 333 p. Cai, Z., Ma, H., and Zhang, L., 2020, Feature selection for airborne LiDAR data filtering—A mutual information method with Parzon window optimization: GIScience & Remote Sensing, v. 57, no. 3, p. 323–337, accessed August 20, 2023, at https://doi.org/​10.1080/​1​5481603.20​ 19.1695406. Drenth, B.J., Reitman, J.J., and Brown, P.J., 2024, Principal facts of regional gravity in the central Upper Peninsula, Michigan, 2022–2023: U.S. Geological Survey data release, accessed November 20, 2024, at https://doi.org/​10.5066/​ P9L6KML9. Elaksher, A.F., 2016, Co-registering satellite images and LIDAR DEMs through straight lines: International Journal of Image and Data Fusion, v. 7, no. 2, p. 103–118, accessed August 20, 2023, at https://doi.org/​10.1080/​1​9479832.20​ 15.1075607. Habib, A., Ghanma, M., Morgan, M., and Al-Ruzouq, R., 2005, Photogrammetric and lidar data registration using linear features: Photogrammetric Engineering and Remote Sensing v. 71, no. 6, p. 699–707, accessed August 20, 2023, at https://doi.org/​10.14358/​PERS.71.6.699. Hinze, W.J., Von Frese, R.R.B., and Saad, A.H., 2012, Gravity and magnetic exploration—Principles, practices, and applications: New York, Cambridge University Press, 528 p.

Hollaus, M., Wagner, W., and Kraus, K., 2005, Airborne laser scanning and usefulness for hydrologic models: Advances in Geosciences, v. 5, p. 57–63, accessed August 20, 2023, at https://doi.org/​10.5194/​adgeo-​5-​57-​2005. Liu, X., 2008, Airborne lidar for DEM generation—Some critical issues: Progress in Physical Geography— Earth and Environment, v. 32, no. 1, p. 31–49, accessed August 20, 2023, at https://doi.org/​10.1177/​ 0309133308089496. Longman, I.M., 1959, Formulas for computing the tidal accelerations due to the moon and the sun: Journal of Geophysical Research, v. 64, no. 12, p. 2351–2355, accessed August 20, 2023, at https://doi.org/​10.1029/​ JZ064i012p02351. Magnin, B.P., and Anderson, E.D., 2024, Gravity data in the Wet Mountains area, southcentral Colorado, 2023: U.S. Geological Survey data release, accessed January 16, 2024, at https://doi.org/​10.5066/​P9KNDYU3. Merrick & Company, 2016, Arkansas River (partial Fremont) in LiDAR download portal: Colorado Hazard Mapping & Risk MAP Portal website, accessed January 24, 2024, at http​s://colora​dohazardma​pping.com/​lidarDownload. Morelli, C., Gantar, C., Honkasalo, T., McConnell, R.K., Tanner, J.G., Szabo, B., Uotila, U., and Whalen, C.T., 1972, The International Gravity Standardization Net 1971 (I.G.S.N. 71): Air Force Cambridge Research Laboratories and European Office of Aerospace Research and Development, Special Publication 4, prepared by International Association of Geodesy, Osservatorio Geofisico Sperimentale [Experimental Geophysical Observatory] under contract AF61 (052) 656, 194 p. [Also available at ht​tps://apps​.dtic.mil/​sti/​pdfs/​ADA006203.pdf.] Murchek, J.T., DeGraff, J.M., and Drenth, B.J., 2025, Gravity data in the Upper Peninsula, Michigan for geophysical profile modeling of the Midcontinent Rift System and associated structures: U.S. Geological Survey data release, accessed March 6, 2025, at https://doi.org/​10.5066/​P14BGVCT. National Geodetic Survey [NGS], 2025, Designation M 65, PID SG0019, State/county MI/Houghton, Country US, USGS quad Laurium: NGS data sheet (datasheet95, ver. 8.12.5.19), accessed May 4, 2025, at https://www.ngs.noaa.gov/cgi-bin/ ds_mark.prl?PidBox=SG0019. National Oceanic and Atmospheric Administration [NOAA], 2012, Lidar 101—An introduction to lidar technology, data, and applications (revised): NOAA Coastal Services Center, Coastal Geospatial Services Division, Coastal Remote Sensing Program report, 72 p., accessed August 20, 2023, at https://oceanservice.noaa.gov/facts/lidar.html.

References Cited  33 National Oceanic and Atmospheric Administration [NOAA], 2023, What is lidar?: NOAA web page, accessed April 13, 2023, at https://oceanservice.noaa.gov/facts/lidar. html

Stoker, J.M., Greenlee, S.K., Gesch, D.B., and Menig, J.C., 2006, CLICK—The new USGS center for lidar information coordination and knowledge: Photogrammetric Engineering and Remote Sensing, v. 72, no. 6, p. 613–616.

Peucker, T.K., Fowler, R.J., Little, J.J., and Mark, D.M., 1976, Digital representation of three-dimensional surfaces by triangulated irregular networks (TIN) (revised): Office of Naval Research, Geography Programs, Technical Report 10, prepared by authors under contract N00014–75–C–0886 {NR 389–171}, 63 p. [Also available at ht​tps://apps​.dtic.mil/​ sti/​pdfs/​ADA094241.pdf.]

Sugarbaker, L.J., Constance, E.W., Heidemann, H.K., Jason, A.L., Lukas, V., Saghy, D.L., and Stoker, J.M., 2014, The 3D Elevation Program initiative—A call for action: U.S. Geological Survey Circular 1399, 35 p., accessed August 20, 2023, at https://doi.org/​10.3133/​cir1399.

Pfeifer, N., and Briese, C., 2007, Geometrical aspects of airborne laser scanning and terrestrial laser scanning in Rönnholm, P., Hyyppä, H., and Hyyppä, J., eds., Proceedings of the ISPRS Workshop on Laser Scanning 2007 and SilviLaser 2007, September 12–14, 2007, Espoo, Finland: International Society of Photogrammetry, Remote Sensing Archives, v. XXXVI, part 3/W52, p. 311–319, accessed August 20, 2023, at ht​tps://www.​isprs.org/​proceedings/​ XXXVI/​3-​W52/​final_​papers/​Pfeifer_​2007_​keynote.pdf. Quantum Spatial, Inc., 2020, 2018 Park, Custer, & Fremont Counties in LiDAR download portal: Colorado Hazard Mapping & Risk MAP Portal website, accessed January 24, 2024, at http​s://colora​dohazardma​pping.com/​ lidarDownload. Reutebuch, S.E., Andersen, H.-E., and McGaughey, R.J., 2005, Light detection and ranging (LIDAR): an emerging tool for multiple resource inventory: Journal of Forestry, v. 103, no. 6, p. 286–292, accessed August 20, 2023, at https://doi.org/​ 10.1093/​jof/​103.6.286. Sanborn Map Company, Inc., 2020, Lidar_Houghton_ Keweenaw_Ontonagon Hydro-Flattened Bare-Earth DEM in MTU GRF DEM Download Tool: Michigan Tech Great Lakes Research Center, Geospatial Research Facility website, accessed January 24, 2024, at https​://geospat​ialresearc​h.mtu.edu/​demindex/​.

U.S. Geological Survey [USGS], 2022, 3D Elevation Program 1-meter resolution digital elevation model, MI 13 County C16 in TNM Download [The National Map Downloader] version 2.0: U.S. Geological Survey website, accessed June 25, 2023, at https://a​pps.nation​almap.gov/​downloader. U.S. Geological Survey [USGS], 2023, 3D Elevation Program 1-meter resolution digital elevation model, WY South Central D20 in TNM Download [The National Map Downloader] version 2.0: U.S. Geological Survey website, accessed June 25, 2023, at https://a​pps.nation​almap.gov/​downloader. U.S. Geological Survey [USGS], 2024, Lidar base specification online [LBS 2024 rev. A.]: U.S. Geological Survey website, accessed January 31, 2024, at h​ttps://www​.usgs.gov/​3DEP/​ lidarspec. U.S. Geological Survey [USGS], [undated], Topographic data quality levels (QLs): U.S. Geological Survey website, accessed April 13, 2023, at h​ttps://www​.usgs.gov/​3d-​ elevation-​program/​topographic-​data-​quality-​levels-​qls. Webster, T.L., and Dias, G., 2006, An automated GIS procedure for comparing GPS and proximal LIDAR elevations: Computers & Geosciences, v. 32, no. 6, p. 713–726, accessed August 20, 2023, at https://doi.org/​10.1016/​ j.cageo.2005.08.009.

For more information about this publication, contact Director Energy and Minerals Mission Area U.S. Geological Survey 12201 Sunrise Valley Drive Reston, VA 20192-0002 For additional information, visit h​ttps://www​.usgs.gov/​mission-​areas/​energy-​and-​minerals Publishing support provided by the Science Publishing Network, Reston Publishing Service Center

Murchek and others—The Feasibility of Using Lidar-Derived Digital Elevation Models for Gravity Data Reduction—OFR 2025–1019, ver. 1.1.

ISSN 2331-1258 (online) https://doi.org/​10.3133/​ofr20251019

Related documents

Record · ID 621510 · SHA-256 f12ddc793bfc0b4d
Conceptio Open Knowledge Archive — every document is proof-bundled with source, license, and retrieval metadata.