Gravity,topography, and crustal evolution of Venus

The gravity anomalies of Venus, although small by comparison with those on Mars and the Moon, are still much larger than those on Earth for large features. On Venus, even the low-degree spherical harmonic terms for Venus' gravity field indicate a close association of broad positive gravity anomalies with major topographic highs. This is striking contrast to the situation on Earth, where the broad regional gravity anomalies show little correlation with continental masses or plate tectonic features, but instead appear to be caused by deep mass anomalies.A method for estimating radial gravity anomalies from line-of-sight acceleration data, their interpolation, and use of iteration for improved radial anomaly estimates is outlined. A preliminary gravity anomaly map of Venus at spacecraft altitude prepared using first estimate values is presented. A profile across the western part of Aphrodite along longitude 85 E was analyzed using time-series techniques. An elastic plate model would require a plate thickness of about 180 to 200 km to match the general amplitude of the observed gravity anomaly (about 33 mgal): a thickness much greater than that found for earth structures and, because of high surface temperatures, unlikely for Venus. An Airy isostatic model convolved with the topography across Aphrodite, however, provides a better match between the predicted and observed gravity anomalies if the nominal crustal thickness is about 70 to 80 km. This thickness is over twice that for continental crust on the earth, and considerably greater than that of the earth's basaltic ocean crust (only 5 km). A different differentiation history for Venus than that of the earth thus is anticipated. High gravity anomalies (+110 mgal) occur over Beta Regio and over the topographic high in eastern Aphrodite; both highs are associated with regions where detected lightning is clustered, and thus the topographic features may be active volcanic constructs. The large gravity anomalies at these two sites of volcanic activity require an explanation different than that indicated for western Aphrodite.     

Gravity studies of Aphrodite Terra on Venus

Aphrodite Terra is the largest highland area on Venus of the size of Africa. It is traversed by the Aphrodite-Beta belt of troughs with a length of 21 000 km. There are two other large belts of troughs on Venus: Themis-Atla, 14 000 km long, and Beta-Phoebe, 8000 km long. In this paper, four gravity profiles across Aphrodite Terra are studied and compared with the morphology.Western Aphrodite and Niobe Planitia to the north seem to be in isostatic equilibrium under the assumption of Airy compensation with a mean crustal thickness of T = 50 km. The graben area in the middle part of Aphrodite Terra shows negative isostatic gravity anomalies indicating deficit masses. The adjacent Atla Regio to the east is regionally Airy compensated with T = 50 km, and the mountains Nokomis, Maat and Ozza Montes are locally undercompensated, i.e. they are associated with surplus masses in the depth. Ulfrun Regio, a hilly terrain just east of Atla Regio is Airy compensated with T = 30 km. These results give a mean crustal thickness around 50 km for Aphrodite Terra. The isostatic disturbed zones in the middle of Aphrodite (grabens) and Atla Regio as well as the undercompensated Beta Regio have been associated with recent volcanism from the observation of the concentrations of electrical discharges in these areas. Atla and Beta Regiones are both located at intersections of the systems of troughs described above.Contribution No. 308, Institut für Geophysik der Universität Kiel, F.R.G.     

Processes of crustal formation and evolution on Venus: An analysis of topography,hypsometry, and crustal thickness variations

