New evidence of microbe origin for ferromanganese nodules from the East Pacific deep sea floor

Using a fluorescence microscope and EPMA, abundant microbe ”bodies“ and clear microbic fluorescent microstructure are determined in the ferromanganese nodules recently collected from the East Pacific deep sea floor. The microbic fluorescent structure shows a close relation to the formation of the ferromanganese nodules. According to their morphological features, the microbes are classified into two types: one is named clumpy microbe, which takes a bar-shaped manganese mineral as a pillar and grows like fasciculate coral, resulting in irregular cauliflorate nodules with rough surfaces; the other is called filamentous microbe, which grows in very thin arcuate and/or concentric circular laminae composed of a microbe layer and a metal (manganese and iron)-rich layer, leading to potato-shaped nodules with relatively smooth surfaces. It also can be seen that the two types of microbes are intergrown together, resulting in nodules complicated in compositions and shapes.     

Biogenesis and binary mineralization of organism and chemism of polymetallic nodules from Pacific Ocean

The biogenic structure and biominerals in the polymetallic nodules reveal their biogenic origin. The experiment of the cation substitution capacity of synthetic 1 nm manganate suggests that chemism participates in the mineralization of the polymetallic nodules. That means the mineralization of the polymetallic nodules is two-component. After discussing the mineralization mechanism, the authors pointed out that the formation of polymetallic modules in deep-sea must meet four prerequisites: oxidative environment, very low sedimentation rate and/or hiatuses, enough nuclei, and suitable for the ultra-microbes-nodule’s builder to live. Project supported by the Deep Sea Geologic Research Project of China Oceanic Mineral Resources Research and Development Association (COMRA)     

Growth history and formation environments of ferromanganese deposits on the Philippine Sea Plate, northwest Pacific Ocean

Abstract   Hydrogenetic ferromanganese crusts are widespread on the floor of the northwestern Pacific Ocean, south and east of the Japanese Islands, despite vigorous tectonic activity, such as subduction and back-arc spreading, since at least the Mid-Paleogene over the Philippine Sea Plate region and nearby. The crusts occur mainly at water depths shallower than 3000 m, but also at greater depths of up to 6000 m. Fine-scale ¹⁰Be/⁹Be dating was undertaken on several 5–10 cm thick hydrogenetic ferromanganese crusts sampled from different geological environments, including inactive submarine volcanoes, tectonic escarpments and abandoned rifts. The results indicate that the crusts have grown at relatively constant rates of 4–7 mm/my without any significant time breaks. These uniform and constant growth rates suggest that the basins have been exposed constantly to oxygenated bottom waters since their formation in the Middle Miocene ( ca 15 Ma) or earlier. Local geological or oceanographic environmental changes might have slowed or increased some of the growth rates resulting in correlation of some internal structures. The Philippine Sea Plate region could have economic potential in areas of thick hydrogenetic ferromanganese crusts over a wide range of water depths.     

Coupled Os–Os enrichments in the Earth’s mantle – core–mantle interaction or recycling of ferromanganese crusts and nodules?

¹⁸⁶Os enrichments in volcanic rocks and peridotite-derived iridosmine grains have been attributed to contributions from Earth’s outer core to the mantle, and apparently constrain the scale of mantle convection and an early timing for inner–outer core segregation more than 3.5 Gyr ago. Here, we highlight that marine ferromanganese crusts and nodules are characterised by high Pt/Os ratios and Pt–Os contents that develop much larger ¹⁸⁶Os excesses over geological time (≥0.2%/Gyr) than those hypothesised for Earth’s outer core (<0.005–0.01%/Gyr). ¹⁸⁷Os/¹⁸⁸Os ratios in ferromanganese crusts are radiogenic due to sequestering of continental Os from seawater. Similarly, ancient ferromanganese materials may have had ¹⁸⁶Os excesses (>0.1%) as a result of high Pt/Os ratios in continental crust, even prior to in-growth of ¹⁸⁶Os after formation due to their high Pt/Os ratios. Past recycling of small amounts of these materials into the Earth’s mantle will produce coupled ¹⁸⁷Os–¹⁸⁶Os excesses and little change in Re and platinum-group-element concentrations, as observed in Hawaiian picrites, and in contrast to the predicted result of outer core addition to the mantle. ¹⁸⁷Os and ¹⁸⁶Os enrichments in the Hawaiian mantle source are potentially consistent with it comprising recycled oceanic lithosphere, pelagic sediments and ferromanganese materials, and questions the notion that Os isotopes can be used to uniquely identify core–mantle interactions and the depth at which mantle sources for volcanism originate.     

