The Earth's surface volcanism exerts first-order controls on the composition of the atmosphere and the climate. On Earth, the majority of surface volcanism occurs at mid-ocean ridges. In this study, based on the dependence of melt fraction on temperature, pressure, and composition, we compute melt production and degassing rate at mid-ocean ridges from three-dimensional global mantle convection models with plate motion history as the surface velocity boundary condition. By incorporating melting in global mantle convection models, we connect deep mantle convection to surface volcanism, with deep and shallow mantle processes internally consistent. We compare two methods to compute melt production: a tracer method and an Eulerian method. Our results show that melt production at mid-ocean ridges is mainly controlled by surface plate motion history, and that changes in plate tectonic motion, including plate reorganizations, may lead to significant deviation of melt production from the expected scaling with seafloor production rate. We also find a good correlation between melt production and degassing rate beneath mid-ocean ridges. The calculated global melt production and CO2 degassing rate at mid-ocean ridges varies by as much as a factor of 3 over the past 200 Myr. We show that mid-ocean ridge melt production and degassing rate would be much larger in the Cretaceous, and reached maximum values at ∼150–120 Ma. Our results raise the possibility that warmer climate in the Cretaceous could be due in part to high magmatic productivity and correspondingly high outgassing rates at mid-ocean ridges during that time.

The lowermost few hundreds of kilometres of the Earth's mantle are elastically anisotropic; seismic velocities vary with direction of propagation and polarization. Observations of strong seismic anisotropy correlate with regions where subducted slab material is expected. In this study, we evaluate the hypothesis that crystal preferred orientation (CPO) in a slab, as it impinges on the core-mantle boundary, is the cause of the observed anisotropy. Next, we determine if fast polarization directions seen by shear waves can be mapped to directions of geodynamic flow. This approach is similar to our previous study performed for a 2-D geodynamic model. In this study, we employ a 3-D geodynamic model with temperaturedependent viscosity and kinematic velocity boundary conditions defined at the surface of the Earth to create a broad downwelling slab. Tracers track the deformation that we assume to be accommodated by dislocation creep. We evaluate the models for the presence of perovskite or post-perovskite and for different main slip systems along which dislocation creep may occur in post-perovskite [(100),(010) and (001)]-resulting in four different mineralogical models of CPO. Combining the crystal pole orientations with single crystal elastic constants results in seismically distinguishable models of seismic anisotropy. The models are evaluated against published seismic observations by analysing different anisotropic components: the radial anisotropy, the splitting for (sub-)vertical phases (i.e. azimuthal anisotropy), and the splitting for subhorizontal phases. The patterns in radial anisotropy confirm our earlier results in 2-D. Observations of radial anisotropy and splitting in subhorizontal phases are mostly consistent with our models of post-perovskite with (010)-slip and (001)-slip. Our model of (001)-slip predicts stronger splitting than for (010)-slip for horizontally propagating phases in all directions. The strongest seismic anisotropy in this model occurs where the slab impinges on the core-mantle boundary. The azimuthal anisotropy pattern for (001)-slip shows fast axis directions at the edges of the slab (sub-)parallel to flow directions, suggesting horizontal flows may be mapped out in the lowermost mantle using seismic observations. © The Authors 2014. Published by Oxford University Press on behalf of The Royal Astronomical Society.

Seismic tomography models reveal two large low shear velocity provinces (LLSVPs) that identify large-scale variations in temperature and composition in the deep mantle. Other characteristics include elevated density, elevated bulk sound speed, and sharp boundaries. We show that properties of LLSVPs can be explained by the presence of small quantities (0.3-3%) of suspended, dense Fe-Ni-S liquid. Trapping of metallic liquid is demonstrated to be likely during the crystallization of a dense basal magma ocean, and retention of such melts is consistent with currently available experimental constraints. Calculated seismic velocities and densities of lower mantle material containing low-abundance metallic liquids match the observed LLSVP properties. Small quantities of metallic liquids trapped at depth provide a natural explanation for primitive noble gas signatures in plume-related magmas. Our model hence provides a mechanism for generating large-scale chemical heterogeneities in Earth's early history and makes clear predictions for future tests of our hypothesis.

As a major greenhouse gas, the continuous increase of carbon dioxide (CO₂) in the atmosphere has caused serious environmental problems, although CO₂ is also an abundant, inexpensive, and nontoxic carbon source. Here, we use metal-organic framework (MOF) with highly ordered hierarchical structure as adsorbent and catalyst for chemical fixation of CO₂ at atmospheric pressure, and the CO₂ can be converted to the formate in excellent yields. Meanwhile, we have successfully integrated highly ordered macroporous and mesoporous structures into MOFs, and the macro-, meso-, and microporous structures have all been presented in one framework. Based on the unique hierarchical pores, high surface area (592 m²/g), and high CO₂ adsorption capacity (49.51 cm³/g), the ordered macroporous-mesoporous MOFs possess high activity for chemical fixation of CO₂ (yield of 77%). These results provide a promising route of chemical CO₂ fixation through MOF materials.