On Venus, present evidence indicates a crust of predominantly basaltic composition and a relatively young average age for the surface (several hundreds of millions of years). Estimates of crustal thickness from several approaches suggest an average crustal thickness of 10–20 km for much of the lowlands and rolling plains and a total volume of crust of about 1 × 10¹⁰ km³, approximately comparable to the present crustal volume of the Earth (1.02 × 10¹⁰ km³). The Earth's oceanic crust is thought to have been recycled at least 10–20 times over Earth history. The near-coincidence in present crustal volumes for the Earth and Venus suggests that either: (1) the presently observed crust of Venus represents the total volume that has accumulated over the history of the planet and that crustal production rates are thus very low, or (2) that crustal production rates are higher and that there is a large volume of missing crust unaccounted for on Venus which may have been lost by processes of crustal recycling.Known processes of crustal formation and thickening (impact-related magma ocean, vertical differentiation, and crustal spreading) are reviewed and are used as a guide to assess regional geologic evidence for the importance of these processes on Venus. Geologic evidence for variations in crustal thickness on Venus (range and frequency distribution of topography, regional slopes, etc.) are outlined. The hypothesis that the topography of Venus could result solely from crustal thickness variations is assessed and tested as an end-member hypothesis. A map of crustal thickness distribution is compiled on the basis of a simple model of Airy isostasy and global Venus topography. An assessment is then made of the significance of crustal thickness variations in explaining the topography of Venus. It is found that the distinctive unimodal hypsometric curve could be explained by: (1) a crust of relatively uniform thickness (most likely 10–20 km thick) comprising over 75% of the surface, (2) local plateaus (tessera) of thickened crust (about 20–30 km) forming less than 15% of the surface, (3) regions of apparent crustal thicknesses of 30–50 km (Beta, Ovda, Thetis, Atla Regiones and Western Ishtar Terra) forming less than 10% of the surface and showing some geologic evidence of crustal thickening processes (these areas can be explained on the basis of geologic observations and gravity data as combinations of thermal effects and crustal thickening), and (4) areas in which Airy isostasy predicts crustal thicknesses in excess of 50 km (the linear orogenic belts of Western Ishtar Terra, less than 1% of the surface).It is concluded that Venus hypsometry can be reasonably explained by a global crust of generally similar thickness with variations in topography being related to (1) crustal thickening processes (orogenic belts and plateau formation) and (2) local variations in the thermal structure (spatially varying thermal expansion in response to spatially varying heat flow). The most likely candidates for the formation and evolution of the crust are vertical differentiation and/or lateral crustal spreading processes. The small average crustal thickness (10–20 km) and the relatively small present crustal volume suggest that if vertical crustal growth processes are the dominant mechanism of crustal growth, than vertical growth has not commonly proceeded to the point where recycling by basal melting or density inversion will occur, and that therefore, rates of crustal production must have been much lower in the past than in recent history. Crustal spreading processes provide a mechanism for crustal formation and evolution that is consistent with observed crustal thicknesses. Crustal spreading processes would be characterized by higher (perhaps more Earth-like) crustal production rates than would characterize vertical differentiation processes, and crust created earlier in the history of Venus and not now observed (missing crust) would be accounted for by loss of crust through recycling processes. Lateral crustal spreading processes for the formation and evolution of the crust of Venus are interpreted to be consistent with many of the observations derived from presently available data. Resurfacing through vertical differentiation processes also clearly occurs, and if it is the major contributor to the total volume of the crust, then very low resurfacing rates are required.Although thermal effects on topography are clearly present and important on both Venus and the Earth, the major difference between the hypsometric curves on Earth (bimodal) and Venus (unimodal) is attributed primarily to the contrast in relative average thickness of the crust between the two terrains on Earth (continental/oceanic; 40/5 km = 35 km, 8:1) and Venus (upland plateaus/lowlands; about 30/15 km = 15 km, 2:1) (35 – 20 km = a difference of 20 km). The Venus unimodal distribution is thus attributed primarily to the large percentage of terrain with relatively uniform crustal thickness, with the skewness toward higher elevations due to the relatively small percentage of crust that is thickened by only about a factor of two. The Earth, in contrast, has a larger percentage of highlands (continents), whose crust is thicker by a factor of eight, on the average, leading to the distinctive bimodal hypsometric curve.Data necessary to firmly establish the dominant type of crustal formation and thickening processes operating and to determine the exact proportion of the topography of Venus that is due to thermal effects versus crustal thickness variations include: (1) global imaging data (to determine the age of the surface, the distribution and age of regions of high heat flux, and evidence for the nature and global distribution of processes of crustal formation and crustal loss), and (2) high resolution global gravity and topography data (to model crustal thickness variations and thermal contributions and to test various hypotheses of crustal growth).'Geology and Tectonics of Venus', special issue edited by Alexander T. Basilevsky (USSR Acad. of Sci. Moscow), James W. Head (Brown University, Providence), Gordon H. Pettengill (MIT, Cambridge, Massachusetts) and R. S. Saunders (J.P.L., Pasadena).     

Deep versus shallow origin of gravity anomalies, topography and volcanism on Earth, Venus and Mars