The water column distributions of thorium isotopes in the western North Pacific

Four vertical profiles of²³⁰Th and²²⁸Th were determined using large volume water samples in the western North Pacific. An almost linear increase of²³⁰Th with depth was observed for all of the profiles for which the unidirectional first order scavenging model was difficult to explain. We developed a model which included a dissolved-particulate transformation as well as parameters of the scavenging model. Application of the model to the vertical distributions of total and the GEOSECS particulate Th isotopes (²³⁰Th and²³⁴Th) yielded the residence time of dissolved Th with respect to adsorption to particles and the turnover time of particulate Th to be 235 days and 57 days, respectively. The Th isotopes appeared to be carried down the water column by fine particles with a mean settling velocity of 1 m/day which continually release Th into sea water as well as pick up Th from the water along their journey to the bottom.For²²⁸Th, a large excess over²³²Th was observed throughout the water column with pronounced high concentrations in surface and bottom waters, suggesting that the²²⁸Th was derived from²²⁸Ra diffused out of sediments. The vertical distributions of²²⁸Th seemed to be significantly influenced by lateral mixing along isopycnals.     

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…     

Organochlorine and butyltin residues in deep-sea organisms collected from the western North Pacific,off-Tohoku,Japan

Organochlorine (OCs) and butyltin (BTs) residues were determined in deep-sea organisms collected from the western North Pacific, off-Tohoku, Japan. Among OCs, concentrations of polychlorinated biphenyls (PCBs) and DDTs (DDTs and its metabolites) were the highest in deep-sea organisms (maximum concentrations of 6700 and 13,000 ng/g lipid wt, respectively). Chlordane compounds (CHLs) were the next most abundant OCs, and hexachlorocyclohexanes (HCHs) and hexachlorobenzene (HCB) were the lowest. BTs were also detected at maximum concentrations of 570 ng/g wet wt. Concentrations of PCBs, CHLs and BTs in deep-sea organisms collected from the western North Pacific, off-Tohoku, were generally lower than those in deep-sea and shallow water organisms from Japanese coastal waters. On the other hand, considerable variations in the concentrations of OCs were found among deep-sea organisms analyzed. Several carnivorous fishes such as snubnosed eel, lanternshark and grenadiers accumulated some OCs such as PCBs, DDTs and CHLs at high concentrations of up to a few μg/g levels. In addition, the residue pattern of OCs and BTs in fishes showed a specific trend according to the sampling depth; higher concentrations of PCBs, DDTs and CHLs and lower concentrations of HCHs, HCB and BTs were found in fishes collected from greater depth (1000 m) compared to those from shallower waters. This trend is consistent with the results of our earlier study on mesopelagic myctophid fishes. Results of this study suggest vertical transport of hydrophobic OCs such as PCBs, DDTs and CHLs and its accumulation in benthic deep-sea organisms.     

An evaluation and improvement of tropical cyclone prediction in the western North Pacific basin from global ensemble forecasts

Forecasts of tropical cyclones(TCs) of the western North Pacific basin during the period of July to August 2018,especially of Rumbia(2018), Ampil(2018) and Jongdari(2018) that made landfall over Shanghai, have opposed great challenges for numerical models and forecasters. The predictive skill of these TCs are analyzed based on ensemble forecasts of ECMWF and NCEP. Results of the overall performance show that ensemble forecasts of ECMWF generally have higher predictive skill of track and intensity forecasts than those of NCEP. Specifically, ensemble forecasts of ECMWF have higher predictive skill of intensity forecasts for Rumbia(2018) and Ampil(2018) than those of NCEP, and both have low predictive skill of intensity forecasts for Jongdari(2018) at peak intensity. To improve the predictive skill of ensemble forecasts for TCs, a method that estimates adaptive weights for members of an ensemble forecast is proposed. The adaptive weights are estimated based on the fit of ensemble priors and posteriors to observations. The performances of ensemble forecasts of ECMWF and NCEP using the adaptive weights are generally improved for track and intensity forecasts. The advantages of the adaptive weights are more prominent for ensemble forecasts of ECMWF than for those of NCEP.     

