Interpretation of Mars southern highlands high amplitude magnetic field with total gradient and fractal source modeling: New insights into the magnetic mystery of Mars

Using model studies, the total gradient (TG) of the Z-component magnetic field is shown to be a useful quantity for delineating sources of satellite-altitude magnetic anomalies; this field is used to constrain the location and lateral boundaries of sources of high amplitude magnetic anomalies of southern highlands of Mars. The TG field suggests two parallel linear and oppositely magnetized sources of 1000 and 1800 km length separated by 1000 km of region of intervening non-parallel sources. The simplest interpretation of the long, linear features is that they are zones of multitudinous crustal scale dikes formed in separate episodes of rifting, and not features associated with the mechanism of seafloor spreading. Forward modeling with uniformly magnetized sources suggests that magnetizations of the order of 10-50 A/m (40 km thickness) over ∼100 km width in the case of the southern source and of 12.5-27.5 A/m (40 km thickness) and ∼200 km width for the northern source are necessary to explain the Z-component amplitudes and features of the TG field. If the crustal magnetization on Mars were to be distributed fractally as on Earth, magnetizations matching the largest amplitude features on Mars may be spatially correlated from a 50-100 km distance range (β ∼ 3) to approaching nearly uniform magnetization (β ∼ 5) values. To keep magnetization intensity as small as possible, the higher end of β values are preferred, whereas, small amplitude anomaly features could be generated from sources with β ∼ 3. Many of the Mars anomaly features could be coalescence effects similar to the coalescence of anomaly features observed on http://icarus.cornell.edu/information/keywords.html.Earth.     

Observation of saturnian stream particles in the interplanetary space

In January 2004 the dust instrument on the Cassini spacecraft detected the first high-velocity grain expelled from Saturn - a so-called stream particle. Prior to Cassini’s arrival at Saturn in July 2004 the instrument registered 801 faint impacts, whose impact signals showed the characteristic features of a high-velocity impact by a tiny grain. The impact rates as well as the directionality of the stream particles clearly correlate with the sector structure of the interplanetary magnetic field (IMF). The Cosmic Dust Analyser (CDA) registered stream particles dominantly during periods when the IMF direction was tangential to the solar wind flow and in the prograde direction. This finding provides clear evidence for a continuous outflow of tiny dust grains with similar properties from the saturnian system. Within the compressed part of co-rotating interaction regions (CIRs) of the IMF, characterized by enhanced magnetic field strength and compressed solar wind plasma, CDA observed impact bursts of faster stream particles. We find that the bursts result from the stream particles being sped up inside the compressed CIR regions. Our analysis of the stream-particle dynamics inside rarefaction regions of the IMF implies that saturnian stream particles have sizes between 2 and 9 nm and exit the saturnian systems closely aligned with the planet’s ring plane with speeds in excess of 70 km s⁻¹.     

Highly reducing conditions during core formation on Mercury: Implications for internal structure and the origin of a magnetic field

The high average density and low surface FeO content of the planet Mercury are shown to be consistent with very low oxygen fugacity during core segregation, in the range 3-6 log units below the iron-wüstite buffer. These low oxygen fugacities, and associated high metal content, are characteristic of high-iron enstatite (EH) and Bencubbinite (CB) chondrites, raising the possibility that such materials may have been important building blocks for this planet. With this idea in mind we have explored the internal structure of a Mercury sized planet of EH or CB bulk composition. Phase equilibria in the silicate mantle have been modeled using the thermodynamic calculator p-MELTS, and these simulations suggest that orthopyroxene will be the dominant mantle phase for both EH and CB compositions, with crystalline SiO₂ being an important minor phase at all pressures. Simulations for both compositions predict a plagioclase-bearing “crust” at low pressure, significant clinopyroxene also being calculated for the CB bulk composition. Concerning the core, comparison with recent high pressure and high temperature experiments relevant to the formation of enstatite meteorites, suggest that the core of Mercury may contain several wt.% silicon, in addition to sulfur. In light of the pressure of the core-mantle boundary on Mercury (∼7 GPa) and the pressure at which the immiscibility gap in the system Fe-S-Si closes (∼15 GPa) we suggest that Mercury’s core may have a complex shell structure comprising: (i) an outer layer of Fe-S liquid, poor in Si; (ii) a middle layer of Fe-Si liquid, poor in S; and (iii) an inner core of solid metal. The distribution of heat-producing elements between mantle and core, and within a layered core have been quantified. Available data for Th and K suggest that these elements will not enter the core in significant amounts. On the other hand, for the case of U both recently published metal/silicate partitioning data, as well as observations of U distribution in enstatite chondrites, suggest that this element behaves as a chalcophile element at low oxygen fugacity. Using these new data we predict that U will be concentrated in the outer layer of the mercurian core. Heat from the decay of U could thus act to maintain this part of Mercury’s core molten, potentially contributing to the origin of Mercury’s magnetic field. This result contrasts with the Earth where the radioactive decay of U represents a negligible contribution to core heating.     

