Formation of lunar mare domes along crustal fractures: Rheologic conditions, dimensions of feeder dikes, and the role of magma evolution

In this study we examine a set of lunar mare domes located in the Hortensius/Milichius/T. Mayer region and in northern Mare Tranquillitatis with respect to their formation along crustal fractures, their rheologic properties, the dimensions of their feeder dikes, and the importance of magma evolution processes during dome formation. Many of these domes display elongated summit vents oriented radially with respect to major impact basins, and several dome locations are also aligned in these preferential directions. Analysis of Clementine UV/VIS and Lunar Prospector gamma ray spectrometer data reveals that the examined mare domes formed from low-Si basaltic lavas of high FeO and low to moderate TiO₂ content. Based on their morphometric properties (diameter, height, volume) obtained by photoclinometric and shape from shading analysis of telescopic CCD images, we derive rheologic quantities (lava viscosity during eruption, effusion rate, duration of the effusion process, magma rise speed) and the dimensions of the feeder dikes. We establish three rheologic groups characterised by specific combinations of rheologic properties and dike dimensions, where the most relevant discriminative parameter is the lava viscosity η. The first group is characterised by and contains the domes with elongated vents in the Milichius/T. Mayer region and two similar domes in northern Mare Tranquillitatis. The second group with comprises the very low aligned domes in northern Mare Tranquillitatis, and the third group with the relatively steep domes near Hortensius and in the T. Mayer region. The inferred dike dimensions in comparison to lunar crustal thickness data indicate that the source regions of the feeder dikes are situated within the upper crust for six of the domes in northern Mare Tranquillitatis, while they are likely to be located in the lower crust and in the upper mantle for the other examined domes. By comparing the time scale of magma ascent with the time scale on which heat is conducted from the magma into the host rock, we find evidence that the importance of magma evolution processes during ascent such as cooling and crystallisation increases with lava viscosity. We conclude that different degrees of evolution of initially fluid basaltic magma are able to explain the broad range of lava viscosities inferred for the examined mare domes. The spectral data reveal that differences in TiO₂ content may additionally account for the systematic difference in lava viscosity between the two examined lunar regions. We show that the described mechanisms are likely to be valid also for other lunar mare domes situated near Cauchy and Arago, regarded for comparison. On the other hand, we find for the Gruithuisen and Mairan highland domes that despite their inferred high lava viscosities of , no significant magma cooling in the dike occurred during ascent, supporting previous findings that the highland domes were formed during a specific phase of non-mare volcanism by highly silicic viscous lavas.     

A lunar cryptomare dome near Mee H and Drebbel F

In this study, we examine the lunar mare dome Mee 1 situated near the craters Mee H and Drebbel F in a region showing evidence of ancient (pre-Orientale) mare volcanism and cryptomare deposits. Regional stratigraphic relations indicate that Mee 1 was formed prior to the Orientale impact at the beginning of the Imbrian period. Based on a combined photoclinometry and shape from shading technique applied to telescopic CCD images of the dome acquired under oblique illumination, we determined a diameter of Mee 1 of 25 km, a height of 250 m, a flank slope of 1.15^°, and a volume of . Based on rheologic modelling of the dome and a viscoelastic model of the feeder dike, we obtained a magma viscosity of , an effusion rate of , a duration of the effusion process of 1.6 years, a magma rise speed of , a width of the feeder dike of 32 m, and a horizontal dike length of 144 km. A comparison of Mee 1 with domes with similar morphometric properties, which are located near Milichius and inside the crater Petavius, reveals strong similarities with respect to the viscosity of the dome-forming magma and the feeder dike geometry, while the effusion rate and magma rise speed of Mee 1 are somewhat higher. The pronounced morphometric differences between Mee 1 and a smaller dome situated close to the crater Doppelmayer and characterised by a similar magma viscosity suggest that the growth of that dome was limited by exhaustion of the magma reservoir, while Mee 1 and the other larger domes display morphometric properties presumably coming closer to the cooling limit. The comparison of the ancient dome Mee 1 with the younger (Eratosthenian) edifices near Milichius and Doppelmayer suggests that the conditions in the upper mantle and the crust favoured high eruption volumes, effusion rates, and magma rise speeds, implying the occurrence of large magma reservoirs preventing the limitation of dome growth by magma exhaustion. On the other hand, we observe similar general morphometric, rheologic, and feeder dike characteristics and, thus, conclude that the formation conditions of lunar mare domes did not change fundamentally during the Imbrian period.     

Lunar intrusive domes: Morphometric analysis and laccolith modelling

