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.     

Resurfacing history of Europa from pole-to-pole geological mapping

We produced geologic maps from two regional mosaics of Galileo images across the leading and trailing hemispheres of Europa in order to investigate the temporal distribution of units in the visible geologic record. Five principal terrain types were identified (plains, bands, ridges, chaos, and crater materials), which are interpreted to result from (1) tectonic fracturing and lineament building, (2) cryovolcanic reworking of surface units, with possible emplacement of sub-surface materials, and (3) impact cratering. The geologic histories of both mapped areas are essentially similar and reflect some common trends: Tectonic resurfacing dominates the early geologic record with the formation of background plains by intricate superposition of lineaments, the opening of wide bands with infilling of inter-plate gaps, and the buildup of ridges and ridge complexes along prominent fractures in the ice. It also appears that lineaments are narrower and more widely spaced with time. The lack of impact craters overprinted by lineaments indicate that the degree of tectonic resurfacing decreased rapidly after ridged plains formation. In contrast, the degree of cryovolcanic resurfacing appears to increase with time, as chaos formation dominates the later parts of the geologic record. These trends, and the transition from tectonic- to cryovolcanic-dominated resurfacing could be attributed to the gradual thickening of Europa's cryosphere during the visible geologic history, that comprises the last 2% or 30-80 Myr of Europa's history: An originally thin, brittle ice shell could be pervasively fractured or melted through by tidal and endogenic processes; the degree of fracturing and plate displacements decreased with time in a thickening shell, and lineaments became narrower and more widely spaced; formation of chaos regions could have occurred where the thickness threshold for solid-state convection was exceeded, and can be aided by preferential tidal heating of more ductile ice. In a long-term context it is not clear at this point whether this inferred thickening trend would reflect a drastic change in the thermal evolution of the satellite, or cyclic or irregular episodes of tectonic and cryovolcanic activity.     

Segmented lineaments on Europa: Implications for the formation of ridge complexes and bright bands

We describe several segmented lineaments on Europa’s surface. These lineaments are extensive, stretching for 100s-1000s of km, and have ridge complex or bright band morphologies. The geometries of the segmented portions of these features are diagnostic of the remote normal and shear stress environment in which they formed and, therefore, constrain ridge complex and bright band formation mechanisms. Analysis of four ridge complexes indicates that they formed in a remote normal stress environment that was tensile and isotropic (or nearly so) and that these lineaments may have formed in a manner more analogous to bands on Europa than to ridges. The stress environment associated with these ridge complexes may also explain the anastomosing nature of their interior morphology. Analysis of two bright bands indicate that one formed in a remote normal stress environment that was tensile and the other was reactivated under a combination of remote tensile normal stress and remote sinistral shear stress. Aspects of the morphologies of these features also indicate that bright bands likely have complex deformation histories that can include multiple episodes of reactivation.     

Pits and uplifts on Europa

A survey of depression and uplift features on Europa, based on Galileo regional mapping images, shows that these features come in a wide range of sizes, with numbers increasing greatly with decreasing size, down to the limits of resolution. Size distributions are similar in the northern leading and southern trailing hemispheres, where they are distinctly different from the southern leading and northern trailing hemispheres, suggesting an oblique, antipodal symmetry pattern, similar to that of chaotic and tectonic terrain. This pattern is suggestive of polar wander. Uplifts are usually polygonal or irregular in shape and rarely are cracked. Patches of chaotic terrain, which we had surveyed earlier, are not included in the current study because their topography is generally unclear, and because there is no a priori known genetic linkage with the pits and uplifts.These results contradict generalizations based on the earlier “pits, spots, and domes” (PSD) taxonomy. Most of the type examples for PSDs were simply patches of chaotic terrain selected from a limited portion of their full size range. The use of the term lenticula to collectively describe PSDs is inconsistent with the IAU definition of lenticula: a small dark spot seen at low resolution. Pits and uplifts do not correlate with lenticulae, although chaos often does. Properties of PSDs that have been widely cited as primary evidence for convective upwelling in thick ice (e.g., that uplifts are generally dome-shaped and often cracked; that pits and domes are regularly spaced; that there is a typical diameter of ∼10 km) were premature and not supported by subsequent data. Most pits and uplifts are less than 10 km across so, if they formed by diapirism or convective upwelling, the sources must have been very shallow, less than 5 km deep. How they actually formed remains unknown.     

Galileo PPR observations of Europa: Hotspot detection limits and surface thermal properties