The relation between gravity anomalies, topography and volcanism can yield important insights about the internal dynamics of planets. From the power spectra of gravity and topography on Earth, Venus and Mars we infer that gravity anomalies have likely predominantly sources below the lithosphere up to about spherical harmonic degree l=30 for Earth, 40 for Venus and 5 for Mars. To interpret the low-degree part of the gravity spectrum in terms of possible sublithospheric density anomalies we derive radial mantle viscosity profiles consistent with mineral physics. For these viscosity profiles we then compute gravity and topography kernels, which indicate how much gravity anomaly and how much topography is caused by a density anomaly at a given depth. With these kernels, we firstly compute an expected gravity-topography ratio. Good agreement with the observed ratio indicates that for Venus, in contrast to Earth and Mars, long-wavelength topography is largely dynamically supported from the sublithospheric mantle. Secondly, we combine an empirical power spectrum of density anomalies inferred from seismic tomography in Earth’s mantle with gravity kernels to model the gravity power spectrum. We find a good match between modeled and observed gravity power spectrum for all three planets, except for 2?l?4 on Venus. Density anomalies in the Venusian mantle for these low degrees thus appear to be very small. We combine gravity kernels and the gravity field to derive radially averaged density anomaly models for the Martian and Venusian mantles. Gravity kernels for l?5 are very small on Venus below ≈800 km depth. Thus our inferences on Venusian mantle density are basically restricted to the upper 800 km. On Mars, gravity anomalies for 2?l?5 may originate from density anomalies anywhere within its mantle. For Mars as for Earth, inferred density anomalies are dominated by l=2 structure, but we cannot infer whether there are features in the lowermost mantle of Mars that correspond to Earth’s Large Low Shear Velocity Provinces (LLSVPs). We find that volcanism on Mars tends to occur primarily in regions above inferred low mantle density, but our model cannot distinguish whether or not there is a Martian analog for the finding that Earth’s Large Igneous Provinces mainly originate above the margins of LLSVPs.     

The stress state of Venusian crust and variations of its thickness: Implication for tectonics and geodynamics

The data obtained for the heights of the relief and the external gravitational field of Venus for spherical harmonics with degree and order up to 18 allow one to start theoretical analysis of the crust-mantle boundary (Venusian Moho) and stress state of the planetary interior. We suppose that Venusian convection is confined by floating massive crust. Apparently the convection in the upper mantle of Venus is separated from that one in the lower mantle and its lateral scale must be essentially smaller than on Earth. So, the convection is reflected to a larger degree of the gravitational field of the planet than for Earth. The spherical harmonic expansion of the topography for Venus correlates with corresponding expansion of the non-equilibrium part of the gravitational potential for n = 3–18. At the same time the relief of Venus is significantly compensated. It is reasonable to suppose that the gravity field for these harmonics is due to crustal thickness variations and, probably, to variations of crustal density. Thus, in the proposed scheme the Moho's relief causes the partial isostatic compensation of the topography.All calculations are carried out for the series of realistic models of Venus taking into consideration an asthenosphere. The asthenosphere is modeled either by a weakened (shear modulus is reduced), or by a liquid inviscid layer. We also suppose that the asthenosphere extends from the base of crust to a depth of 418 km, and the density contrast across the Moho boundary is –0.4 g * cm^–3. If the actual density contrast across the Moho is less than the supposed one by some factor, then one must increase the amplitudes of the roots and inverse roots by the same factor. The results for the Moho's relief and stresses in the crust are presented for the case of the mean thickness of the crust of 50 km, which satisfies the probable upper (connected with phase transitions in waterless basalts) and lower (appearing in the framework of our interpretation) limits.On the whole, the crust-mantle boundary on Venus is evidently smooth, and the stress level in the crust is appreciably smaller than the crustal stresses on the Earth. The strong sensitivity of the stresses character to the parameters of the model of external layers of Venus together with geological data allow us to begin a preliminary investigation of the tectonical structure and geodynamics of the planet.'Geology and Tectonics of Venus', special issue edited by Alexander T. Basilevsky (USSR Acad. of Sci. Moscow), James W. Head (Brown University, Providence), Gordon H. Pettengill (MIT, Cambridge, Massachusetts) and R. S. Saunders (J.P.L., Pasadena).     

Ridge belt tectonics of Ganiki Planitia on Venus

The main major ridge belts of Ganiki Planitia on Venus (Lama, Ahsonnutli and Pandrosos Dorsa) are part of the fan-shaped ridge belt complex along the 200 parallel of longitude. These ridge belts with evidence of crustal shortening support the idea of a large-scale E-W compression. The ridge belt patterns indicate a N-S shear component. These forces are explained by a triangular planitia area which compressed by surrounding terrains. The crustal shortening and ridge belt formation indicates compressional plate movement stresses in the uppermost lithosphere.Three sizes of ridge belt structure are to be found within Ganiki Planitia. (1) The ridge belt spacing of 200–400 km can be used to estimate the depth of the major uppermost homogeneous layer of Venus. There are numerous volcanic coronae, paterae and montes located along the main ridge belts or at their junctions. (2) Mid-size ridge groups or subbelts are to be found within the major ridge belts. These are formed by more local responses to tectonic stresses in the stratified uppermost crust. A wavelength of 40–70 km can be seen as a result of bending of the crustal strata and may relate to its thickness. (3) Small individual ridges are connected with most local stresses, defining places where the surface layers broke along the crests of large ridge belts or mid-scale subbelts. Radial and concentric mare ridge-like structures around coronae indicate that corona formation was effective at a sufficiently close vicinity to fault the surface.     