夏季西北太平洋大气环流异常及其与热带印度洋—太平洋海温变化的关系

本文分析了夏季西北太平洋大气环流异常特征及其与海温变化的关系,发现夏季西北太平洋异常反气旋/气旋(WNPAC/WNPC)是西北太平洋地区对流层中低层存在的重要大气环流异常现象,与东亚—西北太平洋低纬度至高纬度的经向PJ波列及欧亚中高纬度东西纬向波列的变化有关,通过与中高纬度环流变化的联系,对东亚及欧亚中高纬度气候有重要影响.夏季WNPAC/WNPC与热带海温变化的关系存在明显的不对称性,显著的WNPAC一般出现在El Niño衰减年夏季,与前期El Niño成熟年冬季的赤道东太平洋暖海温异常和El Niño衰减年春夏季印度洋海盆尺度的暖海温异常有明显的正相关关系,进一步表明了WNPAC在El Niño事件影响夏季气候中的重要桥梁作用;而夏季显著的WNPC与前期和同期热带海温变化的关系存在明显的不确定性,主要与夏季热带印度洋和赤道中东太平洋之间东暖西冷的热力差异异常引起的孟加拉湾—赤道西太平洋西风异常有关.进一步分析WNPAC/WNPC与海温变化关系不对称的可能原因,发现El Niño和La Niña衰减年夏季热带印度洋和太平洋海温变化所引起的印—太之间海温(热力)差异的一致性特征可能是导致WNPAC/WNPC与海温变化关系不对称的主要原因.     

Interdecadal change of the linkage between the North Atlantic Oscillation and the tropical cyclone frequency over the western North Pacific

The relationship between the North Atlantic Oscillation(NAO) and the tropical cyclone frequency over the western North Pacific(WNPTCF) in summer is investigated by use of observation data. It is found that their linkage appears to have an interdecadal change from weak connection to strong connection. During the period of 1948–1977, the NAO was insignificantly correlated to the WNPTCF. However, during the period of 1980–2009, they were significantly correlated with stronger(weaker) NAO corresponding to more(fewer) tropical cyclones in the western North Pacific. The possible reason for such a different relationship between the NAO and the WNPTCF during the former and latter periods is further analyzed from the perspective of large-scale atmospheric circulations. When the NAO was stronger than normal in the latter period, an anomalous cyclonic circulation prevailed in the lower troposphere of the western North Pacific and the monsoon trough was intensified, concurrent with the eastward-shifting western Pacific subtropical high as well as anomalous low-level convergence and high-level divergence over the western North Pacific. These conditions favor the genesis and development of tropical cyclones, and thus more tropical cyclones appeared over the western North Pacific. In contrast, in the former period, the impact of the NAO on the aforementioned atmospheric circulations became insignificant, thereby weakening its linkage to the WNPTCF. Further study shows that the change of the wave activity flux associated with the NAO during the former and latter periods may account for such an interdecadal shift of the NAO–WNPTCF relationship.     

Relationship between the Asian-Pacific oscillation and the tropical cyclone frequency in the western North Pacific

The relationship between the Asian-Pacific oscillation (APO) and the tropical cyclone frequency over the western North Pacific (WNP) in summer is preliminarily investigated through an analysis of ob- served data. The result has shown clearly that APO is significantly and positively correlated to the tropical cyclone frequency in the WNP. If APO is above (below) the normal in summer, more (less) tropical cyclones will tend to appear in the WNP. The present study also addresses the large-scale at- mospheric general circulation changes underlying the linkage between APO and the WNP tropical cy- clone frequency. It follows that a positive phase of summer APO is concurrent with weakened as well as northward and eastward located western Pacific subtropical high (WPSH), low-level convergence and high-level divergence, and reduced vertical zonal wind shear in the WNP, providing favorable envi- ronment for the tropical cyclone genesis, and thus more tropical cyclones will come into being, and vice versa.     

The relationship between the growth process of the ferromanganese crusts in the pacific seamount and Cenozoic ocean evolvement

