ISAS Methodology

ISAS (In Situ Analysis System) is an optimal interpolation tool developped to synthesize the global dataset of Argo profiles. Complementary data are also taken into account. ISAS is used in research mode at LPO and in operational mode by Coriolis.

ISAS is a data based re-analysis of temperature and salinity fields over the global ocean. It was initially designed to synthetize the temperature and salinity profiles collected by the ARGO program. It has been later extended to accommodate all type of vertical profile as well as time series.  The system aims at monitoring the time evolution of ocean properties for climatological studies and allowing easy computation of climate indices.
ISAS gridded fields are entirely based on in-situ measurements. The methodology and configuration have been conceived to preserve as much as possible the data information content and resolution. ISAS is developed and run in a research laboratory (LPO) in close collaboration with Coriolis, one of ARGO Global Data Assembly Center and unique data provider for the Mercator operational oceanography system.  At the moment the period covered starts in 2002 and only the upper 2000m are considered.

1 - Datasets

ISAS re-analyses are performed iteratively. Le starting point is a global dataset assembled by Coriolis.  For each new re-analysis, delayed mode datasets that have been validated by Pis are downloaded for a specific list of platforms. In case of duplicate, the delayed mode data replaces its equivalent in the global dataset.

1.1       Global datasets

Coriolis is both an ARGO Global Data Assembly Center (ARGO-GDAC) and a data delivery portal for operational oceanography (Mercator).  It produces global datasets for real time or re-analysis uses.

CORA datasets are prepared by Coriolis from the content of their database at a particular time with the aim of operational re-analysis.

Coriolis NRT data are global datasets processed each month from real time data, using an operational version of ISAS-tool. The data set for month m-1 is delivered on the 8th of month m.

1.2       Delayed Mode datasets

The delayed mode processing is either performed by the PIs (who usually deployed the instrument within a scientific project), or by a datacenter (Coriolis for exemple).  Data from a batch of similar platforms are considered at once over a long duration (6 month to the duration of the deployment). The measurements are controlled and may be adjusted if a drift of the sensor has been estimated. Because of this  type of processing the global delayed mode dataset is of inhomogeneous quality and evolves continuously. It is of major importance to keep track of the status and history of the data used.  The main delayed mode datasets used up to now are summarized in Table 1


Data provider

ARGO profilers

NAARC (Coriolis – ARGO-GDAC)

Marine animals

PI (F. Roquet)

Moorings TAO/Pirata/RAMA



Coriolis - GOSUD DM-datasets

Drifting Buoys T-S


Table 1: The different types of platforms providing delayed mode data taken into account in the analysis and the data provider.

1.3       Datasets considered for the ISAS_13 analysis

The analysis of the period 2002-2012 was obtained iteratively, starting from the CORA-2.2 dataset provided by Coriolis. This initial dataset went through additional post-analysis controls. It has been later complemented with the NRT monthly datasets. The characteristics of the datasets are detailed in Table 2.

Datatype selected: CO-01 was the first dataset produced. At that time, the delayed mode processing was not well defined, thus ‘raw data’ type was selected. In all subsequent dataset, ‘adjusted data’ are selected when available.

Pre-processing climatological controls: Since version 5, ISAS-tool includes a preprocessing test that compares each profile to the reference climatology. All profile with data points that fall out of 6 standard deviations are visually checked. For example, the automatic pre-processing of the 2012 data set has detected 19636 outlying profiles, the visual check lead to identify more than 300 platforms as showing some erroneous data. The full time series (for the year) of these platforms have been examined and the bad data, either the full profile, or isolated points, were flagged accordingly.

Post-analysis control: Diagnostic tools have been developed to compare the analyzed anomaly and the data residuals to the a priori statistics. They were applied on the CO-02 dataset after the analysis D2CA1S2, new suspicious data were then eliminated.

