Global mapping of lunar crustal magnetic fields by Lunar Prospector

The Lunar Prospector Electron Reflectometer has obtained the first global map of lunar crustal magnetic fields, revealing that the effects of basin-forming impacts dominate the large-scale distribution of remanent magnetic fields on the Moon. The weakest surface magnetic fields (<0.2 nT) are found within two of the largest and most recent impact basins, Orientale and Imbrium. Conversely, the largest concentrations of strong surface fields (>40 nT) are diametrically opposite to these same basins. This pattern is present though less pronounced for several other post-Nectarian impact basins larger than 500 km in diameter. The reduced strength and clarity of the pattern for older basins may be attributed to: (1) demagnetization from many smaller impacts, which erases antipodal magnetic signatures over time, (2) superposition effects from other large impacts, and (3) variation in the strength of the ambient magnetizing field. The absence of fringing fields stronger than 1 nT around the perimeter of the Imbrium basin or associated with craters within the basin implies that any uniform magnetization of the impact melt must be weaker than ∼¹⁰⁻⁶ G cm³ g⁻¹. This limits the strength of any steady ambient magnetic field to no more than ∼0.1 Oe at the lunar surface while the basin cooled for tens of millions of years following the Imbrium impact 3.8 billion years ago.     

Antipodal effects of lunar basin-forming impacts: Initial 3D simulations and comparisons with observations

3D simulations of basin-scale lunar impacts are carried out to investigate: (a) the origins of strong crustal magnetic fields and unusual terrain observed to occur in regions antipodal to young large basins; and (b) the origin of enhanced magnetic and geochemical anomalies along the northwest periphery of the South Pole-Aitken (SPA) basin. The simulations demonstrate that a basin-forming impact produces a massive, hot, partially ionized cloud of vapor and melt that expands thermally around the Moon, converging near the basin antipode approximately 1 h after the impact for typical impact parameters. In agreement with previous work, analytic calculations of the interaction of this vapor-melt cloud with an initial ambient magnetic field predict a substantial temporary increase in field intensity in the antipodal region. The time of maximum field amplification coincides with a period when impacting ejecta also converge near the antipode. The latter produce antipodal shock stresses within the range of 5-25 GPa where stable shock remanent magnetization (SRM) of lunar soils has been found experimentally to occur. Calculated antipodal ejecta thicknesses are only marginally sufficient to explain the amplitudes of observed magnetic anomalies if mean magnetization intensities are comparable to those produced experimentally. This suggests that pre-existing ejecta materials, which would also contain abundant metallic iron remanence carriers, may be important anomaly sources, a possibility that is consistent with enhanced magnetic anomalies observed peripheral to SPA. The latter anomalies may be produced by amplified secondary ejecta impact shock waves in the thick SPA ejecta mantle occurring near the antipodes of the Imbrium and Serenitatis impacts. Together with converging seismic compressional waves, these antipodal impact shocks may have produced especially deep fracture zones along the northwest edge of SPA near the Imbrium antipode, allowing the ascent of magma with enhanced KREEP concentrations.     

Early martian dynamo generation due to giant impacts

During late-stage planet formation, giant impacts produce localized mantle melt regions within which impactor iron droplets settle to the bottom near a permeability horizon. After accumulation, iron heated by the impact migrates downward to the core through colder, mostly solid mantle. The degree of thermal equilibration and partitioning of viscous heating between impactor iron and silicates depends on the mechanism of iron transport to the core. Simple estimates suggest that, following a giant impact, the temperature difference between iron delivered to the core and the mantle outside the impact heated region can be ∼10³ K. Hot impactor iron mergers with the core where it may be efficiently mixed or remain stratified due to thermal buoyancy. In either case, collisional energy carried to the core by impactor iron helps establish conditions favorable for early core cooling and dynamo generation. In this study, we consider the end-member scenario in which impactor iron forms a layer at the top of the core. Energy transfer from the impactor iron layer to the mantle is sufficient to power a dynamo for up to ∼30 Myr even in the limit of a very viscous mantle and heat flux limited by conduction. Using two-dimensional finite element calculations of mantle convection, we show that large-scale mantle flow driven by the buoyancy of the impact thermal anomaly focuses plumes in the impact region and increases both dynamo strength and duration. Melting within the mantle thermal boundary layer likely leads to formation of a single superplume in the location of the impact anomaly driven upwelling. We suggest that formation of magnetized southern highland crust may be related to spreading and differentiation of an impact melt region during the impact-induced dynamo episode.     

