西江流域环境河水地球化学研究

River is a major component of the global surface water and CO2 cycles. The chemistry of river waters reveals the nature of weathering on a basin-wide scale and helps us understand the exogenic cycles of elements in the continent-river-ocean system. In particular, geochemical investigation of large river gives important information on the biogeochemical cycles of the elements, chemical weathering rates, physical erosion rates and CO2 consumption during the weathering of the rocks within the drainage basin. Its importance has led to a number of detailed geochemical studies on some of the world's large and medium-size river systems. Flowing in the south of China, the Xijiang River is the second largest river in the China with respect to its discharge, after the Yangtze River. Its headwaters drain the YunGui Plateau, where altitude is approximately 2000 meters. Geologically, the carbonate rocks are widely spread in the river drainage basin, which covers an area of about 0.17xl06 km2, i.e., 39% of the whole drainage basin. This study focuses on the chemistry of the Xijiang river system and constitutes the first geochemical investigation into major and trace elements concentrations for both suspended and dissolved loads of this river and its main tributaries, and Sr isotopic composition of the dissolved load is also investigated, in order to determine both chemical weathering and mechanical erosion rates. As compared with the other large rivers of the world, the Xijiang River is characterized by higher major element concentration. The dissolved major cations average 1.17, 0.33, 0.15, and 0.04 mmol I"1 for Ca, Mg, Na, and K, respectively. The total cation concentrations (TZ+) in these rivers vary between 2.2 and 4.4 meq I'1. The high concentration of Ca and Mg, high (Ca+Mg)/(Na+K) ratio (7.9), enormous alkalinity and low dissolved SiO2/HCO3 ratio (0.05) in river waters reveal the importance of carbonate weathering and relatively weak silicate weathering over the river drainage basin. The major elements in river water, such as the alkalis and alkaline-earths, are of different origins: from rain water, silicate weathering, carbonate and evaporite weathering. A mixing model based on mass budget equation is used in this study, which allows the proportions of each element derived from the different source to be calculated. The carbonate weathering is the main source of these elements in the Xijiang drainage basin. The contribution of rainwater, especially for Na, reaches to approximately 50% in some tributaries. Dissolved elemental concentration of the river waters are corrected for rain inputs (mainly oceanic salts), the elemental concentrations derived from the different rock weathering are calculated. As a consequence, silicate, carbonate and total rock weathering rates, together with the consumption rates of atmospheric CO2 by weathering of each of these lithologies have been estimated. They provide specific chemical erosion rates varying between 5.1~17.8 t/km2/yr for silicate, 95.5~157.2 t/km2/yr for carbonate, and 100.6-169.1 t/km2/yr for total rock, respectively. CO2 consumptions by silicate and carbonate weathering approach 13><109and 270.5x10 mol/yr. Mechanical denudation rates deduced from the multi-year average of suspended load concentrations range from 92-874 t/km2/yr. The high denudation rates are mainly attributable to high relief and heavy rainfall, and acid rain is very frequent in the drainage basin, may exceed 50% and the pH value of rainwater may be <4.0, result from SO2 pollution in the atmosphere, results in the dissolution of carbonates and aluminosilicates and hence accelerates the chemical erosion rate. The compositions of minerals and elements of suspended particulate matter are also investigated. The most soluble elements (e.g. Ca, Na, Sr, Mg) are strongly depleted in the suspended phase with respect to upper continent crust, which reflects the high intensity of rock weathering in the drainage basin. Some elements (e.g. Pb, Cu, Co, Cr) show positive anomalies, Pb/Th ratios in suspended matter approach 7 times (Liu Jiang) to 10 times (Nanpan Jiang) the crustal value. The enrichment of these elements in suspended matter reflects the intensity both of anthropogenic pollution and adsorption processes onto particles. The contents of the soluble fraction of rare earth elements (REE) in the river are low, and REE mainly reside in particulate phase. In dissolved phase, the PAAS-normalized distribution patterns show significant HREE enrichment with (La/Yb) SN=0.26~0.94 and Ce depletion with (Ce/Ce*) SN=0.31-0.98, and the most pronounced negative Ce anomalies occur in rivers of high pH. In the suspended phase, the rivers have LREE-enriched patterns relative to PAAS, with (La/Yb) SN=1 -00-1 .40. The results suggest that pH is a major factor controlling both the absolute abundances of REE in solution and the fractionation of REE of dissolved phase. Ce depletion in river waters with high pH values results probably from both preferential removal of Ce onto Fe-Mn oxide coating of particles and CeC^ sedimentation. This process is known to occur in the marine environment and may also occur in high pH rivers. Positive correlations are also observed between La/Yb ratio and DOC, HCO3", PO4", suggesting that colloids and (or) adsorption processes play an important role in the control of these elements.

