共查询到20条相似文献，搜索用时 437 毫秒
1.
Javier Ruiz 《Earth, Moon, and Planets》2010,107(2-4):157-167
Water ice I rheology is a key factor for understanding the thermal and mechanical state of the outer shell of the icy satellites. Ice flow involves several deformation mechanisms (both Newtonian and non-Newtonian), which contribute to different extents depending on the temperature, grain size, and applied stress. In this work I analyze tidally heated and stressed equilibrium convection in the ice shell of Europa by considering a composite viscosity law which includes diffusion creep, basal slip, grain boundary sliding and dislocation creep, and. The calculations take into account the effect of tidal stresses on ice flow and use grain sizes between 0.1 and 100 mm. An Arrhenius-type relation (useful for parameterized convective models) is found then by fitting the calculated viscosity between 170 and 273 K to an exponential regression, which can be expressed in terms of pre-exponential constant and effective activation energy. I obtain convective heat flows between ~40 and ~60 mW m?2, values lower than those usually deduced (~100 mW m?2) from geological indicators of lithospheric thermal state, probably indicating heterogeneous tidal heating. On the other hand, for grain sizes larger than ~0.3 mm the thicknesses of the ice shell and convective sublayer are ~20–30 km and ~5–20 km respectively, values in good agreement with the available information for Europa. So, some fundamental geophysical characteristics of the ice shell of Europa could be arising from the properties of the composite water ice rheology. 相似文献
2.
Two opposing models to explain the geological features observed on Europa’s surface have been proposed. The thin-shell model states that the ice shell is only a few kilometers thick, transfers heat by conduction only, and can become locally thinner until it exposes an underlying ocean on the satellite’s surface. According to the thick-shell model, the ice shell may be several tens of kilometers thick and have a lower convective layer, above which there is a cold stagnant lid that dissipates heat by conduction. Whichever the case, from magnetic data there is strong support for the presence of a layer of salty liquid water under the ice. The present study was performed to examine whether the possibility of convection is theoretically consistent with surface heat flows of ∼100-200 mW m−2, deduced from a thin brittle lithosphere, and with the typical spacing of 15-23 km proposed for the features usually known as lenticulae. It was obtained that under Europa’s ice shell conditions convection could occur and also account for high heat flows due to tidal heating of the convective (nearly isothermal) interior, but only if the dominant water ice rheology is superplastic flow (with activation energy of 49 kJ mol−1; this is the rheology thought dominant in the warm interior of the ice shell). In this case the ice shell would be ∼15-50 km thick. Furthermore, in this scenario explaining the origin of the lenticulae related to convective processes requires ice grain size close to 1 mm and ice thickness around 15-20 km. 相似文献
3.
We present self-consistent, fully coupled two-dimensional (2D) numerical models of thermal evolution and tidal heating to investigate how convection interacts with tidal dissipation under the influence of non-Newtonian grain-size-sensitive creep rheology (plausibly resulting from grain boundary sliding) in Europa’s ice shell. To determine the thermal evolution, we solved the convection equations (using finite-element code ConMan) with the tidal dissipation as a heat source. For a given heterogeneous temperature field at a given time, we determined the tidal dissipation rate throughout the ice shell by solving for the tidal stresses and strains subject to Maxwell viscoelastic rheology (using finite-element code Tekton). In this way, the convection and tidal heating are fully coupled and evolve together. Our simulations show that the tidal dissipation rate can have a strong impact on the onset of thermal convection in Europa’s ice shell under non-Newtonian GSS rheology. By varying the ice grain size (1-10 mm), ice-shell thickness (20-120 km), and tidal-strain amplitude (0-4 × 10−5), we study the interrelationship of convection and conduction regimes in Europa’s ice shell. Under non-Newtonian grain-size-sensitive creep rheology and ice grain size larger than 1 mm, no thermal convection can initiate in Europa’s ice shell (for thicknesses <100 km) without tidal dissipation. However, thermal convection can start in thinner ice shells under the influence of tidal dissipation. The required tidal-strain amplitude for convection to occur decreases as the ice-shell thickness increases. For grain sizes of 1-10 mm, convection can occur in ice shells as thin as 20-40 km with the estimated tidal-strain amplitude of 2 × 10−5 on Europa. 相似文献
4.
Recent investigations into convection in Europa’s ice shell have been based on non-Newtonian (stress-dependent) or Newtonian
(stress-independent) viscosity for water ice. However, despite the wide use of Newtonian convection, experimentally observed
water ice flow is non-Newtonian, and analysis of stability against convection of the ice shell using updated flow laws has
been only performed for non-Newtonian rheologies. Here we use the flow law proposed for diffusion creep to analyze the possibility
of the onset of convection for Newtonian viscosity in relation to the thermal state of Europa. Our findings indicate that
for diffusion creep convection might have started, but that significantly lower heat flows (and equivalently higher shell
thicknesses) and/or grain sizes are required than for superplastic flow, which is the most probable flow mechanism if the
ice shell is convective. 相似文献
5.
We model stagnant-lid convection for water ice I using a multicomponent rheology, combining grain boundary sliding, dislocation and diffusion creep mechanisms. For the superplastic flow-dislocation creep rheology, dislocation creep (n=4) dominates the deformation within the actively convecting sublayer whilst superplastic flow (n=1.8) is the dominant process within the stagnant-lid whilst for the superplastic flow-diffusion creep rheology, superplastic flow is the dominant deformation mechanism within the convecting sublayer while diffusion creep (n=1) is the dominant deformation process in the stagnant-lid. These results suggest deformation in the actively convecting sublayer is likely to be dominated by the mechanism with the largest stress exponent. We also provide heat flux scaling relationships for the superplastic flow, basal slip, dislocation creep-superplastic flow and superplastic flow-diffusion creep rheologies and provide a simple parameterized convection model of an icy satellite thermal evolution. 相似文献
6.
William B. McKinnon 《Icarus》2006,183(2):435-450
It has been argued that the dominant non-Newtonian creep mechanisms of water ice make the ice shell above Callisto's ocean, and by inference all radiogenically heated ice I shells in the outer Solar System, stable against solid-state convective overturn. Conductive heat transport and internal melting (oceans) are therefore predicted to be, or have been, widespread among midsize and larger icy satellites and Kuiper Belt objects. Alternatively, at low stresses (where non-Newtonian viscosities can be arbitrarily large), convective instabilities may arise in the diffusional creep regime for arbitrarily small temperature perturbations. For Callisto, ice viscosities are low enough that convection is expected over most of geologic time above the internal liquid layer for plausible ice grain sizes (?a few mm); the alternative for early Callisto, a conducting shell over a very deep ocean (>450 km), is not compatible with Callisto's present partially differentiated state. Moreover, if convection is occurring today, the stagnant lid would be quite thick (∼100 km) and compatible with the lack of active geology. Nevertheless, Callisto's steady-state heat flow may have fallen below the convective minimum for its ice I shell late in Solar System history. In this case convection ends, the ice shell melts back at its base, and the internal ocean widens considerably. The presence of such an ocean, of order 200 km thick, is compatible with Callisto's moment-of-inertia, but its formation would have caused an ∼0.25% radial expansion. The tectonic effects of such a late, slow expansion are not observed, so convection likely persists in Callisto, possibly subcritically. Ganymede, due to its greater size, rock fraction and full differentiation, has a substantially higher heat flow than Callisto and has not reached this tectonic end state. Titan, if differentiated, and Triton should be more similar to Ganymede in this regard. Pluto, like Callisto, may be near the tipping point for convective shutdown, but uncertainties in its size and rock fraction prevent a more definitive assessment. 相似文献
7.
The history and dynamics of the martian polar deposits (MPD), the largest known water reservoirs on Mars, are of great interest, but estimates of ice grain size are required before detailed modeling can be performed. We clarify the microphysical processes that may control grain size in the MPD. If the MPD are ∼2% dust by mass, the maximum ice grain size is ∼1 mm due to grain boundary pinning by silicate microparticles. Relatively dusty layers in the MPD will have smaller grain sizes. If MPD ice has a very low impurity content and has experienced a significant amount of strain, grains may reach a steady state size of ∼1.5 to 3 mm due to dynamic recrystallization, wherein a steady state grain size is maintained due to the balance of grain growth and destruction during flow. If the near-bed ice in the MPD is warmed close to its melting point and has been extensively sheared, grain sizes at its base may be between 10 and 40 mm, by analogy with warm, dirty, near-bed ice in terrestrial ice sheets. 相似文献
8.
William B. Moore 《Icarus》2006,180(1):141-146
Models of tidal-convective equilibrium for Europa's ice shell are computed using a laboratory-derived composite flow law for ice. Volume diffusion creep is found to dominate the flow law at equilibrium, and thus the thickness of the shell is strongly dependent on the poorly known grain size of the ice. This grain size is, however, constrained to be less than a few millimeters if equilibrium is achieved at the current eccentricity. Europa's ice shell cannot be thinner than about 16 km in equilibrium at present, since tidal dissipation cannot generate enough heat in such a thin shell to balance the heat transport. No conductive equilibria are found; this is likely due to the fact that most of a conductive shell must be cold if temperature gradients are to be large enough to carry the heat. A minimum eccentricity of about 0.0025 (about 1/4 the present value) below which there are no equilibria is also found. 相似文献
9.
Radar observations in the Deuteronilus Mensae region by Mars Reconnaissance Orbiter have constrained the thickness and dust concentration found within mid-latitude ice deposits, providing an opportunity to more accurately estimate the rheology of ice responsible for the formation of lobate debris aprons based on their apparent age of ∼100 Myr. We developed a numerical model simulating ice flow under martian conditions using results from ice deformation experiments, theory of ice grain growth based on terrestrial ice cores, and observational constraints from radar profiles and laser altimetry. By varying the ice grain size, the ice temperature, the subsurface slope, and the initial ice volume we determine the combination of parameters that best reproduce the observed LDA lengths and thicknesses over a period of time comparable to the apparent ages of LDA surfaces (90-300 Myr). We find that an ice temperature of 205 K, an ice grain size of 5 mm, and a flat subsurface slope give reasonable ages for many LDAs in the northern mid-latitudes of Mars. Assuming that the ice grain size is limited by the grain boundary pinning effect of incorporated dust, these results limit the dust volume concentration to less than 4%. However, assuming all LDAs were emplaced by a single event, we find that there is no single combination of grain size, temperature, and subsurface slope which can give realistic ages for all LDAs, suggesting that some or all of these variables are spatially heterogeneous. Based on our model we conclude that the majority of northern mid-latitude LDAs are composed of clean (?4 vol%), coarse (?1 mm) grained ice, but regional differences in either the amount of dust mixed in with the ice, or in the presence of a basal slope below the LDA ice must be invoked. Alternatively, the ice temperature and/or timing of ice deposition may vary significantly between different mid-latitude regions. Either eventuality can be tested with future observations. 相似文献
10.
Coupled thermal-orbital evolution models of Europa and Io are presented. It is assumed that Io, Europa, and Ganymede evolve in the Laplace resonance and that tidal dissipation of orbital energy is an internal heat source for both Io and Europa. While dissipation in Io occurs in the mantle as in the mantle dissipation model of Segatz et al. (1988, Icarus 75, 187), two models for Europa are considered. In the first model dissipation occurs in the silicate mantle while in the second model dissipation occurs in the ice shell. In the latter model, ice shell melting and variations of the shell thickness above an ocean are explicitly included. The rheology of both the ice and the rock is cast in terms of a viscoelastic Maxwell rheology with viscosity and shear modulus depending on the average temperature of the dissipating layer. Heat transfer by convection is calculated using a parameterization for strongly temperature-dependent viscosity convection. Both models are consistent with the present orbital elements of Io, Europa, and Ganymede. It is shown that there may be phases of quasi-steady evolution with large or small dissipation rates (in comparison with radiogenic heating), phases with runaway heating or cooling and oscillatory phases during which the eccentricity and the tidal heating rate will oscillate. Europa's ice thickness varies between roughly 3 and 70 km (dissipation in the silicate layer) or 10 and 60 km (dissipation in the ice layer), suggesting that Europa's ocean existed for geological timescales. The variation in ice thickness, including both convective and purely conductive phases, may be reflected in the formation of different geological surface features on Europa. Both models suggest that at present Europa's ice thickness is several tens of km thick and is increasing, while the eccentricity decreases, implying that the satellites evolve out of resonance. Including lithospheric growth in the models makes it impossible to match the high heat flux constraint for Io. Other heat transfer processes than conduction through the lithosphere must be important for the present Io. 相似文献
11.
The tectonically and cryovolcanically resurfaced terrains of Ganymede attest to the satellite's turbulent geologic history. Yet, the ultimate cause of its geologic violence remains unknown. One plausible scenario suggests that the Galilean satellites passed through one or more Laplace-like resonances before evolving into the current Laplace resonance. Passage through such a resonance can excite Ganymede's eccentricity, leading to tidal dissipation within the ice shell. To evaluate the effects of resonance passage on Ganymede's thermal history we model the coupled orbital-thermal evolution of Ganymede both with and without passage through a Laplace-like resonance. In the absence of tidal dissipation, radiogenic heating alone is capable of creating large internal oceans within Ganymede if the ice grain size is 1 mm or greater. For larger grain sizes, oceans will exist into the present epoch. The inclusion of tidal dissipation significantly alters Ganymede's thermal history, and for some parameters (e.g. ice grain size, tidal Q of Jupiter) a thin ice shell (5 to 20 km) can be maintained throughout the period of resonance passage. The pulse of tidal heating that accompanies Laplace-like resonance capture can cause up to 2.5% volumetric expansion of the satellite and contemporaneous formation of near surface partial melt. The presence of a thin ice shell and high satellite orbital eccentricity would generate moderate diurnal tidal stresses in Ganymede's ice shell. Larger stresses result if the ice shell rotates non-synchronously. The combined effects of satellite expansion, its associated tensile stress, rapid formation of near surface partial melt, and tidal stress due to an eccentric orbit may be responsible for creating Ganymede's unique surface features. 相似文献
12.
Javier Ruiz 《Icarus》2005,177(2):438-446
The heat flow from Europa has profound implications for ice shell thickness and structure, as well as for the existence of an internal ocean, which is strongly suggested by magnetic data. The brittle-ductile transition depth and the effective elastic thickness of the lithosphere are here used to perform heat flow estimations for Europa. Results give preferred heat flow values (for a typical geological strain rate of 10−15 s−1) of 70-110 mW m−2 for a brittle-ductile transition 2 km deep (the usually accepted upper limit for the brittle-ductile transition depth in the ice shell of Europa), 24-35 mW m−2 for an effective elastic thickness of 2.9 km supporting a plateau near the Cilix impact crater, and >130 mW m−2 for effective elastic thicknesses of ?0.4 km proposed for the lithosphere loaded by ridges and domes. These values are clearly higher than those produced by radiogenic heating, thus implying an important role for tidal heating. The ?19-25 km thick ice shell proposed from the analysis of size and depth of impact structures suggests a heat flow of ?30-45 mW m−2 reaching the ice shell base, which in turn would imply an important contribution to the heat flow from tidal heating within the ice shell. Tidally heated convection in the ice shell could be capable to supply ∼100 mW m−2 for superplastic flow, and, at the Cilix crater region, ∼35-50 mW m−2 for dislocation creep, which suggests local variations in the dominant flow mechanism for convection. The very high heat flows maybe related to ridges and domes could be originated by preferential heating at special settings. 相似文献
13.
The proposed past eruption of liquid water on Europa and ongoing eruption of water vapor and ice on Enceladus have led to discussion about the feasibility of cracking a planetary ice shell. We use a boundary element method to model crack penetration in an ice shell subjected to tension and hydrostatic compression. We consider the presence of a region at the base of the ice shell in which the far-field extensional stresses vanish due to viscoelastic relaxation, impeding the penetration of fractures towards a subsurface ocean. The maximum extent of fracture penetration can be limited by hydrostatic pressure or by the presence of the unstressed basal layer, depending on its thickness. Our results indicate that Europa's ice shell is likely to be cracked under 1-3 MPa tension only if it is ?2.5 km thick. Enceladus' ice shell may be completely cracked if it is capable of supporting ∼1-3 MPa tension and is less than 25 km thick. 相似文献
14.
G.B. Hansen 《Icarus》2009,203(2):672-676
The grain size of water ice can be determined from its near-infrared spectrum, which has numerous diagnostic absorption bands of different opacities. Models that have been used to determine water ice grain size from infrared spectra of icy outer Solar System objects have shown discrepancies in modeled grain size of a factor of two or more. Here the single-scattering albedo calculated using the commonly used Hapke model given by Roush [Roush, T.L., 1994. Icarus 108, 243-254] is compared with the exact calculation for spheres from a Mie series. An earlier approximation of single-scattering albedo called the Hapke “slab” model is also used in the comparison. All three models are implemented using the same optical constants for water ice at ∼110 K. Results are displayed for a large range of grain sizes from 1 μm to 1 mm. In general neither Hapke model can mimic the Rayleigh effects from particles sized near the wavelength of light that the Mie model predicts. For 10 μm particles, the slab model matches the Mie calculation quite well, but larger sizes are more discrepant. The Hapke/Roush model grain size needs to be ∼2.5 times larger to mimic the Mie results, and there are additional discrepancies in the continuum levels and band strengths. The Mie calculation for spheres is recommended for analysis of unknown remote sensing measurements, as it can mimic the spectra of oblate, prolate, and hollow particles given by equivalent sphere theories. 相似文献
15.
Ice-shell thickness and ocean depth are calculated for steady state models of tidal dissipation in Europa's ice shell using the present-day values of the orbital elements. The tidal dissipation rate is obtained using a viscoelastic Maxwell rheology for the ice, the viscosity of which has been varied over a wide range, and is found to strongly increase if an (inviscid) internal ocean is present. To determine steady state values, the tidal dissipation rate is equated to the heat-transfer rate through the ice shell calculated from a parameterized model of convective heat transfer or from a thermal conduction model, if the ice layer is found to be stable against convection. Although high dissipation rates and heat fluxes of up to 300 mWm−2 are, in principle, possible for Europa, these values are unrealistic because the states for which they are obtained are thermodynamically unstable. Equilibrium models have surface heat flows around 20 mWm−2 and ice-layer thicknesses around 30 km, which is significantly less than the total thickness of the H2O-layer. These results support models of Europa with ice shells a few tens of kilometers thick and around 100-km-thick subsurface oceans. 相似文献
16.
Data from the recent gravity measurements by the Galileo mission are used to construct wide ranges of interior structure and composition models for the Galilean satellites of Jupiter. These models show that mantle densities of Io and Europa are consistent with an olivine-dominated mineralogy with the ratios of Mg to Fe components depending on mantle temperature for Io and on ice shell thickness for Europa. The mantle density and composition depend relatively little on core composition. The size of the core is largely determined by the core's composition with core radius increasing with the concentration of a light component such as sulfur. For Io, the range of possible core sizes is between 38 and 53% of the satellite's radius. For Europa, there is also a substantial effect of the thickness of the ice layer which is varied between 120 and 170 km on the core size. Core sizes are between 10 and 45% of Europa's radius. The core size of Ganymede ranges between one-quarter and one-third of the surface radius depending on its sulfur content and the thickness of the ice shell. A subset of the Ganymede models is consistent with an olivine-dominated mantle mineralogy. The thickness of the silicate mantle above the core varies between 900 and 1100 km. The outermost ice shell is about 900 km in thickness and is further subdivided by pressure-induced phase transitions into ice I, ice III, ice V, and ice VI layers. Callisto should be differentiated, albeit incompletely. It is proposed that this satellite was never molten at a large scale but differentiated through the convective gradual unmixing of the ice and the metal/rock component. Bulk iron-to-silicon ratios Fe/Si calculated for the inner pair of satellites, Io and Europa, are less than the CI carbonaceous chondrite value of 1.7±0.1, whereas ratios for the outer pair, Ganymede and Callisto, cover a broad range above the chondritic value. Although the ratios are uncertain, in particular for Ganymede and Callisto, the values are sufficiently distinct to suggest a difference in composition between these two pairs of satellites. This may indicate a difference in iron-silicon fractionation during the formation of both classes of satellites in the protojovian nebula. 相似文献
17.
Previous analyses into flexural deformation on the icy satellites of Jupiter and Saturn have assumed static, elastic lithospheres. Viscous creep within the lithosphere, however, can cause evolution over time. Here, we apply a finite-element model that employs a time-dependent elastic–viscous-plastic rheology in order to investigate flexure on icy satellites. Factors that affect this time-dependent response are those that control creep rates; surface temperature, heat flow, and grain size. Our results show that surface temperature is by far the dominant factor. At higher surface temperatures (100–130 K), the evolution of the deformation is such that the thickness of a modeled elastic lithosphere could vary by up to an order of magnitude, depending on the time scale over which the deformation occurred. Because the flexure observed on icy satellites generally indicates transient high heat flow events, our results indicate that the duration of the heat pulse is an important factor. For the icy worlds of Jupiter and Saturn, static models of lithospheric flexure should be used with caution. 相似文献
18.
A whole-moon numerical model of Europa is developed to simulate its thermal history. The thermal evolution covers three phases: (i) an initial, roughly 0.5 Gyr-long period of radiogenic heating and differentiation, (ii) a long period from 0.5 Gyr to 4 Gyr with continuing radiogenic heating but no tidal dissipative heating (TDH), and (iii) a final period covering the last 0.5 Gyr until the present, during which TDH is active. Hydrothermal plumes develop after the initial period of heating and differentiation and transport heat and salt from Europa’s silicate mantle to its ice shell. We find that, even without TDH, vigorous hydrothermal convection in the rocky mantle can sustain flow in an ocean layer throughout Europa’s history. When TDH becomes active, the ice shell melts quickly to a thickness of about 20 km, leaving an ocean 80 km or more deep. Parameterized convection in the ice shell is non-uniform spatially, changes over time, and is tied to the deeper ocean–mantle dynamics. We also find that the dynamics are affected by salt concentrations. An initially non-uniform salt distribution retards plume penetration, but is homogenized over time by turbulent diffusion and time-dependent flow driven by initial thermal gradients. After homogenization, the uniformly distributed salt concentrations are no longer a major factor in controlling plume transport. Salt transport leads to the formation of a heterogeneous brine layer and salt inclusions at the bottom of the ice shell; the presence of salt in the ice shell could strongly influence convection in that layer. 相似文献
19.
Internal layers in ice masses can be detected with ice-penetrating radar. In a flowing ice mass, each horizon represents a past surface that has been subsequently buried by accumulation, and strained by ice flow. These layers retain information about relative spatial patterns of accumulation and ablation (mass balance). Internal layers are necessary to accurately infer mass-balance patterns because the ice-surface shape only weakly reflects spatial variations in mass balance. Additional rate-controlling information, such as the layer age, the ice temperature, or the ice-grain sizes and ice-crystal fabric, can be used to infer the absolute rate of mass balance. To infer mass balance from the shapes of internal layers, we solve an inverse problem. The solution to the inverse problem is the best set or sets of unknown boundary conditions or initial conditions that, when used in our calculation of ice-surface elevation and internal-layer shape, generate appropriate predictions of observations that are available. We also show that internal layers can be used to infer martian paleo-surface topography from a past era of ice flow, even though the topography may have been largely altered by subsequent erosion. We have successfully inferred accumulation rates and surface topography from internal layers in Antarctica. Using synthetic data, we demonstrate the ability of this method to solve the corresponding inverse problem to infer accumulation and ablation rates, as well as the surface topography, for martian ice. If past ice flow has affected the shapes of martian internal layers, this method is necessary to infer the spatial pattern and rate of mass balance. 相似文献
20.
The recently measured dimensionless moment of inertia (MoI) factor for Callisto of 0.3549±0.0042 (Anderson et al., 2001, Icarus, 153, 157-161) poses a problem: its value cannot be explained by a model in which Callisto is completely differentiated into an ice shell above a rock shell and an iron core such as its neighboring satellite Ganymede nor can it be explained by a model of a homogeneous, undifferentiated ice-rock satellite. We show that Callisto may be incompletely differentiated into an outer ice-rock shell in which the volumetric rock concentration is close to the primordial one at the surface and decreases approximately linearly with depth, an ice mantle mostly depleted of rock, and an about 1800 km rock-ice core in which the rock concentration is close to the close-packing limit. The ice-rock shell thickness depends on uncertain rheology parameters and the heat flow and can be roughly 50 to 150 km thick. We show that if Callisto accreted from a mix of metal bearing rock and ice and if the average size of the rocks was of the order of meters to tens of meters, then Callisto may have experienced a gradual, but still incomplete unmixing of the two components. An ocean in Callisto at a depth of 100-200 km is difficult to obtain if the ice is pure H2O and if the ice-rock lithosphere is 100 km or more thick; a water ocean is more plausible for ice contaminated by ammonia, methane or salts; or for pure H2O at a depth of 400-600 km. 相似文献