Base on the Os isotope stratigraphy together with the empirical growth rate models using Co concentrations, the growth ages of the ferromanganese crusts MHD79 and MP3D10 distributed in the seamount of Pacific are confirmed. Through the contrast and research on the previous achievements including ODP Leg 144 and the crusts CD29-2, N5E-06 and N1–15 of the seamount of the Central Pacific, the uniform five growth and growth hiatus periods of them are found, and closely related to the Cenozoic ocean evolvement process. In the Paleocene Carbon Isotope Maximum (PCIM), the rise of the global ocean productivity promoted the growth of the seamount crust; the first growth hiatus (I) of the ferromanganese crust finished. In the Paleocene-Eocene Thermal Maximum (PETM), though the vertical exchange of seawater was weakened, the strong terrestrial chemical weathering led to the input of a great amount of the terrigenous nutrients, which made the bioproductivity rise, so there were no crust hiatuses. During 52–50 Ma, the Early Eocene Optimum Climate (EECO), the two poles were warm, the latitudinal temperature gradient was small, the wind-driven sea circulation and upwelling activity were weak, the terrestrial weathering was also weakened, the open ocean bioproductivity decreased, and the ferromanganese crust had growth hiatus again (II). From early Middle Eocene-Late Eocene, Oligocene, it was a long-term gradually cooling process, the strengthening of the sea circulation and upwelling led to a rise of bioproductivity, and increase of the content of the hydrogenous element Fe, Mn and Co and the biogenous element Cu, Zn, so that was the most favorable stage for the growth of ferromanganese crust (growth periods III and IV) in the studied area. The hiatus III corresponded with the Eocene- Oligocene boundary, is inferred to relate with the global climate transformation, celestial body impact event in the Eocene-Oligocene transition. From the early to the middle Miocene, a large-scale growth hiatus (hiatus period IV) of the ferromanganese crust in the studied area is inferred to relate with temporary warm up climate and ephemeral withdrawal of Antarctic bottom water in the early Miocene. After that, the Antarctic ice sheets extended, the bottom water circumfluence strengthened, the ocean fertility increased, and the once interrupted crust continued to grow in the late Miocene (growth period V). Supported by China Ocean Mineral Resources Research and Development Association “10th Five Year” Topic (Grant No. DY105-01-04-14)     

Carbon isotope stratigraphy of Torinosu‐type limestone in the western Paleo‐Pacific and its implication to paleoceanography in the Late Jurassic and earliest Cretaceous

Carbon isotope stratigraphy of the Late Jurassic and earliest Cretaceous was revealed from Torinosu‐type limestone, which was deposited in a shallow‐marine setting in the western Paleo‐Pacific, in Japan. Two sections were examined; the Nakanosawa section of the late Kimmeridgian to early Tithonian age (Fukushima Prefecture, Northeast Japan), and the Furuichi section of the late Kimmeridgian to early Berriasian age (Ehime Prefecture, Southwest Japan). The age‐model was established using Sr isotope ratio and fossil occurrence. The limestone samples have a low Mn/Sr ratio (mostly <0.5) and lack a distinct correlation between δ¹³C and δ¹⁸O, indicating a low degree of diagenetic alteration. Our composite δ¹³C profile from the two limestone sections shows three stratigraphic correlation points that can be correlated with the profiles of relevant ages from the Alpine Tethyan region: a large‐amplitude fluctuation (the lower upper Kimmeridgian, ～152 Ma), a positive anomaly (above the Kimmeridgian/Tithonian boundary, ～150 Ma), and a negative anomaly (the upper lower Tithonian, ～148 Ma). In addition, we found that δ¹³C values of the Torinosu‐type limestone are ～1‰ lower than the Tethyan values in the late Kimmeridgian. This inter‐regional difference in δ¹³C values is likely to have resulted from a higher productivity and/or an organic burial in the Tethyan region. The difference gradually reduces and disappears in the late Tithonian, where the Tethyan and our δ¹³C records show similar stable values of 1.5–2.0‰. This isotopic homogenization is probably due to changes in the continental distribution and the global ocean circulation, which propagated the ¹³C‐depleted signature from the larger Paleo‐Pacific to the smaller Tethys Ocean during this time.     

