Le comportement du champ magnétique de la Terre, à l'échelle d'un intervalle de polarité, peut être représenté à l'aide de modèles statistiques de type "Giant Gaussian Process" (GGP) permettant de décrire la valeur du champ moyen et la variabilité du champ autour de cette valeur moyenne (ou paléovariation séculaire). Dans un premier temps, le champ magnétique produit par des simulations numériques de géodynamo nous a permis de tester la validité de cette approche. L'étude de la morphologie du champ moyen et de la paléovariation séculaire au sein de ces simulations a montré que les modèles "GGP" peuvent fournir des informations sur le fonctionnement de la géodynamo et sur l'influence des conditions à la limite noyau manteau. Nous avons ensuite développé une nouvelle méthode permettant de tester les modèles "GGP" avec les mesures de directions paléomagnétiques, en prenant en compte les erreurs de mesures. Dans un second temps, les nombreux profils d'anomalies magnétiques marines de surface disponibles dans les bases de données nous ont permis de mettre en évidence de nombreuses anomalies de courte longueur d'onde cohérentes sur de grandes distances (ou "tiny wiggles") et qui sont associées pour la plupart à des variations temporelles de l'intensité du champ géomagnétique. La distribution temporelle de ces tiny wiggles à l'intérieur des chrons et au cours de la période étudiée (83-41 Ma, chrons 33r à 19r) est homogène et suggère qu'ils reflètent simplement l'enregistrement filtré d'un régime de variation séculaire constant, contrastant avec l'évolution de la fréquence des inversions caractérisant la période étudiée.