COMMENTS ON THE NIPALS ALGORITHM

The Non-linear lterative Partial Least Squares(NIPALS)algorithm is used in principal componentanalysis to decompose a data matrix into score vectors and eigenvectors(loading vectors)plus a residualmatrix.N1PALS starts with some guessed starting vector.The principal components obtained by NIPALSdepends on the starting vector;the first principal component could not always be computed.Wold hassuggested a starting vector for NIPALS,but we have found that even if this starting vector is used,thefirst principal component cannot be obtained in all cases.The reason why such a situation occurs isexplained by the power method.A simple modification of the original NIPALS procedure to avoid gettingsmaller eigenvalues is presented.     

Development of a skill assessment tool for the Korea operational oceanographic system

A standard skill assessment(SA) tool was developed and implemented to evaluate the performance of operational forecast models in the Korea operational oceanographic system.The SA tool provided a robust way to assess model skill in the system by comparing predictions and observations,and involved the computation of multiple skill metrics including correlation and error skills.User-and system-based acceptance criteria of skill metrics were applied to determine whether predictions were acceptable for the system.To achieve this,the tool produced a time series comparison plot,a skill score table,and an advanced summarized diagram to effectively demonstrate the multiple skill scores.Moreover,the SA was conducted to evaluate both atmospheric and hydrodynamic forecast variables.For the atmospheric variables,acceptable error criteria were preferable to acceptable correlation criteria over short timescales,since the mean square error overwhelmed the observation variance.Conversely,for the hydrodynamic variables,acceptable root mean square percentage error(e.g.,p erms) criteria were preferable to acceptable error(e.g.,erms) criteria owing to the spatially variable tidal intensity around the Korean Peninsula.Furthermore,the SA indicated that predetermined acceptance error criteria were appropriate to satisfy a target central frequency(fc) for which errors fell within the specified limits(i.e.,the fc equals 70%).     

Global parameter estimation of the Cochlodinium polykrikoides model using bioassay data

Cochlodinium polykrikoides is a notoriously harmful algal species that inflicts severe damage on the aquacultures of the coastal seas of Korea and Japan. Information on their expected movement tracks and boundaries of influence is very useful and important for the effective establishment of a reduction plan. In general, the information is supported by a red-tide(a.k.a algal bloom) model. The performance of the model is highly dependent on the accuracy of parameters, which are the coefficients of functions approximating the biological growth and loss patterns of the C. polykrikoides. These parameters have been estimated using the bioassay data composed of growth-limiting factor and net growth rate value pairs. In the case of the C. polykrikoides, the parameters are different from each other in accordance with the used data because the bioassay data are sufficient compared to the other algal species. The parameters estimated by one specific dataset can be viewed as locally-optimized because they are adjusted only by that dataset. In cases where the other one data set is used, the estimation error might be considerable. In this study, the parameters are estimated by all available data sets without the use of only one specific data set and thus can be considered globally optimized. The cost function for the optimization is defined as the integrated mean squared estimation error, i.e., the difference between the values of the experimental and estimated rates. Based on quantitative error analysis, the root-mean squared errors of the global parameters show smaller values, approximately 25%–50%, than the values of the local parameters. In addition, bias is removed completely in the case of the globally estimated parameters. The parameter sets can be used as the reference default values of a red-tide model because they are optimal and representative. However, additional tuning of the parameters using the in-situ monitoring data is highly required.As opposed to the bioassay data, it is necessary because the bioassay data have limitations in terms of the in-situ coastal conditions.     

Modelling solute leaching during fingered flow by integrating and expanding various theoretical and empirical concepts

Abstract

Wetting front instability (fingered flow) accelerates solute transport through the unsaturated zone to the groundwater table. Whether fingers widen or dissipate close to the groundwater is unclear. Water flow in a two-dimensional artificial capillary fringe below a dry layer exhibiting fingered flow was investigated. The flow diverged strongly in the wet soil, suggesting that fingers dissipate. Expressions for the finger radius in dry and wet soil were combined and adapted to a soil hydraulic property parameterization popular in numerical modelling. The modified equation provided finger radii for soils in humid and arid climates. The fingers in the arid soil were excessively wide. The finger radii were used to model solute transport, assuming fingers dissipated in the subsoil. Modelling was cumbersome for the arid climate. One shower may often be insufficient to trigger fingering in arid regions with short, heavy showers. In soils with shallow groundwater, the diverging subsoil flow determines solute leaching.     

