Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

Seasonal Variation of the East Asian Subtropical Westerly Jet and Its Association with the Heating Field over East Asia

The structure and seasonal variation of the East Asian Subtropical Westerly Jet （EAWJ） and associations with heating fields over East Asia are examined by using NCEP/NCAR reanalysis data. Obvious differences exist in the westerly jet intensity and location in different regions and seasons due to the ocean-land distribution and seasonal thermal contrast, as well as the dynamic and thermodynamic impacts of the Tibetan Plateau. In winter, the EAWJ center is situated over the western Pacific Ocean and the intensity is reduced gradually from east to west over the East Asian region. In summer, the EAWJ center is located over the north of the Tibetan Plateau and the jet intensity is reduced evidently compared with that in winter. The EAWJ seasonal evolution is characterized by the obvious longitudinal inconsistency of the northward migration and in-phase southward retreat of the EAWJ axis. A good correspondence between the seasonal variations of EAWJ and the meridional differences of air temperature （MDT） in the mid-upper troposphere demonstrates that the MDT is the basic reason for the seasonal variation of EAWJ. Correlation analyses indicate that the Kuroshio Current region to the south of Japan and the Tibetan Plateau are the key areas for the variations of the EAWJ intensities in winter and in summer, respectively. The strong sensible and latent heating in the Kuroshio Current region is closely related to the intensification of EAWJ in winter. In summer, strong sensible heating in the Tibetan Plateau corresponds to the EAWJ strengthening and southward shift, while the weak sensible heating in the Tibetan Plateau is consistent with the EAWJ weakening and northward migration.     

Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.     

Individual and combined effects of land use/cover and climate change on Wolf Bay watershed streamflow in southern Alabama

Land use/cover (LULC) and climate change are two main factors affecting watershed hydrology. In this paper, individual and combined impacts of LULC and climate change on hydrologic processes were analysed applying the model Soil and Water Assessment Tool in a coastal Alabama watershed in USA. Temporally and spatially downscaled Global Circulation Model outputs predict a slight increase in precipitation in the study area, which is also projected to experience substantial urban growth in the future. Changes in flow frequency and volume in the 2030s (2016–2040) compared to a baseline period (1984–2008) at daily, monthly and annual time scales were explored. A redistribution of daily streamflow is projected when either climate or LULC change was considered. High flows are predicted to increase, while low flows are expected to decrease. Combined change effect results in a more noticeable and uneven distribution of daily streamflow. Monthly average streamflow and surface runoff are projected to increase in spring and winter, but especially in fall. LULC change does not have a significant effect on monthly average streamflow, but the change affects partitioning of streamflow, causing higher surface runoff and lower baseflow. The combined effect leads to a dramatic increase in monthly average streamflow with a stronger increasing trend in surface runoff and decreasing trend in baseflow. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatially controlled Fe and Si isotope variations: an alternative view on the formation of the Torres del Paine pluton

We present new Fe and Si isotope ratio data for the Torres del Paine igneous complex in southern Chile. The multi-composition pluton consists of an approximately 1 km vertical exposure of homogenous granite overlying a contemporaneous 250-m-thick mafic gabbro suite. This first-of-its-kind spatially dependent Fe and Si isotope investigation of a convergent margin-related pluton aims to understand the nature of granite and silicic igneous rock formation. Results collected by MC-ICP-MS show a trend of increasing δ⁵⁶Fe and δ³⁰Si with increasing silica content as well as a systematic increase in δ⁵⁶Fe away from the mafic base of the pluton. The marginal Torres del Paine granites have heavier Fe isotope signatures (δ⁵⁶Fe = +0.25 ± 0.02 2se) compared to granites found in the interior pluton (δ⁵⁶Fe = +0.17 ± 0.02 2se). Cerro Toro country rock values are isotopically light in both Fe and Si isotopic systems (δ⁵⁶Fe = +0.05 ± 0.02 ‰; δ³⁰Si = ?0.38 ± 0.07 ‰). The variations in the Fe and Si isotopic data cannot be accounted for by local assimilation of the wall rocks, in situ fractional crystallization, late-stage fluid exsolution or some combination of these processes. Instead, we conclude that thermal diffusion or source magma variation is the most likely process producing Fe isotope ratio variations in the Torres del Paine pluton.     

Thermal equation of state of natural tourmaline at high pressure and temperature