Stable isotope paleohydrology and chemostratigraphy of the Albian Wayan Formation from the wedge-top depozone,North American Western Interior Basin

Understanding of the role of atmospheric moisture and heat transport in the climate system of the Cretaceous greenhouse world represents a major challenge in Earth system science.Stable isotopic paleohydrologic data from mid-Cretaceous paleosols in North America,from paleoequatorial to paleoArctic latitudes,have been used to constrain the oxygen isotope mass balance of the Albian hydrologic cycle.Over the range from 40°-50°N paleolatitude,sideritic paleosols predominate,indicating paleoenvironments with positive precipitation-evaporation(P-E)balances.Local exceptions occur on leeward side of the Sevier Orogen,where calcic paleosols in the wedge-top depozone record paleoenvironments with negative P-E balances in the orographic rain shadow.Stratigraphic sections in the Wayan Formation of Idaho(WF)were sampled from the wedge-top depozone.The units consist of stacked m-scale mudstone paleosols separated by m-scale sandstone-siltstone beds.Sections were sampled for organic carbon isotope profiles,and B-horizons from 6 well-developed paleosols were sampled for detrital zircons to determine maximum depositional ages.The first of these from the WF has produced a U-Pb concordia age of 101.0±1.1 Ma.This same WF section has produced a stratigraphic trend of upwardly decreasingδ~(13)C values ranging from-24‰upwards to-27‰VPDB,suggesting correlation to the late Albian C15 C-isotope segment.Pedogenic carbonates from the WF principally consist of micritic calcite,with carbon-oxygen isotope values that array along meteoric calcite lines(MCLs)withδ~(18)O values that range from-9.47‰up to-8.39‰VPDB.At approximately 42°N paleolatitude,these MCL values produce calculated paleoprecipitation values of-8.12‰to-7.04‰VSMOW,a range that is consistent with the estimates produced from other proxies at the same paleolatitudes across North America.These results indicate that despite the orographic rain shadow effect,the processes of meridional atmospheric moisture transport in this locale were similar to those in more humid mid-latitude paleoenvironments elsewhere in the continent.     

New predictors and a new prediction model for the typhoon frequency over western North Pacific

In this paper, the impacts of the atmospheric circulation during boreal winter-spring on the western North Pacific (WNP) typhoon frequency (WNPTF) are studied. Several new factors in winter-spring in- fluencing the typhoon frequency were identified, including the sea ice cover in the North Pacific and the North Pacific oscillation. Based on these results, the multi-linear regression was applied to establishing a new forecast model for the typhoon frequency by using the datasets of 1965―1999. The forecast model shows a high correlation coefficient (0.79) between the model simulated and the actual typhoon frequencies in the period of 1965―1999. The forecast model also exhibits reasonable hindcasts for the typhoon frequencies for the years 2000―2006. Therefore, this work demonstrates that the new pre- dictors are significant for the prediction of the interannual variability of the WNPTF, which could be potentially used in the operational seasonal forecast of the typhoon frequency in the WNP to get a more physically based operational prediction model and higher forecast skill.     

Paleoenvironment evolution of the East Philippine Sea recorded in the new-type ferromanganese crust since the terminal Late Miocene

From systemic research of microstructure, geochemistry, uranium-series and 10Be isotope dating on a new-type deepwater ferromanganese crust from the East Philippine Sea, the paleoenvironment evolu-tion of the target area since the terminal Late Miocene was recovered. The vertical section changes of microstructure and chemical composition are consistent in the studied crust, which indicate three major accretion periods and corresponding paleoenvironment evolution of the crust. The bottom crust zone was formed in the terminal Late Miocene (5.6 Ma) with loose microstructure, high detritus content and high growth rate. Reductions of mineral element content, accretion rate and positive Ce-anomaly degree at 4.6 Ma indicate temporal warming, which went against the crust accretion and finally formed an accretion gap in the terminal Middle Pliocene (2.8―2.7 Ma). The more active Antarctic bottom sea-waters in the Late Pliocene (2.7 Ma) facilitated the fast transfer to the top pure crust zone. Hereafter, with the further apart of volcanic source and the keeping increase of eolian material (1.0 Ma), although surrounding conditions were still favorable, mineral element content still shows an obvious reducing trend. It thereby offers new carrier and data for the unclear paleoceanographic research of the target area since the terminal Late Miocene.     