Divergent plate boundary characteristics and crustal spreading in Aphrodite Terra,Venus: A test of some predictions

It has been proposed that divergence and crustal spreading occur in Western Aphrodite Terra and some adjacent equatorial regions of Venus at rates in the range of a few centimeters per year. If equatorial spreading is common and widespread, then a consequence of this should be: (1) a young average age of the surface of the planet, (2) a trend in age from older terrain in the polar regions to younger terrain toward the equator, and (3) a latitudinal distribution of extensional features in equatorial regions and compressional deformation features in middle to high latitudes. These predictions are tested using published results from Arecibo, VENERA 15/16, and Pioneer Venus data, and it is found that: (1) the northern mid-to-high latitudes are characterized by a young average age, (2) there is a trend in the total number of craters per unit area from high values in the north polar regions to low values toward the equator, and (3) there is evidence for a latitudinal distribution of tectonic features of different types, with extensional features common in equatorial regions and compressional deformation features common in the northern middle to high latitudes. Further tests of these and other predictions can be made using data from the upcoming Magellan mission.     

Tectonic Evolution of Eastern Ishtar Terra,Venus

The complex morphology and topography of Eastern Ishtar Terra have been interpreted as due to tectonic deformation. Models proposed to account for this deformation include: crustal flow through asthenospheric flow and thermal-gravitational sliding; rifting, gravity spreading, and fold belt formation; and horizontal convergence and crustal thickening. In this study we map the detailed structural and topographic fabric of this region in order to explore and test these hypotheses. Eastern Ishtar can be divided into four major provinces: Maxwell Montes/Western Fortuna Tessera, a high plateau and mountain belt dominated by long NNW trending ridges; Central Fortuna Tessera, a low region of orthogonally oriented short WNW trending ridges and long, NNE trending troughs; Eastern Fortuna Tessera, a broad, E-W trending topographic rise characterized by ENE trending troughs and a complex pattern of intersecting ridges; and Northern Fortuna Tessera, a region of steep, NE-facing topographic scarps and ridges that trend WNW. On the basis of structural and topographic relationships, the features within these provinces are found to be inconsistent with a formation through either downslope crustal flow or rifting. We find that the mapped features are most consistent with a formation through convergence, collison, and underthrusting of thickened crustal terranes. These terranes are suggested to have been created through processes of seafloor-type spreading and crustal collision. Based on relationships between the different terranes, several accretional events are proposed in which Eastern Ishtar is produced by the collision of crustal terranes beginning at Lakshmi Planum and extending to the east. This sequence is initiated with the formation of Maxwell Montes and Western Fortuna Tessera during east-west crustal convergence, underthrusting, and stacking. The next step involves the northeast to southwest convergence of a preexisting thick block of tessera in Central Fortuna, which produces shear deformation within Western Fortuna. This northeast to southwest convergence also produces Northern Fortuna Tessera through crustal imbrication, a process recognized along the entire northern boundary of Ishtar Terra. Finally, Laima Tessera converges with Fortuna from the southeast and collides with Eastern Fortuna Tessera producing shear within Eastern Fortuna and the linear convergence zones along the edges of Laima. High resolution images returned by the Magellan spacecraft will enable us to examine the features involved in the proposed production and suturing of crustal terranes.'Geology and Tectonics of Venus', special issue edited by Alexander T. Basilevsky (USSR Acad. of Sci., Moscow), James W. Head (Brown University, Providence), Gordon H. Pettengill (MIT, Cambridge, Massachusetts) and R. S. Saunders (J.P.L., Pasadena).     

Ancient heat flow and crustal thickness at Warrego rise, Thaumasia highlands, Mars: Implications for a stratified crust

Heat flow calculations based on geological and/or geophysical indicators can help to constrain the thickness, and potentially the geochemical stratification, of the martian crust. Here we analyze the Warrego rise region, part of the ancient mountain range referred to as the Thaumasia highlands. This region has a crustal thickness much greater than the martian average, as well as estimations of the depth to the brittle-ductile transition beneath two scarps interpreted to be thrust faults. For the local crustal density (2900 kg m⁻³) favored by our analysis of the flexural state of compensation of the local topography, the crustal thickness is at least 70 and 75 km at the scarp locations. However, for one of the scarp locations our nominal model does not obtain heat flow solutions permitting a homogeneous crust as thick as required. Our results, therefore, suggest that the crust beneath the Warrego rise region is chemically stratified with a heat-producing enriched upper layer thinner than the whole crust. Moreover, if the mantle heat flow (at the time of scarp formation) was higher than 0.3 of the surface heat low, as predicted by thermal history models, then a stratified crust rise seems unavoidable for this region, even if local heat-producing element abundances lower than average or hydrostatic pore pressure are considered. This finding is consistent with a complex geological history, which includes magmatic-driven activity.     