Size and compositional constraints of Ganymede's metallic core for driving an active dynamo

Ganymede has an intrinsic magnetic field which is generally considered to originate from a self-excited dynamo in the metallic core. Driving of the dynamo depends critically on the satellite's thermal state and internal structure. However, the inferred structure based on gravity data alone has a large uncertainty, and this makes the possibility of dynamo activity unclear; variations in core size and composition significantly change the heat capacity and alter the cooling history of the core. The main objectives of this study is to explore the structural conditions for a currently active dynamo in Ganymede using numerical simulations of the thermal history, and to evaluate under which conditions Ganymede can maintain the dynamo activity at present. We have investigated the satellite's thermal history using various core sizes and compositions satisfying the mean density and moment of inertia of Ganymede, and evaluate the temperature and heat flux at the core-mantle boundary (CMB). Based on the following two conditions, we evaluate the possibility of dynamo activity, thereby reducing the uncertainty of the previously inferred interior structure. The first condition is that the temperature at the CMB must exceed the melting point of a metallic core, and the second is that the heat flux through the CMB must exceed the adiabatic temperature gradient. The mantle temperature starts to increase because of the decay of long-lived radiogenic elements in the rocky mantle. After a few Gyr, radiogenic elements are exhausted and temperature starts to decrease. As the rocky mantle cools, the heat flux at the CMB steadily increases. If the temperature and heat flux at the CMB satisfy these conditions simultaneously, we consider the case as capable of driving a dynamo. Finally, we identify the Dynamo Regime, which is the specific range of internal structures capable of driving the dynamo, based on the results of simulations with various structures. If Ganymede's self-sustained magnetic field were maintained by thermal convection, the satellite's metallic core would be relatively large and, in comparison to other terrestrial-type planetary cores, strongly enriched in sulfur. The dynamo activity and the generation of the magnetic field of Ganymede should start from a much later stage, possibly close to the present.     

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.     

The phase diagram of water and the magnetic fields of Uranus and Neptune

The interior of giant planets can give valuable information on formation and evolution processes of planetary systems. However, the interior and evolution of Uranus and Neptune is still largely unknown. In this paper, we compare water-rich three-layer structure models of these planets with predictions of shell structures derived from magnetic field models. Uranus and Neptune have unusual non-dipolar magnetic fields contrary to that of the Earth. Extensive three-dimensional simulations of Stanley and Bloxham (Stanley, S., Bloxham, J. [2004]. Nature 428, 151-153) have indicated that such a magnetic field is generated in a rather thin shell of at most 0.3 planetary radii located below the H/He rich outer envelope and a conducting core that is fluid but stably stratified. Interior models rely on equation of state data for the planetary materials which have usually considerable uncertainties in the high-pressure domain. We present interior models for Uranus and Neptune that are based on ab initio equation of state data for hydrogen, helium, and water as the representative of all heavier elements or ices. Based on a detailed high-pressure phase diagram of water we can specify the region where superionic water should occur in the inner envelope. This superionic region correlates well with the location of the stably-stratified region as found in the dynamo models. Hence we suggest a significant impact of the phase diagram of water on the generation of the magnetic fields in Uranus and Neptune.     

Predicted and observed magnetic signatures of martian (de)magnetized impact craters

