Multi-Agent Modelling of Earth's Dynamics: Towards a Virtual Laboratory of Plate Tectonics

Symposium GEOCEAN en hommage à Jean FrancheteauMACMA (Multi-Agent Convective MAntle) is a new tool developed to simulate plate tectonics and mantle convection in a 2-D cylindrical geometry (Combes et al., 2012)

Modelling the thermal evolution of slow-spreading ridge segments and their off-axis geophysical signature

International audienceSystematic studies conducted between 15°N and 40°N over ridge segments along the slow-spreading Mid-Atlantic Ridge (MAR) have shown that segment characteristics are related to the thermal state of the segments and gradually vary with their length. This paper presents further developments of a 3-D model, based on the presence of a hot zone located under the segment centres (Gac et al. 2003), to (1) quantify the thermal structures and the geophysical signatures of segments of various lengths, considered as representative of the various MAR segments; (2) test if a simple and single model of thermal evolution can account for the characteristics of all segments and (3) explain the past evolution of the segmentation, as is observed off-axis along the MAR. The modelled thermal structure and three simulated geophysical outputs [crustal structure , along-axis variations of the earthquake maximum depth and the mantle Bouguer gravity anomalies (MBA)] are found to be directly controlled by the shape (geometry and dimensions) of the hot zone. A consistent fit between model outputs and along-axis variations of the geophysical observables over the various segments is obtained by varying solely the length of the hot zone. This result shows that segments of different length may in fact constitute the different stages of a single evolution process: the axial geophysical characteristics of the segment would progressively evolve from those of shorter segments to those of longer ones, as the hot zone lengthens along-axis. A subsequent shortening of the segment would result from a simultaneous shortening of the hot zone, segment characteristics reverting back from those of longer segments to those of shorter ones. Three geophysical fields (topography, gravity and magnetic anomalies) are simulated as the results of the thermal evolution of aligned and offset segments the length of which evolves through time. These simulations succeed in fitting observations for the entire range of observed axis offsets between adjacent MAR segments. The segment evolution produces peculiar off-axis isostatic topography and gravity anomaly (MBA), the rhomb-shaped patterns. Our simulations, which model adjacent offset segments having evolved through several cycles of lengthening and shortening, yield a good fit to the isostatic topography and MBA patterns observed in the off-axis region. Finally, the distribution of magnetization depends on the magnetic properties of each type of rocks and on the petrological structure of the lithosphere, which, in turn, results from its thermal structure and evolution. Modelled magnetic anomalies are shown to be in good agreement with off-axis observations along the N21°40' segment (TAMMAR) of the MAR

MACMA : Mantle cooling mechanisms simulated by agents

International audienceMACMA is a new simulating tool based on multiagent systems to build a virtual laboratory in Earth Sciences. Here we study Earth's mantle cooling mechanisms by superposition of analytical and empirical laws accounting for the conservation of mass, energy and momentum, together with the description of plate boundary kinematics

Multiagent simulation of evolutive plate tectonics applied to the thermal evolution of the Earth

International audience[1] The feedback between plate tectonics and mantle convection controls the Earth's thermal evolution via the seafloor age distribution. We therefore designed the MACMA model to simulate time-dependent plate tectonics in a 2D cylindrical geometry with evolutive plate boundaries, based on multiagent systems that express thermal and mechanical interactions. We compute plate velocities using a local force balance and use explicit parameterizations to treat tectonic processes such as trench migration, subduction initiation, continental breakup and plate suturing. These implementations allow the model to update its geometry and thermal state at all times. Our approach has two goals: (1) to test how empirically- and analyticallydetermined rules for surface processes affect mantle and plate dynamics, and (2) to investigate how plate tectonics impact the thermal regime. Our predictions for driving forces, plate velocities and heat flux are in agreement with independent observations. Two time scales arise for the evolution of the heat flux: a linear long-term decrease and high-amplitude short-term fluctuations due to surface tectonics. We also obtain a plausible thermal history, with mantle temperature decreasing by less than 200 K over the last 3 Gyr. In addition, we show that on the long term, mantle viscosity is less thermally influential than tectonic processes such as continental breakup or subduction initiation, because Earth's cooling rate depends mainly on its ability to replace old insulating seafloor by young thin oceanic lithosphere. We infer that simple convective considerations alone cannot account for the nature of mantle heat loss and that tectonic processes dictate the thermal evolution of the Earth