Geology and tectonics of Beta Regio,Guinevere Planitia,Sedna Planitia,and Western Eistla Regio,Venus: Results from Arecibo image data

New radar images (resolution 1.5–2.0 km) obtained from the Arecibo Observatory are used to assess the geology of a portion of the equatorial region of Venus (1 S to 45 N and from 270 eastward to 30). Nine geologic units are mapped on the basis of their radar characteristics and their distribution and correspondences with topography are examined. Plains are the most abundant unit, make up 80%; of the area imaged, and are divided into bright, dark, and mottled. Mottled plains contain abundant lava flows and domes suggesting that volcanism forming plains is a significant process in the equatorial region of Venus. Tesserae are found primarily on Beta Regio and its eastern flank and are interpreted to be locally stratigraphically older units, predating episodes of faulting and plains formation. Isolated regions of tesserae concentrated to the north of Western Eistla Regio are interpreted to predate the formation of plains in this area. The volcanoes Sif Mons, Gula Mons, Sappho, Theia Mons, and Rhea Mons, are found exclusively in highland regions and their deposits are interpreted as contributing only a small percentage to the overall volume of the regional topography. The northern 15 of the image data overlap with Venera 15/16 images making it possible to examine the characteristics of geologic units mapped under various illumination directions and incidence angles. Surface panoramas and geochemical data obtained from Venera landers provide ground truth for map units, evidence that plains are made up of basaltic lava flows, and that linear deformation zones contain abundant blocks and cobbles. On the basis of spatial and temporal relationships between geologic units, the highlands of Beta Regio and Western Eistla Regio are interpreted to have formed in association with areas of mantle upwelling which uplift plains, cause rifting, and in the case of Beta Regio, disrupt a large region of tessera. Zones of linear belt deformation in Beta Regio and Western Eistla Regio are interpreted to be extensional and indicate that at least limited extension has occurred in both regions. The images reveal for the first time that southern Devana Chasma is a region of overlapping rift valleys separated by a distance of 600 km. Linear deformation zones in Guinevere Planitia, separating Beta Regio and Eistla Regio, converge at a region of ovoids forming a discontinuous zone of disruption and completes an equatorial encompassing network of highlands and tectonic features. The similarity between ovoids and coronae suggests a mechanism of formation associated with hotspots or mantle plumes. Analysis of the distribution and density of impact craters suggests a surface age for this part of the planet similar to or slightly less than that determined for the northern high latitudes from Venera 15/16 data (0.3 to 1.5 by) and comparable to that calculated for the southern hemisphere.     

Dynamic topography of the East European craton: Shedding light upon lithospheric structure, composition and mantle dynamics

Most of the East European Craton lacks surface relief; however, the amplitude of topography at the top of the basement exceeds 20 km, the amplitude of topography undulations at the crustal base reaches almost 30 km with an amazing amplitude of ca. 50 km in variation in the thickness of the crystalline crust, and the amplitude of topography variations at the lithosphere–asthenosphere boundary exceeds 200 km. This paper examines the relative contributions of the crust, the subcrustal lithosphere, and the dynamic support of the sublithospheric mantle to maintain surface topography, using regional seismic data on the structure of the crystalline crust and the sedimentary cover, and thermal and large-scale P- and S-wave seismic tomography data on the structure of the lithospheric mantle. For the Precambrian lithosphere, an analysis of Vp/Vs ratio at 100, 150, 200, and 250 km depths does not show any age-dependence, suggesting that while Vp/Vs ratio can be effectively used to outline the cratonic margins, it is not sensitive to compositional variations within the cratonic lithosphere.Statistical analysis of age-dependence of velocity, density, and thermal structure of the continental crust and subcrustal lithosphere in the study area (0–62E, 45–72N) allows to link lithospheric structure with the tectonic evolution of the region since the Archean. Crustal thickness decreases systematically with age from 42–44 km in regions older than 1.6 Ga to 37–40 km in the Paleozoic–Mesoproterozoic structures, and to ca. 31 km in the Meso-Cenozoic regions. However, the isostatic contribution of the crust to the surface topography of the East European Craton is almost independent of age (ca. 4.5 km) due to an interplay of age-dependent crustal and sedimentary thicknesses and lithospheric temperatures.On the contrary, the contribution of the subcrustal lithosphere to the surface topography strongly depends on the age, being slightly positive (+ 0.3 + 0.7 km) for the regions older than 1.6 Ga and negative (− 0.5–1 km) for younger structures. This leads to age-dependent variations in the residual topography, i.e. the topography which cannot be explained by the assumed thermal and density structure of the lithosphere, and which can (at least partly) originate from the dynamic component caused by the mantle flow. Positive dynamic topography at the cratonic margins, which exceeds 2 km in the Norwegian Caledonides and in the Urals, clearly links their on-going uplift with deep mantle processes. Negative residual topography beneath the Archean-Paleoproterozoic cratons (− 1–2 km) indicates either a smaller density deficit (ca. 0.9%) in their subcrustal lithosphere than predicted by global petrologic data on mantle-derived xenoliths or the presence of a strong convective downwelling in the mantle. Such mantle downflows can effectively divert heat from the lithospheric base, leading to a long-term survival of the Archean-Paleoproterozoic lithosphere.     