The current morphology of the martian lithospheric magnetic field results from magnetization and demagnetization processes, both of which shaped the planet. The largest martian impact craters, Hellas, Argyre, Isidis and Utopia, are not associated with intense magnetic fields at spacecraft altitude. This is usually interpreted as locally non- or de-magnetized areas, as large impactors may have reset the magnetization of the pre-impact material. We study the effects of impacts on the magnetic field. First, a careful analysis is performed to compute the impact demagnetization effects. We assume that the pre-impact lithosphere acquired its magnetization while cooling in the presence of a global, centered and mainly dipolar magnetic field, and that the subsequent demagnetization is restricted to the excavation area created by large craters, between 50- and 500-km diameter. Depth-to-diameter ratio of the transient craters is set to 0.1, consistent with observed telluric bodies. Associated magnetic field is computed between 100- and 500-km altitude. For a single-impact event, the maximum magnetic field anomaly associated with a crater located over the magnetic pole is maximum above the crater. A 200-km diameter crater presents a close-to-1-nT magnetic field anomaly at 400-km altitude, while a 100-km diameter crater has a similar signature at 200-km altitude. Second, we statistically study the 400-km altitude Mars Global Surveyor magnetic measurements modelled locally over the visible impact craters. This approach offers a local estimate of the confidence to which the magnetic field can be computed from real measurements. We conclude that currently craters down to a diameter of 200 km can be characterized. There is a slight anti-correlation of −0.23 between magnetic field intensity and impact crater diameters, although we show that this result may be fortuitous. A complete low-altitude magnetic field mapping is needed. New data will allow predicted weak anomalies above craters to be better characterized, and will bring new constraints on the timing of the martian dynamo and on Mars’ evolution.     

The apparent lack of lunar-like swirls on Mercury: Implications for the formation of lunar swirls and for the agent of space weathering

Images returned by the MESSENGER spacecraft from the Mercury flybys have been examined to search for anomalous high-albedo markings similar to lunar swirls. Several features suggested to be swirls on the basis of Mariner 10 imaging (in the craters Handel and Lermontov) are seen in higher-resolution MESSENGER images to lack the characteristic morphology of lunar swirls. Although antipodes of large impact basins on the Moon are correlated with swirls, the antipodes of the large impact basins on Mercury appear to lack unusual albedo markings. The antipodes of Mercury’s Rembrandt, Beethoven, and Tolstoj basins do not have surface textures similar to the “hilly and lineated” terrain found at the Caloris antipode, possibly because these three impacts were too small to produce obvious surface disturbances at their antipodes. Mercury does have a class of unusual high-reflectance features, the bright crater-floor deposits (BCFDs). However, the BCFDs are spectral outliers, not simply optically immature material, which implies the presence of material with an unusual composition or physical state. The BCFDs are thus not analogs to the lunar swirls. We suggest that the lack of lunar-type swirls on Mercury supports models for the formation of lunar swirls that invoke interaction between the solar wind and crustal magnetic anomalies (i.e., the solar-wind standoff model and the electrostatic dust-transport model) rather than those models of swirl formation that relate to cometary impact phenomena. If the solar-wind standoff hypothesis for lunar swirls is correct, it implies that the primary agent responsible for the optical effects of space weathering on the Moon is solar-wind ion bombardment rather than micrometeoroid impact.     

A global model of the internal magnetic field of the Moon based on Lunar Prospector magnetometer observations

A preliminary model of the internal magnetic field of the Moon is developed using a novel, correlative technique on the low-altitude Lunar Prospector magnetic field observations. Subsequent to the removal of a simple model of the external field, an internal dipole model is developed for each pole-to-pole half-orbit. This internal dipole model exploits Lunar Prospector's orbit geometry and incorporates radial and theta vector component data from immediately adjacent passes into the model. These adjacent passes are closely separated in space and time and are thus characteristic of a particular lunar regime (wake, solar wind, magnetotail, magnetosheath) or regimes. Each dipole model thus represents the correlative parts of three adjacent passes, and provides an analytic means of continuing the data to a constant surface of 30 km above the mean lunar radius. The altitude-normalized radial field from the wake and tail regimes is used to build a model in which 99.2% of the 360 by 360 bins covering the lunar surface are filled. This global model of the radial magnetic field is used to construct a degree 178 spherical harmonic model of the field via the Driscoll and Healy sampling theorem. Terms below about degree 150 are robust, and polar regions are considered to be the least reliable. The model resolves additional detail in the low magnetic field regions of the Imbrium and Orientale basins, and also in the four anomaly clusters antipodal to the large lunar basins. The model will be of use in understanding the sources of the internal field, and as a first step in modeling the interaction of the internal field with the solar wind.     