Metamorphic and Tectonic Evolution of the Hwacheon Granulite Complex, Central Korea: Composite P-T Path Resulting from Two Distinct Crustal-Thickening Events

The Hwacheon granulite complex (HGC), occupying the northeasternmargin of the Gyeonggi massif, consists mainly of garnetiferousleucocratic gneiss and leucogranite together with minor kyanite–garnetgneiss, aluminous gneiss, mafic granulite and garnet amphibolite.Mineral assemblages and reaction textures in various rock typesof the HGC document five distinct metamorphic stages: pre- (M₁)and peak (M₂) granulite-facies metamorphism; lower temperature,high (M₃) and low (M₄) pressure upper amphibolite-facies metamorphism;and local retrogression (M₅) producing andalusite-bearing assemblages.Each metamorphic stage can be integrated to give a compositeP–T path consisting of two distinct trajectories, characterizedby clockwise P–T loops at relatively high and low temperatures,respectively. The first P–T trajectory (M₁–M₃) correspondsto a Palaeoproterozoic tectonometamorphic event responsiblefor the formation of the granulite complex at     

RATES OF SLOPE DECLINE,TALUS GROWTH AND CLIFF RETREAT ALONG THE SHOMYO RIVER IN CENTRAL JAPAN: A SPACE‐TIME SUBSTITUTION APPROACH

On the basis of a space–time substitution, rates of cliff retreat and talus development overthe pasttens ofthousands of years were quantified for valley‐side slopes along the Shomyo River with a known rate of waterfall recession. Detailed profiles of the valley‐side slopes were obtained at 19 sites along the river by map reading and field survey, and the morphologic characteristics of the slopes were then measured. By combining the locations of the slopes with the recession rate of Shomyo Falls, it is estimated that the growth rate of the talus slope is 3.9–7.2 mm/yr, the rate of decline of the cliff is 2.0–4.0 × 10^?4 deg/yr, and the retreat rate of the cliff top is 6.2–11.6 mm/yr.     

Tectonometamorphic evolution of the Chuncheon amphibolite, central Gyeonggi massif, South Korea

Abstract The Chuncheon amphibolite, part of the Gubongsan Group which overlies the Yongduri gneiss complex, is interlayered with calc-silicate rock, marble, quartzite, biotite schist and quartzofeldspathic gneiss in the central Gyeonggi massif, South Korea. Metamorphic pressures and temperatures estimated from the amphibolite are 5.5–10.6 kbar and 615–714°C. These P—T conditions are close to those defined by the reaction curve between kyanite and sillimanite, and suggest medium-pressure-type metamorphism of the Chuncheon amphibolite. For two metapelites intercalated with the amphibolite, temperatures are estimated to be 607–699° C, consistent with those obtained from the amphibolite. On the other hand, pressures estimated from these metapelites are significantly different, 4–6 kbar and 9–13 kbar, when rim and core compositions of garnet are, respectively, used. These P—T estimates obtained from the amphibolite and metapelite suggest a nearly isothermal decompression of 3–7 kbar during denudation. Rapid decompression is likely on the basis of the results of mineral chemistry, phase equilibria and geothermobarometer. Moreover, in conjunction with the occurrence of kyanite in the adjacent Gyeonggi gneiss complex, P—T estimates of the Chuncheon amphibolite and metapelite suggest a clockwise P—T—t path. This evolutionary path may be related to the amalgamation of continents during the late Proterozoic event which corresponds to the Jinningian orogeny in the Qinling belt of China.     

Transition from the Zeolite to Prehnite-Pumpellyite Facies in the Karmutsen Metabasites, Vancouver Island, British Columbia

