The great thickness debate: Ice shell thickness models for Europa and comparisons with estimates based on flexure at ridges

Estimates of the thickness of the ice shell of Europa range from <1 to >30 km. The higher values are generally assumed to be estimates of the entire ice shell thickness, which may include a lower ductile layer of ice, whereas many of the smaller thickness estimates are based on analyses that only consider that portion of the ice layer that behaves elastically at a particular strain rate. One example of the latter is flexure analysis, in which the elastic ice layer is modeled as a plate or sphere that is flexed under the weight of a surface load. We present calculations based on flexure analysis in which we model the elastic ice layer as flexing under a line-load caused by ridges. We use precisely located, parallel flanking cracks as indicators of the location of greatest tensile stress induced by flexure. Our elastic thickness results are spatially variable: ∼500-2200 m (two sites) and ∼200-1000 m (one site). Thorough analysis of Europan flexure studies performed by various researchers shows that the type of model selected causes the greatest variability in the thickness results, followed by the choice of Young's modulus, which is poorly constrained for the Europan ice shell. Comparing our results to those of previously published flexure analyses for Europa, we infer spatial variability in the elastic ice thickness (at the time of load emplacement), with smooth bands having the thinnest elastic ice thickness of all areas studied. Because analysis of flexure-induced fracturing can only reveal the elastic thickness at the time of load emplacement, calculated thickness variability between features having different ages may also reflect a temporal variability in the thickness of Europa's ice shell.     

Effects of plasticity on convection in an ice shell: Implications for Europa

Europa's surface exhibits numerous pits, uplifts, and disrupted chaos terrains that have been suggested to result from convection in the ice shell. To test this hypothesis, we present numerical simulations of convection in an ice shell including the effects of plasticity, which provides a simple continuum representation for brittle or semibrittle deformation along discrete fractures. Plastic deformation occurs when stresses reach a specified yield stress; at lower stresses, the fluid flow follows a Newtonian, temperature-dependent viscosity. Four distinct modes of behavior can occur. For yield stresses exceeding ∼1 bar, plastic effects are negligible and stagnant-lid convection, with no surface motion and minimal topography, results. At intermediate yield stresses, a stagnant lid forms but deforms plastically, leading to surface velocities up to several millimeters per year. Slightly smaller yield stresses allow episodic, catastrophic overturns of the upper conductive lid, with (transient) stagnant lids forming in between overturn events. The smallest yield stresses allow continual recycling of the upper lid, with simultaneous, gradual ascent of warm ice to the surface and descent of cold, near-surface ice into the interior. The exact yield stresses over which each regime occurs depend on the ice-shell thickness, melting-temperature viscosity, and activation energy for viscous creep. To form hummocky matrix and translate chaos plates by several kilometers, substantial surface strain must accompany chaos formation, and the three plasticity-dominated convection modes described here can provide such deformation. Our simulations suggest that, if yield stresses of ∼0.2-1 bar are relevant to Europa, then convection in Europa's ice shell can produce chaos-like structures at the surface. However, our simulations have difficulty explaining Europa's numerous pits and uplifts. When plasticity forces the upper lid to participate in the convection, dynamic topography of ∼50-100-m amplitude results, but the topographic structures generally have diameters of 30-100 km, an order of magnitude wider than typical pits and uplifts. None of our simulations produced isolated pits or uplifts of any diameter.     

The chemical nature of Europa surface material and the relation to a subsurface ocean

The surface composition of Europa is of special interest due to the information it might provide regarding the presence of a subsurface ocean. One source of this information is the infrared reflectance spectrum. Certain surface regions of Europa exhibit distorted H₂O vibrational overtone bands in the 1.5 and 2.0 μm region, as measured by the Galileo mission Near Infrared Mapping Spectrometer (NIMS). These bands are clearly the result of highly concentrated solvated contaminants. However, two interpretations of their identity have been presented. One emphasizes hydrated salt minerals and the other sulfuric acid, although each does not specifically rule out some of the other. It has been pointed out that accurate chemical identification of the surface composition must depend on integrating spectral data with geochemical models, and information on the tenuous atmosphere sputtered from the surface. It is also extremely important to apply detailed chemistry when interpreting the spectral data, including knowledge of mineral dissolution chemistry and the subsequent optical signatures of ion solvation in low-temperature ice. We present studies of flash frozen acid and salt mixtures as Europa surface analogs and demonstrate that solvated protons, metal cations and inorganic anions all influence the spectra and must all, collectively, be considered when assigning Europa spectral features. These laboratory data show best correlation with NIMS Europa spectra for multi-component mixtures of sodium and magnesium bearing sulfate salts mixed with sulfuric acid. The data provide a concentration upper bound of 50-mol% for MgSO₄ and 40-mol% for Na₂SO₄. This newly reported higher sodium and proton content is consistent with low-temperature aqueous differentiation and hydrothermal processing of carbonaceous chondrite-forming materials during the formation and early evolution of Europa.     