Precession of the lunar core

Goldreich (Goldreich, P. [1967]. J. Geophys. Res. 72, 3135) showed that a lunar core of low viscosity would not precess with the mantle. We show that this is also the case for much of lunar history. But when the Moon was close to the Earth, the Moon’s core was forced to follow closely the precessing mantle, in that the rotation axis of the core remained nearly aligned with the symmetry axis of the mantle. The transition from locked to unlocked core precession occurred between 26.0 and 29.0 Earth radii, thus it is likely that the lunar core did not follow the mantle during the Cassini transition. Dwyer and Stevenson (Dwyer, C.A., Stevenson, D.J. [2005]. An Early Nutation-Driven Lunar Dynamo. AGU Fall Meeting Abstracts GP42A-06) suggested that the lunar dynamo needs mechanical stirring to power it. The stirring is caused by the lack of locked precession of the lunar core. So, we do not expect a lunar dynamo powered by mechanical stirring when the Moon was closer to the Earth than 26.0-29.0 Earth radii. A lunar dynamo powered by mechanical stirring might have been strongest near the Cassini transition.     

A 3D hybrid simulation study of the electromagnetic field distributions in the lunar wake

A three-dimensional hybrid code is used to study the electromagnetic disturbances in the solar wind that arise due to the absorption effect of the Moon. Due to the nearly insulating nature of the Moon, interplanetary magnetic fields (IMFs) can move through the interior without hindrance. However, the near-vacuum created in the wake region due to the lunar absorption effect will lead to enhancement of the strength of the magnetic field by a factor of about 1.4 in the middle of the lunar wake and lead to depletions at two sides. The situations arising from different orientations of the interplanetary magnetic fields relative to the radial direction are compared. Asymmetries of the inward diffusions both along and perpendicular to the field lines are also observed. The electric field formed from the plasma convection could reach a magnitude of 0.2–0.8 mV/m at the border of the wake. The role of the electric field on the inward accelerations is important to the geometry of the lunar wake.     

Numerical dynamo models of Uranus' and Neptune's magnetic fields

The non-axisymmetric, non-dipolar magnetic fields of Uranus and Neptune are markedly different from the axially-dipolar dominated fields of the other planets in our Solar System with active dynamos. Stanley and Bloxham [Stanley, S., Bloxham, J., 2004. Nature 428, 151-153] used numerical modeling to demonstrate that Uranus' and Neptune's unusual fields could be the result of a different convective region geometry in these planets. They found that a numerical dynamo operating in a thin shell surrounding a stably-stratified fluid interior produces magnetic field morphologies similar to those of Uranus and Neptune. This geometry for the convective region was initially proposed by Hubbard et al. [Hubbard, W.B., Podolak, M., Stevenson, D.J., 1995. In: Cruickshank, D. (Ed.), Neptune and Triton. Univ. of Arizona Press, Tucson, pp. 109-138] to explain both the magnetic field morphology as well as the low intrinsic heat flows from these planets. Here we examine the influence of varying the stable layer radius in numerical models and compare the results to thin shell models surrounding solid inner cores. We find that a limited range of stable-layer shell thicknesses exist in which Uranus/Neptune-like field morphologies result. This allows us to put constraints on the size of the convective layers in Uranus and Neptune.     

Magnetic anomalies near Apollinaris Patera and the Medusae Fossae Formation in Lucus Planum, Mars