Synchrotron-based in situ angle-dispersive X-ray diffraction experiments were conducted on a natural uvite-dominated tourmaline sample by using an external-heating diamond anvil cell at simultaneously high pressures and temperatures up to 18 GPa and 723 K, respectively. The angle-dispersive X-ray diffraction data reveal no indication of a structural phase transition over the P–T range of the current experiment in this study. The pressure–volume–temperature data were fitted by the high-temperature Birch–Murnaghan equation of state. Isothermal bulk modulus of K ₀ = 96.6 (9) GPa, pressure derivative of the bulk modulus of \(K_{0}^{\prime } = 12.5 \;(4)\), thermal expansion coefficient of α ₀ = 4.39 (27) × 10^?5 K^?1 and temperature derivative of the bulk modulus (?K/?T) _P  = ?0.009 (6) GPa K^?1 were obtained. The axial thermoelastic properties were also obtained with K _a0 = 139 (2) GPa, \(K_{a0}^{\prime }\) = 11.5 (7) and α _a0 = 1.00 (11) × 10^?5 K^?1 for the a-axis, and K _c0 = 59 (1) GPa, \(K_{c0}^{\prime }\) = 11.4 (5) and α _c0 = 2.41 (24) × 10^?5 K^?1 for the c-axis. Both of axial compression and thermal expansion exhibit large anisotropic behavior. Thermoelastic parameters of tourmaline in this study were also compared with that of the other two ring silicates of beryl and cordierite.     

Biogenic volatile compound emissions from a temperate forest,China: model simulation

Emissions of biogenic volatile organic compounds (BVOC) were measured using a relaxed eddy accumulation (REA) technique on an above-canopy tower in a temperate forest (Changbai Mountain, Jilin province, China) during the 2010 and 2011 summer seasons. Solar global radiation and photosynthetically active radiation (PAR) were also measured. Based on PAR energy dynamic balance, an empirical BVOC emission and PAR transfer model was developed that includes the processes of BVOC emissions and PAR transfer above the canopy level, including PAR absorption and consumption, and scattering by gases, liquids, and particles (GLPs). Simulated emissions of isoprene and monoterpenes were in agreement with observations. The averages of the relative estimator biases for the flux were 39.3 % for isoprene, and 27.1 % for monoterpenes in the 2010 and 2011 growing seasons, with NMSE (normalized mean square error) values of 0.133 and 0.101, respectively. The observed and simulated mean diurnal variations of isoprene and monoterpenes in the 2010 and 2011 growing seasons were evaluated for the validation of the empirical model. Under observed atmospheric conditions, the sensitivity analysis showed that emissions of isoprene and monoterpenes were more sensitive to changes in PAR than to water vapor content or to the magnitude of the scattering factor. The emissions of isoprene and monoterpenes in the 2010 and 2011 growing seasons (from June to September) were estimated using this empirical model along with hourly observational data, with mean hourly emissions of 1.71 and 1.55 mg m^?2 h^?1 for isoprene, and 0.48 and 0.47 mg m^?2 h^?1 for monoterpenes in 2010 and 2011, respectively. As formaldehyde (HCHO) is considered as the main oxidation product of isoprene and monoterpenes, it is necessary to investigate the link between HCHO and BVOC emissions. GOME-2 HCHO vertical column densities (VCDs) can be used to estimate BVOC emission fluxes in the Changbai Mountain temperate forest.     

闽西北5月雨量特多年份环流特征及其成因分析

利用闽西北三明市的11个县、市1961—2014年5月降雨量资料、常规气象观测资料和NCEP 2.5×2.5°再分析资料等对闽西北历史上5月雨量特多年的降雨量时空分布特征和大气环流特征及其成因进行了分析,结果表明:5月区域性、全区性雨量特多年的环流异常特征是500hPa欧亚中高纬度自西向东环流呈"+-+"的波列分布,鄂霍次克海北侧高度正距平区的稳定少动,起到阻挡冷空气东移的作用,冷空气沿乌拉尔山西侧高度正距平中心前的偏北气流不断南下,使萨彦岭一带高度负距平中心加强,其底部不断分裂小槽东移影响闽西北。925hPa华南北部维持一条稳定切变线;对于2014年5月特例的环流形势分析结果可知,500hPa东亚大槽明显偏东,引导冷空气南下与强大副高西北侧西南暖湿气流交绥在闽西北上空。925hPa闽西北处于气旋式辐合区内。北支槽、南支槽均非常活跃,东移影响闽西北,是导致该地5月区域性降雨量特多的主要原因。