Effects of convection instability due to incompatibility between ocean dynamics and surface forcings

The study demonstrates that an incompatibility between a surface temperature climatology and a given ocean model, into which the climatology is assimilated via Haney restoration, can cause model ocean climate drift and interdecadal oscillations when the ocean is switched to a weaker restoration. This is made using an idealized Atlantic Ocean model driven by thermal and wind forcing only. Initially, the temperature climatology is forcefully assimilated into the model, and an implied heat flux field is diagnosed. During this stage any incompatibility is suppressed. The restoring boundary condition is then switched to a new forcing consisting of a part of the diagnosed flux and a part of the restoring forcing in such a way that at the moment of the switching the heat flux is identical to that prior to the switching. Under this new forcing condition, the incompatibility becomes manifest, causing changes in convection patterns, and producing drift and interdecadal oscillations. The mechanisms are described.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

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.     

A Modeling Study of the Effects of Anomalous Snow Cover over the Tibetan Plateau upon the South Asian Summer Monsoon

The effect of anomalous snow cover over the Tibetan Plateau upon the South Asian summer monsoon is investigated by numerical simulations using the NCAR regional climate model (RegCM2) into which gravity wave drag has been introduced. The simulations adopt relatively realistic snow mass forcings based on Scanning Multi-channel Microwave Radiometer (SNINIR) pentad snow depth data. The physical mechanism and spatial structure of the sensitivity of the South Asian early summer monsoon to snow cover anomaly over the Tibetan Plateau are revealed. The main results are summarized as follows. The heavier than normal snow cover over the Plateau can obviously reduce the shortwave radiation absorbed by surface through the albedo effect, which is compensated by weaker upward sensible heat flux associated with colder surface temperature, whereas the effects of snow melting and evaporation are relatively smaller.The anomalies of surface heat fluxes can last until June and become unobvious in July. The decrease of the Plateau surface temperature caused by heavier snow cover reaches its maximum value from late April to early May. The atmospheric cooling in the mid-upper troposphere over the Plateau and its surrounding areas is most obvious in May and can keep a fairly strong intensity in June. In contrast, there is warming to the south of the Plateau in the mid-lower troposphere from April to June with a maximum value in May.The heavier snow cover over the Plateau can reduce the intensity of the South Asian summer monsoon and rainfall to some extent, but this influence is only obvious in early summer and almost disappears in later stages.     

近50年中国地面气候变化基本特征

采用国家基准气候站和基本气象站的地面资料,系统地分析了中国大陆地区1951年以来近地表主要气候要素演化的时间和空间特征。结果表明,中国近50 a来年平均地表气温变暖幅度约为1.1℃,增温速率接近0.22℃/(10 a),比全球或半球同期平均增温速率明显偏高。地表气温增暖主要发生在最近的20余年,其季节和空间特征与前人分析结论基本一致。降水量变化趋势对所取时间段和区域范围敏感。1951年以来全国平均降水量变化趋势不明显,但1956年以来略有增加。降水变化的空间特征明显而相对稳定,东北北部、包括长江中下游的东南部地区和西部广大地区降水增加,而华北地区以及东北东南部和西北东部地区降水明显减少。分析还发现,近50a来全国平均的日照时数、平均风速、水面蒸发等气候要素均呈显著下降趋势,但积雪地带的最大积雪深度却有所增加。中国日照时间和水面蒸发量变化的空间特征很相似,减少最明显的地区均发生在华北和华东,新疆次之。影响中国年代以上尺度气候变化的因子错综复杂,人类活动引起的大气中温室气体浓度增高可能在一定程度上影响了中国近50 a来的气候,但考虑到尚存的不确定性,目前仍不能给出明确结论。中国东部大部分地区日照时间和水面蒸发量减少可能均起源于人为排放的气溶胶影响,平均风速减弱也有利于水面蒸发量下降,而在西部地区云量和降水量的变化可能更重要。     

The eastward subtropical countercurrent on isopycnal surface in the western North Pacific

The classic Sverdrup theory suggests that the water movement in the central subtropical gyre of North Pa-cific be slowly westward or southwestward.In the late sixties of the20th century,the existence of a peculiar eastward narrow flow between20°N and25°N in spring was theoretically predicted.It was named the Subtropical Countercurrent(STCC),although direct observational evidences were not yet sufficient to con-firm whether or not such eastward flow between20°N and25°N was a persistent…     

Geochemical and Sr-Nd isotopic study of alkaline syenites in Liangtun-Kuangdonggou, Liaoning Province, China： Evidence for enriched mantle before 1.86 Ga and implications

The 1.86 Ga Liangtun-Kuangdonggou complex (LKC) is one of the oldest alkaline syenite bodies so far discovered in China. This syenite suite has elevated contents of total alkali (K2O Na2O), with an average of 10.50%, and a mean Rittmann Index (σ) of 6.48. The intrusions have slightly higher concentrations of K2O than those of Na2O on a weight percent basis, indicating the rocks belong to potassium-rich alkaline syenite series. Total rare-earth element concentrations (∑REE ) of the rocks are relatively high, ranging from 324×10 -6 to 1314×10 -6, with a mean value of 666×10 -6. The REE patterns are subparallel and rightward steep with (La/Yb)N >33, showing mild negative to positive Eu anomalies (δEu: 0.63-1.15). All samples exhibit strong LILE and LREE enrichments and TNT (Nb, Ta, Ti) and P depletions in multi-element spidergrams. On the εSr(t)-εNd(t) correlation diagram, most analytical data points plot within the enriched mantle field with low ( 87Sr/86Sr)i ratios (0.7045-0.7051) and negative εNd(t) values (-3.72--3.97), falling among those kimberlites from Fuxian County, Liaoning Provinve, from Mengyin County, Shandong Province and the Ⅱ-type kimberlites from South Africa. These characteristics imply that the LKC-rocks may have the same source as the above-mentioned kimberlites, i.e., they have close connections to the materials derived from enriched mantle reservoirs, further revealing that the upper mantle beneath the northeastern part of the North China Plate had been highly enriched before 1.86 Ga. Geodynamically, the LKC-rocks were formed in a within-plate environment with close genetic connections to rift-related alkaline magma activities possibly controlled by ancient mantle plumes.     

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.     

Target setting and allocation of carbon emissions abatement based on DEA and closest target: an application to 20 APEC economies

Regarded as an effective method for treating the global warming problem, carbon emissions abatement (CEA) allocation has become a hot research topic and has drawn great attention recently. However, the traditional CEA allocation methods generally set efficient targets for the decision-making units (DMUs) using the farthest targets, which neglects the DMUs’ unwillingness to maximize (minimize) some of their inputs (outputs). In addition, the total CEA level is usually subjectively determined without any consideration of the current carbon emission situations of the DMUs. To surmount these deficiencies, we incorporate data envelopment analysis and its closest target technique into the CEA allocation problem. Firstly, a two-stage approach is proposed for setting the optimal total CEA level for the DMUs. Then, another two-stage approach is given for allocating the identified optimal total CEA among the DMUs. Our approach provides more flexibility when setting new input and output targets for the DMUs in CEA allocation. Finally, the proposed approaches are applied for CEA target setting and allocation for 20 Asia-Pacific Economic Cooperation economies.