The nature of strong martian crustal field sources is investigated by mapping and modeling of Mars Global Surveyor magnetometer data near Apollinaris Patera, a previously proposed volcanic source, supplemented by large-scale correlative studies. Regional mapping yields evidence for positive correlations of orbital anomalies with both Apollinaris Patera and Lucus Planum, a nearby probable extrusive pyroclastic flow deposit that is mapped as part of the Medusae Fossae Formation. Iterative forward modeling of the Apollinaris Patera magnetic anomaly assuming a source model consisting of one or more uniformly magnetized near-surface disks indicates that the source is centered approximately on the construct with a scale size several times larger and comparable to that of the Apollinaris Patera free-air gravity anomaly. A significantly lower rms deviation is obtained using a two-disk model that favors a concentration of magnetization near the construct itself. Estimates for the dipole moment per unit area of the Lucus Planum source together with maximum thicknesses of ∼3 km based on topographic and radar sounding data lead to an estimated minimum magnetization intensity of ∼50 A/m within the pyroclastic deposits. Intensities of this magnitude are similar to those obtained experimentally for Fe-rich Mars analog basalts that cooled in an oxidizing (high fO₂) environment in the presence of a strong (?10 μT) surface field. Further evidence for the need for an oxidizing environment is provided by a broad spatial correlation of the locations of phyllosilicate exposures identified to date using Mars Express OMEGA data with areas containing strong crustal magnetic fields and valley networks in the Noachian-aged southern highlands. This indicates that the presence of liquid water, which is a major crustal oxidant, was an important factor in the formation of strong magnetic sources. The evidence discussed here for magnetic sources associated with relatively young volcanic units suggests that a martian dynamo existed during the late Noachian/early Hesperian, after the last major basin-forming impacts and the formation of the northern lowlands.     

Modeling of major martian magnetic anomalies: Further evidence for polar reorientations during the Noachian

Maps of the vector components of the Mars crustal magnetic field are constructed at the mapping altitude (360 to 410 km) using a selected set of data obtained with the Mars Global Surveyor magnetometer during 2780 orbits of the planet in 1999. Forward modeling calculations are then applied to six relatively strong and isolated, dominantly dipolar, magnetic anomalies for the primary purpose of estimating bulk directions of magnetization. Assuming that the magnetizing field was a (dipolar) core dynamo field centered in the planet, paleomagnetic pole positions are calculated for the six primary source bodies together with that for a seventh anomaly analyzed earlier. In agreement with several previous studies, it is found that six of the seven pole positions are clustered in what is now the northern lowlands in a region centered northwest of Olympus Mons (mean pole position: 34°±10° N, 202°±58° E). Assuming that the dynamo dipole moment vector was approximately parallel to the rotation axis, the modeling results therefore suggest a major reorientation of Mars relative to its rotation axis after magnetization was acquired. Such a reorientation may have been stimulated by internal mass redistributions associated with the formation of the northern lowlands and Tharsis, for example. A comparison of the mean paleo (magnetic) equator to the global distribution of crustal fields shows that magnetic anomalies tend to occur at low paleolatitudes. The same appears to be true for the Noachian-aged valley networks, which exhibit a broad spatial correlation with the magnetic anomalies. A possible interpretation is that the formation of magnetic anomalies and the valley networks was favored in the tropics where melting of water ice and snow was a stronger source of both surface valley erosion and groundwater recharge during the earliest history of the planet. This would be consistent with models in which hydrothermal alteration of crustal rocks played a role in producing the unusually strong martian magnetic anomalies.     

East-west trending magnetic anomalies in the Southern Hemisphere of Mars: Modeling analysis and interpretation

Maps of the vector components of the martian crustal magnetic field over the strongly magnetized Terra Cimmeria/Sirenum region are constructed using Mars Global Surveyor magnetometer data. Although pronounced east-west trending anomalies are present on the radial and north field component maps at the mapping altitude (∼360-380 km), these trends are much less prominent at the lower aerobraking altitude (∼90-150 km). Comparisons with similar maps produced using artificial data at the aerobraking altitude indicate that elongated sources in this region may have maximum lengths along the martian surface of ∼500 km and maximum aspect ratios of ∼2. Iterative forward modeling of several relatively isolated anomalies in the mapped region yields paleomagnetic pole positions consistent with those estimated in previous studies of other anomalies using mapping phase and science phasing orbit data. On this basis, it is inferred that sources in the studied region are most probably magnetized primarily in northward or southward directions. Using this additional constraint, iterative forward modeling is then applied to determine a magnetization distribution that is consistent with data at both the aerobraking altitude and the mapping altitude. The model magnetization distribution, which includes 41 discrete sources, again indicates no highly elongated sources. An examination of surface geology in the region as well as a consideration of the global distribution of anomalies suggests that magmatic intrusions (e.g., subsurface dike swarms), cooling in the presence of water, are the most likely sources of the magnetic anomalies.     

