Clay mineral distributions in the southern Yellow Sea and their significance

To better understand the characteristics of the clay minerals in the southern Yellow Sea, the X-ray quantitative determinations have been carried out for the surface samples obtained from the Yellow Sea. With newly compiled clay mineral synoptic maps, the depositional processes were described for four main clay minerals (illite, chlorite, kaolinite and smectite). The analysis shows that most clay minerals are of terrigenous source with the Huanghe River acting as the major sediment supplier. Besides, the source of muddy sediments in the Yellow Sea was also discussed. As for the central Yellow Sea mud (CYSM), the sediments in its northern part mainly come from the Huanghe River, and those in the rest are of multi-origin. Very similarly, a large amount of sediments in the northern part of the southeastern Yellow Sea Mud (SEYSM) derive from the Keum River and Yeong-san River, while those in the southern part are of multi-origin.

A quasi-synoptic interpretation of water mass distribution and circulation in the western North Pacific II: Circulation

Using the data of conductivity-temperature-depth (CTD) intensive observations conducted during Oct.-Nov. 2005, this study provides the first three-dimension quasi-synoptic description of the circulation in the western North Pacific. Several novel phenomena are revealed, especially in the deep ocean where earlier observations were very sparse. During the observations, the North Equatorial Current (NEC) splits at about 12A degrees N near the sea surface. This bifurcation shifts northward with depth, reaching about 20A degrees N at 1 000 m, and then remains nearly unchanged to as deep as 2 000 m. The Luzon Undercurrent (LUC), emerging below the Kuroshio from about 21A degrees N, intensifies southward, with its upper boundary surfacing around 12A degrees N. From there, part of the LUC separates from the coast, while the rest continues southward to join the Mindanao Current (MC). The MC extends to 2 000 m near the coast, and appears to be closely related to the subsurface cyclonic eddies which overlap low-salinity water from the North Pacific. The Mindanao Undercurrent (MUC), carrying waters from the South Pacific, shifts eastward upon approaching the Mindanao coast and eventually becomes part of the eastward undercurrent between 10A degrees N and 12A degrees N at 130A degrees E. In the upper 2 000 dbar, the total westward transport across 130A degrees E between 7.5A degrees N and 18A degrees N reaches 65.4 Sv (1 Sv = 10(-6) m(3)s(-1)), the northward transport across 18A degrees N from Luzon coast to 130A degrees E is up to 35.0 Sv, and the southward transport across 7.5A degrees N from Mindanao coast to 130A degrees E is 27.9 Sv.

Mean sea surface heights of the South and East China Seas from ocean circulation model and geodetic leveling

The mean sea surface heights (sea surface topography) of the South China, East China, Yellow and Bohai Seas are derived from an ocean general circulation model and surface air pressure. The circulation model covers the global oceans, with fine grid (1/6degrees) covering the East Asian marginal seas and coarse grid (31) covering the rest part of the global oceans. The result shows that the China 1985 National Altitude Datum is 24.7 cm above the me-an sea surface height of the world oceans. The mean sea surface in the coastal ocean adjacent to China is higher in the south than in the north. Intercomparison of the model results with the geodetic leveling measurements at 28 coastal tidal stations shows a standard deviation of 4.8 cm and a fitting coefficient of 95.3%. After correction through linear regression, the standard deviation is reduced to 4.5 cm. This indicates that the accuracy of model results is sufficient for practical application. Based on the model results, the mean sea surface heights for the study area with a resolution of 1/6 degree are given. This result also links the mean sea levels at islands with those on the mainland coast and gives the mean sea surface heights at tidal stations in the Taiwan Island, the Dongsha Islands, the Yisha Islands and the Nansha Islands relative to the China 1985 National Altitude Datum.

Analysis of the sea-level variability along the Chinese coast and estimation of the impact of a CO2-perturbed atmospheric circulation

