GNSS Ocean Reflected Signals Per Høeg DTU Space Technical University of Denmark
Content Experimental setup Instrument Measurements and observations Spectral characteristics, analysis and retrieval method 2
Experiment Altitude: > 3000 meters Year: 2004 60 3
GNSS Ocean Reflected Signals GNSS ocean reflected signals describe the height and the roughness of the ocean The characteristics of the reflected signal depend on the scattering properties of the sea surface and the footprint of the reflection zone The footprint size and shape in turn depends on the signal incidence angle and the relative velocities of transmitter and receiver Scattering properties of the sea surface (the roughness parameter) relates to the ocean wave characteristics 60 4
Motivation Studies of Sea Surface Height Processes GEROS on the ISS Columbus Module 5
Ocean Sea Surface Height Processes Oceanic observations carry signals of a wide range of related processes. The observed fingerprints of these processes have temporal time scales from 1 hour to thousands of years and spatial scales from ten to tens thousands of kilometres. The figure illustrates the spatial and temporal scales for these processes and indicates phenomena, which can be investigated with GEROS- ISS data complementary to and distinct from the planned NASA SWOT mission and ESA and NASA radar altimetry missions (Revised from Zuffada et al., 2005). 6
Power Spectrum of Sea Surface Heights (SSH) 500 km The black SSH power spectrum is for reference based on Jason altimeter observations (pass 132). The red curve gives the error spectrum of the NASA SWOT mission. The solid black line is the expected spectral continuation. The intersection of the spectral signal with the noise floor at 10km determines the resolving capabilities for the SWOT instrument (JPL, NASA, 2009). 7
Elevation (deg) Elevation (deg) Experiment Azimuth-Elevation (2004-10-4 19:34:0 p:\p-roft\m-functions\almanacs\yuma267.txt) 30 26 19 10 4 16 25 25 20 Elevation (2004-10-5 0:0:0 p:\p-roft\m-functions\almanacs\yuma267.txt) 30 24 15 26 25 31 10 19 20 14 5 17 16 15 7 0 4-5 -10 1 21 27 22 8 23 28 15 13 6 20 5 10 5 10 6 30 18 21 15 25 22 1 23 13 8 27 28 24-175 -170-165 -160-155 -150-145 -140-135 -130-125 -120 Azimuth relative North (deg) 0 20 31 Azimuth (deg) -5 14 7 5 17 18 30-10 0 5 10 15 20 Time UTC (hours) Time of day (hours) 8
Experiment Antenna Front-End Command and Data Storage Receiver 9
Tangent height (km) Doppler [km/s] Instrument CODE NCO / GEN LOOP FILTER Dot Prod. Discrim. Front End Code Phase E/P/L E/P/L Punctual I/Q 1,10ms OL: 1 khz Meas. Packaging Code Phase CCSDS Data Carrier NCO Carrier Phase CARRIER NCO CL OL LOOP FILTER arctan Discrim. 30 25 20 15 10 5 6.9 6.85 6.8 6.75 6.7 6.65 6.6 6.55 6.5 6.45 0 20 40 60 80 100 120 Time [s] ON-BOARD DOPPLER MODEL Closed Loop tracking: Carrier Phase is phase locked to the received signal. Open Loop tracking (raw mode): Carrier Phase is measured relative on-board Doppler model 0 0 200 400 600 800 1000 1200 Atmospheric Doppler (Hz) 10
Amplitude Q Amplitude Amplitude [-] Instrument (Navigation Data De-modulation) 0.2 0.1 0.3 0.2 0.1 20 ms 20 ms 0 24.74 24.76 24.78 24.8 24.82 24.84 24.86 24.88 24.9 24.92 Time [s] 4 x 104 200 3 100 Phase [ l ] 0 24.76 24.78 24.8 24.82 24.84 24.86 24.88 24.9 24.92 1 0.5 0-0.5-1 24.76 24.78 24.8 24.82 24.84 24.86 24.88 24.9 24.92 Time [s] 2 1 0-50 0 50 Frequency [Hz] 0-100 -200-200 0 200 I 11
Measurements (De-trended Power Spectra) Background signal frequency drift are removed using a least squares fit parabola to the main phase. 12
Measurements (Power Spectra at Different Elevation Angles) Power spectra as function of frequency difference (f-f 0 ) from the main signal peak f 0. Clock related spectral slope of f -2. No clock correction has been applied. 13
Measurements (Trend Analysis) Variation of the slope in the frequency region 1-5 Hz of the power spectrum, as function of elevation angle. The red straight line is the linear least squares fit to the curve in the elevation angle interval, 0-5 degrees. 14
Spectral Analysis Direct signal h/sin cos2 Ocean reflected signal h/sin h Path length difference: h / sin cos2 15
Spectral Analysis 16
Pressure (mbar) Temperature (deg C) Wind speed (m/s) Meteorological Conditions 12 10 8 6 4 2 0 2 4 6 8 10 12 14 Day of October 2004 29 28 27 26 25 24 2 4 6 8 10 12 14 Day of October 2004 1016 1014 1012 1010 2 4 6 8 10 12 14 Day of October 2004 17
Particle Filtering Method The State Model Tracking the position of a reflection point can be viewed as estimation of the state of the dynamical system based on sets of measurements. We assume that the state (parameterized by its position, height and velocity) is a firstorder Markov process of the form (time-varying stochastic process): F G W k k k 1 k k 1 K is the reflection point state vector at the time k θ k = p k W k is the unknown system noise input v k Here, the state vector consists of the 3-dimensional position and velocity vectors The discrete time matrices: 2 t IM t IM IM Fk and Gk 2 OM I M t IM 18
Particle Filtering Method The Observation Model The system has M received signals and a single GNSS transmitter. The time-variant m th received signal in the k th snapshot is given as: y ( t) x ( t) z ( t) ( t) m, k m, k m, k m, k The signal component: x t t u t t e ( j2 f ( t)) c m, k ( ) ( ) ( ( )) m, k m, k m, k, () t 2 mk is the complex zero-mean Gaussian noise process with spectral height mk, () t f c denotes the complex amplitude of the wave propagation path carrier frequency, ut () the transmitted signal the delay of the signal received from the m th received signal at time t in the k mk, () t th observation window 19
Particle Filtering Method The Particle Filtering Algorithm A sequential Bayesian estimation method is used to estimate recursively the conditional probability density function p y y 1: k 1: k and denote respectively, a sequence of observations and state vectors, from the 1 st to the k th measurement cycle The proposed particle filter technique makes use of a fixed number of particles, where each particle is associated with a state vector Steps for each new observation: 1) Predict the states of particles and calculate the weights 2) Re-sampling 3) Estimate the power density function i k i p k i vk i k i k 20
Sea Surface Reflection Zone 200 100 0-100 -200-200 -100 0 100 100-200 -100 0 100 100 Horizontal probability density functions at the reflection location (in meters) 21
Conclusions The spectral variances are driven by the atmospheric physical conditions, sea surface roughness and wind dynamics. Spectral noise characteristics follows a power distribution related to the clock noise and the receiver amplitude estimation. Estimated spectral variances link directly to the turbulence structure function constant of the measured atmosphere region. 22
Thank you! 23