The Galileo photopolarimeter–radiometer (PPR) made over 100 observations of Europa’s surface temperature. We have used these data to constrain a diurnal thermal model and, thus, map the thermal inertia and bolometric albedo over 20% of the surface. We find an increased thermal inertia at mid-latitudes that is widespread in longitude and does not appear to correlate with geology, albedo, or other observables. Our derived thermophysical properties can be used to predict volatile stability across the surface over the course of a day and in planning of infrared instruments on future missions. Furthermore, while observations in the thermal infrared can and have been used to find endogenic activity, no such activity was detected at Europa. We have calculated the detection limits of these PPR observations and find that 100 km² hotspots with temperatures of 116–1200 K could exist undetected on the surface, depending on the location.     

The evil twin of Agenor: tectonic convergence on Europa

The dilemma of the surface-area budget on Europa is resolved by identification of sites of crustal convergence, which have balanced the continual and common creation of new surface along dilational bands and pull-aparts. Convergence bands are characterized by a distinctive, albeit subdued, morphology. The prominent, unusual lineament Agenor is one of several examples. We also find diametrically opposite Agenor a similar bright linear feature surrounded by markings that allow reconstruction, which shows it to be a convergence feature. Until recently, identification of convergence sites was difficult because these features are subtle and do not exhibit structures (like the Himalayas or plate subduction) familiar from convergence of thick solid crusts on terrestrial planets.     

Crack azimuths on Europa: the G1 lineament sequence revisited

The tectonic sequence in the anti-jovian area covered by regional mapping images from Galileo's orbit E15 is determined from a study of cross-cutting relationships among lineament features. The sequence is used to test earlier results from orbit G1, based on lower resolution images, which appeared to display a progressive change in azimuthal orientation over about 90° in a clockwise sense. Such a progression is consistent with expected stress variations that would accompany plausible non-synchronous rotation. The more recent data provide a more complete record than the G1 data did. We find that to fit the sequence into a continual clockwise change of orientation would require at least 1000° (>5 cycles) of azimuthal rotation. If due to non-synchronous rotation of Europa, this result implies that we are seeing back further into the tectonic record than the G1 results had suggested. The three sets of orientations found by Geissler et al. now appear to have been spaced out over several cycles, not during a fraction of one cycle. While our more complete sequence of lineament formation is consistent with non-synchronous rotation, a statistical test shows that it cannot be construed as independent evidence. Other lines of evidence do support non-synchronous rotation, but azimuths of crack sequences do not show it, probably because only a couple of cracks form in a given region in any given non-synchronous rotation period.     

Crack azimuths on Europa: time sequence in the southern leading face

The formation sequence of prominent ridges and other tectonic lineaments on the southern portion of the leading hemisphere of Europa is determined from cross-cutting relationships. These selected features formed fairly recently relative to most of the surface; older lineaments no longer retain clear evidence of cross-cutting. If we assume that this sequence represents the order of formation of cracks that underlie each lineament, and that the orientation of each crack was determined by tidal stress whose azimuth varies monotonically counter-clockwise with time, the azimuth must have rotated more than 740°, which would correspond to the change in tidal stress over two periods of nonsynchronous rotation (relative to the direction of Jupiter). However, that interpretation is not necessarily compelling, because the observed orientations of cross-cutting lineaments are not densely spaced over these cycles; in fact, the sequence would fit nearly as well into an arbitrary model with rotation in the opposite sense from that predicted by theory. This tectonic record may have formed over many more rotational cycles, such that typically only a few cracks form per cycle, which would be consistent with evidence from considerations of cycloidal crack patterns. Sets of cracks that cluster near certain azimuth orientations appear to be parts of globe-encircling lineament systems and may result from other effects, perhaps polar wander that occurred rapidly relative to nonsynchronous rotation.     

Spectral comparison of heavily hydrated salts with disrupted terrains on Europa

Hydrated magnesium sulfate salts have been proposed as major components of the disrupted, reddish terrains on the surface of Europa. This is based on near-infrared reflectance spectra which contain distorted and asymmetric water absorption features typical of moderately hydrated materials such as hexahydrite (MgSO₄⋅6H₂O) and epsomite (MgSO₄⋅7H₂O). Hydrated magnesium sulfates having many waters of hydration could produce improved spectral matches. Here we present cryogenic laboratory spectra of highly hydrated sulfur-bearing salts, including hexahydrite, epsomite, bloedite (Na₂Mg(SO₄)₂⋅4H₂O), mirabilite (Na₂SO₄⋅10H₂O), sodium sulfide nonahydrate (Na₂S⋅9H₂O), supersaturated MgSO₄, NaHCO₃, and Na₂SO₄ brines, and magnesium sulfate dodecahydrate (MgSO₄⋅12H₂O). All have been measured under conditions of pressure and temperature appropriate to the surface environment of Europa. Novel methods for preparation, verification and analysis of MgSO₄⋅12H₂O, which is not stable at standard temperature and pressure (STP), are described. At 100 K, all of these materials exhibit distorted and asymmetric absorption features similar to those in the Europa observations, as well as several weaker, narrow absorptions having widths ranging from 15 to 80 nm. While the agreement with Galileo NIMS observations of dark terrains on Europa is indeed better for highly hydrated salts than for salts of lower hydration states, we conclude that none of these materials alone can account for all of the observed spectral character. As previously suggested, Europa's reddish material appears to be a complex mixture of sulfate hydrates and other materials.     

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.     

