ssha in the Solomon Sea

by olvac last modified Mar 26, 2015 04:45 PM


Bathymetry and along track ssha data

  • Saral/Altika (S/A) along track in the Solomon Sea

      The ssha data (1Hz) are collected from CTOH.

bathy_plus_tracks  

  • How is ssha computed?

     ssha = altitude - range - correction - MSSH_CNES11 

 

 

Procedure of reprocessing of S/A data:

    0. Elimination de données dont la valeur absolue >30cm

    1. Elimination des données localisées sur une bathy <100m

    2. Elimination des données dont la rms > 2.8*sigma (=rms)

   

Les figures suivantes présentent la ssha issue de l'application des procedures 0, 1 et 2 aux données brutes de ssha de S/A. En orange la ssha du premier cycle et en rouge les autres cycles (en trait plein rouge indique le  dernier cycle). Rappelons qu'il ya 16 cycles traités. En noir la bathymétrie (modèle) masqué à 100m.

 

 

Track 0776 ;    Track 0318 ;     Track 0490 ;      Track 0862 ;      Track 0948 ;     Track 0032 ;     Track 0576 ; Track 0247 ;    

Track 0705 ;     Track 0075 ;     Track 0161 ;     Track 0619 ;      Track 0533 ;     Track 0991 ;     Track 0447 

 

 

Quelques traces avec l'élimination à 4*sigma

Track 0776 ;    Track 0318 ;   Track 0075

 

 

Tableau récapitulatif

a. Tableau présentant les plot des traces de tous les cycles (16 cycles) avec un code couleur par cycle en considération de l'élimination à 4*sigma et 2.8*sigma. Sur ces plots (pour chaque trace) montrent deux sous figures présentant 8 cycles chacune (sauf la trace 0032 qui présente 15 cycle au lieu de 16 car il n'y a pas cette trace au 1er cycle dans la mer des Salomon).

Sigma

Tracks

C1-8/9-16

4*sigma 0776 0318 0490 0862 0948 0032 0576 0247 0705 0075 0161 0619 0533 0991 0447
2.8*sigma 0776 0318 0490 0862 0948 0032 0576 0247 0705 0075 0161 0619 0533 0991 0447

 

 

 

 b. Tableau présentant les plot des traces de tous les cycles (16 cycles) avec un code couleur par cycle en considération de l'élimination à 4*sigma et 2.8*sigma. La différence avec les figure précédentes du tableau a est que les 16 cycles sont mis sur des figures séparées (et non tous à la fois). Donc chaque plot contient 16 sous figures (sauf la trace 0032 qui présente 15) présentant 1 cycle chacune. C'est pour mieux suivre la variabilité temporelle.

 

Sigma

Tracks

C1-16

4*sigma 0776 0318 0490 0862 0948 0032 0576 0247 0705 0075 0161 0619 0533 0991 0447
2.8*sigma 0776 0318 0490 0862 0948 0032 0576 0247 0705 0075 0161 0619 0533 0991 0447

 

 

 

