Holocene and Eemian climatic optima in the Korean Peninsula based on textural and carbon isotopic records from the stalagmite of the Daeya Cave,South Korea

Textural and stable isotopic records from the absolute-dated stalagmite of the Daeya Cave (DY-1) provide new insights into the climatic evolution of the Korean Peninsula during the Holocene and Eemian climatic optima. The stalagmite yielded ages of 8572 ± 227 to 5907 ± 158 and 1,23,456 ± 535 to 1,19,837 ± 1089 years, which coincide with the Holocene and Eemian climatic optima, respectively. The stalagmite’s δ¹³C record closely resembles previously reported Chinese speleothem δ¹⁸O data. Thus it can be suggested that textural and geochemical results of the DY-1 reflect East Asian monsoon intensity, which is forced by summer insolation patterns in the northern hemisphere. Lighter carbon isotopic compositions, well-developed fibrous calcite crystals, and their relatively faster growth rate in the stalagmite sample are interpreted to reflect the warmest and wettest climate conditions of the Holocene and Eemian interglacials. Both climatic optima took place when insolation was decreasing from its maximum level, temperature in Greenland was highest, and sea level approached its maximum level. These climatic optima also coincide with decreasing Antarctic temperatures. Compared the DY-1 data to other proxies, it is suggested that the Holocene and Eemian climatic optima developed through a balance among boreal insolation, monsoon intensity, and sea level (also continental ice volume), which are the main climatic forcing factors in the northern hemisphere. These trends also follow the bi-polar seesaw mechanism as previously described.     

A geochronological and petrological study of anatectic paragneiss and associated granite dykes from the Day Nui Con Voi metamorphic core complex,North Vietnam: constraints on the timing of metamorphism within the Red River shear zone

The Red River shear zone (RRSZ) is a major left‐lateral strike‐slip shear zone, containing a ductilely deformed metamorphic core bounded by brittle strike‐slip and normal faults, which stretches for >1000 km from Tibet through Yunnan and North Vietnam to the South China Sea. The RRSZ exposes four high‐grade metamorphic core complexes along its length. Various lithologies from the southernmost core complex, the Day Nui Con Voi (DNCV), North Vietnam, provide new constraints on the tectonic and metamorphic evolution of this region prior to and following the initial India–Asia collision. Analysis of a weakly deformed anatectic paragneiss using P–T pseudosections constructed in the MnO–Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (MnNCKFMASHTO) system provides prograde, peak and retrograde metamorphic conditions, and in situ U–Th–Pb geochronology of metamorphic monazite yields texturally controlled age constraints. Tertiary metamorphism and deformation, overprinting earlier Triassic metamorphism associated with the Indosinian orogeny and possible Cretaceous metamorphism, are characterized by peak metamorphic conditions of ~805 °C and ~8.5 kbar between c. 38 and 34 Ma. Exhumation occurred along a steep retrograde P–T path with final melt crystallizing at the solidus at ≥~5.5 kbar at ~790 °C. Further exhumation at ~640–700 °C and ~4–5 kbar at c. 31 Ma occurred at subsolidus conditions. U–Pb geochronological analysis of monazite from a strongly deformed pre‐kinematic granite dyke from the flank of the DNCV provides further evidence for exhumation at this time. Magmatic grains suggest initial emplacement at 66.0 ± 1.0 Ma prior to the India–Asia collision, whereas grains with metamorphic characteristics indicate later growth at 30.6 ± 0.4 Ma. Monazite grains from a cross‐cutting post‐kinematic dyke within the core of the DNCV antiform provide a minimum age constraint of 25.2 ± 1.4 Ma for the termination of fabric development. A separate and significant episode of monazite growth at c. 83–69 Ma is suggested to be the result of fluid‐assisted recrystallization following the emplacement of magmatic units.     

Modelling combined open channel flow by artificial neural networks

Combined open channel flow is encountered in many hydraulic engineering structures and processes, such as irrigation ditches and wastewater treatment facilities. Extensive experimental studies have conducted to investigate combined flow characteristics. Nevertheless, there is no simple relationship that can fully describe the velocity profiles in a turbulent flow. The artificial neural network (ANN) has great computational capability for solving various complex problems, such as function approximation. The main objective of this study is to evaluate the applicability of the ANN for simulating velocity profiles, velocity contours and estimating the discharges accordingly. The velocity profiles measured by an acoustic doppler velocimeter in the open channel of the Chihtan purification plant, Taipei, with different discharges at fixed measuring section and different depths are presented. The total number of data sets is 640 and the data sets are split into two subsets, i.e. training and validation sets. The backpropagation algorithm is used to construct the neural network. The results demonstrate that the velocity profiles can be modelled by the ANN, and the ANN constructed can nicely fit the velocity profiles and can precisely predict the discharges for the conditions investigated. Copyright © 2005 John Wiley & Sons, Ltd.     