Production of noble gases near the surface of Europa and the prospects for in situ chronology

The age of the surface of Europa is probably tens of Myr or less, but is poorly constrained. Two different geochronology schemes could potentially be applied to near-surface samples to provide far more precise ages. First, the surface salts apparently contain enough potassium to make potassium-argon dating feasible. Second, the bombardment of the surface with both galactic cosmic rays and protons trapped within the jovian magnetosphere will cause nuclear reactions that can lead to measurable buildups of cosmogenic noble gases, which can be used to determine both cosmic-ray exposure ages and erosion, deposition, or mixing rates for surface modification. The major differences between Europa's salt-rich ice and the rocks (meteorites, lunar samples and terrestrial rocks) in which cosmogenic noble gases are normally measured are that the abundance of target elements for nuclear reactions creating neon and argon are lower (because of the high water content), but neutron-capture reactions, particularly ³⁵Cl(n,γ)³⁶Cl(β⁻)³⁶Ar, are enhanced because of the thermalizing effects of the water. As well as presenting calculations of noble gas production near the surface of Europa, we also show that the measurements required to determine ages are within reach of technology in the near future, if an instrument can be landed on the surface of Europa.     

Origin and evolution of Castalia Macula, an anomalous young depression on Europa

Europa's Castalia Macula region was comprehensively imaged by the Galileo spacecraft on several orbits, at both local and regional resolutions and with different illumination geometries. Using these datasets we have mapped and identified the different geological units within the Castalia area, and derived digital elevation models (DEMs) of the topography within most of the Castalia Macula region. Using these data sets in combination allows us to map the geology and topography of this area in greater detail than perhaps any other site on Europa. Castalia Macula consists of unusually dark and reddish material, most of which is confined to a broad topographic depression 350 m deep located between two large uplifted domes 900 and 750 m high, to the north and south, respectively. The preservation of topography at the bottom of Castalia Macula indicates that dark material initially filled the depression to a certain depth but was subsequently removed via drainage, resulting in a dark stain up to the original equipotential surface. Superposition and topographic relationships suggest that the Castalia Macula plains deposit formed prior to uplift of both domes, and at least two distinct episodes of chaos formation have occurred near and on top of the northern dome. It appears that Castalia Macula is comparatively young and was active relatively recently, therefore it could provide an ideal place to sample material that has recently been erupted from the subsurface, and may have been in communication with Europa's ocean. These factors combine to make Castalia Macula a very attractive site for a future Europa lander.     

Neutral particle release from Europa’s surface

In this paper, we look at space weathering processes on the icy surface of Jupiter’s moon Europa. The heavy energetic ions of the jovian plasma (H⁺, O⁺, S⁺, C⁺) can erode the surface of Europa via ion sputtering (IS), ejecting up to 1000 H₂O molecules per ion. UV photons impinging the Europa’s surface can also result in neutral atom release via photon-stimulated desorption (PSD) and chemical change (photolysis). In this work, we study the efficiency of the IS and PSD processes for ejecting water molecules, simulating the resulting neutral H₂O density. We also estimate the contribution to the total neutral atom release by the Ion Backscattering (IBS) process. Moreover, we estimate the possibility of detecting the sputtered high energy atoms, in order to distinguish the action of the IS process from other surface release mechanisms. Our main results are: (1) The most significant sputtered-particle flux and the largest contribution to the neutral H₂O density come from the incident S⁺ ions; (2) the H₂O density produced via PSD is lower than that due to sputtering by ∼1.5 orders of magnitude; (3) in the energy range below 1 keV, the IBS can be considered negligible for the production of neutrals, whereas in the higher energy range it becomes the dominant neutral emission mechanism; (4) the total sputtering rate for Europa is 2.0 × 10²⁷ H₂O s⁻¹; and (5) the fraction of escaping H₂O via IS is 22% of the total sputtered population, while the escape fraction for H₂O produced by PSD is 30% of the total PSD population. Since the PSD exosphere is lower than the IS one, the major agent for Europa’s surface total net erosion is IS on both the non-illuminated and illuminated side. Lastly, the exospheric neutral density, estimated from the Galileo electron density measurements appears to be higher than that calculated for H₂O alone; this favors the scenario of the presence of O₂ produced by radiolysis and photolysis.     