c. Tableau présentant les plots des traces de la donnée brute (ie. qui n'ont subit aucun pré-filtrage)

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

 

d. Carte représentant la RMS de la SLA issu des produits AVISO de la période 1993-2006.

      i.  RMS AVISO sans trace S/A

      ii. RMS AVISO avec trace S/A

e. Tableau présentant les plots des traces de la donnée brute (ie. qui n'ont subit aucun pré-filtrage) masquée avecune bathy à 100m de profondeur. Le trait en rouge bourdeau désigne la RMS de la SLA issu des produits AVISO de la période 1993-2006.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

f. Tableau présentant les données de traces avec eliminations des points < 100m de profondeur (bathy) + 2.8*rms (rms issus de aviso 1993-2006).

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

g. Tableau présentant les données de traces avec eliminations des points < 100m de profondeur (bathy) + 4*sigma*rms (aviso)

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

 ----------------------------------------------------------------------------------------------------------------------12/12/2014

h. Tableau présentant les données de traces avec eliminations des points < 100m de profondeur (bathy) + 4*sigma*rms (aviso) + filtre median (nombre de point considéré: 15).

           h.i. La premiere colonne représente la SSHA filtré et la deuxième colonne représente la différence entre la SSHA (bathy 100m + 4 sigma)  et celle filtrée par le filtre median.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

          h.ii. Même plot que ceux dans le tableau h.i. sauf qu'il n'y a pas la deuxième colonne et aussi chaque cycle est dessiné  individuellement. Courbes en gras (couleur foncée) correspondent aux données filtré par le filtre median et celle en trait fin les données non filtrées.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

I. Tableau presque identique à celui en h sauf qu'avant d'appliquer le filtre median, les valeurs masquées (dues au masque appliqué par la bathy à 100m et 4*sigma) ont été supprimées afin d'avoir un champ de données continue. Les valeurs masquées sont remises après, une fois le filtre médian calculé.

           I.i. La premiere colonne représente la SSHA filtré et la deuxième colonne représente la différence entre la SSHA (bathy 100m + 4 sigma)  et celle filtrée par le filtre median.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

            I.ii. Même plot que I.i. sauf qu'il n'y a pas la deuxième colonne et aussi chaque cycle est dessiné individuellement. Courbes en gras (couleur foncée) correspondent aux données filtré par le filtre median et celle en trait fin les données non filtrées.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

-----------------------------------------------------------------------------------------------------------------15/12/2014

j. Tableau présentant le signal brut (trait fin) et le signal final (gras)  (masque 100m bathy + 4*rms AVISO + masque abs((x(i)-xm(i))) < alpha * AMLD(i) ); avec s(i) le signal ayant été masqué par  100m bathy + 4*rms AVISO ;

     xm(i) le signal filtré par le filtre médian (= medfiltre(x(i), 15).

     alpha = 10

     AMLD(i) est le filtre médian de (abs(x(i)-xm(i),m): c'est l'équivalent de l'écart-type sur m points, mais sans les spikes (m=51)

                   j.a. Figures 8 par 8

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

                jb. Même plot que le tableau j.a sauf qu'il a les plots chaque cycle présenté individuellement

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

k. Même plot que le tableau j. sauf qu'il n'y a pas les plots du signal brut. Il n'y a que le signal final traité (voir l'explication dans le titre du tableau j.).

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

l. Les vitesses géostrophiques calculées à partir des signaux finaux obtenus par la méthode décrite en j.

  • Tableau présentant les vitesses géostrophiques cycle par cycle (pour chaque trace)

 

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

 

---------------------------------------------------------------------------------------------------------------------------------------------------------

NEW PART: RE-CALCUL de SSHA DES DONNES SARAL

B. Recalcul de la SSHA et SSH en tenant compte que des corrections géophysiques

Formule: ssha = alt - range - iono_corr_gim - model_dry_tropo_corr - rad_wet_tropo_corr - sea_state_bias - solid_earth_tide - ocean_tide_sol1 - pole_tide - inv_bar_corr - hf_fluctuations_corr - mean_sea_surface

SSH = ssha + mean_sea_surface

 

 i.

Trackscycles
0776 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
905 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
318 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
361 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

 

ii. SSHA recalculée et Ug (=sla+mdt) issues de mssh PM Envisat

 Tracks77690536131886294849032576247
705 161 619 75 533 991 447 232 690

 

iii. Ug  = sla+mdt (en utilisant PM Envisat dans le calcul de la sla)

Procedure of reprocessing of S/A data:

    0. Elimination de données dont la valeur absolue >30cm

    1. Elimination des données localisées sur une bathy <100m

    2. Elimination des données dont la rms > 2.8*sigma (=rms)

   

Les figures suivantes présentent la ssha issue de l'application des procedures 0, 1 et 2 aux données brutes de ssha de S/A. En orange la ssha du premier cycle et en rouge les autres cycles (en trait plein rouge indique le  dernier cycle). Rappelons qu'il ya 16 cycles traités. En noir la bathymétrie (modèle) masqué à 100m.

 

 

Track 0776 ;    Track 0318 ;     Track 0490 ;      Track 0862 ;      Track 0948 ;     Track 0032 ;     Track 0576 ; Track 0247 ;    

Track 0705 ;     Track 0075 ;     Track 0161 ;     Track 0619 ;      Track 0533 ;     Track 0991 ;     Track 0447 

 

 

Quelques traces avec l'élimination à 4*sigma

Track 0776 ;    Track 0318 ;   Track 0075

 

 

Tableau récapitulatif

a. Tableau présentant les plot des traces de tous les cycles (16 cycles) avec un code couleur par cycle en considération de l'élimination à 4*sigma et 2.8*sigma. Sur ces plots (pour chaque trace) montrent deux sous figures présentant 8 cycles chacune (sauf la trace 0032 qui présente 15 cycle au lieu de 16 car il n'y a pas cette trace au 1er cycle dans la mer des Salomon).

Sigma

Tracks

C1-8/9-16

4*sigma 0776 0318 0490 0862 0948 0032 0576 0247 0705 0075 0161 0619 0533 0991 0447
2.8*sigma 0776 0318 0490 0862 0948 0032 0576 0247 0705 0075 0161 0619 0533 0991 0447

 

 

 

 b. Tableau présentant les plot des traces de tous les cycles (16 cycles) avec un code couleur par cycle en considération de l'élimination à 4*sigma et 2.8*sigma. La différence avec les figure précédentes du tableau a est que les 16 cycles sont mis sur des figures séparées (et non tous à la fois). Donc chaque plot contient 16 sous figures (sauf la trace 0032 qui présente 15) présentant 1 cycle chacune. C'est pour mieux suivre la variabilité temporelle.

 

Sigma

Tracks

C1-16

4*sigma 0776 0318 0490 0862 0948 0032 0576 0247 0705 0075 0161 0619 0533 0991 0447
2.8*sigma 0776 0318 0490 0862 0948 0032 0576 0247 0705 0075 0161 0619 0533 0991 0447

 

 

 

c. Tableau présentant les plots des traces de la donnée brute (ie. qui n'ont subit aucun pré-filtrage)

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

 

d. Carte représentant la RMS de la SLA issu des produits AVISO de la période 1993-2006.

      i.  RMS AVISO sans trace S/A

      ii. RMS AVISO avec trace S/A

e. Tableau présentant les plots des traces de la donnée brute (ie. qui n'ont subit aucun pré-filtrage) masquée avecune bathy à 100m de profondeur. Le trait en rouge bourdeau désigne la RMS de la SLA issu des produits AVISO de la période 1993-2006.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

f. Tableau présentant les données de traces avec eliminations des points < 100m de profondeur (bathy) + 2.8*rms (rms issus de aviso 1993-2006).

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

g. Tableau présentant les données de traces avec eliminations des points < 100m de profondeur (bathy) + 4*sigma*rms (aviso)

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

 ----------------------------------------------------------------------------------------------------------------------12/12/2014

h. Tableau présentant les données de traces avec eliminations des points < 100m de profondeur (bathy) + 4*sigma*rms (aviso) + filtre median (nombre de point considéré: 15).

           h.i. La premiere colonne représente la SSHA filtré et la deuxième colonne représente la différence entre la SSHA (bathy 100m + 4 sigma)  et celle filtrée par le filtre median.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

          h.ii. Même plot que ceux dans le tableau h.i. sauf qu'il n'y a pas la deuxième colonne et aussi chaque cycle est dessiné  individuellement. Courbes en gras (couleur foncée) correspondent aux données filtré par le filtre median et celle en trait fin les données non filtrées.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

I. Tableau presque identique à celui en h sauf qu'avant d'appliquer le filtre median, les valeurs masquées (dues au masque appliqué par la bathy à 100m et 4*sigma) ont été supprimées afin d'avoir un champ de données continue. Les valeurs masquées sont remises après, une fois le filtre médian calculé.

           I.i. La premiere colonne représente la SSHA filtré et la deuxième colonne représente la différence entre la SSHA (bathy 100m + 4 sigma)  et celle filtrée par le filtre median.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

            I.ii. Même plot que I.i. sauf qu'il n'y a pas la deuxième colonne et aussi chaque cycle est dessiné individuellement. Courbes en gras (couleur foncée) correspondent aux données filtré par le filtre median et celle en trait fin les données non filtrées.

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

-----------------------------------------------------------------------------------------------------------------15/12/2014

j. Tableau présentant le signal brut (trait fin) et le signal final (gras)  (masque 100m bathy + 4*rms AVISO + masque abs((x(i)-xm(i))) < alpha * AMLD(i) ); avec s(i) le signal ayant été masqué par  100m bathy + 4*rms AVISO ;

     xm(i) le signal filtré par le filtre médian (= medfiltre(x(i), 15).

     alpha = 10

     AMLD(i) est le filtre médian de (abs(x(i)-xm(i),m): c'est l'équivalent de l'écart-type sur m points, mais sans les spikes (m=51)

                   j.a. Figures 8 par 8

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

                jb. Même plot que le tableau j.a sauf qu'il a les plots chaque cycle présenté individuellement

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

k. Même plot que le tableau j. sauf qu'il n'y a pas les plots du signal brut. Il n'y a que le signal final traité (voir l'explication dans le titre du tableau j.).

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

l. Les vitesses géostrophiques calculées à partir des signaux finaux obtenus par la méthode décrite en j.

  • Tableau présentant les vitesses géostrophiques cycle par cycle (pour chaque trace)

 

Tracks07760318049008620948003205760247
Tracks 0705 0075 0161 0619 0533 0991 0447

 

 

 

---------------------------------------------------------------------------------------------------------------------------------------------------------

NEW PART: RE-CALCUL de SSHA DES DONNES SARAL                                                      --------------27/01/215

B. Recalcul de la SSHA et SSH en tenant compte que des corrections géophysiques

Formule: ssha = alt - range - iono_corr_gim - model_dry_tropo_corr - rad_wet_tropo_corr - sea_state_bias - solid_earth_tide - ocean_tide_sol1 - pole_tide - inv_bar_corr - hf_fluctuations_corr - mean_sea_surface

SSH = ssha + MDT

 

 i.

Trackscycles
0776 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
905 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
318 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
361 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

 

ii. SSHA recalculée en utilisant les corrections geophysiques et la mssh de CLS11 qui a intégré ~20 ans de données. Ensuite cette SSHA a été filtré par le filtre median evoqué plus haut issues de mssh PM Envisat

    ii.a Signal de SSHA recalculé  (line fine) et SSHA recalculé + filtre median (gras).

 Tracks77690536131886294849032576247
705 161 619 75 533 991 447 232 690

        ii.b  Signal SSHA recalculé + filtre median (deux figures de la première colonne ) et SSH (mdt+SSHA) calculé à partir de ce signal final

 Tracks77690536131886294849032576247
705 161 619 75 533 991 447 232 690

   ----------------------------------------------------------------------------------------------------------------------------------------------------------28/01/2015

    iii.  Vitesse geostrophique déduite de la SSH calculée précédemment :  Ug  = ~gradient(SSHA+mdt)

          iii.a .

 Tracks77690536131886294849032576247
705 161 619 75 533 991 447 232 690

 

------------------------------------------------------------------------------------------------------------------------------------------------------------13/03/2015

          iiib . Ug calculée à partir du lissage des signaux finaux de SSH (issues de iib. ): filtrage à ~100km en utilisant le filtre SPZ (ferret)

  1.           Interpolation sur une grille régulière (~3km) +  Application du filtre median à 5point
  2.           Filtrage spz (ferret) à ~100km de la SSH avant le calcul de la geostrophie

 

Traces : 776
filtre 50km

ssh&ugeo

ugeo_filt

filtre 100km

ssh&ugeo

ugeo_filt

 

------------------------------------------------------------------------------------------------------------------------------------------------------------17/03/2015

C.a. Hovmuller of SSHA (iib) of  some tracks for all available cycles (c1..20)

Tracks 776 518 905 361 490 991

 

C.b. Carte de vitesse geostrophique de certaines traces SRL (filtre 100km)

Document Actions

logo cnes logo IRD Logo université de Toulouse Logo université Paul Sabatier Logo CNRS
Logo bibliothèque OBS Logo Observatoire Midi Pyrénées