Magnetic data were postprocessed following the methodology:
1)Initial Processing in MagMap2000
-Download base files from INTERMAGNET observatory (
http://www.intermagnet.org/apps/dl_data_prel_e.php)
-Despike base files and visually inspect for dropouts. Data were only lightly edited and may require additional noise-editing and/or heading corrections to make fully interpretable.
-Apply time-offset to DIURNAL (base) readings for GMT to local PDT + longitude difference from base to survey = -7hrs + 0.15hrs =-24,660sec **NOTE: 0.0666667hr/deg longitude correction
-Import .survey files, despike
-Export files as Geosoft .xyz
2)Import to Oasis-Montaj
3)Calculated DGRF for G882 data and FRN base station = DGRF_G882 and DGRF_FRN. IGRF 13 Generation model was calculated per Alken and others (2021),
https://doi.org/10.1186/s40623-020-01288-x. *NOTE: for survey dates 10/11/19 and 10/12/19, DGRF calculated in Oasis montaj (OM) was incorrect. DGRF values for these two days were corrected by applying a DC shift to the OM values based on DGRF value calculated using the following online calculator (
http://geomag.bgs.ac.uk/data_service/models_compass/igrf_calc.html).
4)Back calculated FRN Base Magnetometer measurements because they are NOT exported from MagMap2000
G882-DIURNAL_1 = FRN_base
5)Calculated DGRF value at base station FRN = 48138nT
6)Calculated Total Field (Residual) Anomaly (Resid_anom) using the following equation:
A(x) = [F(x,t) - I(x)] - [FB(t) - IB]
where:
A(x)is total field anomaly at position x
F(x,t) is the G882 reading at time t and position x
I(x) is IGRF/DGRF at the G882
FB(t) is the base station reading at time t
IB is the IGRF/DGRF at the base station
Projection is UTM z10N, wgs84.