This study examines a set of lunar domes with very low flank slopes which differ in several respects from the frequently occurring lunar effusive domes. Some of these domes are exceptionally large, and most of them are associated with faults or linear rilles of presumably tensional origin. Accordingly, they might be interpreted as surface manifestations of laccolithic intrusions formed by flexure-induced vertical uplift of the lunar crust (or, alternatively, as low effusive edifices due to lava mantling of highland terrain, or kipukas, or structural features). All of them are situated near the borders of mare regions or in regions characterised by extensive effusive volcanic activity. Clementine multispectral UVVIS imagery indicates that they do not preferentially occur in specific types of mare basalt. Our determination of their morphometric properties, involving a combined photoclinometry and shape from shading technique applied to telescopic CCD images acquired at oblique illumination, reveals large dome diameters between 10 and more than 30 km, flank slopes below 0.9°, and volumes ranging from 0.5 to 50 km³. We establish three morphometric classes. The first class, In1, comprises large domes with diameters above 25 km and flank slopes of 0.2°-0.6°, class In2 is made up by smaller and slightly steeper domes with diameters of 10-15 km and flank slopes between 0.4° and 0.9°, and domes of class In3 have diameters of 13-20 km and flank slopes below 0.3°. While the morphometric properties of several candidate intrusive domes overlap with those of some classes of effusive domes, we show that a possible distinction criterion are the characteristic elongated outlines of the candidate intrusive domes. We examine how they differ from typical effusive domes of classes 5 and 6 defined by Head and Gifford [Head, J.W., Gifford, A., 1980. Lunar mare domes: classification and modes of origin. Moon Planets 22, 235-257], and show that they are likely no highland kipukas due to the absence of spectral contrast to their surrounding. These considerations serve as a motivation for an analysis of the candidate intrusive domes in terms of the laccolith model by Kerr and Pollard [Kerr, A.D., Pollard, D.D., 1998. Toward more realistic formulations for the analysis of laccoliths. J. Struct. Geol. 20(12), 1783-1793], to estimate the geophysical parameters, especially the intrusion depth and the magma pressure, which would result from the observed morphometric properties. Accordingly, domes of class In1 are characterised by intrusion depths of 2.3-3.5 km and magma pressures between 18 and 29 MPa. For the smaller and steeper domes of class In2 the magma intruded to shallow depths between 0.4 and 1.0 km while the inferred magma pressures range from 3 to 8 MPa. Class In3 domes are similar to those of class In1 with intrusion depths of 1.8-2.7 km and magma pressures of 15-23 MPa. As an extraordinary feature, we describe in some detail the concentric crater Archytas G associated with the intrusive dome Ar1 and discuss possible modes of origin. In comparison to the candidate intrusive domes, terrestrial laccoliths tend to be smaller, but it remains unclear if this observation is merely a selection effect due to the limited resolution of our telescopic CCD images. An elongated outline is common to many terrestrial laccoliths and the putative lunar laccoliths, while the thickness values measured for terrestrial laccoliths are typically higher than those inferred for lunar laccoliths, but the typical intrusion depths are comparable.     

Lunar mare domes: Classification and modes of origin

A classification of over 200 lunar mare domes shows that they have two major modes of occurrence: (1) low, flat, generally circular structures with convex shapes, slopes less than about 5°, and displaying summit craters, and (2) irregular structures often adjacent to highland regions and rarely containing summit craters. On the basis of morphologic and morphometric similarities, the first mode of occurrence appears to be analogous to small terrestrial shield volcanoes, and to represent primary volcanic constructs, while the second class of domes appears to result from secondary volcanic effects (flooding of highland material to produce kipukas and draping of lavas to produce irregular dome-like topography).Domes comparable to small shield volcanoes generally range from 3–17 km in diameter and up to several hundred meters in height and occur predominantly in groupings in the lunar equatorial region in northeast Tranquillitatis (Cauchy area), between Kepler and Copernicus (Hortensius area), and in the Marius Hills. In the Marius Hills, domes generally lack summit craters and have a rough surface texture formed in part by superposed cones and steep-sided flows. Elsewhere, domes representing volcanic sources are smooth-surfaced and usually contain a summit crater. These features are similar in general morphology to small terrestrial lava shields. They are generally intermediate in volume, slope, and height between small shields of terrestrial basaltic plains (such as the Snake River Plains) and larger Icelandic shields. Summit craters on lunar domes are considerably larger than craters on terrestrial shields of comparable diameters, apparently due to a combination of factors, including vent enlargement during extrusion, possibly higher lunar extrusion rates, different amounts of collapse, and impact erosion.Most vent-related domes appear to be associated with, and are thus approximately the same age as, surrounding lava plains, although relationships in specific areas have not yet been established. On the basis of age ranges of mare deposits established by Apollo samples, mare vent-related domes formed over an approximately one billion year period starting about 3.7 b.y. ago. Extrusion rates were apparently relatively low compared to the very high values characteristic of flows associated with major lunar sinuous rilles and terrestrial flood basalts, but may have been relatively high compared to similar terrestrial shields. Large shield volcanoes equivalent to the terrestrial Hawaiian-type or to the martian edifices such as Olympus Mons, do not occur on the Moon. Lack of these features may be due to the low viscosities and high effusion rates typical of many lunar eruptions and the lack of continuous eruptions from single sources.Paper presented at the European Workshop on Planetary Sciences, organised by the Laboratorio di Astrofisica Spaziale di Frascati, and held between April 23–27, 1979, at the Accademia Nazionale del Lincei in Rome, Italy.     