The analysis retained as ISAS_13 is based on the combination of the global dataset CO-06 and the delayed mode dataset DM-01. It corresponds to the D7CA2S0 experiment.


  • 2002-2008: CORA-2.2 – Coriolis re-analysis global dataset
  • Raw data selected


  • 2002-2008: CORA-2.2 – Coriolis re-analysis global dataset
  • 2009-2011: Coriolis Near Real Time global dataset
  • Adjusted data selected if available


  • 2002-2008:
    • CORA-2.2 – Coriolis re-analysis global dataset
    • Post-analysis QC performed
  • 2009-2011: Coriolis Near Real Time global dataset
    • Post-analysis QC performed
    • 2010 month 6,7,8,9: reprocessed
  • 2012: Coriolis Near Real Time global dataset
  • Adjusted data selected if available


  • DM-ARGO-Ovide 2002-2011 (44 floats)
  • DM-ARGO-GoodHope 2002-2011 (94 floats)
  • DM-Elephant Seals 2002-2011 (26 animals)
  • DM-TAO 2012 (test)

Table 2: Characteristics of the datasets used in the analysis.

2       Configuration:

2.1      Grid: 

  • Horizontal grid: 0.5° Mercator (frozen north of au 66.5°N)
  • Vertical grid: 152 levels from 0 to 2000 m. Levels are:  0, 3, and 5 m then every 5m down to 100m, every 10m down to 800m, every 20m down to 2000m.
    • The reference field is a mean annual cycle defined monthly.
    • The variance around this cycle is assumed identical for all months. It was computed as the mean of the monthly variance.

2.2     Mean field and variances

As the data set, the reference fields have been computed iteratively.

Climatology name

Characteristics of the analysis used


NODC climatology interpolated on the ISAS grid

All NODC data before 2004


Monthly means of the analysis D1CW5S1

Data: CO-1, reference WOA05, period averaged : 2002-2008


Monthly means of the analysis D2CA1S2

Data: CO-2, reference ARV09, period averaged:  2004-2010

Table 3: Reference fields

2.3      ISAS_13 parameters

Horizontal covariances are approximated as the sum of 2 gaussian functions defined by their scales L1 et L2. These scales can be anisotropic (Lix and Liy are defined independently).

Experiment name



CO-6 + DM-1

Reference fields


Period analyzed



W1 = 1,  L1time = 30 days

L1x = L1y = 300 km,  


W2 = 2, L2time = 30 days

1.5*dx <L2x  = 6 * Rossby Radius < 600 km          

1.5*dy <L2y  = 6 * Rossby Radius < 300 km


Wur = 6

Table 4: Parameters that define the optimal interpolation for the analysis ISAS_13.

3 References

Initial description of the ISAS method :

Gaillard, F.,  E. Autret, V.Thierry, P. Galaup, C. Coatanoan, and T. Loubrieu , 2009 : Quality control of large Argo data sets. JAOT, Vol. 26, No. 2.  337–351.

Validation of the global analysis (arv09) (ISAS version used : isas-V5) :

von Schuckmann, K., F. Gaillard, and P.-Y. Le Traon (2009), Global hydrographic variability patterns during 2003-2008, J. Geophys. Res., 114, C09007, doi:10.1029/2008JC005237.

Description of the ISAS method V6 :

Gaillard Fabienne (2012). ISAS-Tool Version 6 : Method and configuration.

Validation of the global analysis 2004-2010 (arv11)  (ISAS version used : isas-V6) :

 Kolodziejczyk, N. and F.  Gaillard (2012). Observation of spiciness interannual variability in the Pacific pycnocline. Journal of Geophysical Research ( JGR ) - Oceans, 117(C12018). Publisher's official version ,  Open Access version :

Kolodziejczyk Nicolas, Gaillard Fabienne (2013). Variability of the Heat and Salt Budget in the Subtropical Southeastern Pacific Mixed Layer between 2004 and 2010: Spice Injection Mechanism. Journal Of Physical Oceanography, 43(9), 1880-1898. Publisher's official version : , Open Access version :