Numerical modeling of endogenic thermal anomalies on Europa

A variety of recent resurfacing features have been observed on Europa, which may produce thermal anomalies detectable by a future mission. However, the likelihood of such a detection depends on their size and lifetimes. The results of this numerical study suggest that the lifetime of a thermal anomaly associated with the emplacement of 100 m of water onto the surface of Europa is several hundred years, and ∼10 years for 10 m of water. If warm ice is emplaced on the surface instead of liquid water, these lifetimes decrease by up to a factor of two. Exploration of model parameters indicates that a thin insulating surface layer can double thermal anomaly lifetimes, anomalies emplaced at a latitude of 80° can remain detectable nearly a factor of two longer than those at equatorial latitudes, and anomalies on the night side can remain detectable for up to ∼20% longer than those on the day side. High temperatures are very short-lived as the surface ice cools very rapidly to below 200 K due to sublimation cooling. Assuming steady-state resurfacing, the number of detectable thermal anomalies associated with the emplacement of 100 m of water would be on the order of 10 if the typical resurfacing area is 15 km². If recent resurfacing is dominated by chaos regions with typical areas of 100 to 1000 km² and lifetimes of 1000 to 4000 years, the number of detectable thermal anomalies would be on the order of 1 to 10.     

Distribution of hydrate on Europa: Further evidence for sulfuric acid hydrate

Sulfuric acid hydrate has been proposed as an important species on Europa's surface, the acid being produced by radiolysis of surficial sulfur compounds. We investigated the spectral properties of disordered and crystalline forms of sulfuric acid and suggest that the hydration properties of Europa's hypothesized sulfuric acid lie between two end members: liquid sulfuric acid and its higher crystalline hydrates. The spectra of these end members are similar except for spectral shifts at the band edges. We measured the optical constants of sulfuric acid octahydrate and used these with simple radiative transfer calculations to fit Europa spectra obtained by Galileo's Near Infrared Mapping Spectrometer (NIMS). The global distribution of the hydrate that we associate here with hydrated sulfuric acid shows a strong trailing-side enhancement with a maximum fractional hydrate abundance of 90% by volume, corresponding to a sulfur atom to water molecule ratio of 10%. The hydrate concentration spatially correlates with the ultraviolet and visible absorption of the surface and with the sulfur dioxide concentration. The asymmetric global distribution is consistent with Iogenic plasma ion implantation as the source of the sulfur, possibly modified by electron irradiation and sputtering effects. The variegated distribution also correlates with geologic forms. A high spatial resolution image shows resolved lineae with less hydrate appearing within the lineae than in nearby crustal material. The low concentration of hydrated material in these lineae argues against their conveying sulfurous material to the surface from the putative ocean.     

Model constraints on the opening rates of bands on Europa

A mid-ocean-ridge spreading analog is used to constrain the opening rates and brittle-ductile transition depths for two axisymmetric ridged bands on Europa. Estimates of brittle-ductile transition depth based on the morphologies of Yelland and Ino Lineae are combined with a conductive cooling model based on a mid-ocean ridge analog to estimate the opening rates and active lifetimes of the bands. This model limits local strain rates to ∼¹⁰⁻¹⁵-¹⁰⁻¹² s⁻¹, opening rates to 0.2-40 mm yr⁻¹, and active lifetimes of the bands to 0.1-30 Myr. If the observed structures in the outer portions of ridged bands are indeed normal faults, the estimated range for the tensile strength of ice on Europa is 0.4-2 MPa, consistent with nonsynchronous rotation as the dominant driving mechanism for band opening.     

A description of surface features in north Tyrrhena Terra, Mars: Evidence for extension and lava flooding

