Transitions in thermal convection with strongly temperature-dependent viscosity in a wide box

Numerical models are systematically presented for time-dependent thermal convection of Newtonian fluid with strongly temperature-dependent viscosity in a two-dimensional rectangular box of aspect ratio 3 at various values of the Rayleigh number Ra_b defined with viscosity at the bottom boundary up to 1.6×10⁸ and the viscosity contrast across the box r_η up to 10⁸. We found that there are two different series of bifurcations that take place as r_η increases. One series of bifurcations causes changes in the behavior of the thermal boundary layer along the surface boundary from small-viscosity-contrast (SVC) mode, through transitional (TR) mode, to stagnant-lid (ST) mode, or from SVC mode directly to ST mode, depending on Ra_b. Another series of bifurcations causes changes in the aspect ratio of convection cells; convection with an elongated cell can take place at moderate r_η (10³–10^5.5 at Ra_b=6×10⁶), while only convection of aspect ratio close to 1 takes place at small r_η and large r_η. The parameter range of r_η and Ra_b for elongated-cell convection overlaps the parameter range for SVC and ST modes and include the entire parameter range for TR mode. In the elongated-ST regime, the lid of highly viscous fluid along the top boundary is not literally ‘stagnant’ but can horizontally move at a velocity high enough to induce a convection cell with aspect ratio much larger than 1.     

Contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments

We examined the contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments by conducting detailed hydrochemical observations in five catchments that ranged from zero to second order. End‐member mixing analysis (EMMA) was performed to identify the geographical sources of stream water. Throughfall, hillslope groundwater, shallow bedrock groundwater, and deep bedrock groundwater were identified as end members. The contribution of each end member to storm runoff differed among the catchments because of the differing quantities of riparian groundwater, which was recharged by the bedrock groundwater prior to rainfall events. Among the five catchments, the contribution of throughfall was highest during both baseflow and storm flow in a zero‐order catchment with little contribution from the bedrock groundwater to the riparian reservoir. In zero‐order catchments with some contribution from bedrock groundwater, stream water was dominated by shallow bedrock groundwater during baseflow, but it was significantly influenced by hillslope groundwater during storms. In the first‐order catchment, stream water was dominated by shallow bedrock groundwater during storms as well as baseflow periods. In the second‐order catchment, deeper bedrock groundwater than that found in the zero‐order and first‐order catchments contributed to stream water in all periods, except during large storm events. These results suggest that bedrock groundwater influences the upscaling of storm‐runoff generation processes by affecting the linkages of geomorphic units such as hillslopes, riparian zones, and stream channels. Our results highlight the need for a three‐dimensional approach that considers bedrock groundwater flow when studying the upscaling of storm‐runoff generation processes. Copyright © 2014 John Wiley & Sons, Ltd.     

Experimental study of reaction between perovskite and molten iron to 146 GPa and implications for chemically distinct buoyant layer at the top of the core

Partitioning of oxygen and silicon between molten iron and (Mg,Fe)SiO₃ perovskite was investigated by a combination of laser-heated diamond-anvil cell (LHDAC) and analytical transmission electron microscope (TEM) to 146 GPa and 3,500 K. The chemical compositions of co-existing quenched molten iron and perovskite were determined quantitatively with energy-dispersive X-ray spectrometry (EDS) and electron energy loss spectroscopy (EELS). The results demonstrate that the quenched liquid iron in contact with perovskite contained substantial amounts of oxygen and silicon at such high pressure and temperature (P–T). The chemical equilibrium between perovskite, ferropericlase, and molten iron at the P–T conditions of the core–mantle boundary (CMB) was calculated in Mg–Fe–Si–O system from these experimental results and previous data on partitioning of oxygen between molten iron and ferropericlase. We found that molten iron should include oxygen and silicon more than required to account for the core density deficit (<10%) when co-existing with both perovskite and ferropericlase at the CMB. This suggests that the very bottom of the mantle may consist of either one of perovskite or ferropericlase. Alternatively, it is also possible that the bulk outer core liquid is not in direct contact with the mantle. Seismological observations of a small P-wave velocity reduction in the topmost core suggest the presence of chemically-distinct buoyant liquid layer. Such layer physically separates the mantle from the bulk outer core liquid, hindering the chemical reaction between them.     

Verification of two-dimensional non-linear analysis of sand-structure system by examining the results of the shaking-table test

The verification of a two-dimensional non-linear analysis based on a hybrid constitutive model, which consists of the Ramberg-Osgood model extended to two-dimensional problems and a dilatancy model, is discussed through a comparison of the shaking-table test results for a one-storey structure standing on a dry sandy deposit with those for the same on a saturated sandy deposit. Since the relationships G vs. γ and h vs. γ for shear strains of 10^?5?10^?3 can be obtained accurately by an inversion analysis of the resonance curve of the sandy deposit the proposed non-linear method can represent well the observed non-linear response of the dry or the saturated sandy deposit including the structure. This method, however, should be applied carefully to assess the response of a soil-structure system when the three-dimensional interaction affects significantly the response of the model. An equivalent linear analysis using a sway-rocking model is applied to simulate the non-linear ground motion including the three-dimensional interaction, and it is found that the sway-rocking model can also represent well the non-linear response of the system.     

Molecular dynamics study of the structures and bulk moduli of crystals in the system CaO-MgO-Al2O3-SiO2

Molecular dynamics (MD) simulations have been used to calculate the structures and bulk moduli of crystals in the system CaO-MgO-Al₂O₃-SiO₂ (CMAS) using an interatomic potential model (CMAS94), which is composed of pairwise additive Coulomb, van der Waals, and repulsive interactions. The crystals studied, total of 27, include oxides, Mg meta- and ortho-silicates, Al garnets, and various Ca or Al bearing silicates, with the coordination number of cations ranging 6 to 12 for Ca, 4 to 12 for Mg, 4 to 6 for Al, and 4 and 6 for Si. In spite of the simplicity of the CMAS94 potential and the diversity of the structural types treated, MD simulations are quite satisfactory in reproducing well the observed structural data, including the crystal symmetries, lattice parameters, and average and individual nearest neighbour Ca-O, Mg-O, Al-O, and Si-O distances. In addition MD simulated bulk moduli of crystals in the CMAS system compare well with the observed values.     

Hydrological control of dissolved organic carbon dynamics in a forested headwater catchment,Kiryu Experimental Watershed,Japan

To evaluate the influence of hydrological processes on dissolved organic carbon (DOC) dynamics in a forested headwater catchment, DOC concentration was observed along the flow path from rainfall to stream water via throughfall, soil water, groundwater, and spring water for 4 years, and DOC flux through the catchment was calculated. The spatial and temporal variations in DOC concentration and flux were compared with physical hydrological observations and the mean residence time of water. In the upslope soil layer, DOC concentrations were not significantly correlated with water fluxes, suggesting that DOC concentrations were not strictly controlled by water fluxes. In the upslope perennial groundwater, DOC concentration was affected by the change in the amount of microbial degradation of DOC produced by changes in the mean residence time of water. In stream water, the temporal variation in DOC concentration was usually affected by changes in DOC concentration of the inflow component via vertical infiltration from above the perennial groundwater. During dry periods, however, the component from inflow via vertical infiltration was negligible and DOC in the upslope perennial groundwater became the major component of stream water DOC. The temporal variation in stream water DOC concentration during baseflow was affected by rainfall patterns over several preceding months. Therefore, records of rainfall over several preceding months are one of the most important factors for predicting changes in DOC concentration on a catchment scale. Copyright © 2007 John Wiley & Sons, Ltd.