Photometric anomalies of the lunar surface studied with SMART-1 AMIE data

We present new results from the mapping of lunar photometric function parameters using images acquired by the spacecraft SMART-1 (European Space Agency). The source data for selected lunar areas imaged by the AMIE camera of SMART-1 and the data processing are described. We interpret the behavior of photometric function in terms of lunar regolith properties. Our study reveals photometric anomalies on both small (sub-kilometer) and large (tens of kilometers) scales. We found the regolith mesoscale roughness of lunar swirls to be similar in Mare Marginis, Mare Ingenii, and the surrounding terrains. Unique photometric properties related to peculiarities of the millimeter-scale regolith structure for the Reiner Gamma swirl are confirmed. We identified several impact craters of subkilometer sizes as the source of photometric anomalies created by an increase in mesoscale roughness within the proximal crater ejecta zones. The extended ray systems reveal differences in the photometric properties between proximal and distant ejecta blankets. Basaltic lava flows within Mare Imbrium and Oceanus Procellarum indicate higher regolith porosity for the redder soils due to differences in the chemical composition of lavas.     

An improved crustal magnetic field map of Mars from electron reflectometry: Highland volcano magmatic history and the end of the martian dynamo

We apply improved kinetic modeling of electron transport in the martian thermosphere to fit pitch angle distributions measured by the Mars Global Surveyor (MGS) Magnetometer/Electron Reflectometer (MAG/ER), together with appropriate filtering, binning, averaging and error correction techniques, to create the most reliable ER global map to date of crustal magnetic field magnitude at 185 km altitude, with twice the spatial resolution and considerably higher sensitivity to crustal fields than global maps of magnetic field components produced with MAG data alone. This map compares favorably to sparsely sampled dayside MAG data taken at similar altitudes, insofar as a direct comparison is meaningful. Using this map, we present two case studies. The first compares the magnetic signatures of two highland volcanoes, concluding that the comparatively greater thermal demagnetization at Syrtis Major compared with Tyrrhena Patera is likely due to a higher ratio of intruded to extruded magmas. The second uses the map along with topographic data to compare the magnetic signatures and crater retention ages of the demagnetized Hellas impact basin and magnetized Ladon impact basin. From this comparison, we determine that the martian global dynamo magnetic field went from substantial to very weak or nonexistent in the absolute model age time interval 4.15±0.05 to 4.07±0.05 Ga ago.     

Permeability of JSC-1A: A lunar soil simulant

The permeability of lunar soil simulant, JSC-1A, is measured over a range of bulk densities from 1550 to 2000 kg m⁻³. The corresponding viscous flow permeability is 1 × 10⁻¹² m² to 6.1 × 10⁻¹² m² for this bulk density range. Implications of these values on the contamination of regolith by rockets, on barrier/enhancement to bulk flow of ice, and on cratering are discussed. Although the particle size and shape distribution of the JSC-1A are extremely wide, the permeability measurements agree surprisingly well with the Carman-Kozeny equation. The results provide evidence that the Carman-Kozeny model could be applicable to other naturally occurring soils if effective soil properties are considered.     

The role of magnetic buoyancy in a Babcock-Leighton type solar dynamo

We study the effects of incorporating magnetic buoyancy in a model of the solar dynamo&amp;#x2014;which draws inspiration from the Babcock-Leighton idea of surface processes generating the poloidal field. We present our main results here.     

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.