The economics of the marine sector: general equilibrium and ecological analysis

Recent transfers of potential marine wealth to coastal states, resulting from the Third United Nations Convention on the Law of the Sea, have enhanced our awareness that the marine sector might ultimately prove to be a new frontier. The objective of this paper is to initiate work toward the development of a systematic appraoch to the economics of this sector. To this end, the marine sector is assumed to be the fourth factor of production, taking its place with the traditional trio of land, capital and labor. The productivity of this sector is described and analyzed with the aid of a simple general equilibrium model.     

Time-delay control of structures

Time-delay causes unsynchronized application of the control forces which may not only degrade the performance of the control system but also even induce instability to the dynamic system. Time-delay control algorithm is developed in this paper to solve this practical problem. The control system is first formulated in discrete-time form. In the presence of time-delay, the motion equation of the discrete-time control system remains a difference equation which can be transformed into first-order difference equation by augmenting the state variables. Optimal time-delay control algorithm is derived based on the augmented system. The time-delay control forces are simply generated from the time-delay states multiplied by the constant feedback gain. Numerical simulation is illustrated to verify the feasibility of the proposed control algorithm. Since time-delay effect is incorporated in the mathematical model for the structural control system throughout the derivation of the proposed algorithm, system performance and dynamic stability are guaranteed.     

Intercomparison of210Pb measurements at GEOSECS station 500 in the northeast Pacific

During reoccupation of the GEOSECS-I test station in May, 1979, more than eighty 30-liter Niskin samples were collected in profile, many as replicates, for²¹⁰Pb intercomparison measurements by the WHOI, SIO and Yale groups. In addition to the inter-laboratory comparisons, the SIO group also carried out extensive experiments to test the effect of sample scavenging method. Pb equilibration time (storage effect), and filtration process on the measured²¹⁰Pb results.The intercomparison measurements indicate that there is a general agreement between the various sets of data. The sample set which allows a direct comparison at the same depth was available in most cases only between two of the three groups. The direct paired comparison shows that (1) the WHOI data are systematically 3% lower than the SIO data; (2) there are no systematic differences observed between the SIO and Yale data although the scatter is rather large; (3) the Yale data are systematically higher than the WHOI data by about 8%.The SIO experiments show that (1) the two scavenging methods employed (Fe(OH)₃ and Co-APDC co-precipitation) yield identical²¹⁰Pb results; (2) variation of Pb carrier equilibration time or of storage time has no discernible effect; (3) the filtration apparatus and procedure employed at this station do not result in²¹⁰Pb loss or contamination.The²¹⁰Pb profile structure and absolute concentration measured earlier at the same location (GOGO-II test station and GEOSECS station 347) agree with those of station 500 within 10%. The present profile shows a minimum²¹⁰Pb concentration around 500 m depth, marking the penetration depth of the flux of excess²¹⁰Pb from the atmosphere. There is a mild mid-depth maximum around 2500–3000 m. The²¹⁰Pb/²²⁶Ra activity ratio decreases monotonically from about 1 at the²¹⁰Pb minimum to about 0.5 near the bottom. The particulate²¹⁰Pb profile shows a systematic increase from the subsurface water to the bottom water by a factor of 5. This feature has been observed in many GEOSECS particulate²¹⁰Pb profiles.     

210Pb in the Pacific: the GEOSECS measurements of particulate and dissolved concentrations