Duration and extent of lunar volcanism: Comparison of 3D convection models to mare basalt ages

It is widely accepted that lunar volcanism started before the emplacement of the mare fills ( b.p.) and lasted for probably more than 3.0 Ga. While the early volcanic activity is relatively easy to understand from a thermal point of view, the late stages of volcanism are harder to explain, because a relatively small body like the Earth's Moon is expected to cool rapidly and any molten layer in the interior should solidify rather quickly. We present several thermal evolution models, in which we varied the boundary conditions at the model surface in order to evaluate the influence on the extent and lifetime of a molten layer in the lunar interior. To investigate the influence of a top insulating layer we used a fully three-dimensional spherical shell convection code for the modelling of the lunar thermal history. In all our models, a partial melt zone formed nearly immediately after the simulation started (early in lunar history), consistent with the identification of lunar cryptomare and early mare basalt volcanism on the Moon. Due to the characteristic thickening of the Moon's lithosphere the melt zone solidified from above. This suggests that the source regions of volcanic rock material proceeded to increasing depth with time. The rapid growth of a massive lithosphere kept the Moon's interior warm and prevented the melt zone from fast freezing. The lifetimes of the melt zones derived from our models are consistent with basalt ages obtained from crater chronology. We conclude that an insulating megaregolith layer is sufficient to prevent the interior from fast cooling, allowing for the thermal regime necessary for the production and eruption of young lava flows in Oceanus Procellarum.     

Lava tubes and channels in the lunar Marius Hills

The Marius Hills region, a volcanic plateau in Oceanus Procellarum, contains numerous rilles, rille-like structures, and chains of elongate craters. Most of these structures characteristically: (1) originate on or near irregular shaped craters associated with features previously interpreted as volcanic domes, (2) trend downslope onto Plateau Plains, (3) generally taper in width and become shallower, (4) are often discountinuous, (5) occupy the center, or apparent crest of a broad ridge, (6) may contain cut-off branches and distributary structures, and (7) may have local reversals in longitudinal slope. Structures having these characteristics are interpreted to be lava channels or partly collapsed lava tubes. Terrestrial lava tubes form exclusively, and commonly, in fluid basalt flows. Recent evidence indicates that viscosities of lunar mare basalt lava flows were conducive for lava tube formation.Terrestrial analogs are offered for structures described in the Marius Hills. The analogs are comparable in qualitative and quantitative geomorphic aspects, excluding that of width. Scaling consideration of lunar reduced gravity accounts for increased width of the lunar structures. Linear and curvilinear rilles trending along equal elevations are interpreted to result from fracturing or faulting.     

Individual lava flow thicknesses in Oceanus Procellarum and Mare Serenitatis determined from Clementine multispectral data

We use multispectral reflectance data from the lunar Clementine mission to investigate the impact ejecta deposits of simple craters in two separate lunar mare basalt regions, one in Oceanus Procellarum and one in Mare Serenitatis. Over 100 impact craters are studied, and for a number of these we observe differences between the TiO₂ (and FeO) contents of their ejecta deposits and the lava flow units in which they are located. We demonstrate that, in the majority of cases, these differences cannot plausibly be attributed to uncorrected maturity effects. These observations, coupled with morphometric crater relationships that provide maximum crater excavation depths, allow the investigation of sub-surface lava flow stratigraphy. We provide estimated average thicknesses for a number of lava flow units in the two study regions, ranging from ∼80 m to ∼600 m. In the case of the Serenitatis study area, our results are consistent with the presence of sub-surface horizons inferred from recent radar sounding measurements from the JAXA Kaguya spacecraft. The average lava flow thicknesses we obtain are used to make estimates of the average flux of volcanic material in these regions. These are in broad agreement with previous studies, suggesting that the variation in mare basalt types we observe with depth is similar to the lateral variations identified at the surface.     