The relationship between monthly sea-level data measured at stations located along the Chinese coast and concurrent large-scale atmospheric forcing in the period 1960-1990 is examined. It is found that sea-level varies quite coherently along the whole coast, despite the geographical extension of the station set. A canonical correlation analysis between sea-level and sea-level pressure (SLP) indicates that a great part of the sea-level variability can be explained by the action of the wind stress on the ocean surface. The relationship between sea-level and sea-level pressure is analyzed separately for the summer and winter half-years. In winter, one factor affecting sea-level variability at all stations is the SLP contrast between the continent and the Pacific Ocean, hence the intensity of the winter Monsoon circulation. Another factor that affects coherently all stations is the intensity of the zonal circulation at mid-latitudes. In the summer half year, on the other hand, the influence of SLP on sea-level is spatially less coherent: the stations in the Yellow Sea are affected by a more localized circulation anomaly pattern, whereas the rest of the stations is more directly connected to the intensity of the zonal circulation. Based on this analysis, statistical models (different for summer and winter) to hindcast coastal sealevel anomalies from the large-scale SLP field are formulated. These models have been tested by fitting their internal parameters in a test period and reproducing reasonably the sea-level evolution in an independent period. These statistical models are also used to estimate the contribution of the changes of the atmospheric circulation on sea-level along the Chinese coast in an altered climate. For this purpose the ouput of 150 year-long experiment with the coupled ocean-atmosphere model ECHAM1-LSG has been analyzed, in which the atmospheric concentration of greenhouse gases was continuously increased from 1940 until 2090, according to the Scenario A projection of the Intergovermental Panel on Climate Change. In this experiment the meridional (zonal) circulation relevant for sea-level tends to become weaker (stronger) in the winter half year and stronger (weaker) in summer. The estimated contribution of this atmospheric circulation changes to coastal sea-level is of the order of a few centimeters at the end of the integration, being in winter negative in the Yellow Sea and positive in the China Sea with opposite signs in the summer half-year.

Effects of diet and body size on phosphorus utilization of Liza haematocheila T. & S.

A 21-d laboratory experiment was conducted to study, the phosphorus (P) utilization of two different diets by redlip mullet Liza haematocheila T. & S. Sand-filtered water in salinity 30 and temperature 25 degrees C was used. Twenty-nine fish individuals were divided into three groups: 11 to group 1 (G1) fed on diet 1, 11 to group 2 (G2) fed on diet 2, and 7 to contrast group. Diet 1 was a commercial feed, more valuable in nutrition than diet 2 that similar to natural detritus. The results show the intake phosphorus (IP) of G1 was significantly higher than that of G2, and both increased linearly with body size at a certain amount of diet. The retention phosphorus (RP) in fish of G1 was lower than G2. The relationship between retention phosphorus and body size was positive and stronger in G2. Significant difference in faecal phosphorus (FP) was found between G1 and G2. Body size significantly impacted the excretion phosphorus (EP) in G1 but G2. The loss of intake phosphor-us in G1 was 10.83-20.27 mg per g fish weight gain, higher than that in G2 for 6.63-9.56. Of the phosphor-us, about 10% was allocated into growth, 50% in faeces, and the rest lost in excretion. The main part of phosphorus was lost in faeces but excretion. The phosphorus budget of the fish could be described as 100IP = 7.40RP + 47.39FP + 36.63EP (Diet 1) or 100IP = 11.93RP + 56.64FP + 21.76EP (Diet 2).

First record of marine crab, Eucrate alcocki SerSne, in SerSne and Lohavanijaya, 1973 (Crustacea, Decapoda, Brachyura, Euryplacidae) from India

This paper reports the occurrence of a marine brachyuran crab species Eucrate alcocki SerSne, in SerSne et al., 1973, of the family Euryplacidae Stimpson, 1871, for first time from India, based on a male specimen from Parangipettai fish landing centre in Bay of Bengal, Southeast Coast of India. Although morphologically corresponding with what is currently defined as E. alcocki, the color pattern of the carapace of the present specimen is rather different from that of the Chinese material-only the anterior fifth of the carapace is marked with scattered red spots, the rest of the surface is yellowish, with four unusually shaped red blotches which almost look like Sanskrit characters.