Crustal recycling, mantle dehydration, and the thermal evolution of Mars

We have reinvestigated the coupled thermal and crustal evolution of Mars taking new laboratory data concerning the flow behavior of iron-rich olivine into account. The low mantle viscosities associated with the relatively higher iron content of the martian mantle as well as the observed high concentrations of heat producing elements in a crust with a reduced thermal conductivity were found to promote phases of crustal recycling in many models. As crustal recycling is incompatible with an early separation of geochemical reservoirs, models were required to show no episodes of crustal recycling. Furthermore, admissible models were required to reproduce the martian crust formation history, to allow for the formation of partial melt under present day mantle conditions and to reproduce the measured concentrations of potassium and thorium on the martian surface. Taking dehydration stiffening of the mantle viscosity by the extraction of water from the mantle into account, we found that admissible models have low initial upper mantle temperatures around 1650 K, preferably a primordial crustal thickness of 30 km, and an initially wet mantle rheology. The crust formation process on Mars would then be driven by the extraction of a primordial crust after core formation, cooling the mantle to temperatures close to the peridotite solidus. According to this scenario, the second stage of global crust formation took place over a more extended period of time, waning at around 3500 Myr b.p., and was driven by heat produced by the decay of radioactive elements. Present-day volcanism would then be driven by mantle plumes originating at the core-mantle boundary under regions of locally thickened, thermally insulating crust. Water extraction from the mantle was found to be relatively efficient and close to 40% of the total inventory was lost from the mantle in most models. Assuming an initial mantle water content of 100 ppm and that 10% of the extracted water is supplied to the surface, this amount is equivalent to a 14 m thick global surface layer, suggesting that volcanic outgassing of H₂O could have significantly influenced the early martian climate and increased the planet’s habitability.     

Rayleigh-Taylor instability as a mechanism for corona formation on Venus

In this study we explore the idea that coronae have formed on Venus as a result of gravitational (Rayleigh-Taylor) instability of the lithosphere. The lithosphere is represented by a system of stratified homogeneous viscous layers (low-density crust over high density mantle, over lower density layer beneath the lithosphere). A small harmonic perturbation imposed on the base of the lithosphere is observed to result in gravitational instability under the constraint of assumed axisymmetry. Topography develops with time under the influence of dynamic stress associated with downwelling or upwelling, and spatially variable crustal thickening or thinning. Topography may therefore be elevated or depressed above a mantle downwelling, but the computed gravity anomaly is always negative above a mantle downwelling in a homogeneous asthenosphere. The ratio of peak gravity to topography anomaly depends primarily on the ratio of crust to lithospheric viscosity. Average observed ratios are well resolved for two groups of coronae (∼40 mgal km⁻¹), consistent with models in which the crust is perhaps 5 times stronger than the lithosphere. Group 3a (rim surrounding elevated central region) coronae are inferred to arise from a central upwelling model, whereas Group 8 (depression) coronae are inferred to arise from central downwelling. Observed average coronae radii are consistent with a lithospheric thickness of only 50 km. An upper low-density crustal layer is 10-20 km thick, as inferred from the amplitude of gravity and topography anomalies.     

Diurnal CO variations in the Venus mesophere from CO microwave spectra