Low-temperature magnetic properties of iron-bearing sulfides and their contribution to magnetism of cometary bodies

In this study we present a review of low-temperature magnetic properties of alabandite (Fe, Mn)S, daubreelite FeCr₂S₄, pyrrhotite Fe_1−xS and troilite FeS updated with new experimental data. The results indicate that besides FeNi alloys mainly daubreelite with its Curie temperature T_C ∼ 150 K and strong induced and remanent magnetizations may be a significant magnetic mineral in cold environments and may complement that of FeNi or even dominate magnetic properties of sulfide rich bodies at temperatures below T_C.Comets are known to contain iron-bearing sulfides within dusty fraction and their surfaces are subject to temperature variations in the range of 100-200 K down to the depth of several meters while the cometary interior is thermally stable at several tens of Kelvin which is within the temperature range where alabandite, daubreelite or troilite are “magnetic”. Thus not only FeNi alloys, but also sulfides have to be considered in the interpretation of magnetic data from cometary objects such as will be delivered by Rosetta mission. Modeling indicates that magnetic interactions between cometary nucleus containing iron-bearing sulfides and interplanetary magnetic field would be difficult, but not impossible, to detect from orbit. Rosetta’s Philae lander present on the surface would provide more reliable signal.     

Lithospheric drift on early Mars: Evidence in the magnetic field

The crustal magnetic anomalies on Mars may represent hot spot tracks resulting from lithospheric drift on ancient Mars. As evidence, an analysis of lineation patterns derived from the ΔBr magnetic map is presented. The ΔBr map, largely free of external magnetic field effects, allows excellent detail of the magnetic anomaly pattern, particularly in areas of Mars where the field is relatively weak. Using cluster analysis, we show that the elongated anomalies in the martian magnetic field form concentric small circles (parallels of latitude) about two distinct north pole locations. If these pole locations represent ancient spin axes, then tidal force on the early lithosphere by former satellites in retrograde orbits may have pulled the lithosphere in an east-west direction over hot mantle plumes. With an active martian core dynamo, this may have resulted in the observed magnetic anomaly pattern of concentric small circles. As further evidence, we observe that, of the 15 martian giant impact basins that were possibly formed while the core dynamo was active, seven lie along the equators of our two proposed paleopoles. We also find that four other re-magnetized giant impact basins lie along a great circle about the mean magnetic paleopole of Mars. These 11 impact basins, likely the result of fallen retrograde satellite fragments, indicate that Mars once had moons large enough to cause tidal drag on the early martian lithosphere. The results of this study suggest that the magnetic signatures of this tidal interaction have been preserved to the present day.     

Alfven Wave Reflection model of field-aligned currents at Mercury

An Alfven Wave Reflection (AWR) model is proposed that provides closure for strong field-aligned currents (FACs) driven by the magnetopause reconnection in the magnetospheres of planets having no significant ionospheric and surface electrical conductance. The model is based on properties of the Alfven waves, generated at high altitudes and reflected from the low-conductivity surface of the planet. When magnetospheric convection is very slow, the incident and reflected Alfven waves propagate along approximately the same path. In this case, the net field-aligned currents will be small. However, as the convection speed increases, the reflected wave is displaced relatively to the incident wave so that the incident and reflected waves no longer compensate each other. In this case, the net field-aligned current may be large despite the lack of significant ionospheric and surface conductivity. Our estimate shows that for typical solar wind conditions at Mercury, the magnitude of Region 1-type FACs in Mercury’s magnetosphere may reach hundreds of kilo-Amperes. This AWR model of field-aligned currents may provide a solution to the long-standing problem of the closure of FACs in the Mercury’s magnetosphere.     

