Preliminary study on assessment of nutrient transport in the Taihu Basin based on SWAT modeling

With the Taihu Basin as a study area, using the spatially distributed and mechanism-based SWAT model, preliminary simulations of nutrient transport in the Taihu Basin during the period of 1995:_2002 has been carried out. The topography, soil, meteorology and land use with industrial point pollution discharge, the loss of agricultural fertilizers, urban sewerage, and livestock drainages were all considered in the boundary conditions of the simulations. The model was calibrated and validated against water quality monitoring data from 2001 to 2002. The results show that the annual total productions of nitrogen (TN) and phosphorus (TP) into Lake Taihu are 40000t and 2000t respectively. Nutrient from the Huxi Region is a major resource for Lake Taihu. The non-point source (surface source) pollution is the main form of catchment sources of nutrients into Lake Taihu, occupied TN 53% and TP 56% respectively. TN and TP nutrients from industrial point pollution discharge are 30% and 16%, and sewerage in both forms of point source and non-point source are TN 31% and TP 47%. Both the loss of agricultural fertilizers and livestock drainages from the catchment should be paid more attention as an important nutrient source. The results also show that SWAT is an effective model for the simulation of temporally and spatially nutrient changes and for the assessment of the trends in a catchment scale.     

An Improved Approach for Parameterizing Surface-Layer Turbulent Transfer Coefficients in Numerical Models

Based on classic iterative computation results, new equations to calculate the surface turbulent transfer coefficients are proposed, which allow for large ratios of the momentum and heat roughness lengths. Compared to the Launiainen scheme, our proposed scheme generates results closer to classical iterative computations. Under unstable stratification, the relative error in the Launiainen scheme increases linearly with increasing instability, even exceeding 15%, while the relative error of the present scheme is always less than 8.5%. Under stable stratification, the Launiainen scheme uses two equations, one for 0 < Ri _B ≤ 0.08 and another for 0.08 < Ri _B ≤ 0.2, and does not consider the condition that Ri _B > 0.2, while its relative errors in the region 0 < Ri _B ≤ 0.2 exceed 31 and 24% for momentum and heat transfer coefficients, respectively. In contrast, the present scheme uses only one equation for 0 < Ri _B ≤ 0.2 and another equation for Ri _B > 0.2, and the relative error of the present scheme is always less than 14%.     

Late Quaternary clay minerals off Middle Vietnam in the western South China Sea: Implications for source analysis and East Asian monsoon evolution

High-resolution clay mineral records combined with oxygen isotopic stratigraphy over the past 450 ka during late Quaternary from Core MD05-2901 off Middle Vietnam in the western South China Sea are re-ported to reconstruct a history of East Asian monsoon evolution. Variations in Illite, chlorite, and kaolinite contents indicate a strong glacial-interglacial cyclicity, while changes in smectite content present a higher frequency cyclicity. The provenance analysis indicates a mixture of individual clay minerals from various sources surrounding the South China Sea. Smectite derived mainly from the Sunda shelf and its major source area of the Indonesian islands. Illite and chlorite originated mainly from the Mekong and Red rivers. Kaolinite was provided mainly by the Pearl River. Spectral analysis of the kaolin-ite/(illite chlorite) ratio displays a strong eccentricity period of 100 ka, implying the ice sheet-forced win-ter monsoon evolution; whereas higher frequency changes in the smectite content show an ice sheet-forced obliquity period of 41 ka, and precession periods of 23 and 19 ka and a semi-precession period of 13 ka as well, implying the tropical-forced summer monsoon evolution. The winter monsoon evolution is generally in coherence with the glacial-interglacial cyclicity, with intensified winter monsoon winds during glacials and weakened winter monsoon winds during interglacials; whereas the summer monsoon evolution provides an almost linear response to the summer insolation of low latitude in the Northern Hemisphere, with strengthened summer monsoon during higher insolation and weakened summer monsoon during lower insolation. The result suggests that the high-latitude ice sheet and low-latitude tropical factor could drive the late Quaternary evolution of East Asian winter and summer monsoons, respectively, implying their diplex and self-contained forcing mechanism.     

Late Carboniferous porphyry copper mineralization at La Voluntad,Neuquén,Argentina: Constraints from Re–Os molybdenite dating

The La Voluntad porphyry Cu–Mo deposit in Neuquén, Argentina, is one of several poorly known porphyry-type deposits of Paleozoic to Early Jurassic age in the central and southern Andes. Mineralization at La Voluntad is related to a tonalite porphyry from the Chachil Plutonic Complex that intruded metasedimentary units of the Piedra Santa Complex. Five new Re–Os molybdenite ages from four samples representing three different vein types (i.e., quartz–molybdenite, quartz–sericite–molybdenite and quartz–sericite–molybdenite ± chalcopyrite–pyrite) are identical within error and were formed between ~312 to ~316 Ma. Rhenium and Os concentrations range between 34 to 183 ppm and 112 to 599 ppb, respectively. The new Re–Os ages indicate that the main mineralization event at La Voluntad, associated to sericitic alteration, was emplaced during a time span of 1.7 ± 3.2 Ma and that the deposit is Carboniferous in age, not Permian as previously thought. La Voluntad is the oldest porphyry copper deposit so far recognized in the Andes and indicates the presence of an active magmatic arc, with associated porphyry style mineralization, at the proto-Pacific margin of Gondwana during the Early Pennsylvanian.     

Well placement optimization subject to realistic field development constraints