Microwave spectra of carbon monoxide (¹²CO) in the mesosphere of Venus were measured in December 1978, May and December 1980, and January, September, and November 1982. These spectra are analyzed to provide mixing profiles of CO in the Venus mesosphere and best constrain the mixing profile of CO between ～ 100 and 80 km altitude. From the January 1982 measurement (which, of all our spectra, best constrains the abundance of CO below 80 km altitude) we find an upper limit for the CO mixing ratio below 80 km altitude that is two to three times smaller than the stratospheric (～65 km) value of 4.5 ± 1.0 × 10^?5 determined by P. Connes, J. Connes, L.D. Kaplan, and W. S. Benedict (1968, Astrophys. J.152, 731–743) in 1967, indicating a possible long-term change in the lower atmospheric concentration of CO. Intercomparison among the individual CO profiles derived from our spectra indicates considerable short-term temporal and/or spatial variation in the profile of CO mixing in the Venus mesosphere above 80 km. A more complete comparison with previously published CO microwave spectra from a number of authors specifies the basic diurnal nature of mesospheric CO variability. CO abundance above ～ 95 km in the Venus atmosphere shows approximately a factor of 2–4 enhancement on the nightside relative to the dayside of Venus. Peak nightside CO abundance above ～95 km occurs very near to the antisolar point on Venus (local time of peak CO abundance above ～95 km occurs at 0.6_?0.6^+0.7 hr after midnight on Venus), strongly suggesting that retrograde zonal flow is substantially reduced at an altitude of 100 km in the Venus mesosphere. In contrast, CO abundances between 80 and 90 km altitude show a maximum that is shifted from the antisolar point toward the morningside of Venus (local time of peak CO abundance between 80 and 90 km occurs at 8.5 ± 1.0 hr past midnight on Venus). The magnitude of the diurnal variation of CO abundance between 80 and 90 km is again, approximately a factor of 2–4. Disk-averaged spectra of Venus do not determine the exact form for the diurnal distribution of CO in the Venus mesosphere as indicated by comparison of synthetic spectra, based upon model distributions, and the measured spectra. However, the offset in phase for the diurnal variation for the >95 km and 80–90-km-altitude regions requires an asymmetric (in solar zenith angle) distribution.     

Constraints on mantle plumes on Venus: Implications for volatile history

The analysis of Venus’ gravity field and topography suggests the presence of a small number of deep mantle plumes (～9). This study predicts the number of plumes formed at the core–mantle boundary, their characteristics, and the production of partial melt from adiabatic decompression. Numerical simulations are performed using a 3D spherical code that includes large viscosity variations and internal heating. This study investigates the effect of several parameters including the core–mantle boundary temperature, the amount of internal heating, and the mantle viscosity. The smallest number of plumes is achieved when no internal heating is present. However, scaling Earth’s radiogenic heating to Venus suggests a value of ～16 TW. Cases with internal heating produce more realistic lid thickness and partial melting, but produce either too many plumes or no plumes if a high mantle temperature precludes the formation of a hot thermal boundary layer. Mantle viscosity must be reduced to at least 10²⁰ Pa s in order to include significant internal heating and still produce hot plumes. In all cases that predict melting, melting occurs throughout the upper mantle. Only cases with high core temperature (>1700 K) produce dry melting. Over time the upper mantle may have lost significant volatiles. Depending on the water content of the lower mantle, deep plumes may contribute to present-day atmospheric water via volcanic outgassing. Assuming 50 ppm water in mantle, 10 plumes with a buoyancy flux of 500 kg/s continuously erupting for 4 myr will outgas an amount of water on the order of that in the lower atmosphere. A higher level of internal heating than achieved to date, as well as relatively low mantle viscosity, may be required to achieve simulations with ～10 plumes and a thinner lid. Alternatively, if the mantle is heating up due to the stagnant lid, the effect is equivalent to having lower rates of internal heating. A temperature increase of 110 K/byr is equivalent to ?13 TW. This value along with the internal heating of 3 TW used in this study may represent the approximate heat budget of Venus’ mantle.     

Towards self-consistent modeling of the martian dichotomy: The influence of one-ridge convection on crustal thickness distribution

In order to find an explanation for the origin of the martian crustal dichotomy, a number of recent papers have examined the effect of layered viscosity on the evolution of a degree-1 mantle convection, e.g. Roberts and Zhong [Roberts, J.H., Zhong, S., 2006. J. Geophys. Res. 111. E06013] and Yoshida and Kageyama [Yoshida, M., Kageyama, A., 2006. J. Geophys. Res. 111, doi:10.1029/2005JB003905. B03412]. It was found that a mid-mantle viscosity jump, combined with highly temperature- and depth-dependent rheology, are effective in developing a degree-1 convection within a short timescale. Such a layered viscosity profile could be justified by martian mineralogy. However, the effect of a degree-1 convective planform on the crustal thickness distribution has not yet been demonstrated. It is not obvious whether a thinner crust, due to sublithospheric erosion and crustal thinning, or a thicker crust, due to enhanced crustal production, would form above the hemisphere of mantle upwelling. Also, the general shape of the dichotomy, which is not strictly hemispherical, has not yet been fully investigated. Here we propose a model of the crustal patterns produced by numerical simulations of martian mantle convection, using the finite-volume multigrid code StagYY [Tackley, P.J., 2008. Phys. Earth Planet. Int. 107, 7-18, doi:10.1016/j.pepi.2008.08.005] A self-consistent treatment of melting, crustal formation and chemical differentiation has been added to models of three-dimensional thermal convection. This allows us to obtain global maps of the crustal thickness distribution as it evolves with time. The obtained results demonstrate that it is indeed possible to form a crustal dichotomy as a consequence of near degree-1 mantle convection early in Mars' history. We find that some of the observed patterns show intriguing first order similarities to the elliptical shape of the martian dichotomy. In all models, the region of thick crust is located over the region of mantle upwelling, which itself is a ridge-like structure spread over roughly one half of the planet, a planform we describe as “one-ridge convection.”     