Phylogeography of the Qinghai-Tibetan Plateau endemic Juniperus przewalskii (Cupressaceae) inferred from chloroplast DNA sequence variation

The vegetation of the northeast Qinghai-Tibetan Plateau is dominated by alpine meadow and desert-steppe with sparse forests scattered within it. To obtain a better understanding of the phylogeography of one constituent species of the forests in this region, we examined chloroplast trnT-trnF and trnS-trnG sequence variation within Juniperus przewalskii, a key endemic tree species. Sequence data were obtained from 392 trees in 20 populations covering the entire distribution range of the species. Six cpDNA haplotypes were identified. Significant population subdivision was detected (G(ST) = 0.772, N-ST = 0.834), suggesting low levels of recurrent gene flow among populations and significant phylogeographic structure (N-ST > G(ST), P < 0.05). Eight of the nine disjunct populations surveyed on the high-elevation northeast plateau were fixed for a single haplotype (A), while the remaining, more westerly population, contained the same haplotype at high frequency together with two low frequency haplotypes (C and F). In contrast, most populations that occurred at lower altitudes at the plateau edge were fixed or nearly fixed for one of two haplotypes, A or E. However, two plateau edge populations had haplotype compositions different from the rest. In one, four haplotypes (A, B, D and E) were present at approximately equivalent frequencies, which might reflect a larger refugium in the area of this population during the last glacial period. Phylogenetic analysis indicated that the most widely distributed haplotype A is not ancestral to other haplotypes. The contrasting phylogeographic structures of the haplotype-rich plateau edge area and the almost haplotype-uniform plateau platform region indicate that the plateau platform was recolonized by J. przewalskii during the most recent postglacial period. This is supported by the findings of a nested clade analysis, which inferred that postglacial range expansion from the plateau edge followed by recent fragmentation is largely responsible for the present-day spatial distribution of cpDNA haplotypes within the species.

森林资源模型仿真研究

A dynamic model is built up by means of system dynamics method for the development offorest resources in Liaoning province.Prediction and strategy research are made for the fo-rest ecology and timber production in 2000 in Liaoning.The impactions of various designson the development of forest resources are understood by adjusting different parameters andadding disturbances to the model.Various parameters are determined according to the dataprovided by Land Dept.and Forestry Dept.of Liaoning province,the data of provincialforest statistics in 1978 is taken as the original values for running model,which is perfor-med on IBM PC/X7 computer in DYNAMO language of system dynamics.

大洋有孔虫微量元素研究

The most novel aspect of this thesis is the combination analysis of the boron isotopes and trace elements. What’s more, it also provides a reliable analytical technique, which is suitable for both boron isotopes and trace elements. Al/Ca values can be used to monitor the clay removal during the sample preparation. It is found that when Al/Ca>100 mol/mol, the measured boron isotopic compositions are always several permil lower than those properly cleaned. B/Ca ratios can be used to calculate the exact boron loaded for each sample. Otherwise, too much loading will lead to too long time for the whole analytical sequence, and too less loading might incur serious blank problem. One other benefit besides those discussed above is that the combination analysis of boron isotopes and trace elements on the same sample allows reconstruction of the marine carbonate system and atmospheric pCO2 without assumption of the other parameter. In the marine carbonate system, with the seawater pH from the foraminiferal 11B, one has to make an assumption on the other variable to obtain the rest four variables. A series studies found that U/Ca and B/Ca are potential proxies for seawater [CO32-]and [HCO3-], respectively. Since they are measured on the same sample with boron isotopes, hence, there is no spatial or temporal ambiguity in the incorporation of the two controlling parameters. With 11B and U/Ca, the reconstructed atmospheric pCO2 variations match the atmospheric pCO2 record from the Vostok ice core within ±20 ppm. The incorporations of U and B into foraminiferal carbonates are controlled by the overall growth rate of individual foraminifers and other possible factors. The reliable application of these proxies still require further calibrations. In a similar fashion, the combination analysis of boron isotopes and Mg/Ca also has great advantages. Mg/Ca has been proved to be a reliable proxy for the surface seawater temperature. With the combination analysis, one can determine the phase between changes in atmospheric pCO2 and surface seawater temperature, thus explore the cause and mechanism of the changes in atmospheric pCO2. .

