La boîte à outils Gibbs-SeaWater (GSW) est une bibliothèque dédiée au calcul des propriétés physico-chiimques de l'eau de mer, basée sur la norme TEOS-10 (Thermodynamic Equation of Seawater - 2010). Son implémentation dans Matlab permet aux océanographes et aux ingénieurs de traiter et d'analyser de manière standardisée des données océanographiques complexes. Le code source fournit un ensemble de fonctions pour calculer la densité, la salinité, la vitesse du son, l'atténuation lumineuse et modéliser les interactions océan-atmosphère.
Présentation de la boîte à outils GSW
La boîte à outils GSW offre un cadre unifié pour le traitement des données marines, basé sur un état de référence thermodynamique. Elle est utilisée dans des domaines variés comme la climatologie, l'ingénierie sous-marine, la pêche et la biologie marine pour transformer des mesures brutes en informations exploitables sur l'état de l'océan.
Fonctionnalités principales
Les modules incluent le calcul des propriétés thermodynamiques (température, salinité, pression), la détermination de la vitesse du son et du modèle d'atténuation de la lumière, ainsi que l'estimation de la solubilité des gaz et du pH. L'accès aux fonctions se fait via des appels simples, et la documentation intégrée facilite leur utilisation.
Calcul des propriétés physiques de l'eau de mer
Relation entre température et densité
La température influence directement la densité de l'eau de mer. Une augmentation de température provoque une expansion de l'eau, réduisant ainsi sa densité. Ce mécanisme est fondamental pour la circulation thermohaline. Le calcul de la densité à partir de la température, de la salinité et de la pression peut être effectué comme suit :
% Paramètres d'entrée
T = 15; % Température en degrés Celsius
SP = 35; % Salinité pratique (PSS-78)
p = 0; % Pression relative en dbar (surface)
% Calcul de la densité in situ (rho)
rho_in_situ = gsw_rho_t_exact(SP, T, p);
fprintf('La densité calculée est de %.4f kg/m^3.\n', rho_in_situ);
Impact de la salinité sur la conductivité
La salinité de l'eau de mer est étroitement liée à sa conductivité électrique. Des instruments de mesure (CTD) utilisent cette relation pour déduire la salinité. La conversion s'effectue via des fonctions dédiées :
% Valeur de conductivité mesurée
C = 4.2; % Conductivité en S/m
% Déduction de la salinité à partir de la conductivité
SP_calc = gsw_SP_from_C(C, T, p);
fprintf('La salinité estimée est de %.4f g/kg.\n', SP_calc);
Calcul des propriétés chimiques
Calcul du pH
Le pH de l'eau de mer, généralement compris entre 7.5 et 8.4, est un indicateur clé de l'acidification océanique. Son calcul repose sur le logarithme négatif de l'activité des ions hydrogène et nécessite la connaissance du système carbonaté. L'équation de base est pH = -log[H⁺], mais des itérations avec des constantes d'équilibre sont souvent requises pour une précision océanographique.
Solubilité des gaz
La teneur en gaz dissous, en particulier l'oxygène et le dioxyde de carbone, est critique pour les écosystèmes marins. La loi de Henry gouverne leur solubilité, qui dépend de la température, de la salinité et de la pression partielle du gaz. Un modèle simplifié pour l'oxygène dissous peut s'écrire :
% Données d'entrée pour le calcul
kH_O2 = 1.265e-3; % Constante de Henry pour O2, en mol/(L·atm)
pO2 = 0.21; % Pression partielle d'O2 dans l'atmosphère, en atm
% Concentration molaire en oxygène dissous
concentration_O2 = kH_O2 * pO2;
fprintf('Concentration en O2 dissous : %.5f mol/L.\n', concentration_O2);
Des facteurs de correction pour la température et la salinité sont appliqués dans les modèles complets pour obtenir la saturation réelle.
Vitesse du son et atténuation lumineuse
Calcul de la vitesse du son
La vitesse du son dans l'eau de mer est fonction de la température, de la salinité et de la pression. Sa connaissance est essentielle pour les relevés acoustiques (sonars). L'influence combinée de ces paramètres peut être explorée par une analyse vectorielle :
% Définition des plages de paramètres
vect_T = 0:5:35; % Températures de 0 à 35 °C
vect_S = 30:5:40; % Salinités de 30 à 40 g/kg
% Création de grilles de calcul
[GRILLE_T, GRILLE_S] = ndgrid(vect_T, vect_S);
% Calcul vectoriel de la vitesse du son à pression atmosphérique
vitesse_son = gsw_sound_speed(GRILLE_T, GRILLE_S, 0);
% Affichage des résultats (matrice 8x3)
disp(vitesse_son);
Modélisation de l'atténuation de la lumière
La lumière s'atténue exponentiellement avec la profondeur selon la loi de Beer-Lambert : I(z) = I₀ · exp(-Kd · z), où Kd est le coefficient d'atténuation diffuse. Ce modèle est fondamental pour la biogénie et la vision sous-marine.
% Paramètres du modèle
I0 = 1; % Intensité de référence en surface
Kd = 0.05; % Coefficient d'atténuation (m⁻¹)
z = 0:0.1:10; % Profondeur de 0 à 10 m
% Calcul de l'intensité résiduelle
intensite = I0 * exp(-Kd .* z);
% Tracé du profil
figure;
plot(z, intensite, 'LineWidth', 1.5);
xlabel('Profondeur (m)');
ylabel('Intensité relative');
title('Profil d''atténuation lumineuse dans l''eau');
grid on;
Modélisation des interactions océan-atmosphère
Échange de chaleur à l'interface
Les flux de chaleur à la surface de la mer (rayonnement solaire, émission infrarouge, chaleur sensible et latente) régulent le climat global. Le calcul de ces flux intégrant des variables atmosphériques (température, humidité, vent) permet de simuler les rétroactions climatiques.
Processus d'évaporation et de précipitation
Le taux d'évaporation, influencé par la température de surface, l'humidité relative et le vent, affecte la salinité et le cycle hydrologique. Les précipitations, en apportant de l'eau douce, diluent la salinité de surface. Un modèle d'évaporation simplifié pourrait être implémenté pour des études de sensibilité.