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.     

The diel vertical migration of the sound-scattering layer in the Yellow Sea Bottom Cold Water of the southeastern Yellow Sea：focus on its relationship with a temperature structure

Using the hydroacoustic method with a 200 kHz scientific echo sounding system,the diel vertical migration(DVM) of the sound-scattering layer(SSL) in the Yellow Sea Bottom Cold Water(YSBCW) of the southeastern Yellow Sea was studied in April(spring) and August(summer) of 2010 and 2011.For each survey,13–27 hours of acoustic data were continuously collected at a stationary station.The acoustic volume scattering strength(Sv) data were analyzed with temperature profile data.In the spring of both 2010 and 2011,the SSL clearly showed the vertical migration throughout the entire water column,moving from the surface layer at night to near the bottom during the day.Conductivity,temperature,and depth data indicated that the entire water column was well mixed with low temperature of about 8 C.However,the SSL showed different patterns in the summers of 2010 and 2011.In the summer of 2010(28 C at the surface),the SSL migrated to near the bottom during the day,but there were two SSLs above and below the thermocline at depth of 10–30 m at night.In the summer of 2011(20 C at the surface),the SSL extended throughout the entire water column at night,possibly owing to an abrupt change in sea weather conditions caused by the passage of a Typhoon Muifa over the study area.It was concluded that the DVM patterns in summer in the YSBCW area may be greatly influenced by a strengthened or weakened thermocline.     

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.