Mercury’s magnetospheric magnetic field after the first two MESSENGER flybys

The “paraboloid” model of Mercury’s magnetospheric magnetic field is used to determine the best-fit magnetospheric current system and internal dipole parameters from magnetic field measurements taken during the first and second MESSENGER flybys of Mercury on 14 January and 6 October 2008. Together with magnetic field measurements taken during the Mariner 10 flybys on 29 March 1974 and 16 March 1975, there exist three low-latitude traversals separated in longitude and one high-latitude encounter. From our model formulation and fitting procedure a Mercury dipole moment of 196 nT ·  (where R_M is Mercury’s radius) was determined. The dipole is offset from Mercury’s center by 405 km in the northward direction. The dipole inclination to Mercury’s rotation axis is relatively small, ∼4°, with an eastern longitude of 193° for the dipole northern pole. Our model is based on the a priori assumption that the dipole position and the moment orientation and strength do not change in time. The root mean square (rms) deviation between the Mariner 10 and MESSENGER magnetic field measurements and the predictions of our model for all four flybys is 10.7 nT. For each magnetic field component the rms residual is ∼6 nT or about 1.5% of the maximum measured magnetic field, ∼400 nT. This level of agreement is possible only because the magnetospheric current system parameters have been determined separately for each flyby. The magnetospheric stand-off distance, the distance from the planet’s center to the inner edge of the tail current sheet, the tail lobe magnetic flux, and the displacement of the tail current sheet relative to the Mercury solar-magnetospheric equatorial plane have been determined independently for each flyby. The magnetic flux in the tail lobes varied from 3.8 to 5.9 MWb; the subsolar magnetopause stand-off distance from 1.28 to 1.43 R_M; and the distance to the inner edge of the current sheet from 1.23 to 1.32 R_M. The differences in the current systems between the first and second MESSENGER flybys are attributed to the effects of strong magnetic reconnection driven by southward interplanetary magnetic field during the latter flyby.     

The analytic properties of the flapping current sheets in the earth magnetotail

The current sheet in Earth’s magnetotail often flaps, and the flapping waves could be induced propagating towards the dawn and dusk flanks, which could make the current sheet dynamic. To explore the dynamic characteristics of current sheet associated with the flapping motion holistically and provide reasonable physical interpretations, detailed direct calculation and analysis have been applied to one approximate analytic model of magnetic field in the flapping current sheet. The main results from the model demonstrate: (1) the magnetic fluctuation amplitude is attenuated from the center of current sheet to the lobe regions; The larger wave amplitude would induce the larger magnetic amplitude; (2) the curvature of magnetic field lines (MFLs), with maximum at the center of current sheet, is only dependent on the displacement Z along the south-north direction from the center of current sheet, regardless of the tilt of current sheet; (3) the half-thickness of neutral sheet, h, the minimum curvature radius of MFLs, R_cmin, and the tilt angle of current sheet, δ, satisfies h=R_cmin cos δ; (4) the gradient of magnetic strength forms a double-peak profile, and the peak value would be more intense if the local current sheet is more tilted; (5) current density j and its j_y, j_z components reach the extremum at the center of CS. j and j_z would be more intense if the local current sheet is more tilted, but it is not the case for j_y; and (6) the field-aligned component of current density mainly appears in the neutral sheet, and the sign of it would change alternatively as the flapping waves passing by. To check the validity of the model, one simulation on the virtual measurements has been made, and the results are in well consistence with actual observations of Cluster.     

Comparison of the Characteristics of Magnetic Clouds and Magnetic Cloud-Like Structures for the Events of 1995 – 2003

