V. Theoretical Fundamentals of Inertial Gravimetry

Total Page:16

File Type:pdf, Size:1020Kb

V. Theoretical Fundamentals of Inertial Gravimetry National Geodetic Survey V. Theoretical Fundamentals of Inertial Gravimetry • Basic gravimetry equation • Essential IMU data processing for gravimetry • Kalman filter approaches • Rudimentary error analysis – spectral window • Instrumentation Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.1 Inertial Gravimetry • Vector gravimetry, instead of scalar gravimetry − determine three components of gravity in the n-frame • Use precision accelerometer triad, instead of gravimeter − OTF (off-the-shelf) units designed for inertial navigation rather than gravimetry • Usually consider strapdown mechanization, instead of stabilized platform • Need precision gyroscopes to minimize effect of orientation error on horizontal components • Two documented approaches of data processing − either integrate accelerometer data or differentiate GPS data Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.2 Direct Method • Recall strapdown mechanization Gravitational Vector n= n i − ib in n-frame gCCib() xa accelerometers transformation from GPS inertial frame to n-frame transformation from body gyros frame to inertial frame xi – kinematic accelerations obtained from GPS-derived positions, x, in i-frame. ab – inertial accelerations measured by accelerometers in body frame tk b b δ vk − accelerometer data typically are δ va= ()t dt ⇒≈ak (e.g.) k ∫ δt tk−1 better approximations can be formulated ◦ where δ tt = kk − t − 1 is the data interval, e.g., δt = 1/50 s Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.3 i Determination of C b (1) • Let e be the unit vector about which a rotation ζ ζ by the angle, ζ, rotates the b-frame to the i-frame 3b cos()ζ 2 q1 b eζ bsin()ζ 2 q2 e = c → quaternion vector: q = = ζ ζ csin() 2 q3 d d sin()ζ 2 q4 2b qqqq2222+−−22()() qqqq + qqqq− 1 2 3 4 23 14 24 13 i = −2222 −+− + 1b Cb 22()qqqq23 14 qqqq 1 2 3 4() qqqq 34 12 2222 22()()qqqq24+ 13 qqqq34−12 qqqq 1−−+ 2 3 4 d 1 • Quaternions satisfy the differential equation: qq= A dt 2 0 ωωω123 ω1 −−ω10 ωω 32 b − where A = and ω =ω = gyro data ib 2 −−ωω230 ω 1 ω3 −−ωω32 ω 10 − this is a linear D.E. with no singularities Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.4 i Determination of C b (2) • Solution to D.E., if A is assumed constant, 1123 1 qqˆˆ=+++IkΘΘ Θ +− , = 1, 2, k2 kk 8 48 k k1 tk δ = b − using gyro data, typically given as θωk∫ ibdt tk −1 tk − where = ΘAk ∫ dt tk −1 − note: A is assumed constant only in the solution to the D.E., not in using the gyro data (model error, not data error in Θ) − q ˆ 0 is given by an initialization procedure − solution is a second-order algorithm, neglecting terms of order δt3 − higher-order algorithms are easily developed (Jekeli 2000) Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.5 Determination of δgn (1) • Kalman filter approach to minimizing estimation errors − estimate IMU systematic errors and gravity disturbance vector n − formulate in i-frame and assume negligible error in Ci − system state updates (observations) are differences in accelerations yag=+−()ii x i iiiiii i i i ii =+()aaδδδ +−()() gg −+ xx gx =γ + ΩΩie ie =iδ b −×ψ i ib −−δδ i i CCbba a xg δ ai − δai includes accelerometer errors and orientation errors − over-script, ~, denotes indicated (measured) quantity Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.6 Determination of δgn (2) • System state vector, εD d − orientation errors, ψ i; ψωi= −C ibδ dt b ib − IMU biases, scale-factor errors, assumed as random constants − gravity disturbance components, modeled, e.g., as second-order i in Gauss-Markov processes in n-frame, with δδgg= Cn ββ00 2 00 2 NN ddn nn2 δg=−20000 β δg−+ βδgw dt 2 Edt Egδ ββ2 00 DD00 ◦ wδg is a white noise vector process with appropriately selected variances ◦ βN , βE , βD are parameters appropriately selected to model the correlation time of the processes = 2.146 βNED,, • System dynamics equation d εε=FG + w where wD is a vector of white noise processes dt D DD D D Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.7 Determination of δgn (3) • Alternatively, omit gravity disturbance model − observations assume normal gravitation is correct yag=+−()ii x i i iii i i i i ii =()()aagx +δδ +− + x gx=γ + ΩΩie ie =iδδ b −×ψ i ib − i CCbba ax δ ai − optimal estimates of IMU systematic errors by Kalman filter yield yˆ − gravity disturbance estimates: δ gi ≈− yyˆ ◦ assumes residual IMU systematic errors are small and white noise can be filtered − successfully applied technique (Kwon and Jekeli, 2001) Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.8 Indirect Method • Recall inertial navigation equations in n-frame n dv nb n n n n =Cba −+()ΩΩie in vg + dt dxe = Cenv dt n − integrate (i.e., get IMU navigation solution) and solve for g n using a model and GNSS tracking data − analogous to traditional satellite tracking methods to determine global gravitational field, except gravity model is linear stochastic process, not spherical harmonic model − navigation solution from OTF INS should not be integrated with GNSS! ◦ IMU and GNSS must be treated as separate sensors, just like in scalar gravimetry ◦ use raw accelerometer and gyro data to obtain free-inertial navigation solution - one could pre-process IMU/GNSS data to solve for IMU systematic errors, neglecting gravity Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.9 Determination of δgn (1) • Solve for gravity disturbance treated as an error state (among many others) in the linear perturbation of navigation equations − typical error states collected in state vector, εΙ , include: ◦ position errors, velocity errors, orientation errors ◦ IMU systematic errors (biases, etc.) ◦ gravity disturbance components d εε=FG + w where wI is a vector of white noise processes dt I II I I − integration is done numerically (e.g., using linear finite differences) εεIk()()()()t= ΦG Ikk tt, −−11 Ik t+ Ik tw Ik() t where Φ = state transition matrix − observations are differences, IMU-indicated minus GNSS positions, treated as updates to the corresponding system states y()()()()tk= H tt kεv Ik+ t kwhere v is a vector of discrete white noise processes Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.10 Determination of δgn (2) • Gravity disturbance model − stochastic process; e.g., second-order Gauss-Markov process, ββ00 2 00 2 NN ddn nn2 δg=−20000 β δg −+ βδgw dt 2 Edt Egδ ββ2 00 DD00 ◦ wδg is a white noise vector process with appropriately selected variances ◦ βN , βE , βD are parameters appropriately selected to model the correlation time of the processes = 2.146 βNED,, • Kalman filter/smoother estimate, is optimal in the sense of minimum mean square error − theoretically the gravity model is an approximation since the gravity field is not a linear, finite-dimensional, set of independent along-track signals as required/modeled in the system state formalism − successful estimation depends on stochastic separability of gravity disturbance from accelerometer errors and coupled gyro errors Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.11 Rudimentary Error analysis • Assume that n- and i-frames coincide (approx. valid for < 1 hour) negligible n nn n n n nb n b nn g≈−⇒ xa δδgx=+−− ΨCbb a Cδ aΓ δ p gravitation errors in errors in errors in gravitational position errors kinematic sensor inertial gradients errors acceleration orientation acceleration • δg PSD is obtained from models of IMU and position error PSDs and PSD of vehicle acceleration ΦΦΦΦΦ=+++ δg1 δψ x1 32 a ψ 23 aa δ1 ΦΦΦΦΦ=+++ δg3 δψ x3 21 a ψ 12 aa δ3 Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.12 PSD Models Gravitational Field Accelerations* MathCad: accpsdfit.mcd MathCad: degvarfit4.mcd EGM96 dashed-line: psd model Φ a ] 2 2 1'x1' ∆g /(cy/m) /Hz] 2 2 ) 2 model Φ a3 [(m/s ) [mgal f ( g δ Φ Φ a1 frequency [cy/m] frequency [Hz] * data from Twin Otter aircraft – Accelerometers: white noise (±25 mGal) – Gyros: rate bias (± 0.003°/hr) plus white noise (±0.06 °/hr/√Hz) approximated by simple PSD ° – Orientation: initial bias (0.005 ) models − Position: white noise (±0.1 m) − Aircraft speed = 250 km/hr; altitude = 1000 m Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C. Jekeli, OSU 5.13 PSDs of Gravity Errors, δg1, δg3, Versus Along-Track PSDs of Signals, g1, g3 0.01 3 1 10 Spectral Window: 4 S 1 10 g3 Φ δ x1,3 5 S 1 10 g1 signal-to-noise ratio > 1 ]/Hz ]/Hz 4 /s 6 2 1 10 g : [m 1 7 1 10 Φ δg1 8 1 10 g3 : + Φδg 9 3 1 10 4 3 1 10 1 10 0.01 0.1 Frequency [Hz] • Signal-to-noise ratio > 1 over larger bandwidth for g3 • Significant error source is orientation bias, especially for g1, g2 • Signal PSD’s move to right and down with increased velocity Airborne Gravity for Geodesy Summer School , 23-27 May 2016 Theoretical Fundamentals of Inertial Gravimetry, C.
Recommended publications
  • Wireless Accelerometer - G-Force Max & Avg
    The Leader in Low-Cost, Remote Monitoring Solutions Wireless Accelerometer - G-Force Max & Avg General Description Monnit Sensor Core Specifications The RF Wireless Accelerometer is a digital, low power, • Wireless Range: 250 - 300 ft. (non line-of-sight / low profile, capacitive sensor that is able to measure indoors through walls, ceilings & floors) * acceleration on three axes. Four different • Communication: RF 900, 920, 868 and 433 MHz accelerometer types are available from Monnit. • Power: Replaceable batteries (optimized for long battery life - Line-power (AA version) and solar Free iMonnit basic online wireless (Industrial version) options available sensor monitoring and notification • Battery Life (at 1 hour heartbeat setting): ** system to configure sensors, view data AA battery > 4-8 years and set alerts via SMS text and email. Coin Cell > 2-3 years. Industrial > 4-8 years Principle of Operation * Actual range may vary depending on environment. ** Battery life is determined by sensor reporting Accelerometer samples at 800 Hz over a 10 second frequency and other variables. period, and reports the measured MAXIMUM value for each axis in g-force and the AVERAGE measured g-force on each axis over the same period, for all three axes. (Only available in the AA version.) This sensor reports in every 10 seconds with this data. Other sampling periods can be configured, down to one second and up to 10 minutes*. The data reported is useful for tracking periodic motion. Sensor data is displayed as max and average. Example: • Max X: 0.125 Max Y: 1.012 Max Z: 0.015 • Avg X: 0.110 Avg Y: 1.005 Avg Z: 0.007 * Customer cannot configure sampling period on their own.
    [Show full text]
  • Measurement of the Earth Tides with a MEMS Gravimeter
    1 Measurement of the Earth tides with a MEMS gravimeter ∗ y ∗ y ∗ ∗ ∗ 2 R. P. Middlemiss, A. Samarelli, D. J. Paul, J. Hough, S. Rowan and G. D. Hammond 3 The ability to measure tiny variations in the local gravitational acceleration allows { amongst 4 other applications { the detection of hidden hydrocarbon reserves, magma build-up before volcanic 5 eruptions, and subterranean tunnels. Several technologies are available that achieve the sensi- p 1 6 tivities required for such applications (tens of µGal= Hz): free-fall gravimeters , spring-based 2; 3 4 5 7 gravimeters , superconducting gravimeters , and atom interferometers . All of these devices can 6 8 observe the Earth Tides ; the elastic deformation of the Earth's crust as a result of tidal forces. 9 This is a universally predictable gravitational signal that requires both high sensitivity and high 10 stability over timescales of several days to measure. All present gravimeters, however, have limita- 11 tions of excessive cost (> $100 k) and high mass (>8 kg). We have built a microelectromechanical p 12 system (MEMS) gravimeter with a sensitivity of 40 µGal= Hz in a package size of only a few 13 cubic centimetres. We demonstrate the remarkable stability and sensitivity of our device with a 14 measurement of the Earth tides. Such a measurement has never been undertaken with a MEMS 15 device, and proves the long term stability of our instrument compared to any other MEMS device, 16 making it the first MEMS accelerometer to transition from seismometer to gravimeter. This heralds 17 a transformative step in MEMS accelerometer technology.
    [Show full text]
  • A GUIDE to USING FETS for SENSOR APPLICATIONS by Ron Quan
    Three Decades of Quality Through Innovation A GUIDE TO USING FETS FOR SENSOR APPLICATIONS By Ron Quan Linear Integrated Systems • 4042 Clipper Court • Fremont, CA 94538 • Tel: 510 490-9160 • Fax: 510 353-0261 • Email: [email protected] A GUIDE TO USING FETS FOR SENSOR APPLICATIONS many discrete FETs have input capacitances of less than 5 pF. Also, there are few low noise FET input op amps Linear Systems that have equivalent input noise voltages density of less provides a variety of FETs (Field Effect Transistors) than 4 nV/ 퐻푧. However, there are a number of suitable for use in low noise amplifier applications for discrete FETs rated at ≤ 2 nV/ 퐻푧 in terms of equivalent photo diodes, accelerometers, transducers, and other Input noise voltage density. types of sensors. For those op amps that are rated as low noise, normally In particular, low noise JFETs exhibit low input gate the input stages use bipolar transistors that generate currents that are desirable when working with high much greater noise currents at the input terminals than impedance devices at the input or with high value FETs. These noise currents flowing into high impedances feedback resistors (e.g., ≥1MΩ). Operational amplifiers form added (random) noise voltages that are often (op amps) with bipolar transistor input stages have much greater than the equivalent input noise. much higher input noise currents than FETs. One advantage of using discrete FETs is that an op amp In general, many op amps have a combination of higher that is not rated as low noise in terms of input current noise and input capacitance when compared to some can be converted into an amplifier with low input discrete FETs.
    [Show full text]
  • Preliminary Estimation and Validation of Polar Motion Excitation from Different Types of the GRACE and GRACE Follow-On Missions Data
    remote sensing Article Preliminary Estimation and Validation of Polar Motion Excitation from Different Types of the GRACE and GRACE Follow-On Missions Data Justyna Sliwi´ ´nska 1,* , Małgorzata Wi ´nska 2 and Jolanta Nastula 1 1 Space Research Centre, Polish Academy of Sciences, 00-716 Warsaw, Poland; [email protected] 2 Faculty of Civil Engineering, Warsaw University of Technology, 00-637 Warsaw, Poland; [email protected] * Correspondence: [email protected] Received: 18 September 2020; Accepted: 22 October 2020; Published: 23 October 2020 Abstract: The Gravity Recovery and Climate Experiment (GRACE) mission has provided global observations of temporal variations in the gravity field resulting from mass redistribution at the surface and within the Earth for the period 2002–2017. Although GRACE satellites are not able to realistically detect the second zonal parameter (DC20) of geopotential associated with the flattening of the Earth, they can accurately determine variations in degree-2 order-1 (DC21, DS21) coefficients that are proportional to variations in polar motion. Therefore, GRACE measurements are commonly exploited to interpret polar motion changes due to variations in the global mass redistribution, especially in the continental hydrosphere and cryosphere. Such impacts are usually examined by computing the so-called hydrological polar motion excitation (HAM) and cryospheric polar motion excitation (CAM), often analyzed together as HAM/CAM. The great success of the GRACE mission and the scientific robustness of its data contributed to the launch of its successor, GRACE Follow-On (GRACE-FO), which began in May 2018 and continues to the present. This study presents the first estimates of HAM/CAM computed from GRACE-FO data provided by three data centers: Center for Space Research (CSR), Jet Propulsion Laboratory (JPL), and GeoForschungsZentrum (GFZ).
    [Show full text]
  • Utilising Accelerometer and Gyroscope in Smartphone to Detect Incidents on a Test Track for Cars
    Utilising accelerometer and gyroscope in smartphone to detect incidents on a test track for cars Carl-Johan Holst Data- och systemvetenskap, kandidat 2017 Luleå tekniska universitet Institutionen för system- och rymdteknik LULEÅ UNIVERSITY OF TECHNOLOGY BACHELOR THESIS Utilising accelerometer and gyroscope in smartphone to detect incidents on a test track for cars Author: Examiner: Carl-Johan HOLST Patrik HOLMLUND [email protected] [email protected] Supervisor: Jörgen STENBERG-ÖFJÄLL [email protected] Computer and space technology Campus Skellefteå June 4, 2017 ii Abstract Utilising accelerometer and gyroscope in smartphone to detect incidents on a test track for cars Every smartphone today includes an accelerometer. An accelerometer works by de- tecting acceleration affecting the device, meaning it can be used to identify incidents such as collisions at a relatively high speed where large spikes of acceleration often occur. A gyroscope on the other hand is not as common as the accelerometer but it does exists in most newer phones. Gyroscopes can detect rotations around an arbitrary axis and as such can be used to detect critical rotations. This thesis work will present an algorithm for utilising the accelerometer and gy- roscope in a smartphone to detect incidents occurring on a test track for cars. Sammanfattning Utilising accelerometer and gyroscope in smartphone to detect incidents on a test track for cars Alla smarta telefoner innehåller idag en accelerometer. En accelerometer analyserar acceleration som påverkar enheten, vilket innebär att den kan användas för att de- tektera incidenter så som kollisioner vid relativt höga hastigheter där stora spikar av acceleration vanligtvis påträffas.
    [Show full text]
  • Detection of Caves by Gravimetry
    Detection of Caves by Gravimetry By HAnlUi'DO J. Cmco1) lVi/h plates 18 (1)-21 (4) Illtroduction A growing interest in locating caves - largely among non-speleolo- gists - has developed within the last decade, arising from industrial or military needs, such as: (1) analyzing subsUl'face characteristics for building sites or highway projects in karst areas; (2) locating shallow caves under airport runways constructed on karst terrain covered by a thin residual soil; and (3) finding strategic shelters of tactical significance. As a result, geologists and geophysicists have been experimenting with the possibility of applying standard geophysical methods toward void detection at shallow depths. Pioneering work along this line was accomplishecl by the U.S. Geological Survey illilitary Geology teams dUl'ing World War II on Okinawan airfields. Nicol (1951) reported that the residual soil covel' of these runways frequently indicated subsi- dence due to the collapse of the rooves of caves in an underlying coralline-limestone formation (partially detected by seismic methods). In spite of the wide application of geophysics to exploration, not much has been published regarding subsUl'face interpretation of ground conditions within the upper 50 feet of the earth's surface. Recently, however, Homberg (1962) and Colley (1962) did report some encoUl'ag- ing data using the gravity technique for void detection. This led to the present field study into the practical means of how this complex method can be simplified, and to a use-and-limitations appraisal of gravimetric techniques for speleologic research. Principles all(1 Correctiolls The fundamentals of gravimetry are based on the fact that natUl'al 01' artificial voids within the earth's sUl'face - which are filled with ail' 3 (negligible density) 01' water (density about 1 gmjcm ) - have a remark- able density contrast with the sUl'roun<ling rocks (density 2.0 to 1) 4609 Keswick Hoad, Baltimore 10, Maryland, U.S.A.
    [Show full text]
  • Airborne Geoid Determination
    LETTER Earth Planets Space, 52, 863–866, 2000 Airborne geoid determination R. Forsberg1, A. Olesen1, L. Bastos2, A. Gidskehaug3,U.Meyer4∗, and L. Timmen5 1KMS, Geodynamics Department, Rentemestervej 8, 2400 Copenhagen NV, Denmark 2Astronomical Observatory, University of Porto, Portugal 3Institute of Solid Earth Physics, University of Bergen, Norway 4Alfred Wegener Institute, Bremerhaven, Germany 5Geo Forschungs Zentrum, Potsdam, Germany (Received January 17, 2000; Revised August 18, 2000; Accepted August 18, 2000) Airborne geoid mapping techniques may provide the opportunity to improve the geoid over vast areas of the Earth, such as polar areas, tropical jungles and mountainous areas, and provide an accurate “seam-less” geoid model across most coastal regions. Determination of the geoid by airborne methods relies on the development of airborne gravimetry, which in turn is dependent on developments in kinematic GPS. Routine accuracy of airborne gravimetry are now at the 2 mGal level, which may translate into 5–10 cm geoid accuracy on regional scales. The error behaviour of airborne gravimetry is well-suited for geoid determination, with high-frequency survey and downward continuation noise being offset by the low-pass gravity to geoid filtering operation. In the paper the basic principles of airborne geoid determination are outlined, and examples of results of recent airborne gravity and geoid surveys in the North Sea and Greenland are given. 1. Introduction the sea-surface (H) by airborne altimetry. This allows— Precise geoid determination has in recent years been facil- in principle—the determination of the dynamic sea-surface itated through the progress in airborne gravimetry. The first topography (ζ) through the equation large-scale aerogravity experiment was the airborne gravity survey of Greenland 1991–92 (Brozena, 1991).
    [Show full text]
  • Basic Principles of Inertial Navigation
    Basic Principles of Inertial Navigation Seminar on inertial navigation systems Tampere University of Technology 1 The five basic forms of navigation • Pilotage, which essentially relies on recognizing landmarks to know where you are. It is older than human kind. • Dead reckoning, which relies on knowing where you started from plus some form of heading information and some estimate of speed. • Celestial navigation, using time and the angles between local vertical and known celestial objects (e.g., sun, moon, or stars). • Radio navigation, which relies on radio‐frequency sources with known locations (including GNSS satellites, LORAN‐C, Omega, Tacan, US Army Position Location and Reporting System…) • Inertial navigation, which relies on knowing your initial position, velocity, and attitude and thereafter measuring your attitude rates and accelerations. The operation of inertial navigation systems (INS) depends upon Newton’s laws of classical mechanics. It is the only form of navigation that does not rely on external references. • These forms of navigation can be used in combination as well. The subject of our seminar is the fifth form of navigation – inertial navigation. 2 A few definitions • Inertia is the property of bodies to maintain constant translational and rotational velocity, unless disturbed by forces or torques, respectively (Newton’s first law of motion). • An inertial reference frame is a coordinate frame in which Newton’s laws of motion are valid. Inertial reference frames are neither rotating nor accelerating. • Inertial sensors measure rotation rate and acceleration, both of which are vector‐ valued variables. • Gyroscopes are sensors for measuring rotation: rate gyroscopes measure rotation rate, and integrating gyroscopes (also called whole‐angle gyroscopes) measure rotation angle.
    [Show full text]
  • Full Auto-Calibration of a Smartphone on Board a Vehicle Using IMU and GPS Embedded Sensors
    Full auto-calibration of a smartphone on board a vehicle using IMU and GPS embedded sensors Javier Almaz´an, Luis M. Bergasa, J. Javier Yebes, Rafael Barea and Roberto Arroyo Abstract| Nowadays, smartphones are widely used Smartphones provide two kinds of measurements. The in the world, and generally, they are equipped with first ones are relative to the world, such as GPS and many sensors. In this paper we study how powerful magnetometers. The second ones are relative to the the low-cost embedded IMU and GPS could become for Intelligent Vehicles. The information given by device, such as accelerometers and gyroscopes. Knowing accelerometer and gyroscope is useful if the relations the relation between the smartphone reference system between the smartphone reference system, the vehicle and the world reference system is very important in order reference system and the world reference system are to reference the second kind of measurements globally [3]. known. Commonly, the magnetometer sensor is used In other words, working with this kind of measurements to determine the orientation of the smartphone, but its main drawback is the high influence of electro- involves knowing the pose of the smartphone in the magnetic interference. In view of this, we propose a world. novel automatic method to calibrate a smartphone Intelligent Vehicles can be helped by using in-vehicle on board a vehicle using its embedded IMU and smartphones to measure some driving indicators. On the GPS, based on longitudinal vehicle acceleration. To one hand, using smartphones requires no extra hardware the best of our knowledge, this is the first attempt to estimate the yaw angle of a smartphone relative to a mounted in the vehicle, offering cheap and standard vehicle in every case, even on non-zero slope roads.
    [Show full text]
  • Resolving Geophysical Signals by Terrestrial Gravimetry: a Time Domain Assessment of the Correction-Induced Uncertainty
    Originally published as: Mikolaj, M., Reich, M., Güntner, A. (2019): Resolving Geophysical Signals by Terrestrial Gravimetry: A Time Domain Assessment of the Correction-Induced Uncertainty. - Journal of Geophysical Research, 124, 2, pp. 2153—2165. DOI: http://doi.org/10.1029/2018JB016682 RESEARCH ARTICLE Resolving Geophysical Signals by Terrestrial Gravimetry: 10.1029/2018JB016682 A Time Domain Assessment of the Correction-Induced Key Points: Uncertainty • Global-scale uncertainty assessment of tidal, oceanic, large-scale 1 1 1,2 hydrological, and atmospheric M. Mikolaj , M. Reich , and A. Güntner corrections for terrestrial gravimetry • Resolving subtle gravity signals in the 1Section Hydrology, GFZ German Research Centre for Geosciences, Potsdam, Germany, 2Institute of Environmental order of few nanometers per square Science and Geography, University of Potsdam, Potsdam, Germany second is challenged by the statistical uncertainty of correction models • Uncertainty computed for selected periods varies significantly Abstract Terrestrial gravimetry is increasingly used to monitor mass transport processes in with latitude and altitude of the geophysics boosted by the ongoing technological development of instruments. Resolving a particular gravimeter phenomenon of interest, however, requires a set of gravity corrections of which the uncertainties have not been addressed up to now. In this study, we quantify the time domain uncertainty of tide, global Supporting Information: atmospheric, large-scale hydrological, and nontidal ocean loading corrections. The uncertainty is • Supporting Information S1 assessed by comparing the majority of available global models for a suite of sites worldwide. The average uncertainty expressed as root-mean-square error equals 5.1 nm/s2, discounting local hydrology or air Correspondence to: pressure. The correction-induced uncertainty of gravity changes over various time periods of interest M.
    [Show full text]
  • Macroscopic Superpositions and Gravimetry with Quantum Magnetomechanics Received: 07 July 2016 Mattias T
    www.nature.com/scientificreports OPEN Macroscopic superpositions and gravimetry with quantum magnetomechanics Received: 07 July 2016 Mattias T. Johnsson, Gavin K. Brennen & Jason Twamley Accepted: 24 October 2016 Precision measurements of gravity can provide tests of fundamental physics and are of broad Published: 21 November 2016 practical interest for metrology. We propose a scheme for absolute gravimetry using a quantum magnetomechanical system consisting of a magnetically trapped superconducting resonator whose motion is controlled and measured by a nearby RF-SQUID or flux qubit. By driving the mechanical massive resonator to be in a macroscopic superposition of two different heights our we predict that our interferometry protocol could, subject to systematic errors, achieve a gravimetric sensitivity of Δg/g ~ 2.2 × 10−10 Hz−1/2, with a spatial resolution of a few nanometres. This sensitivity and spatial resolution exceeds the precision of current state of the art atom-interferometric and corner-cube gravimeters by more than an order of magnitude, and unlike classical superconducting interferometers produces an absolute rather than relative measurement of gravity. In addition, our scheme takes measurements at ~10 kHz, a region where the ambient vibrational noise spectrum is heavily suppressed compared the ~10 Hz region relevant for current cold atom gravimeters. Gravimetry is the measurement of the local acceleration due to gravity on a test body. As gravity is an infinite range force that cannot be shielded against, it provides a useful probe for a variety of effects over a wide range of distances. Given a theory of gravity, gravimetry provides information about local mass distribution or, given a known, local mass distribution, it provides information that can be used test theories of gravity.
    [Show full text]
  • Pendulum MEMS Gravimeters for Semi-Absolute Gravimetry
    EGU2020-1244 https://doi.org/10.5194/egusphere-egu2020-1244 EGU General Assembly 2020 © Author(s) 2021. This work is distributed under the Creative Commons Attribution 4.0 License. Pendulum MEMS gravimeters for semi-absolute gravimetry Richard Middlemiss1, Giles Hammond2, Richard Walker2, and Abhinav Prasad2 1University of Glasgow, James Watt School of Engineering, Glasgow, United Kingdom of Great Britain and Northern Ireland ([email protected]) 2University of Glasgow, Institute for Gravitational Research, Glasgow, United Kingdom of Great Britain and Northern Ireland By measuring tiny variations in the Earth’s gravitational acceleration, g, one can infer density variations beneath the ground. Since magmatic systems contain rock of differing density, changes in gravity over time can tell us when/where magma is moving. Traditional gravity sensors (gravimeters) were costly and heavy, but with the advent of the technology used to make mobile phone accelerometers (MEMS – Microelectromechanical-systems), this is changing. At Glasgow University we have already developed the first MEMS gravity sensor and we are now working with several other European institutions to make a network of gravity sensors around Mt Etna – NEWTON-g. It will be the first multi-pixel gravity imager – enabling unprecedented resolution of Etna’s plumbing system. While this work is ongoing, a second generation of MEMS gravity sensor is now under development. The first-generation sensor comprises a mass on a spring, which moves in response to changing values of g. This, however, can only ever be used to measure changes in gravity, which means it can be difficult to tell the difference between a geophysical signal and instrumental drift.
    [Show full text]