We report here on particulate and dissolved²¹⁰Pb profiles at 16 stations, and on total²¹⁰Pb profiles at 3 stations, all occupied during the Pacific GEOSECS expedition. Comparison with measurements at Yale on GEOSECS library samples indicates that during separation of particulate lead from dissolved lead, our filtered water samples suffered some loss of²¹⁰Pb in the filtration system; this effect appears to have reduced the dissolved²¹⁰Pb activities by ～ 20% in stations where the water was filtered. However, for these first Pacific data on the²¹⁰Pb distribution between the two phases, this effect does not significantly interfere with our recognition of the major features of both particulate and dissolved²¹⁰Pb distributions.The dissolved²¹⁰Pb profiles in general vary geographically, following the²²⁶Ra profiles. In deep water,²²⁶Ra increases northward and eastward from the southwest Pacific, from ～ 22dpm/100kg, to over 40 dpm/100 kg in the northeast Pacific. Our dissolved²¹⁰Pb profiles show a similar increase in deep water, varying from about 10 to 20 dpm/100 kg along this line, and are commonly characterized by a mid-depth maximum. This²¹⁰Pb maximum reflects the mid-depth²²⁶Ra maximum of the Pacific Deep Water observed along the western boundary current.In surface water at low latitudes there is a significant²¹⁰Pb flux from the atmosphere, which produces a²¹⁰Pb/²²⁶Ra activity ratio generally greater than unity. This flux penetrates as deep as 600 m, as indicated by an “induced”²¹⁰Pb minimum caused by the surface maximum. The surface water²¹⁰Pb excess decreases toward high southern latitudes and vanishes in the Circumpolar region.The particulate²¹⁰Pb profiles show a general increase with depth, from ～ 0.3dpm/100kg in subsurface water to ～ 1.5dpm/100kg in bottom water, with or without a mid-depth maximum that reflects the²²⁶Ra or dissolved²¹⁰Pb maximum. The particulate²¹⁰Pb normally comprises about 2% of the total²¹⁰Pb in subsurface water, and this fraction increases to about 10% near the bottom. As the filtration loss is not taken into account, the fraction of particulate²¹⁰Pb quoted here is an upper limit. Since the particulate matter concentrations are quite uniform in the water column below a few hundred meters, the²¹⁰Pb activity of the particulate matter also increases with depth. The particulate matter has a²¹⁰Pb concentration of ～ 100dpm/g in subsurface water, but the concentration increases to ～ 500dpm/g or more toward the bottom. This indicates that there is a cumulative adsorption of Pb onto the suspended particles as they are sinking through the water column.     

Mapping landslide susceptibility from small datasets: A case study in the Pays de Herve (E Belgium)

A landslide susceptibility map is proposed for the Pays de Herve (E Belgium), where large landslides affect Cretaceous clay outcrop areas. Based on a Bayesian approach, this GIS-supported probabilistic map identifies the areas most susceptible to deep landslides. The database is comprised of the source areas of ten pre-existing landslides (i.e. a sample of 154 grid cells) and of six environmental data layers, namely lithology, proximity to active faults, slope angle and aspect, elevation and distance to the nearest valley-floor. A 30-m-resolution DEM from the Belgian National Geographical Institute is used for the analysis. Owing to the small size of the sample, a special cross-validation procedure of the susceptibility map is performed, which uses in an iterative way each of the landslides to test the predictive power of the map derived from the other landslides. Four different sets of variables are used to produce four susceptibility maps, whose prediction curves are compared. While the prediction rates associated with the models not involving the “proximity to active fault” criterion are comparable to those of the models considering this variable, strong weaknesses inherent in the fault data on which the latter rely suggest that the final susceptibility map should be based on a model that excludes any reference to fault. This highlights the difference between a triggering factor and determining factors, and in the same time broadens the scope of the produced map. A single reactivated slide is also used to test the possibility of predicting future reactivation of existing landslides in the area. Finally, the need for geomorphological control over the mathematical treatment is underlined in order to obtain realistic prediction maps.     

Experimental and numerical study of soil-reinforcement effects on the low-strain stiffness and bearing capacity of shallow foundations

This paper presents results of meticulous laboratory testing and numerical simulations on the effect of reinforcement on the low-strain stiffness and bearing capacity of shallow foundations on dry sand. The effect of the location and the number of reinforcement layers is studied in the laboratory, whereas numerical simulations are used to study the reinforcement-foundation interaction. Laboratory tests show an increase of 100, 200, and 275% not only in bearing capacity but also in low-strain stiffness (linear load–displacement behaviour) of a square foundation when one, two, and three layers of reinforcement are used, respectively. The specimen preparation technique is found to be crucial for the repeatability and reliability of the laboratory results (less than 5% variability). Numerical simulations demonstrate that if reinforcements are placed up to a depth of one footing width (B) below the foundation, better re-distribution of the load to deeper layers is achieved, thus reducing the stresses and strains underneath the foundation. Numerical simulations and experimental results clearly identify a critical zone between 0.3 and 0.5B, where maximum benefits not only on the bearing capacity but also on the low-strain stiffness of the foundation are obtained. Therefore, soil reinforcement can also be used to reduce low-strain vibrations of foundations.