Using nine years of solar wind plasma and magnetic field data from the Wind mission, we investigated the characteristics of both magnetic clouds (MCs) and magnetic cloud-like structures (MCLs) during 1995 – 2003. A MCL structure is an event that is identified by an automatic scheme (Lepping, Wu, and Berdichevsky, Ann. Geophys. 23, 2687, 2005) with the same criteria as for a MC, but it is not usually identifiable as a flux rope by using the MC (Burlaga et al., J. Geophys. Res. 86, 6673, 1981) fitting model developed by Lepping, Jones, and Burlaga (Geophys. Res. Lett. 95(11), 957, 1990). The average occurrence rate is 9.5 for MCs and 13.6 for MCLs per year for the overall period of interest, and there were 82 MCs and 122 MCLs identified during this period. The characteristics of MCs and MCL structures are as follows: (1) The average duration, Δt, of MCs is 21.1 h, which is 40% longer than that for MCLs (Δt=15 h); (2) the average (minimum B _z found in MC/MCL measured in geocentric solar ecliptic coordinates) is −10.2 nT for MCs and −6 nT for MCLs; (3) the average Dst_min (minimum Dst caused by MCs/MCLs) is −82 nT for MCs and −37 nT for MCLs; (4) the average solar wind velocity is 453 km s⁻¹ for MCs and 413 km s⁻¹ for MCLs; (5) the average thermal speed is 24.6 km s⁻¹ for MCs and 27.7 km s⁻¹ for MCLs; (6) the average magnetic field intensity is 12.7 nT for MCs and 9.8 nT for MCLs; (7) the average solar wind density is 9.4 cm⁻³ for MCs and 6.3 cm⁻³ for MCLs; and (8) a MC is one of the most important interplanetary structures capable of causing severe geomagnetic storms. The longer duration, more intense magnetic field and higher solar wind speed of MCs, compared to those properties of the MCLs, are very likely the major reasons for MCs generally causing more severe geomagnetic storms than MCLs. But the fact that a MC is an important interplanetary structure with respect to geomagnetic storms is not new (e.g., Zhang and Burlaga, J. Geophys. Res. 93, 2511, 1988; Bothmer, ESA SP-535, 419, 2003).     

The evolution of rifting on the volcanic margin of the Pelotas Basin and the contextualization of the Paraná–Etendeka LIP in the separation of Gondwana in the South Atlantic

The Pelotas Basin is the classical example of a volcanic passive margin displaying large wedges of seaward-dipping reflectors (SDR). The SDR fill entirely its rifts throughout the basin, characterizing the abundant syn-rift magmatism (133–113 Ma). The Paraná–Etendeka Large Igneous Province (LIP), adjacent to west, constituted the pre-rift magmatism (134–132 Ma). The interpretation of ultra-deep seismic lines showed a very different geology from the adjacent Santos, Campos and Espírito Santo Basins, which constitute examples of magma-poor passive margins. Besides displaying rifts totally filled by volcanic rocks, diverse continental crustal domains were defined in the Pelotas Basin, such as an outer domain, probably constituted by highly stretched and permeated continental igneous crust, and a highly reflective lower crust probably reflecting underplating.The analysis of rifting in this portion of the South Atlantic is based on seismic interpretation and on the distribution of regional linear magnetic anomalies. The lateral accretion of SDR to the east towards the future site of the breakup and the temporal relationship between their rift and sag geometries allows the reconstitution of the evolution of rifting in the basin. Breakup propagated from south to north in three stages (130–127.5; 127.5–125; 125–113 Ma) physically separated by oceanic fracture zones (FZ). The width of the stretched, thinned and heavily intruded continental crust also showed a three-stage increase in the same direction and at the same FZ. Consequently, the Continental-Oceanic Boundary (COB) shows three marked shifts, from west to east, from south to north, resulting into rift to margin segmentation. Rifting also propagated from west to east, in the direction of the final breakup, in each of the three segments defined. The importance of the Paraná–Etendeka LIP upon the overall history of rupturing and breakup of Western Gondwanaland seems to have been restricted in time and in space only to the Pelotas Basin.