We studied north Tyrrhena Terra, an approximately 39,000 km² area, located in the transition region straddling the Amenthes and Mare Tyrrhenum Mars Chart quadrangles 14 and 22, respectively. The study area comprises ancient terrains with infilled craters, ridges and valleys. Interpretation of orbiter data of ancient terrains is inherently difficult, but valuable information can be obtained using multiple datasets and analyzing various geological features. Using data from the High Resolution Stereo Camera on board Mars Express, complemented by Mars Global Surveyor MOLA DEM and MOC Narrow Angle datasets, we observed and interpreted surface morphologies at a scale suitable for geologic investigation. Morphometric examination of a 31 km diameter large impact crater indicated that tectonism and volcanism were responsible for its morphologic modification. Small impact crater depth/diameter relationships indicated that smooth surfaces of valleys are composed of highly consolidated material. Surface cracks and lobate fronts further suggested that the rocks are volcanic. Examination of tectonic features revealed that in the study area: a dominant NW-SE fabric is related to a ridge/bench-scarp-valley repetition consistent with synthetic and antithetic normal faulting; a NNW-SSE lineament represents the surface expression of normal faulting post-dating all other tectonic features. A weak NE-SW fabric is observable as small sublinear depressions, and at the contact between units internal to one large crater. One 20 km diameter crater in the study area was interpreted to be a caldera, infilled by thick volcanic rock layers. Identification of wrinkle ridges further indicated that thick layered lava flows infilled the main depressions of the study area. The available evidence suggests that the study area underwent multiple episodes of extension and volcanism.     

Cycloid crack sequences on Europa: Relationship to stress history and constraints on growth mechanics based on cusp angles

The original model developed to explain cycloidal cracks on Europa interprets cycloids as tensile fractures that grow in a curved path in response to the constantly rotating diurnal tidal stress field. Cusps form when a new cycloid crack segment propagates at an angle to the first in response to a rotation of the principal tidal stress orientation during a period of no crack growth. A recent revised model states that a cycloid cusp forms through the creation of a secondary fracture called a tailcrack at the tip of an existing cycloid segment during shearing motion induced by the rotating tidal stress field. As the tailcrack propagates away from the cusp, it becomes the next cycloid segment in the chain. The qualitative tailcrack model uniquely accounts for the normal and shear stresses that mechanically must resolve onto the tip of an existing cycloid segment at the instant of cusp formation. In this work, we provide a quantitative framework and test of the hitherto purely conceptual tailcrack model. We first present a relative age sequence inferred from geologic mapping of multiply cross-cutting cycloids in Europa's trailing hemisphere and place this into the context of the global stress history. The age sequence requires a cumulative minimum of 630° of shell reorientation due to nonsynchronous rotation to account for the observed range of orientations of cycloids of different ages. We determined the back-rotated longitudes of formation of two cycloid chain examples and used mathematical modeling of europan tidal stresses to show that the tailcrack model for cusp formation is not only viable, but places constraints on the overall development of a cycloid chain by controlling the timing of cusp development within Europa's orbit. For all cusps analyzed, the exact ratio of resolved shear to normal stress required to form the cusp angles by a process of tailcracking, as governed by the principles of linear elastic fracture mechanics, is produced at the tip of a shearing cycloid segment during Europa's orbit. Cusp formation occurs after the point in the orbit at which the maximum tensile principal tidal stress occurs, implying that tensile tidal stresses are not directly responsible for cusp development. Instead, cusps develop when a tailcrack forms at the tip of a cycloid segment in response to the highly perturbed stress field induced during concomitant opening and shearing at the tip of the cycloid segment.     

Secondary and sesquinary craters on Europa

We address impact cratering on Io and Europa, with the emphasis on the origin of small craters on Europa as secondary to the primary impacts of comets on Io, Europa, and Ganymede. In passing we also address the origin of secondary craters generated by Zunil, a well-studied impact crater on Mars that is a plausible analog to impact craters on Io. At nominal impact rates, and taking volcanic resurfacing into account, we find that there should be 1.3 impact craters on Io, equally likely to be of any diameter between 100 m and 20 km. The corresponding model age of Europa's surface is between 60 and 100 Ma. This range of ages does not include a factor three uncertainty stemming from the uncertain sizes and numbers of comets. The mass of basaltic impact ejecta from Io to reach Europa is found to meet or exceed the micrometeoroid flux as a source of rock-forming elements to Europa's ice crust. To describe impact ejecta in more detail we adapt models for impact-generated spalls and Grady-Kipp fragments originally developed by Melosh. Our model successfully reproduces the observed size-number distributions of small craters on both Mars and Europa. However, the model predicts that a significant fraction of the 200-500 m diameter craters on Europa are not traditional secondary craters but are instead sesquinary craters caused by impact ejecta from Io that had gone into orbit about Jupiter. This prediction is not supported by observation, which implies that high speed spalls usually break up into smaller fragments that make smaller sesquinary craters. Iogenic basalts are also interesting because they provide stratigraphic horizons on Europa that in principle could be used to track historic motions of the ice, and they provide materials suitable to radiometric dating of Europa's surface.