Tungsten isotopes in ferroan anorthosites: Implications for the age of the Moon and lifetime of its magma ocean

New W isotope data for ferroan anorthosites 60025 and 62255 and low-Ti mare basalt 15555 show that these samples, contrary to previous reports [Lee, D.C., et al., 1997. Science 278, 1098-1103; Lee, D.C., et al., 2002. Earth Planet. Sci. Lett. 198, 267-274], have a W isotope composition that is indistinguishable from KREEP and other mare basalts. This requires crust extraction on the Moon later than ∼60 Myr after CAI formation, consistent with ¹⁴⁷Sm-¹⁴³Nd ages for ferroan anorthosites. The absence of ¹⁸²Hf-induced ¹⁸²W variations in the Moon is consistent with formation of the Moon at after CAI formation that has been inferred based on the indistinguishable ¹⁸²W/¹⁸⁴W ratios of the bulk Moon and the bulk silicate Earth. The uncertainties on the age of the Moon and the age of the oldest lunar samples currently hamper a precise determination of the duration of magma ocean solidification and are consistent with both an almost immediate crystallization and a more protracted timescale of ∼100 Myr.     

Rümker hills: a lunar volcanic dome complex

The Rümker Hills, a volcanic dome-flow complex in the northern Oceanus Procellarum, is characterized by overlapping plains-forming units with lobate scarps, volcanic domes, a 60 km ring, and a scarp which separates the plateau from surrounding mare materials. Plains-forming units are interpreted as fluid volcanic flows, and domes as viscous extrusions. One dome may be a stratovolcano. The ring system is discordant with regional structural trends and probably has a local origin. The Rümker Hills is the closest lunar analog to the large martian shield structures revealed on the Mariner 9 photographs of Mars.     

Meridional variations of temperature, C2H2 and C2H6 abundances in Saturn's stratosphere at southern summer solstice

Measurements of the vertical and latitudinal variations of temperature and C₂H₂ and C₂H₆ abundances in the stratosphere of Saturn can be used as stringent constraints on seasonal climate models, photochemical models, and dynamics. The summertime photochemical loss timescale for C₂H₆ in Saturn's middle and lower stratosphere (∼40-10,000 years, depending on altitude and latitude) is much greater than the atmospheric transport timescale; ethane observations may therefore be used to trace stratospheric dynamics. The shorter chemical lifetime for C₂H₂ (∼1-7 years depending on altitude and latitude) makes the acetylene abundance less sensitive to transport effects and more sensitive to insolation and seasonal effects. To obtain information on the temperature and hydrocarbon abundance distributions in Saturn's stratosphere, high-resolution spectral observations were obtained on September 13-14, 2002 UT at NASA's IRTF using the mid-infrared TEXES grating spectrograph. At the time of the observations, Saturn was at a LS≈270°, corresponding to Saturn's southern summer solstice. The observed spectra exhibit a strong increase in the strength of methane emission at 1230 cm⁻¹ with increasing southern latitude. Line-by-line radiative transfer calculations indicate that a temperature increase in the stratosphere of ≈10 K from the equator to the south pole between 10 and 0.01 mbar is implied. Similar observations of acetylene and ethane were also recorded. We find the 1.16 mbar mixing ratio of C₂H₂ at −1° and −83° planetocentric latitude to be and , respectively. The C₂H₂ mixing ratio at 0.12 mbar is found to be at −1° planetocentric latitude and at −83° planetocentric latitude. The 2.3 mbar mixing ratio of C₂H₆ inferred from the data is and at −1° and −83° planetocentric latitude, respectively. Further observations, creating a time baseline, will be required to completely resolve the question of how much the latitudinal variations of C₂H₂ and C₂H₆ are affected by seasonal forcing and/or stratospheric circulation.     

New estimates for Io eruption temperatures: Implications for the interior

The initial interpretation of Galileo data from Jupiter's moon, Io, suggested eruption temperatures . Tidal heating models have difficulties explaining Io's prodigious heat flow if the mantle is , although we suggest that temperatures up to may be possible. In general, Io eruption temperatures have been overestimated because the incorrect thermal model has been applied. Much of the thermal emission from high-temperature hot spots comes from lava fountains but lava flow models were utilized. We apply a new lava fountain model to the highest reported eruption temperature, the SSI observation of the 1997 eruption at Pillan. This resets the lower temperature limit for the eruption from 1600 to . Additionally, viscous heating of the magma may have increased eruption temperature by as a result of the strong compressive stresses in the ionian lithosphere. While further work is needed, it appears that the discrepancy between observations and interior models is largely resolved.     