Organohalogen compounds in deep-sea fishes from the western North Pacific, off-Tohoku, Japan: Contamination status and bioaccumulation profiles

Twelve species of deep-sea fishes collected in 2005 from the western North Pacific, off-Tohoku, Japan were analyzed for organohalogen compounds. Among the compounds analyzed, concentrations of DDTs and PCBs (up to 23,000 and 12,400 ng/g lipid wt, respectively) were the highest. The present study is the foremost to report the occurrence of brominated flame retardants such as PBDEs and HBCDs in deep-sea organisms from the North Pacific region. Significant positive correlations found between δ¹⁵N (‰) and PCBs, DDTs and PBDEs suggest the high biomagnification potential of these contaminants in food web. The large variation in δ¹³C (‰) values observed between the species indicate multiple sources of carbon in the food web and specific accumulation of hydrophobic organohalogen compounds in benthic dwelling carnivore species like snubnosed eel. The results obtained in this study highlight the usefulness of deep-sea fishes as sentinel species to monitor the deep-sea environment.     

春季Hadley环流异常对夏季西北太平洋热带气旋频数影响的数值模拟试验

观测事实揭示出春季Hadley环流与夏季西北太平洋热带气旋频数之间存在显著的负相关关系.由春季Hadley环流异常引起的西北太平洋地区夏季纬向风垂直切变、大气辐合辐散等的异常变化是这一关系存在的内在原因.本文通过数值试验对这一关系的真实性进行了验证,即利用中国科学院大气物理研究所发展的9层大气环流模式（IAP9L-AGCM）模拟了春季Hadley环流异常偏强情景,并分析了该情景下影响西北太平洋热带气旋生成的环境场的响应.结果表明,在春季Hadley环流偏强情景下,夏季西北太平洋地区纬向风垂直切变幅度加大,低空大气异常辐散,高空大气异常辐合,东亚夏季风减弱,这种环流背景不利于热带气旋生成和发展,因此,西北太平洋热带气旋频数异常偏少.数值模拟结果与已有的诊断结果相吻合,进而证实了春季Hadley环流与夏季西北太平洋热带气旋频数负相关关系的存在.因此,春季Hadley环流信号可以用于西北太平洋热带气旋活动的气候预测.     

Tropical cyclones and multiscale climate variability: The active western North Pacific Typhoon season of 2018

The 2018 typhoon season in the western North Pacific(WNP) was highly active, with 26 named tropical cyclones(TCs) from June to November, which exceeded the climatological mean(22) and was the second busiest season over the past twenty years. More TCs formed in the eastern region of the WNP and the northern region of the South China Sea(SCS). More TCs took the northeast quadrant in the WNP, recurving from northwestward to northward and causing heavy damages in China's Mainland(69.73 billion yuan) in 2018. Multiscale climate variability is conducive to an active season via an enhanced monsoon trough and a weakened subtropical high in the WNP. The large-scale backgrounds in 2018 showed a favorable environment for TCs established by a developing central Pacific(CP) El Ni?o and positive Pacific meridional mode(PMM)episode on interannual timescales. The tropical central Pacific(TCP) SST forcing exhibits primary control on TCs in the WNP and large-scale circulations, which are insensitive to the PMM. During CP El Ni?o years, anomalous convection associated with the TCP warming leads to significantly increased anomalous cyclonic circulation in the WNP because of a Gill-type Rossby wave response. As a result, the weakened subtropical high and enhanced monsoon trough shift eastward and northward, which favor TC genesis and development. Although such increased TC activity in 2018 might be slightly suppressed by interdecadal climate variability, it was mostly attributed to the favorable interannual background. In addition, high-frequency climate signals,such as intraseasonal oscillations(ISOs) and synoptic-scale disturbances(SSDs), interacted with the enhanced monsoon trough and strongly modulated regional TC genesis and development in 2018.     

Evolution mechanism of the western Pacific subtropical high

Ever since Charney et al.[1] studied the multiple equilibrium states in atmosphere with highly truncated spectral method in 1979, many Chinese researchers, such as Li Maicun et al. (1983)[2], Liu Chongjian et al. (1983)[3], Miu Jinhai et al. (1985)[4] and…