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.     

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.     

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.     

Relationship between positive cerium anomaly and adsorbed water in Antarctic lunar meteorites

Abstract— Concentrations of adsorbed water in single mineral grains of Antarctic lunar meteorites were determined with micro infrared (IR) spectroscopy. A relationship was found between the mineral's ability to adsorb water and the extent of Ce anomaly in rare earth element (REE) patterns precisely determined by the isotope dilution method using a thermal ionization mass spectrometer. Asuka (A) 881757, a lunar meteorite from the mare basalt without Ce anomaly, showed no trace of IR absorption due to adsorbed water. On the contrary, Yamato (Y) 791197-109, Y-86032-98, Y-86032-95, Y-791197-115 and Y-82192-55A from the lunar highland exhibiting positive Ce anomaly showed IR absorption due to adsorbed water in some of their minerals. The detected water would be of terrestrial origin, because it was not structurally bound and easy to exchange judging from the spectral band shape. The contrast in concentration of adsorbed water between the lunar highland and the mare basalt is derived from a difference in the density of microfractures in mineral grains. Average concentrations of adsorbed water in the lunar highland meteorites were 3.8 mg/cm³ for pyroxene and olivine, and 1.7 mg/cm³ for plagioclase. This contrast between minerals is noteworthy because it has been known that Ce anomaly of pyroxene and olivine is larger than that of plagioclase both for Antarctic lunar meteorites and some lunar rocks. Furthermore, more adsorbed water was detected for minerals in meteorites that exhibit larger Ce anomaly. The present observations demonstrate that the extent of Ce anomaly correlated with the concentration of adsorbed water, which suggests that active mineral surface resulting in adsorption of water could be a trace of interaction forming Ce anomaly. Terrestrial weathering on Antarctica and REE fractionation on the Moon are discussed for possible origins of the Ce anomaly.     

Denudation history of the Kiso Range,central Japan,and its tectonic implications: Constraints from low‐temperature thermochronology

Fission‐track (FT) and (U–Th–Sm)/He (He) analyses are used to constrain the denudation pattern and history of the Kiso Range, a Japanese fault‐block mountain range which has been uplifted since ca 0.8 Ma. We obtained nine zircon FT ages ranging 59.3–42.1 Ma, 18 apatite FT ages ranging 81.9–2.3 Ma, and 13 apatite He ages ranging 36.7–2.2 Ma. The apatite FT and He ages are divided into an older group comparable to the zircon FT age range and a younger group of <18 Ma. The younger ages are interpreted as a reflection of uplift of the Kiso Range because they were obtained only to the east of the Seinaiji‐touge Fault, and the event age estimated from apatite FT data is consistent with the timing of the onset of the Kiso Range uplift. On the basis of the distribution of the younger ages, we propose westward tilting uplift of the Kiso Range between the boundary fault of the Inadani Fault Zone and Seinaiji‐touge Fault, which implies a model of bedrock uplift that is intermediate between two existing models: a pop‐up model in which the Kiso Range is squeezed upward between the two faults and a tilted uplift model which assumes that the Kiso Range is uplifted and tilted to the west by the Inadani Fault Zone. The original land surface before the onset of uplift/denudation of the Kiso Range is estimated to have been uplifted to an elevation of 2700–4900 m. We estimated denudation rates at 1.3–4.0 mm/y and maximum bedrock uplift rates at 3.4–6.1 mm/y since ca 0.8 Ma. The Seinaiji‐touge fault is interpreted as a back thrust of the west‐dipping Inadani Fault Zone. The older group of apatite FT and He ages is interpreted to reflect long‐term peneplanation with a probable denudation rate of <0.1 mm/y.     

An Experimentally Constrained Petrogenetic Grid in the Silica-Saturated Portion of the System KFMASH at High Temperatures and Pressures

Experiments in the quartz-saturated part of the system KFMASHunder fO₂ conditions of the haematite–magnetite bufferand using bulk compositions with X_Mg of 0·81, 0·72,0·53 define the stability limits of several mineral assemblageswithin the P–T field 9–12 kbar, 850–1100°C.The stability limits of the mineral assemblages orthopyroxene+ spinel + cordierite ± sapphirine, orthopyroxene + garnet+ sapphirine, sapphirine + cordierite + orthopyroxene and garnet+ orthopyroxene + spinel have been delineated on the basis ofP–T and T–X pseudosections. Sapphirine did not appearin the bulk composition of X_Mg = 0·53. A partial petrogeneticgrid applicable to high Mg–Al granulites metamorphosedat high fO₂, developed in our earlier work, was extended tohigher pressures. The experimental results were successfullyapplied to several high-grade terranes to estimate P–Tconditions and retrograde P–T trajectories. KEY WORDS: KFMASH equilibria; experimental petrogenetic grid at high fO₂     

POST-VOLCANIC EROSION RATES OF SHOMYO FALLS IN TATEYAMA, CENTRAL JAPAN

Post‐eruptive fluvial erosion of welded pyroclastic flow deposits often depends on the recession of waterfalls because of their rapid erosion involved. We examine the recession rate of Shomyo Falls, which consists of Pleistocene welded pyroclastic flow deposits in Tateyama, north‐central Japan. The mean recession rate of the waterfall obtained from lithological and topographical evidences is 0.08‐0.15 m/a for 100000 a. However, the recession rate estimated by means of an empirical equation comprising physical parameters of erosive force and bedrock resistance is 0.006‐0.011 m/a with small uncertainties. The discrepancy between the geology‐based and equation‐based recession rates indicates that some factors, not taken account of in the equation, significantly influence the recession rate. We suggest that a factor in the rapid erosion of the waterfall is a large amount of transported sediment acting as abrasive material, which is supplied from high mountains in the watershed above the waterfall.