太阳辐射对热层和电离层变化性的影响

Solar ultraviolet (UV) radiation at wavelengths less than 400 nm is an important source of energy for aeronomic processes throughout the solar system. Solar UV photons are absorbed in planetary atmospheres, as well as throughout the heliosphere, via photodissociation of molecules, photoionization of molecules and atoms, and photoexcitation toexcitation including resonance scattering. In this paper, the solar irradiances data measured by TIMED SEE, as well as the solar proxies such as F10.7 and Mg II, thermosphere neutral density of CHAMP measurements and topside ionospheric plasmas densities from DMSP, are used to analyze solar irradiance effects on the variabilities of the thermosphere and the ionosphere. First, thermosphere densities near 410 km altitude are analyzed for solar irradiance variability effects during the period 2002-2004. Correlations between the densities and the solar irradiances for different spectral lines and wavelength ranges reveal significantly different characteristics. The density correlates remarkably well with all the selected solar irradiances except the lower chromospheric O I (130.4 nm) emission. Among the chosen solar proxies, the Mg II core-to-wing ratio index, EUV (30-120 nm) and F10.7 show the highest correlations with the density for short-term (< ~27 days) variations. For both long- (> ~27 days) and short-term variations, linear correlation coefficients exhibit a decreasing trend from low latitudes towards high latitudes. The density variability can be effectively modeled (capturing 71% of the variance) using multiple solar irradiance indices, including F10.7, SEUV (the EUV 30-120 nm index), and SFUV (the FUV 120-193 nm index), in which a lag time of 1 day was used for both F10.7 and SEUV, and 5 days for SFUV. In our regression formulation SEUV has the largest contribution to the density variation (40%), with the F10.7 having the next largest contribution (32%) and SFUV accounting for the rest (28%). Furthermore, a pronounced period of about 27.2 days (mean period of the Sun's rotation) is present in both density and solar irradiance data of 2003 and 2004, and a pronounced period of about 54.4 days (doubled period of the solar rotation) is also revealed in 2004. However, soft X-ray and FUV irradiances did not present a pronounced 54.4 day period in 2004, in spite of their high correlation with the densities. The Ap index also shows 54-day periodicities in 2004, and magnetic activity, together with solar irradiance, affects the 54-day variation in density significantly. In addition, NRLMSISE00, DTM-2000 and JB2006 model predictions are compared with density measurements from CHAMP to assess their accuracy, and the results show that these models underestimate the response of the thermosphere to variations induced by solar rotation. Next, the equatorial topside ionospheric plasmas densities Ni are analyzed for solar irradiance variability effects during the period 2002-2005. Linear correlations between Ni and the solar irradiances for different wavelength ranges reveal significantly different characteristics. XUV (0-35 nm) and EUV (115-130 nm) show higher correlation with Ni for the long-term variations, whereas EUV (35-115 nm) show higher correlation for the short-term variations. Moreover, partial correlation analysis shows that the long-term variations of Ni are affected by both XUV (0-35 nm) and EUV (35-115 nm), whereas XUV (0-35 nm) play a more important role; the short-term variations of Ni are mostly affected by EUV (35-115 nm). Furthermore, a pronounced period of about 27 days is present in both Ni and solar irradiance data of 2003 and 2004, and a pronounced period of about 54 days is also revealed in 2004. Finally, prompted by previous studies that have suggested solar EUV radiation as a means of driving the semiannual variation, we investigate the intra-annual variation in thermosphere neutral density near 400 km during 2002-2005. The intra-annual variation, commonly referred to as the ‘semiannual variation’, is characterized by significant latitude structure, hemispheric asymmetries, and inter-annual variability. The magnitude of the maximum yearly difference, from the yearly minimum to the yearly maximum, varies by as much as 60% from year to year, and the phases of the minima and maxima also change by 20-40 days from year to year. Each annual harmonic of the intra-annual variation, namely, annual, semiannual, ter-annual and quatra-annual, exhibits a decreasing trend from 2002 through 2005 that is correlated with the decline in solar activity. In addition, some variations in these harmonics are correlated with geomagnetic activity, as represented by the daily mean value of Kp. Recent empirical models of the thermosphere are found to be deficient in capturing most of the latitude dependencies discovered in our data. In addition, the solar flux and geomagnetic activity proxies that we have employed do not capture some latitude and inter-annual variations detected in our data. It is possible that these variations are partly due to other effects, such as seasonal-latitudinal variations in turbopause altitude (and hence O/N2 composition) and ionosphere coupling processes that remain to be discovered in the context of influencing the intra-annual variations depicted here. Our results provide a new dataset to challenge and validate thermosphere-ionosphere general circulation models that seek to delineate the thermosphere intra-annual variation and to understand the various competing mechanisms that may contribute to its existence and variability. We furthermore suggest that the term “intra-annual” variation be adopted to describe the variability in thermosphere and ionosphere parameters that is well-captured through a superposition of annual, semiannual, ter-annual, and quatra-annual harmonic terms, and that “semiannual’ be used strictly in reference to a pure 6-monthly sinusoidal variation. Moreover, we propose the term “intra-seasonal” to refer to those shorter-term variations that arise as residuals from the above Fourier representation.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

A型、D型人格影响冠心病预后的心理生理机制

Coronary heart disease （CHD）is a common cardiovascular disease in the elderly, is also a typical psychosomatic disease. Personality factors are very important in many psychological factors impacting on the prognosis of patients with CHD. The most influential personality factors to CHD are Type A and Type D personality. The previous research has shown that although Type A personality increased the prevalence of CHD, it cannot predict the development and prognosis after diagnosis. In contradict, Type D personality can predict prognosis. There is still no clinic-based or theory-based answer to the question: Why Type A personality cannot predict the outcome while Type D personality could predict the prognosis independently. The current research conducted a systematic investigation to the above question, which included one comparison study between CHD patients and control group, and four studies on reaction experiment and answered the question: why Type A personality cannot predict whereas Type D personality could effectively predict prognosis of CHD. The findings of the current research were: Type A and Type D personality influence CHD prognosis through different psychological mechanisms: both dimensions of Type D personality have direct influence on social support, whereas neither dimensions of Type A personality related to social support, directly of indirectly. Negative affection component of Type D personality significantly related to anxiety and depression, Social repression significantly related to anxiety but not depression. Both dimensions of Type A personality significantly related to anxiety but not depression. Neither under rest or diaphragmatic breathing conditions, Type A personality had no significant influence on vestibular autonomic reaction among healthy young males. Neither Type A nor Type D personality had significant influence on vestibular autonomic reaction among old CHD patients under rest condition. Type D personality predicted lower sympathetic excitation under rest condition, and lower cardiac vagal tone under diaphragmatic breathing condition among healthy young males. When actively reacted to stimuli (math calculation) under rest condition, Type A personality increased sympathetic excitation among healthy young males. When actively reacted to stimuli (math calculation) under diaphragmatic breathing condition, Type A personality increased cardiac vagal tone among the same group of subjects. When actively reacted to stimuli under neither condition, Type D personality showed no significant influence on vestibular autonomic reaction among young males. When passively reacted to stimuli under neither condition, Type A personality showed no significant influence on vestibular autonomic reaction among young males. When passively stimulated followed rest, Type D personality increased sympathetic excitation and decreased cardiac vagal tone among young males. When passively stimulated followed diaphragmatic breathing, Typed showed no significant influence on vestibular autonomic reaction among young males. The above results indicated that Type A and Type D personalities had different psychological mechanisms to the outcome of CHD treatment: neither dimensions of Type A personality had direct or indirect effects on social support; both dimensions of Type D personality had direct and indirect effects on social support. Negative affection component of Type D personality significantly related to anxiety and depression, Social repression significantly related to anxiety but not depression. Both dimensions of Type A personality significantly related to anxiety but not depression. Social support positively related to the outcome after CHD treatment. The biological mechanisms of Type A and Type B personality to CHD prognosis differed in the following ways: Type A personality increased sympathetic excitation when actively stimulated, but had no influence when passively stimulated among young male subjects. When passively stimulated after rest, Type D personality predicted high sympathetic excitation and low cardiac vagal tone among young males, but not vestibular autonomic reaction among young males. Key words: Type A personality, Type D personality, Coronary Heart Disease (CHD), Prognosis, Psychobiological Mechanisms

相同知识背景下技能之间迁移机制的研究

Since the middle of 1980's, the mechanisms of transfer of training between cognitive subskills rest on the same body of declarative knowledge has been highly concerned. The dominant theory is theory of common element (Singley & Anderson, 1989) which predict that there will be little or no transfer between subskills within the same domain when knowledge is used in different ways, even though the subskills might rest on a common body of declarative knowledge. This idea is termed as "principle of use specificity of knowledge" (Anderson, 1987). Although this principle has gained some empirical evidence from different domains such as elementary geometry (Neves & Anderson, 1981) and computer programming (McKendree & Anderson, 1987), it is challenged by some research (Pennington et al., 1991; 1995) in which substantially larger amounts of transfer of training was found between substills that rest on a shared declarative knowledge but share little procedures (production rules). Pennington et al. (1995) provided evidence that this larger amounts of transfer are due to the elaboration of declarative knowledge. Our research provide a test of these two different explanation, by considering transfer between two subskills within the domain of elementary geometry and elementary algebra respectively, and the inference of learning method ("learning from examples" and "learning from declarative-text") and subject ability (high, middle, low) on the amounts of transfer. Within the domain of elementary geometry, the two subskills of generating proofs" (GP) and "explaining proofs" (EP) which are rest on the declarative knowledge of "theorems on the characters of parallelogram" share little procedures. Within the domain of elementary algebra, the two subskills of "calculation" (C) and "simplification" (S) which are rest on the declarative knowledge of "multiplication of radical" share some more procedures. The results demonstrate that: 1. Within the domain of elementary geometry, although little transfer was found between the two subskills of GP and EP within the total subjects, different results occurred when considering the factor of subject's ability. Within the high level subjects, significant positive transfer was found from EP to GP, while little transfer was found on the opposite direction (i. e. from GP to EP). Within the low level subjects, significant positive transfer was found from EP to GP, while significant negative transfer was found on the opposite direction. For the middle level subject, little transfer was found between the two subskills. 2. Within the domain of elementary algebra, significant positive transfer was found from S to C, while significant negative transfer was found on the opposite direction (i. e. from C to S), when considering the total subjects. The same pattern of transfer occurred within the middle level subjects and low level subject. Within the high level subjects, no transfer was found between the two subskills. 3. Within theses two domains, different learning methods yield little influence on transfer of training between subskills. Apparently, these results can not be attributed to either common procedures or elaboration of declarative knowledge. A kind of synthetic inspection is essential to construct a reasonable explanation of these results which should take into account the following three elements: (1) relations between the procedures of subskills; (2) elaboration of declarative knowledge; (3) elaboration of procedural knowledge. 排Excluding the factor of subject, transfer of training between subskills can be predicted and explained by analyzing the relations between the procedures of two subskills. However, when considering some certain subjects, the explanation of transfer of training between subskills must include subjects' elaboration of declarative knowledge and procedural knowledge, especially the influence of the elaboration on performing the other subskill. The fact that different learning methods yield little influence on transfer of training between subskills can be explained by the fact that these two methods did not effect the level of declarative knowledge. Protocol analysis provided evidence to support these hypothesis. From this research, we conclude that in order to expound the mechanisms of transfer of training between cognitive subskills rest on the same body of declarative knowledge, three elements must be considered synthetically which include: (1) relations between the procedures of subskills; (2) elaboration of declarative knowledge; (3) elaboration of procedural knowledge.