Spatiotemporal variations of internal P-loading and the related mechanisms in the large shallow Lake Chaohu

Spatiotemporal variations of P species and adsorption behavior in water column, interstitial water, and sediments were investigated in the large shallow eutrophic Lake Chaohu. Orthophosphate (Ortho-P) and total phosphorus (TP) concentrations were significantly higher in the western part than in the eastern part of the lake, due to different nutrient inputs from the surrounding rivers. Moreover, particulate phosphorus (PP) concentration was in a similar spatial pattern to Ortho-P and TIP concentrations, and also showed significantly positive correlation with the biomass of Microcystis, indicating more uptake and store of phosphorus by Microcystis than by other algae. Increase of pH and intensive utilization of P by phytoplankton were the main factors promoting P (especially Fe-P) release from the sediment to interstitial water during the cyanobacterial blooms in Lake Chaohu. Spatial dynamics in TP concentration, P species and adsorption behavior of the sediment, coupled with the statistical analyses, suggested that the spatial heterogeneity of P contents in the sediment was influenced by various factors, e.g. human activities, soil geochemistry and mineral composition. In spite of similar TP contents in the sediments, increase in proportion of Fe-P concentration in the sediment may result in a high risk of P release.

Burst Detection in Turbulent Channel Flows Based on Large Eddy Simulation Databases

Reliable turbulent channel flow databases at several Reynolds numbers have been established by large eddy simulation (LES), with two of them validated by comparing with typical direct numerical simulation (DNS) results. Furthermore, the statistics, such as velocity profile, turbulent intensities and shear stress, were obtained as well as the temporal and spatial structure of turbulent bursts. Based on the LES databases available, the conditional sampling methods are used to detect the structures of burst events. A method to deterimine the grouping parameter from the probability distribution function (pdf) curve of the time separation between ejection events is proposed to avoid the errors in detected results. And thus, the dependence of average burst period on thresholds is considerably weakened. Meanwhile, the average burst-to-bed area ratios are detected. It is found that the Reynolds number exhibits little effect on the burst period and burst-to-bed area ratio.

Increasing Critical Sensitivity of the Load/Unload Response Ratio before Large Earthquakes with Identified Stress Acclumulation Pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.

State Vector: A New Approach to Prediction of the Failure of Brittle Heterogeneous Media and Large Earthquakes

The concept of state vector stems from statistical physics, where it is usually used to describe activity patterns of a physical field in its manner of coarsegrain. In this paper, we propose an approach by which the state vector was applied to describe quantitatively the damage evolution of the brittle heterogeneous systems, and some interesting results are presented, i.e., prior to the macro-fracture of rock specimens and occurrence of a strong earthquake, evolutions of the four relevant scalars time series derived from the state vectors changed anomalously. As retrospective studies, some prominent large earthquakes occurred in the Chinese Mainland (e.g., the M 7.4 Haicheng earthquake on February 4, 1975, and the M 7.8 Tangshan earthquake on July 28, 1976, etc) were investigated. Results show considerable promise that the time-dependent state vectors could serve as a kind of precursor to predict earthquakes.

Generation of large-scale vortex dislocations in a three-dimensional wake-type flow

Numerical study of three-dimensional evolution of wake-type flow and vortex dislocations is performed by using a compact finite diffenence-Fourier spectral method to solve 3-D incompressible Navier-Stokes equations. A local spanwise nonuniformity in momentum defect is imposed on the incoming wake-type flow. The present numerical results have shown that the flow instability leads to three-dimensional vortex streets, whose frequency, phase as well as the strength vary with the span caused by the local nonuniformity. The vortex dislocations are generated in the nonuniform region and the large-scale chain-like vortex linkage structures in the dislocations are shown. The generation and the characteristics of the vortex dislocations are described in detail.

Large-Eddy-Simulation of 3-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids

The 3-dimensiqnal incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation (LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinu

Effects of Subgrid-Scale Modeling on Time Correlations in Large Eddy Simulation

The effects of the unresolved subgrid-scale (SGS) motions on the energy balance of the resolved scales in large eddy simulation (LES) have been investigated actively because modeling the energy transfer between the resolved and unresolved scales is crucial to constructing accurate SGS models. But the subgrid scales not only modify the energy balance, they also contribute to temporal decorrelation of the resolved scales. The importance of this effect in applications including the predictability problem and the evaluation of sound radiation by turbulent flows motivates the present study of the effect of SGS modeling on turbulent time correlations. This paper compares the two-point, two-time Eulerian velocity correlation in isotropic homogeneous turbulence evaluated by direct numerical simulation (DNS) with the correlations evaluated by LES using a standard spectral eddy viscosity. It proves convenient to express the two-point correlations in terms of spatial Fourier decomposition of the velocity field. The LES fields are more coherent than the DNS fields: their time correlations decay more slowly at all resolved scales of motion and both their integral scales and microscales are larger than those of the DNS field. Filtering alone is not responsible for this effect: in the Fourier representation, the time correlations of the filtered DNS field are identical to those of the DNS field itself. The possibility of modeling the decorrelating effects of the unresolved scales of motion by including a random force in the model is briefly discussed. The results could have applications to the problem of computing sound sources in isotropic homogeneous turbulence by LES

Carbon transition efficiency and process cost in high-rate, large-area deposition of diamond films by DC arc plasma jet

A new DC plasma torch in which are jet states and deposition parameters can be regulated over a wide range has been built. It showed advantages in producing stable plasma conditions at a small gas flow rate. Plasma jets with and without magnetically rotated arcs could be generated. With straight are jet deposition, diamond films could be formed at a rate of 39 mu m/h on Mo substrates of Phi 25 mm, and the conversion rate of carbon in CH4 to diamond was less than 3%. Under magnetically rotated conditions, diamond films could be deposited uniformly in a range of Phi 40 mm at 30 mu m/h, with a quite low total gas flow rate and high carbon conversion rate of over 11%. Mechanisms of rapid and uniform deposition of diamond films with low gas consumption and high carbon transition efficiency are discussed.

Large Eddy Simulation of Complex Turbulent Flows: Physical Aspects and Research Trends

In the current paper, we have primarily addressed one powerful simulation tool developed during the last decades-Large Eddy Simulation (LES), which is most suitable for unsteady three-dimensional complex turbulent flows in industry and natural environment. The main point in LES is that the large-scale motion is resolved while the small-scale motion is modeled or, in geophysical terminology, parameterized. With a view to devising a subgrid-scale(SGS) model of high quality, we have highlighted analyzing physical aspects in scale interaction and-energy transfer such as dissipation, backscatter, local and non-local interaction, anisotropy and resolution requirement. They are the factors responsible for where the advantages and disadvantages in existing SGS models come from. A case study on LES of turbulence in vegetative canopy is presented to illustrate that LES model is more based on physical arguments. Then, varieties of challenging complex turbulent flows in both industry and geophysical fields in the near future-are presented. In conclusion; we may say with confidence that new century shall see the flourish in the research of turbulence with the aid of LES combined with other approaches.

Large strain field near a crack tip in a rubber sheet

The distribution of stress-strain near a crack tip in a rubber sheet is investigated by employing the constitutive relation given by Gao (1997). It is shown that the crack tip field is composed of two shrinking sectors and one expanding sector. The stress state near the crack tip is in uniaxial tension. The analytical solutions are obtained for both expanding and shrinking sectors.

Space experimental device on Marangoni drop migrations of large Reynolds numbers

The space experimental device for testing the Marangoni drop migrations has been discussed in the present paper. The experiment is one of the spaceship projects of China. In comparison with similar devices, it has the ability of completing all the scientific experiments by both auto controlling and telescience methods. It not only can perform drop migration experiments of large Reynolds numbers but also has an equi-thick interferential system.

Teracluster LSSC-II - Its Designing Principles and Applications in Large Scale Numerical Simulations

The teracluster LSSC-II installed at the State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences is one of the most powerful PC clusters in China. It has a peek performance of 2Tflops. With a Linpack performance of 1.04Tflops, it is ranked at the 43rd place in the 20th TOP500 List (November 2002), 51st place in the 21st TOP500 List (June 2003), and the 82nd place in the 22nd TOP500 List (November 2003) with a new Linpack performance of 1.3Tflops. In this paper, we present some design principles of this cluster, as well as its applications in some largescale numerical simulations.

A Lagrangian for von Karman Equations of Large Deflection Problem of Thin Circular Plate

By the semi-inverse method proposed by He, a Lagrangian is established for the large deflection problem of thin circular plate. Ritz method is used to obtain an approximate analytical solution of the problem. First order approximate solution is obtained, which is similar to those in open literature. By Mathematica a more accurate solution can be deduced.