This work considers the well placement problem in reservoir management and field development optimization. In particular, it emphasizes embedding realistic and practical constraints into a mathematical optimization formulation. Such constraints are a prerequisite for the wider use of mathematical optimization techniques in well placement problems, since constraints are a way to incorporate reservoir engineering knowledge into the problem formulation. There are important design limitations that are used by the field development team when treating the well placement problem, and these limitations need to be articulated and eventually formalized within the problem before conducting the search for optimal well placements. In addition, these design limitations may be explicit or implicit. In this work, various design limitations pertaining to well locations have been developed in close collaboration with a field operator on the Norwegian Continental Shelf. Moreover, this work focuses on developing constraint-handling capability to enforce these various considerations during optimization. In particular, the Particle Swarm Optimization (PSO) algorithm is applied to optimize for the well locations, and various practical well placement constraints are incorporated into the PSO algorithm using two different constraint-handling techniques: a decoder procedure and the penalty method. The decoder procedure maps the feasible search space onto a cube and has the advantage of not requiring parameter tuning. The penalty method converts the constrained optimization problem into an unconstrained one by introducing an additional term, which is called a penalty function, to the objective function. In contrast to the penalty method, only feasible solutions are evaluated in the decoder method. Through numerical simulations, a comparison between the penalty method and the decoder technique is performed for two cases. We show that the decoder technique can easily be implemented for the well placement problem, and furthermore, that it performs better than the penalty method in most of the cases.     

Broad search for unstable resonant orbits in the planar circular restricted three-body problem

Unstable resonant orbits in the circular restricted three-body problem have increasingly been used for trajectory design using optimization and invariant manifold techniques. In this study, several methods for computing these unstable resonant orbits are explored including grid searches, flyby maps, and continuation. Families of orbits are computed focusing on orbits with multiple loops near the secondary in the Jupiter–Europa system, and their characteristics are explored. Different parameters such as period and stability are examined for each set of resonant orbits, and the continuation of several specific orbits is explored in more detail.     

Non-Linear Theory and Power-Law Models for Information Integration and Mineral Resources Quantitative Assessments

Singular physical or chemical processes may result in anomalous amounts of energy release or mass accumulation that, generally, are confined to narrow intervals in space or time. Singularity is a property of different types of non-linear natural processes including cloud formation, rainfall, hurricanes, flooding, landslides, earthquakes, wildfires, and mineralization. The end products of these non-linear processes can be modeled as fractals or multifractals. Hydrothermal processes in the Earth’s crust can result in ore deposits characterized by high concentrations of metals with fractal or multifractal properties. Here we show that the non-linear properties of the end products of singular mineralization processes can be applied for prediction of undiscovered mineral deposits and for quantitative mineral resource assessment, whether for mineral exploration or for regional, national and global planning for mineral resource utilization. In addition to the general theory and framework for the non-linear mineral resources assessment, this paper focuses on several power-law models proposed for characterizing non-linear properties of mineralization and for geoinformation extraction and integration. The theories, methods, and computer system discussed in this paper were validated using a case study dealing with hydrothermal Au mineral potential in southern Nova Scotia, Canada.     

A Note on Reviving the Goddard Satellite-Based Surface Turbulent Fluxes （GSSTF） Dataset

Accurate sea surface flux measurements are crucial for understanding the global water and energy cycles. The oceanic evaporation, which is a major component of the global oceanic fresh water flux, is useful for predicting oceanic circulation and transport. The global Goddard Satellite-based Surface Turbulent Fluxes Version-2 (GSSTF2; July 1987--December 2000) dateset that was officially released in 2001 has been widely used by scientific community for global energy and water cycle research, and regional and short period data analyses. We have recently been funded by NASA to resume processing the GSSTF dataset with an objective of continually producing a uniform dataset of sea surface turbulent fluxes, derived from remote sensing data. The dataset is to be reprocessed and brought up-to-date (GSSTF2b) using improved input datasets such as a recently upgraded NCEP/DOE sea surface temperature reanalysis, and an upgraded surface wind and microwave brightness temperature V6 dataset (Version 6) from the Special Sensor Microwave Imager (SSM/I) produced by Remote Sensing Systems (RSS). A second new product (GSSTF3) is further proposed with a finer temporal (12-h) and spatial (0.25^ox0.25^o) resolution. GSSTF2b (July 1987--December 2008) and GSSTF3 (July 1999--December 2009) will be released for the research community to use by late 2009 and early 2011, respectively.     

Errors and correction of precipitation measurements in China

In order to discover the range of various errors in Chinese precipitation measurements and seek a correction method, 30 precipitation evaluation stations were set up countrywide before 1993. All the stations are reference stations in China. To seek a correction method for wind-induced error, a precipitation correction instrument called the ＂horizontal precipitation gauge＂ was devised beforehand. Field intercomparison observations regarding 29,000 precipitation events have been conducted using one pit gauge, two elevated operational gauges and one horizontal gauge at the above 30 stations. The range of precipitation measurement errors in China is obtained by analysis of intercomparison measurement results. The distribution of random errors and systematic errors in precipitation measurements are studied in this paper. A correction method, especially for wind-induced errors, is developed. The results prove that a correlation of power function exists between the precipitation amount caught by the horizontal gauge and the absolute difference of observations implemented by the operational gauge and pit gauge. The correlation coefficient is 0.99. For operational observations, precipitation correction can be carried out only by parallel observation with a horizontal precipitation gauge. The precipitation accuracy after correction approaches that of the pit gauge. The correction method developed is simple and feasible.