Axial magnetic anomalies over slow-spreading ridge segments: insights from numerical 3-D thermal and physical modelling

International audienceThe axial magnetic anomaly amplitude along Mid-Atlantic Ridge segments is systematically twice as high at segment ends compared with segment centres. Various processes have been proposed to account for such observations, either directly or indirectly related to the thermal structure of the segments: (1) shallower Curie isotherm at segment centres, (2) higher Fe-Ti content at segment ends, (3) serpentinized peridotites at segment ends or (4) a combination of these processes. In this paper the contribution of each of these processes to the axial magnetic anomaly amplitude is quantitatively evaluated by achieving a 3-D numerical modelling of the magnetization distribution and a magnetic anomaly over a medium-sized, 50 km long segment. The magnetization distribution depends on the thermal structure and thermal evolution of the lithosphere. The thermal structure is calculated considering the presence of a permanent hot zone beneath the segment centre. The 'best-fitting' thermal structure is determined by adjusting the parameters (shape, size, depth, etc.) of this hot zone, to fit the modelled geophysical outputs (Mantle Bouguer anomaly, maximum earthquake depths and crustal thickness) to the observations. Both the thermoremanent magnetization, acquired during the thermal evolution, and the induced magnetization, which depends on the present thermal structure, are modelled. The resulting magnetic anomalies are then computed and compared with the observed ones. This modelling exercise suggests that, in the case of aligned and slightly offset segments, a combination of higher Fe-Ti content and the presence of serpentinized peridotites at segment ends will produce the observed higher axial magnetic anomaly amplitudes over the segment ends. In the case of greater offsets, the presence of serpentinized peridotites at segment ends is sufficient to account for the observations

Seismicity and accretion processes along the Mid-Atlantic Ridge south of the Azores using data from the MARCHE autonomous hydrophone array

The seismicity of the South Atlantic Ocean has been recorded by the MARCHE network of 4 autonomous underwater hydrophones (AUH) moored within the SOFAR channel on the flanks of the Mid-Atlantic Ridge (MAR). The instruments were deployed south of the Azores Plateau between 32° and 39°N from July 2005 to August 2008. The low attenuation properties of the SOFAR channel for earthquake T-wave propagation result in a detection threshold reduction from a magnitude completeness level (Mc) of ~4.3 for MAR events recorded by the land-based seismic networks to Mc=2.1 using this hydrophone array. A spatio-temporal analysis has been performed among the 5600 events recorded inside the MARCHE array. Most events are distributed along the ridge between lat. 39°N on the Azores Platform and the Rainbow (36°N) segment. In the hydrophone catalogue, acoustic magnitude (Source Level, SL) is used as a measure of earthquake size. The source level above which the data set is complete is SLc=205 dB. We look for seismic swarms using the cluster software of the SEISAN package. The criterion used are a minimum SL of 210 to detect a possible mainshock, and a radius of 30 km and a time window of 40 days after this mainshock (Cevatoglu, 2010, Goslin et al., 2012). 7 swarms with more than 15 events are identified using this approach between 32°et 39°N of latitude. The maximum number of earthquake in a swarm is 57 events. This result differs from the study of Simao et al. (2010) as we processed a further year of data and selected sequences with fewer events. Looking at the distribution of the SL as a function of time after the mainshock, we discuss the possible mechanism of these earthquakes : tectonic events with a "mainshock-aftershock" distribution fitting a modified Omori law or volcanic events showing more constant SL values. We also present the geophysical setting of these 7 swarms, using gravity, bathymetry, and available local geological data. This study illustrates the potential of hydrophone data to monitor segment-scale ridges processes in the vicinity of the Lucky Strike seafloor observatory (lat. 37°20'N), the Azores node of the EMSO (European Multidiciplinary Subsea Observatory) system