The upper Triassic Karmutsen metabasites from northeast VancouverIsland, B.C., are thermally metamorphosed by the intrusion ofthe Coast Range Batholith. The amygdaloidal metabasites developedin the outer portion of the contact aureole show a progressivemetamorphism from zeolite to prehnite-pumpellyite facies. Thesize of an equilibrium domain is extremely small for these metabasites,and the individual amygdule assemblages are assumed to be inequilibrium. Two major calcite-free assemblages (+chlorite+quartz)are characteristic: (i) laumontite+pumpellyite+epidote in thezeolite facies and (ii) prehnite+pumpellyite+epidote in theprehnite-pumpellyite facies. The assemblages and compositionsof Ca-Al silicates are chemographically and theoretically interpretedon the basis of the predicted P-T grid for the model basalticsystem, CaO-MgO-A1₂O₃-Fe₂O₃-SiO₂-H₂O. The results indicate:(1) local equilibrium has been approached in mineral assemblagesand compositions; (2) the X_Fe³⁺ values in the coexisting Ca-Alsilicates decrease from epidote, through pumpellyite to prehnite;(3) with increasing metamorphic grade, the Fe³⁺ contents ofepidotes in reaction assemblages decrease in the zeolite facies,then increase in the prehnite-pumpellyite facies rocks. Suchvariations in the assemblages and mineral compositions are controlledby a sequence of continuous and discontinuous reactions, andallow delineation of T-X_Fe³⁺ relations at constant pressure.The transition from the zeolite to prehnite-pumpellyite faciesof the Karmutsen metabasites is defined by a discontinuous reaction:0·18 laumontite+pumpellyite+0·15 quartz = 1·31prehnite+ 0·78 epidote+0·2 chlorite+ 1·72H₂O, where the X_Fe³⁺ values of prehnite, pumpellyite and epidoteare 0·03, 0·10 and 0·18, respectively.These values together with available thermodynamic data andour preliminary experimental data are used to calculate theP-T condition for the discontinuous reaction as P = 1·1±0·5 kb and T = 190±30°C. The effectsof pressure on the upper stability of the zeolite facies assemblagesare discussed utilizing T-X_Fe³⁺ diagrams. The stability of thelaumontite-bearing assemblages for the zeolite facies metamorphismof basaltic rocks may be defined by either continuous or discontinuousreactions depending on the imposed metamorphic field gradient.Hence, the zeolite and prehnite-pumpellyite facies transitionboundary is multivariant.     

夏季东亚高空急流与天气尺度波动的气候特征之间的联系

利用NCEP/NCAR逐日资料(1979~2003年),分析了东亚夏季高空急流和天气尺度波动气候特征之间的联系。通过分析急流与天气尺度波动活动量Q(用于描述某一时间段内天气尺度波动活动强弱程度)、热量涡动输送(V′T′)的关系,发现:在气候平均意义下,天气尺度波动活动量在6、8月较大,这意味着在这两个月中,天气尺度波动较为剧烈。这和急流的强度变化是相对应的,即强的天气尺度波动活动量对应着强的急流,弱的天气尺度波动活动量对应着弱的急流。热量涡动输送也表现出与天气尺度波动活动量类似的特征。年际变化上,夏季平均和6、7、8各月平均的热量涡动输送和急流强度都存在显著的正相关。对强、弱急流年份条件下热量涡动输送合成分析的结果表明:热量涡动输送偏南时,对应着强的急流;热量涡动输送偏北时,对应着弱的急流。以上的研究结果说明,在气候平均意义和年际变化上夏季东亚高空急流和天气尺度波动之间有着密切的联系。     

High-resolution circulation forecasting of the Maenggol Channel, south coast of Korea

The Maenggol Channel and Uldolmok Strait, located on the south-west coast of Korea, have notably strong and complex currents due to tidal effects and to local geological factors. In these areas, electric power has been generated using strong tidal currents, the speed of which is more than 3 m/s during spring tides. The region also provides a shortcut for navigation. These tidal conditions are therefore sometimes useful, but may also cause terrible accidents or severe economic damage, in the absence of accurate information regarding ocean conditions. In April 2014, the passenger ferry MV Sewol capsized in the Maenggol Channel, with 295 passengers killed and 9 still missing. While this was unquestionably a man-made disaster, strong currents were one of the contributing causes. It was also difficult to conduct scuba diving rescue operations given strong current speeds, and accurate prediction of the time when the tide would turn was thus critically needed. In this research, we used the high-resolution coastal circulation forecasting system of KOOS (Korea Operational Oceanographic System) for analysis and simulation of strong tidal currents in such areas with many small islands, using measurements and modeling from this research area. For accurate prediction of tidal currents, small grid size-modeling was needed, and in this study, we identified a suitable grid size that offers efficiency as well as accuracy.