Surficial properties in Melas Chasma, Mars, from Mars Odyssey THEMIS data

We examine the nature of the surface layer in a small area of the Melas Chasma region as determined from high-resolution thermal and visible Mars Odyssey Thermal Emission Imaging System (THEMIS) data as well as how our conclusions compare to past analyses. At THEMIS resolution, the thermal structure is dominated by local control and all significant thermal variations can be linked to morphology. Thus, THEMIS provides us with detailed images that contain thermophysical information as well, allowing us to create a surficial geologic map intended to reflect the surface structure of the region. Eight units have been defined: (i) blanketed plateaus with thermally distinct craters and fractures, (ii) blanketed canyon walls with rocky edges, (iii) indurated and/or rocky canyon wall slide material partially covered by aeolian material, (iv) an anomalous wall region with fluvial-like depressions partially filled with particulate material, (v) indurated and/or rocky ridged and non-ridged canyon floor landslide material mingled with aeolian material, (vi) sand sheets, (vii) indurated and/or rocky rounded blocks intermingled with small areas of aeolian material, and (viii) transverse dunes. The THEMIS thermal data support conclusions from previous studies but also reveal much more structure than was seen in the past. We have found that all significant thermal variations in this region can be linked to morphology but all morphological variations cannot be linked to significant thermal variations. THEMIS visible images provide an intermediate resolution that bridges the gap between MOC and Viking and allow for a more meaningful interpretation of the geologic context of a region. Surfaces indicate that landslides were an important geologic process long ago, shaping the canyon walls and floor, while aeolian processes have subsequently altered the surface layer in many locations and may still be active.     

Properties of the Hermean regolith: IV. Photometric parameters of Mercury and the Moon contrasted with Hapke modelling

A comparison of the photometric properties of Mercury and the Moon is performed, based on their integral phase curves and disk-resolved image data of Mercury obtained with the Swedish Vacuum Solar Telescope. Proper absolute calibration of integral V-band magnitude observations reveals that the near-side of the Moon is 10-15% brighter than average Mercury, and 0-5% brighter for the “bolometric” wavelength range 400-1000 nm. As shown, this is supported by recent estimates of their geometric albedos. Hapke photometric parameters of their surfaces are derived from identical approaches, allowing a contrasting study between their surface properties to be performed. Compared to the average near-side Moon, Mercury has a slightly lower single-scattering albedo, an opposition surge with smaller width and of marginally smaller amplitude, and a somewhat smoother surface with similar porosity. The width of the lobes of the single-particle scattering function are smaller for Mercury, and the backward scattering anisotropy is stronger. In terms of the double Henyey-Greenstein b-c parameter plot, the scattering properties of an average particle on Mercury is closer to the properties of lunar maria than highlands, indicating a higher density of internal scatterers than that of lunar particles. The photometric roughness of Mercury is well constrained by the recent study of Mallama et al. (2002, Icarus 155, 253-264) to a value of about 8°, suggesting that the surfaces sampled by the highest phase angle observations (Borealis, Susei, and Sobkou Planitia) are lunar mare-like in their textural properties. However, Mariner 10 disk brightness profiles obtained at intermediate phase angles indicate a surface roughness of about twice this value. The photometric parameters of the Moon are more difficult to constrain due to limited phase angle coverage, but the best Hapke fits are provided by rather small surface roughnesses. Better-calibrated, multiple-wavelength observations of the integral and disk-resolved brightnesses of both bodies, and obtained at higher phase angle values in the case of the Moon, are urgently needed to arrive at a more consistent picture of the contrasting light scattering properties of their surfaces.     

Iapetus dark and bright material:: giving compositional interpretation some latitude

Narrowband reflectance spectra (0.53-1.0 μm) of Iapetus' leading and trailing sides were obtained in 2000 to test the presence of an absorption feature located near 0.67 μm seen in reflectance spectra of Iapetus' dark material and Hyperion's surface material. No feature was observed. The difference in reflectance across the UV/VIS/NIR spectral region, and the dependence of the presence or absence of this absorption feature on angular separation from the apex of Iapetus in its orbit, phase angle, and heliocentric distance (affecting temperature), were examined. A trend of increased reddening, and the presence of the absorption feature, correlate with an angular separation from the apex of ? approximately 10°. Spectral information is lost when the contribution of the bright water ice signal to the reflectance spectrum increases sufficiently. In order to optimize compositional studies of Iapetus, we encourage future ground-based and space-based spectral observations to maximize the concentration of dark material in the instrumental field of view.     

Surface changes on Io during the Galileo mission