The orbit, mass, and albedo of transneptunian binary (66652) 1999 RZ253

We have observed (66652) 1999 RZ₂₅₃ with the Hubble Space Telescope at seven separate epochs and have fit an orbit to the observed relative positions of this binary. Two orbital solutions have been identified that differ primarily in the inclination of the orbit plane. The best fit corresponds to an orbital period, days, semimajor axis a=4660±170 km and orbital eccentricity e=0.460±0.013 corresponding to a system mass m=3.7±0.4×¹⁰¹⁸ kg. For a density of the albedo at 477 nm is p₄₇₇=0.12±0.01, significantly higher than has been commonly assumed for objects in the Kuiper belt. Multicolor, multiepoch photometry shows this pair to have colors typical for the Kuiper belt with a spectral gradient of 0.35 per 100 nm in the range between 475 and 775 nm. Photometric variations at the four epochs we observed were as large as 12±3% but the sampling is insufficient to confirm the existence of a lightcurve.     

Spectral unmixing for mineral identification in pancam images of soils in Gusev crater, Mars

The objective of this work is to propose an automated unmixing technique for the analysis of 11-channel Mars Exploration Rover Panoramic Camera (MER/Pancam) spectra. Our approach is to provide a screening tool for identifying unique/distinct reflectance spectra. We demonstrate the utility of this unmixing technique in a study of the mineralogy of the bright salty soils exposed by the rover wheels in images of Gusev crater regions known as Paso Robles (Sols 400,426), Arad (Sol 721), and Tyrone (Sol 790). The unmixing algorithm is based on a novel derivation of the Nonnegative Matrix Factorization technique and includes added features that preclude the adverse effects of low abundance materials that would otherwise skew the unmixing. In order to create a full 11-channel spectrum out of the left and right eye stereo pairs, we also developed a new registration procedure that includes rectification and disparity calculation of the images. We identified two classes of endmember spectra for the bright soils imaged on Sols 426 and 790. One of these endmember classes is also observed for soils imaged on Sols 400 and 721 and has a unique spectral shape that is distinct from most iron oxide, sulfate and silicate spectra and differs from typical martian surface spectra. Instead, its unique spectral character resembles the spectral shape of the ferric sulfate minerals fibroferrite (Fe³⁺(SO₄)(OH) · 5H₂O) and ferricopiapite and the phosphate mineral ferristrunzite . The other endmember class is less consistent with specific minerals and is likely a mixture of altered volcanic material and some bright salts. Further analyses of data from Sols 400 and 790 using an anomaly detection algorithm as a tool for detecting low abundance materials additionally suggests the identification of the sulfate mineral paracoquimbite (Fe₂(SO₄)₃ · 9H₂O). This spectral study of Pancam images of the bright S- and P-enriched soils of Gusev crater identifies specific ferric sulfate and ferric phosphate minerals that are consistent with the unique spectral properties observed here in the 0.4-1 μm range.     

Lunar surface agglutinates: Mapping composition anomalies

From the Clementine UVVIS imagery of the lunar surface, the abundance of agglutinates in the lunar regolith and their composition in terms of FeO and Al₂O₃ oxides have been predicted. Data on the spectral, chemical, and mineralogic measurements of about 30 lunar soil samples from the Lunar Samples Characterization Consortium (LSCC) collection were used. The fulfilled prognosis confirms that the mare agglutinates are enriched in Al₂O₃ and depleted of FeO, while the highland agglutinates are depleted of Al₂O₃ and enriched in FeO. This behavior can be caused by the global transport of the lunar surface material induced by cosmogenic factors.     

Identification,distribution and significance of lunar volcanic domes

Over 300 previously unrecognized volcanic domes were identified on Lunar Orbiter photographs using the following criteria: (1) the recognition of land forms on the Moon similar in morphology to terrestrial volcanic domes, (2) structural control, (3) geomorphic discordance, and (4) the recognition of land forms modified by dome-like swellings. Many terrestrial volcanic domes are similar in morphology to lunar domes. This analogy suggests that some lunar hills are in fact extrusive volcanic domes. Many of the domes identified in this paper seem to be related to basins and craters, and with the exception of local tectonic grid control few domes are related to any observable Moon-wide pattern. Domes are not uniquely found on maria. Dome formation probably spans a wide range of lunar time and activity in areas where domes are located may be continuing to the present as revealed by the close correlation of dome distribution with the distribution of lunar transient events. The overall morphology of a lunar dome is a poor indicator of the composition of the rock that forms the dome. Contribution No. 33 Planetary Geology Group, University of New Mexico.