1. Participants à l’étude
- Sélection des participants
- Méthode d’échantillonnage : Un échantillonnage en grappe de 30 membres présélectionnés et alternatifs de l’équipe d’expédition intérieure antarctique a été pris comme sujets d’étude.
- Sujets d’étude : inclure tous les candidats présélectionnés et les membres de réserve de l’équipe chinoise d’expédition intérieure antarctique. Excluez les personnes incapables de suivre le programme d’entraînement pendant la période d’entraînement (par exemple, celles qui n’ont pas complété le plan d’entraînement en raison d’une maladie, d’un inconfort physique ou d’autres raisons).
- Collecte de données
- Documentation démographique : Enregistrez des données sociodémographiques complètes, incluant l’âge, le niveau d’éducation et l’expérience antérieure lors d’expéditions polaires (Tableau 1).
- Veillez à ce que des mesures strictes d’anonymisation soient appliquées pour toutes les données psychométriques, notamment : l’attribution d’identifiants anonymisés uniques (par exemple, EXP##) pour remplacer les identifiants personnels, la suppression des identifiants directs (par exemple, noms, adresses) des ensembles de données, le stockage de données identifiables (par exemple, les formulaires de consentement) dans des fichiers chiffrés protégés par mot de passe, séparés des données de recherche.
| n=30 | N | % |
| Ethnicité | | |
| Chinois Han | 28 | 93.3 |
| Minorité ethnique (tibétaine et mandchou) | 2 | 6.67 |
| Plus haut niveau d’éducation | | |
| Lycée et en dessous | 1 | 3.3 |
| Diplôme universitaire ou licence | 11 | 36.7 |
| Master | 8 | 26.7 |
| Doctorat | 10 | 33.3 |
| Enfant unique | | |
| Oui | 16 | 53.3 |
| Non | 14 | 46.7 |
| État civil | | |
| Célibataire | 15 | 50 |
| Mariés | 15 | 50 |
| Nombre d’expéditions scientifiques participantes | | |
| Non participé | 24 | 80 |
| 1 fois | 2 | 6.7 |
| 2 fois ou plus | 4 | 13.3 |
| Âge : 32,72±7,73 | | |
Tableau 1 : Données démographiques. (A) Ce tableau présente la composition ethnique, le niveau d’éducation, le statut d’enfant unique, le statut marital et la fréquence de participation à la recherche pour l’échantillon, la dernière ligne indiquant les données d’âge (Moyenne ± SD). (B) N : Nombre d’échantillons ; % : Pourcentage par rapport à la taille totale de l’échantillon (n=30), les données montrent la moyenne ± SD.
2. Préparation d’instruments psychologiques et d’évaluation du sommeil
- Échelles psychométriques standardisées : Réalisez une évaluation par questionnaire administrée par un chercheur dans un format d’entretien standardisé afin d’assurer la compréhension et de minimiser les biais d’auto-déclaration. Guide verbalement les participants à travers chaque point, avec des clarifications fournies si nécessaire.
- Inventaire des Big Five - 44 (BFI-44)14
- Sélectionnez le questionnaire BFI-44 de 44 items pour évaluer cinq dimensions de la personnalité. Demandez aux participants de noter les énoncés à l’aide d’une échelle de Likert à 5 points (1 = Fortement en désaccord avec 5 = Fortement d’accord).
- Calculez les scores de domaine pour Extraversion, Agréabilité, Névrosisme, Ouverture et Conscience. Des scores plus élevés indiquent une expression plus forte des traits correspondants. L’échelle comprend 28 éléments à clé positive.
- Échelles de stress dépressif, anxiété - 21 (DASS-21)15
- Administrez le DASS-21 de 21 items avec une échelle Likert à 4 points (1= Ne s’appliquait pas du tout à 4=S’appliquait beaucoup à moi). Enregistrer des scores en sous-échelle pour Dépression, Anxiété et Stress. Des scores élevés reflètent une plus grande gravité des troubles émotionnels.
- Indice de qualité du sommeil de Pittsburgh (PSQI)10
- Distribuez des questionnaires PSQI de 18 items aux participants. Enseignez la fin de la lecture en 10 minutes en utilisant un score de 0 à 3 par domaine : qualité subjective du sommeil, latence du sommeil, troubles du sommeil, utilisation de médicaments hypnotiques, dysfonctionnement diurne. Calculez les scores globaux (échelle de 0 à 21). Des scores totaux plus élevés au PSQI indiquent une qualité globale de sommeil inférieure.
- Surveillance du sommeil
REMARQUE : L’appareil utilisé ici pour l’acte du poignet était ActiGraph wGT3X-BT16.
- Initialiser les appareils à l’aide du logiciel avec des horodatages standardisés.
- Lancez ActiLife v6.15.0. Cliquez sur Appareil > initialiser > ensemble. Réglez la fréquence d’échantillonnage : 60 Hz, Heure de début : YYYY-MM-DD HH :MM :SS (synchronisé avec l’horloge atomique), ID du sujet : EXP## (format anonymisé).
- Fixez un point de contrôle pour vérifier que le journal d’initialisation affiche « Initialisation réussie ».
- Calibrez le logiciel en utilisant les protocoles du fabricant après l’avoir placé sur le poignet non dominant. Insérez l’appareil dans la station d’accueil USB. Cliquez sur Calibrer puis suivez le protocole de test du vecteur gravitationnel à l’écran.
- Collectez des données d’activité 24h/24 à un taux d’échantillonnage de 60 Hz. Effectuez un point de contrôle visuel et confirmez que la LED d’état verte clignote tous les 15 secondes (fonctionnement normal).
- Effectuer l’acquisition de données. Après le déploiement, insérez l’appareil dans la station d’accueil. Sélectionnez ActiLife. Cliquez sur Télécharger > format RAW .gt3x > Processus utilisant l’algorithme Cole-Kripke. Utilisez l’algorithme Cole-Kripke pour la classification sommeil-veille. Valider les paramètres de sommeil dérivés de l’appareil via les journaux de sommeil desparticipants 17.
- Effectuer un échantillonnage avec un enregistrement continu 24/8 à 60 Hz (assurez l’espace entre la valeur/unité : 60 Hz)
3. Évaluation cognitive et acquisition de fNIRS
- Lettre N-retourtâche 8,9
- Avec le logiciel E-Prime 2.0, les lettres majuscules sont présentées aléatoirement au centre d’un écran LCD de 16 pouces. Les participants doivent porter des jugements : appuyer sur la touche F si la lettre actuelle correspond à celle présentée aux N positions précédentes ou appuyer sur la touche J si elle ne correspond pas (N=1, 1-retour ; N=2, 2-arrière ; N=3, 3-arrière).
- Enregistrer les métriques de précision (ACC) et de temps de réaction (RT) pour 50 essais par niveau de difficulté.
- Dans les paramètres de sortie E-Prime, configurez le DataFile.Name comme « EXP##Nback.edat2 ». Pour le temps d’apparition du stimulus, sélectionnez les variables : ACC (Correct = 1 / Incorrect = 0) et RT (en millisecondes). Cliquez sur Exporter > format Excel.
- Cliquez sur Analyser > modèle linéaire général > mesures répétéesFacteurs dans le sujet : niveau N(1,2,3). Fixez des points de contrôle aux données brutes : 50 essais/niveau = 150 essais au total.
- Pour la phase d’entraînement, dans E-Prime, cliquez sur BlockProc > Propriétés de pratique. Essais maximaux définis : 10, écran de retour : Correct ! (vert) / Incorrect (rouge), et Avancement : Appuyez sur la touche « Q » pour répéter la phase d’entraînement, appuyez sur la touche « P » pour passer à l’expérience formelle. Un point de contrôle est fixé, limitant la phase d’entraînement à un maximum de 3 répétitions (Figure 1).

Figure 1 : Paradigme de la lettre N-dos. Ce schéma illustre le paradigme de tâche à la lettre n-retour sur trois niveaux de difficulté (1-retour, 2-arrière et 3-retour). Veuillez cliquer ici pour voir une version agrandie de cette figurine.
- Acquisition de données fNIRS
REMARQUE : Le NirScan-500C, Huichuang, Danyang a été utilisé ici.
- Effectuez la configuration optode. Positionnez un couvre-tête fNIRS à 48 canaux selon le système international 10-20.
- Régler l’arrangement spatial avec la distance source-détecteur à 30 mm, la couverture : cortex préfrontal (PFC) et régions temporales bilatérales, et la matrice d’optodes : 15 optodes sources + 16 optodes détecteurs (Tableau 2, Tableau 3, Figure 2).
- Paramètres d’acquisition du signal
- Définir les paramètres d’acquisition du signal à l’aide du logiciel NirScan-500C. Sélectionnez la lumière proche infrarouge à double longueur d’onde : réglez les longueurs d’onde à 730 nm ± 5 nm et 850 nm ± 5 nm.
- Configurez le taux d’échantillonnage : réglez à 11 Hz en cliquant sur Paramètres > Taux d’acquisition. Activez la synchronisation tâche-événement : Activez les déclencheurs matériels en sélectionnant Synchroniser > déclencheur externe.
- Enregistrer les variations brutes d’intensité optique pour la conversion ultérieure en concentrations d’hémoglobine oxygénée/désoxygénée en utilisant la loi modifiée de Beer-Lambert. Vérifiez les paramètres dans le journal logiciel (par exemple, longueurs d’onde : 730/850 nm, débit : 11 Hz).

Figure 2 : Disposition des canaux, position et configuration de la sonde fNIRS. (A) Les cercles rouges représentent les sources lumineuses, les carrés bleus les détecteurs, et les zones numérotées représentent les paires sources-détecteurs (canaux) les plus proches pour mesurer les activités cérébrales. (B) Positions coenregistrées des optodes sur un atlas cérébral standard. La position anatomique de chaque canal de l’atlas cérébral est rapportée en détail dans le tableau 1. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
| Chaîne | Zone de Brodmann | Zone anatomique | Pourcentage de couverture |
| CH1 | 21 | Gyrus temporel moyen | 97.44% |
| CH2 | 22 | Gyrus temporal supérieur | 43.21% |
| CH3 | 38 | Aire temporopolaire | 51.84% |
| CH4 | 47 | Gyrus préfrontal inférieur | 100.00% |
| CH5 | 47 | Gyrus préfrontal inférieur | 56.35% |
| CH6 | 10 | Aire frontopolaire | 62.75% |
| CH7 | 10 | Aire frontopolaire | 51.96% |
| CH8 | 10 | Aire frontopolaire | 100.00% |
| CH9 | 11 | Aire orbitofrontale | 59.66% |
| CH10 | 11 | Aire orbitofrontale | 46.67% |
| CH11 | 10 | Aire frontopolaire | 98.66% |
| CH12 | 38 | Aire temporopolaire | 57.25% |
| CH13 | 21 | Gyrus temporel moyen | 98.63% |
| CH14 | 22 | Gyrus temporal supérieur | 43.46% |
| CH15 | 21 | Gyrus temporel moyen | 92.79% |
| CH16 | 21 | Gyrus temporel moyen | 75.85% |
| CH17 | 22 | Gyrus temporal supérieur | 62.67% |
| CH18 | 43 | Zone subcentrale | 36.79% |
| CH19 | 44 | les pars opercularis_ partie de la superficie de Broca | 57.05% |
| CH20 | 6 | Cortex moteur prémoteur et moteur supplémentaire | 73.08% |
| CH21 | 10 | Aire frontopolaire | 49.10% |
| CH22 | 46 | Cortex préfrontal dorsolatéral | 78.62% |
| CH23 | 10 | Aire frontopolaire | 93.20% |
| CH24 | 9,46 | Cortex préfrontal dorsolatéral | 77.95% |
| CH25 | 10 | Aire frontopolaire | 87.63% |
| CH26 | 10 | Aire frontopolaire | 100.00% |
| CH27 | 10 | Aire frontopolaire | 100.00% |
| CH28 | 10 | Aire frontopolaire | 99.65% |
| CH29 | 47 | Gyrus préfrontal inférieur | 58.95% |
| CH30 | 10 | Aire frontopolaire | 92.46% |
| CH31 | 45 | Pars triangularis Aire de Broca | 64.19% |
| CH32 | 46 | Cortex préfrontal dorsolatéral | 86.52% |
| CH33 | 21 | Gyrus temporel moyen | 56.73% |
| CH34 | 22 | Gyrus temporal supérieur | 53.70% |
| CH35 | 42 | Cortex de l’Association primaire et auditive | 64.69% |
| CH36 | 43 | Zone subcentrale | 46.23% |
| CH37 | 9,46 | Cortex préfrontal dorsolatéral | 60.54% |
| CH38 | 9 | Cortex préfrontal dorsolatéral | 52.09% |
| CH39 | 9,46 | Cortex préfrontal dorsolatéral | 92.77% |
| CH40 | 10 | Aire frontopolaire | 89.88% |
| CH41 | 9,46 | Cortex préfrontal dorsolatéral | 91.08% |
| CH42 | 9 | Cortex préfrontal dorsolatéral | 84.34% |
| CH43 | 10 | Aire frontopolaire | 92.46% |
| CH44 | 9 | Cortex préfrontal dorsolatéral | 83.67% |
| CH45 | 9,46 | Cortex préfrontal dorsolatéral | 90.37% |
| CH46 | 45 | Pars triangularis Aire de Broca | 39.27% |
| CH47 | 9,46 | Cortex préfrontal dorsolatéral | 98.16% |
| CH48 | 6 | Cortex moteur prémoteur et moteur supplémentaire | 80.97% |
Tableau 2 : Répartition des régions cérébrales du canal de mesure fNIRS après étalonnage par locateur 3D. (A) Caractéristiques de distribution spatiale de 48 canaux fNIRS calibrés par un système de positionnement 3D à travers les aires de Brodmann et les régions anatomiques cérébrales correspondantes. Le pourcentage de couverture indique la proportion effective d’acquisition du signal à partir des régions cérébrales ciblées. (B) Canal : numéro de canal optique de mesure fNIRS (CH1-CH48) ; Aire de Brodmann : nombre d’aire de Brodmann (par exemple, 22 = Cortex d’association auditive) ; Domaine anatomique : terminologie standard selon Terminologia Anatomica (1998) ; Pourcentage de couverture : Pourcentage de couverture efficace des paires optodes pour la région cérébrale ciblée.
| Zone anatomique | Chaîne |
| L-aPFC | 11,25,27,28,30,43 |
| R-aPFC | 6,7,8,21,23,26,40 |
| L-DLPFC | 32,44,45,47 |
| R-DLPFC | 22,24,37,38,39,41,42 |
Tableau 3 : Zones d’intérêt régional. (A) Afficher la distribution de la configuration des canaux pour des régions cérébrales fonctionnelles spécifiques dans le système de mesure fNIRS, y compris le cortex préfrontal antérieur (aPFC) et le cortex préfrontal dorsolatéral (DLPFC) dans les hémisphères gauche et droit. (B) Abréviations : L = hémisphère gauche ; R = hémisphère droit ; aPFC = cortex préfrontal antérieur ; DLPFC = Cortex préfrontal dorsolatéral.
4. Acquisition de données
- Mise en œuvre de l’évaluation psychologique
REMARQUE : Les membres de l’équipe ont quitté le niveau de la mer (Shanghai, 0 m d’altitude) et ont effectué un vol commercial (environ 8 h) vers Lhassa, située à 3 700 m d’altitude, où ils ont mené une formation de recherche scientifique de 72 heures (3 jours). Par la suite, l’équipe a transféré un transport terrestre (environ 2 h) à Yangbajain, à 4 300 m d’altitude, pour poursuivre une formation scientifique supplémentaire de 72 h (3 jours). La collecte des données a été effectuée immédiatement à l’arrivée de l’équipe à chacune des trois altitudes mentionnées précédemment (Figure 3).
- Mener une administration standardisée : lancer des sessions de présentation de groupe de pré-évaluation. Administrez des questionnaires avec une contrainte de 30 minutes.
- Atténuer les biais : lire les instructions standardisées issues de protocoles scriptés. Éliminez les formulations suggestives en suivant un texte pré-approuvé.
- Synchronisez les données multimodales : Permettre la synchronisation entre les tâches cognitives et l’acquisition de fNIRS. Verrouillez le temps des stimuli pour les déclencheurs fNIRS en sélectionnant E-Prime > Paramètres de déclenchement > Synchronisez avec fNIRS.
- Assurez-vous que toutes les sessions respectent des scripts et des contraintes de temps identiques.
- Surveillance des paramètres de sommeil
- Mettre en œuvre le protocole d’actigraphie continue : collecter des données d’actigraphie du poignet de 24 h/8h. Enregistrer simultanément un journal de sommeil.
- Assurez l’intégrité des données : effectuez des contrôles quotidiens de la batterie de 07h00 à 08h00 en cliquant sur Statut > Appareil > Batterie. Validez la continuité des données en cliquant sur ActiLife > Outils > Vérification de l’intégrité des données.

Figure 3 : Organigramme d’évaluation psychologique. (A) Ce diagramme illustre un protocole d’évaluation psychométrique en trois étapes à travers des expositions graduées à altitude : basse altitude (0-2 m s.m.) : mesures de base à Shanghai (jours 1-2), altitude moyenne (3650 m s.m.) : mesures de phase d’acclimatation à Lhassa (jours 3-5), haute altitude (4300 m s.m.) : mesures d’exposition élevée à Yangbajing (jours 6-8). (B) Modules d’évaluation identiques mis en place à chaque étape : Tâche cognitive : lettre N-retour (3 niveaux de difficulté) ; Échelles psychologiques : échelle dépression-anxiété et stress (DASS-21), Pittsburgh Sleep Quality Index (PSQI) ; Surveillance du sommeil : subjective (journal du sommeil), objective (actigraphie). Veuillez cliquer ici pour voir une version agrandie de cette figurine.
5. Analyse des données et statistiques
- Pour le prétraitement des données dans un tableur, appliquez des critères d’exclusion pour supprimer les ensembles de données incomplets avec >10 % de réponses manquantes en sélectionnant le filtre Données > > Valeurs manquantes. Filtrez les données cognitives pour écarter les valeurs aberrantes RT (<150 ms ou> moyenne ±2,5 SD) en utilisant Data > Sort & Filter > Advanced, et calculez le HitRT (en millisecondes) et le rejRT (en millisecondes). Gérez les données manquantes pour implémenter l’imputation multiple dans SPSS 29.0 en cliquant sur Transformer > Imputation multiple (5 itérations).
- Traitement du signal fNIRS
REMARQUE : Cela a été fait avec NirSpark v1.8.8, Huichuang, Danyang.
- Convertir les signaux bruts en densité optique (OD). Cliquez sur fichier > Importer Raw > sélectionnez fichier .nirs. Sélectionnez Convertir en densité optique et appliquez l’algorithme : OD = -log10(I/I 0). Vérifiez la plage OD 0,1-1,2 UA dans la vue > métriques de canal. Rejetez les canaux avec OD >1.5 (mauvais contact) en sélectionnant Outils > Rejeter les canaux.
- Calculez les changements hémodynamiques. Naviguez dans Analyse > paramètres hémodynamiques >. Régler la loi en Bière-Lambert modifiée, DPF à 66, et mode image à adult_head_template. Effectuer le calcul comme [HbO, HbR] = mBLL(OD, [730,850], [6,6]) et HbT = HbO + HbR. Confirmer que la valeur du DPF correspond à l’âge des participants (par exemple, 6,0 pour 25 à 45 ans).
- Réduisez le bruit des signaux en cliquant sur Prétraitement > Filtre > Paramètres. Sélectionnez Bandpass Butterworth du 4e ordre et réglez la plage à 0-0,2 Hz. Effectuez la correction de mouvement via la méthode de détection automatique (std_thr=6, amp_thr=0,5) et d’interpolation spline, puis réglez le seuil de rejet des artefacts à SD >±6. Vérifiez l’amélioration du SNR >10 dBin en cliquant sur Voir > Comparaison avant/post.
- Quantifiez l’activation comme suit. Définissez la référence comme la moyenne HbO pendant toute la période post-stimulus de la tâche. Calculez la réponse à la tâche en HbO moyen pendant 2 à 20 secondes après le stimulus. Générez des moyennes de ROI en sélectionnant les masques d’analyse > ROI > les modèles anatomiques.
- Calcul de l’efficacité du sommeil
- Importez et prétraitez les données en cliquant sur Fichier > Importer > Sélectionnez fichier .gt3x. Réglez l’Époque à 60 secondes et validez le temps d’usure en sélectionnant Outils > Validation du Temps d’Usure > algorithme Choi. Pour effectuer l’actigraphie et la synchronisation du journal du sommeil, cliquez sur Outils > Sommeil > Détection des périodes de sommeil. Ensuite, calibrez les données d’acte par rapport au journal du sommeil. Assurez-vous que les horodatages du journal de sommeil soient dans ± 5 minutes des données d’acte.
- Extraire les paramètres comme suit. Naviguez vers Rapports > Rapport de sommeil. Paramètres de sélection : Latence, TIB, TST, WASO. Calculer l’efficacité du sommeil = (TST / TIB) x 100 en sélectionnant Options > Formules. Exportez les résultats en cliquant sur Exporter > Excel .xlsx > Nom du fichier : EXP##_SleepEfficiency.xlsx. Confirmez que l’actogramme affiche ≥ 3 cycles de sommeil/éveil clairs.
- Tests statistiques
- Réalisez des statistiques descriptives. Déclarez les données normalement distribuées comme moyenne ± SD et les données non normales comme médiane [IQR].
- Exécutez des tests paramétriques. Effectuez des tests t à un échantillon en cliquant sur Analyser > Comparer les moyennes > test T à un seul échantillon. Effectuez RM-ANOVA en cliquant sur Analyser > Modèle linéaire général > Mesures répétées. Appliquez la correction Greenhouse-Geisser pour les violations de sphéricité.
- Appliquer des alternatives non paramétriques si nécessaire. Utilisez des tests non paramétriques sur des échantillons apparentés pour les tests de rang signé Wilcoxon. Le contrôle covarie en sélectionnant Analyser > Modèle linéaire général > ANCOVA. Fixez le seuil de signification en appliquant α=0,05 (à deux queues) avec correction de Bonferroni pour les comparaisons multiples.