Divergent evolution of Earth and Venus: Influence of degassing,tectonics, and magnetic fields

Knowledge of the earliest evolution of Earth and Venus is extremely limited, but it is obvious from their dramatic contrasts today that at some point in their evolution conditions on the two planets diverged. In this paper we develop a geophysical systems box model that simulates the flux of carbon through the mantle, atmosphere, ocean, and seafloor, and the degassing of water from the mantle. Volatile fluxes, including loss to space, are functions of local volatile concentration, degassing efficiency, tectonic plate speed, and magnetic field intensity. Numerical results are presented that demonstrate the equilibration to a steady state carbon cycle, where carbon and water are distributed among mantle, atmosphere, ocean, and crustal reservoirs, similar to present-day Earth. These stable models reach steady state after several hundred million years by maintaining a negative feedback between atmospheric temperature, carbon dioxide weathering, and surface tectonics. At the orbit of Venus, an otherwise similar model evolves to a runaway greenhouse with all volatiles in the atmosphere. The influence of magnetic field intensity on atmospheric escape is demonstrated in Venus models where either a strong magnetic field helps the atmosphere to retain about 60 bars of water vapor after 4.5 Gyr, or the lack of a magnetic field allows for the loss of all atmospheric water to space in about 1 Gyr. The relative influences of plate speed and degassing rate on the weathering rate and greenhouse stability are demonstrated, and a stable to runaway regime diagram is presented. In conclusion, we propose that a stable climate-tectonic-carbon cycle is part of a larger coupled geophysical system where a moderate surface climate provides a stabilizing feedback for maintaining surface tectonics, the thermal cooling of the deep interior, magnetic field generation, and the shielding of the atmosphere over billion year time scales.     

Models of the internal structure of Venus

A survey is made of the physics of the interiors of Venus. The introduction explains the main concepts used in the construction of models of Venus and the history of the question; observational data are gathered and analyzed. The method of constructing the models of the planet is explained and earth-like models of Venus and parametrically simple PVM models are discussed. Within the compass of a physical nodel of Venus, the thermodynamics of the mantle and core is constructed and questions are discussed concerning the heat conduction, temperature distribution in the lithosphere and the thermal flux from the interior of Venus, the electrical conduction and mechanical quality, and large-scale steady stresses in the mantle of Venus. A rheological model of the crust and mantle is constructed. In conclusion, the question as to the distribution of radioactivity and convection in the interior of the planet is discussed.     

Tectonics of Venus: A review

The article presents a new tectonic scheme of Venus and gives the following interpretation of the planet's main structural units: (1) plains — areas of flood volcanism over stretched crust; (2) dome-like uplifts — areas of uplifting and volcanic activity above the mantle hot-spots; (3) coronae —former dome-like uplifts, partially subsided and diffused by gravity; (4) ridge belts — fold zones; (5) tesserae — fragments of ductile compression and shortening of crust; (6) supercoronae — coronae formed in the course of further evolution and relaxation of Beta-type uplifts. Ishtar Terra is considered to be a fragment of an ancient tessera paleocontinent, on the edge of which the Lakshmi supercorona is superimposed. Aphrodite Terra is considered as a belt of mantle hot-spot structures (dome-like uplifts, coronae, supercoronae, volcanoes, rifts).Three types of planetary belts have been distinguished on Venus: uplifted 'weakened' belts with an abundance of mantle hot-spot structures; a northern fan of ridge belts; and belts of low basalt plains. The center of the planetary system of uplifted weakened belts is situated in Atla Regio.The present tectonic structure of Venus is inferred to have formed during two stages of evolution characterized by different tectonic regimes. Stage I is a regime of soft ductile plates (formation of tessera uplifts and volcanic plains). Stage II is a formation of 'weakened' uplifted planetary belts, various tectonic regimes of mantle hot-spots, and plains-forming volcanism.'Geology and Tectonics of Venus', special issue edited by Alexander T. Basilevsky (USSR Acad. of Sci. Moscow), James W. Head (Brown University, Providence), Gordon H. Pettengill (MIT, Cambridge, Massachusetts) and R. S. Saunders (J.P.L., Pasadena).