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%.     

Assessment of PPP integer ambiguity resolution using GPS,GLONASS and BeiDou (IGSO,MEO) constellations

Although integer ambiguity resolution (IAR) can improve positioning accuracy considerably and shorten the convergence time of precise point positioning (PPP), it requires an initialization time of over 30 min. With the full operation of GLONASS globally and BDS in the Asia–Pacific region, it is necessary to assess the PPP–IAR performance by simultaneous fixing of GPS, GLONASS, and BDS ambiguities. This study proposed a GPS + GLONASS + BDS combined PPP–IAR strategy and processed PPP–IAR kinematically and statically using one week of data collected at 20 static stations. The undifferenced wide- and narrow-lane fractional cycle biases for GPS, GLONASS, and BDS were estimated using a regional network, and undifferenced PPP ambiguity resolution was performed to assess the contribution of multi-GNSSs. Generally, over 99% of a posteriori residuals of wide-lane ambiguities were within ±0.25 cycles for both GPS and BDS, while the value was 91.5% for GLONASS. Over 96% of narrow-lane residuals were within ±0.15 cycles for GPS, GLONASS, and BDS. For kinematic PPP with a 10-min observation time, only 16.2% of all cases could be fixed with GPS alone. However, adding GLONASS improved the percentage considerably to 75.9%, and it reached 90.0% when using GPS + GLONASS + BDS. Not all epochs could be fixed with a correct set of ambiguities; therefore, we defined the ratio of the number of epochs with correctly fixed ambiguities to the number of all fixed epochs as the correct fixing rate (CFR). Because partial ambiguity fixing was used, when more than five ambiguities were fixed correctly, we considered the epoch correctly fixed. For the small ratio criteria of 2.0, the CFR improved considerably from 51.7% for GPS alone, to 98.3% when using GPS + GLONASS + BDS combined solutions.     

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.     

On the period of the periodic orbits of the restricted three body problem

We will show that the period T of a closed orbit of the planar circular restricted three body problem (viewed on rotating coordinates) depends on the region it encloses. Roughly speaking, we show that, \(2 T=k\pi +\int _\Omega g\) where k is an integer, \(\Omega \) is the region enclosed by the periodic orbit and \(g:{\mathbb {R}}^2\rightarrow {\mathbb {R}}\) is a function that only depends on the constant C known as the Jacobian constant; it does not depend on \(\Omega \). This theorem has a Keplerian flavor in the sense that it relates the period with the space “swept” by the orbit. As an application we prove that there is a neighborhood around \(L_4\) such that every periodic solution contained in this neighborhood must move clockwise. The same result holds true for \(L_5\).     

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.     

Recurrent mesoproterozoic continental magmatism in South-Central Norway

We report U–Pb dates and Lu–Hf isotope data, obtained by LAM-ICPMS, for zircons from metamorphic rocks of the Setesdalen valley, situated in the Telemark block south of the classic Telemark region of southern Norway. The samples include infracrustal rocks from the metamorphic basement, metaigneous rocks and metasediments from the Byglandsfjorden supracrustal cover sequence, and metaigneous rocks which intruded the whole succession. The main crustal evolution took place from 1,550–1,020 Ma, beginning with the emplacement of juvenile tonalitic melts; the contribution of older crustal material increased with time. Around 1,320 Ma, further addition of juvenile material occurred, involving both mafic and felsic melts, metamorphism and deformation. Acid magmas with high FeO*/MgO were intruded at 1,215 Ma, coinciding with underplating elsewhere in South Norway. The period starting at 1,215 Ma is represented by supracrustal rocks, principally metarhyolites with minor mafic material and immature sediments of the Byglandsfjorden Group. The crust generation processes ended with the intrusion of diorites and granodiorites at 1,030 Ma, late in the Sveconorwegian orogeny. Regional processes of metamorphism and deformation (around 1,290 and 1,000 Ma) can be related to the assembly of Rodinia. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.