A careful survey of Galileo SSI global monitoring images revealed more than 80 apparent surface changes that took place on Io during the 5 year period of observation, ranging from giant plume deposits to subtle changes in the color or albedo of patera surfaces. Explosive volcanic activity was discovered at four previously unrecognized centers: an unnamed patera to the south of Karei that produced a Pele-sized red ring, a patera to the west of Zal that produced a small circular bright deposit, a large orange ring detected near the north pole of Io, and a small bright ring near Io's south pole. Only a handful of Io's many active volcanoes produced large scale explosive eruptions, and several of these erupted repeatedly, leaving at least 83% of Io's surface unaltered throughout the Galileo mission. Most of the hot spots detected from SSI, NIMS and ground-based thermal observations caused no noticeable surface changes greater than 10 km in extent over the five year period. Surface changes were found at every location where active plumes were identified, including Acala which was never seen in sunlight and was only detected through auroral emissions during eclipse. Two types of plumes are distinguished on the basis of the size and color of their deposits, confirming post-Voyager suggestions by McEwen and Soderblom [Icarus 55 (1983) 191]. Smaller plumes produce near-circular rings typically 150-200 km in radius that are white or yellow in color unless contaminated with silicates, and frequently coat their surroundings with frosts of fine-grained SO₂. The larger plumes are much less numerous, limited to a half dozen examples, and produce oval, orange or red, sulfur-rich rings with maximum radii in the north-south direction that are typically in the range from 500 to 550 km. Both types of plumes can be either episodic or quasi-continuous over a five year period. Repeated eruptions of the smaller SO₂-rich plumes likely contribute significantly to Io's resurfacing rate, whereas dust ejection is likely dominated by the tenuous giant plumes. Both types of plume deposits fade on time-scales of months to years through burial and alteration. Episodic seepages of SO₂ at Haemus Montes, Zal Montes, Dorian Montes, and the plateau to the north of Pillan Patera may have been triggered by activity at nearby volcanic centers.     

The nucleus of Comet Borrelly: a study of morphology and surface brightness

Stereo images obtained during the DS1 flyby were analyzed to derive a topographic model for the nucleus of Comet 19P/Borrelly for morphologic and photometric studies. The elongated nucleus has an overall concave shape, resembling a peanut, with the lower end tilted towards the camera. The bimodal character of surface-slopes and curvatures support the idea that the nucleus is a gravitational aggregate, consisting of two fragments in contact. Our photometric modeling suggests that topographic shading effects on Borrelly's surface are very minor (<10%) at the given resolution of the terrain model. Instead, albedo effects are thought to dominate Borrelly's large variations in surface brightness. With 90% of the visible surface having single scattering albedos between 0.008 and 0.024, Borrelly is confirmed to be among the darkest of the known Solar System objects. Photometrically corrected images emphasize that the nucleus has distinct, contiguous terrains covered with either bright or dark, smooth or mottled materials. Also, mapping of the changes in surface brightness with phase angle suggests that terrain roughness at subpixel scale is not uniform over the nucleus. High surface roughness is noted in particular near the transition between the upper and lower end of the nucleus, as well as near the presumed source region of Borrelly's main jets. Borrelly's surface is complex and characterized by distinct types of materials that have different compositional and/or physical properties.     

Photometry and models of eight near-Earth asteroids

We present new observations and models of the shapes and rotational states of the eight near-Earth Asteroids (1580) Betulia, (1627) Ivar, (1980) Tezcatlipoca, (2100) Ra-Shalom, (3199) Nefertiti, (3908) Nyx, (4957) Brucemurray, and (5587) 1990 SB. We also outline some of their solar phase curves, corrected to common reference geometry with the models. Some of the targets may feature sizable global nonconvexities, but the observable solar phase angles were not sufficiently high for confirming these. None is likely to have a very densely cratered surface. We discuss the role of the intermediate topographic scale range in photometry, and surmise that this scale range is less important than large or small scale lengths.     

Investigating the Vesta-vestoid-HED connection

We present spectral observations of Minor Planet 4 Vesta, of five V-type asteroids which are physically near Vesta, and of two V-type NEAs. We use these spectra to determine the presence or absence of a weak feature at 506.5 nm which is indicative of the presence of spin-forbidden Fe²⁺ in sixfold coordination. As with our earlier observations [Cochran and Vilas, Icarus 134 (1998) 207-212], we find this feature at all observed rotational phases of Vesta and again see the trend that spectra at longitudes between 240° and 360° have a smaller 506.5 nm feature equivalent width than spectra obtained at other longitudes. Additionally, we searched for this feature in V-class main-belt and NEA asteroids and positively detected the feature in main-belt Asteroid 2579 Spartacus and possibly in 3376 Armandhammer. The other objects lacked the feature. Our results are compared with previous observations of this feature by Vilas et al. [Icarus 147 (2000) 119-128]. The spatial distribution of the bodies as a function of the presence of this feature was investigated. We discuss the implication of the presence of this feature and the depth of the 0.9 μm pyroxene band for the scenario that pieces of Vesta were transported, via the 3:1 and ν₆ resonances, to the NEAs, and thence to inclusion in our meteorite collections as HED meteorites.