Distinct groundwater recharge sources and geochemical evolution of two adjacent sub-basins in the lower Shule River Basin,northwest China

Based on analysis of groundwater hydrogeochemical and isotopic data, this study aims to identify the recharge sources and understand geochemical evolution of groundwater along the downstream section of the Shule River, northwest China, including two sub-basins. Groundwater samples from the Tashi sub-basin show markedly depleted stable isotopes compared to those in the Guazhou sub-basin. This difference suggests that groundwater in the Tashi sub-basin mainly originates from meltwater in the Qilian Mountains, while the groundwater in the Guazhou sub-basin may be recharged by seepage of the Shule River water. During the groundwater flow process in the Tashi sub-basin, minerals within the aquifer material (e.g., halite, calcite, dolomite, gypsum) dissolve in groundwater. Mineral dissolution leads to strongly linear relationships between Na⁺ and Cl^? and between Mg²⁺+ Ca²⁺ and SO₄ ^2??+?HCO₃ ^?, with stoichiometry ratios of approximately 1:1 in both cases. The ion-exchange reaction plays a dominant role in hydrogeochemical evolution of groundwater in the Guazhou sub-basin and causes a good linear relationship between (Mg²⁺+ Ca²⁺)–(SO₄ ^2??+?HCO₃ ^?) and (Na⁺+ K⁺)–Cl^? with a slope of ?0.89 and also results in positive chloroalkaline indices CAI 1 and CAI 2. The scientific results have implications for groundwater management in the downstream section of Shule River. As an important irrigation district in Hexi Corridor, groundwater in the Guazhou sub-basin should be used sustainably and rationally because its recharge source is not as abundant as expected. It is recommended that the surface water should be used efficiently and routinely, while groundwater exploitation should be limited as much as possible.     

The radio counterpart of the likely TeV binary HESS J0632+057

The few known γ-ray binary systems are all associated with variable radio and X-ray emission. The TeV source HESS J0632+057, apparently associated with the Be star MWC 148, is plausibly a new member of this class. Following the identification of a variable X-ray counterpart to the TeV source we conducted Giant Metrewave Radio Telescope (GMRT) and Very Large Array (VLA) observations in 2008 June–September to search for the radio counterpart of this object. A point-like radio source at the position of the star is detected in both 1280-MHz GMRT and 5-GHz VLA observations, with an average spectral index, α, of ∼0.6. In the VLA data there is significant flux variability on ∼month time-scales around the mean flux density of ≈0.3 mJy. These radio properties (and the overall spectral energy distribution) are consistent with an interpretation of HESS J0632+057 as a lower power analogue of the established γ-ray binary systems.     

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.     

Applying essentially non‐oscillatory interpolation to controlled‐source electromagnetic modelling*

Modelling and inversion of controlled‐source electromagnetic (CSEM) fields requires accurate interpolation of modelled results near strong resistivity contrasts. There, simple linear interpolation may produce large errors, whereas higher‐order interpolation may lead to oscillatory behaviour in the interpolated result. We propose to use the essentially non‐oscillatory, piecewise polynomial interpolation scheme designed for piecewise smooth functions that contains discontinuities in the function itself or in its first or higher derivatives. The scheme uses a non‐linear adaptive algorithm to select a set of interpolation points that represent the smoothest part of the function among the sets of neighbouring points. We present numerical examples to demonstrate the usefulness of the scheme. The first example shows that the essentially non‐oscillatory interpolation (ENO) scheme better captures an isolated discontinuity. In the second example, we consider the case of sampling the electric field computed by a finite‐volume CSEM code at a receiver location. In this example, the ENO interpolation performs quite well. However, the overall error is dominated by the discretization error. The other examples consider the comparison between sampling with essentially non‐oscillatory interpolation and existing interpolation schemes. In these examples, essentially non‐oscillatory interpolation provides more accurate results than standard interpolation, especially near discontinuities.     

Geocenter variations derived from a combined processing of LEO- and ground-based GPS observations

GNSS observations provided by the global tracking network of the International GNSS Service (IGS, Dow et al. in J Geod 83(3):191–198, 2009) play an important role in the realization of a unique terrestrial reference frame that is accurate enough to allow a detailed monitoring of the Earth’s system. Combining these ground-based data with GPS observations tracked by high-quality dual-frequency receivers on-board low earth orbiters (LEOs) is a promising way to further improve the realization of the terrestrial reference frame and the estimation of geocenter coordinates, GPS satellite orbits and Earth rotation parameters. To assess the scope of the improvement on the geocenter coordinates, we processed a network of 53 globally distributed and stable IGS stations together with four LEOs (GRACE-A, GRACE-B, OSTM/Jason-2 and GOCE) over a time interval of 3 years (2010–2012). To ensure fully consistent solutions, the zero-difference phase observations of the ground stations and LEOs were processed in a common least-squares adjustment, estimating all the relevant parameters such as GPS and LEO orbits, station coordinates, Earth rotation parameters and geocenter motion. We present the significant impact of the individual LEO and a combination of all four LEOs on the geocenter coordinates. The formal errors are reduced by around 20% due to the inclusion of one LEO into the ground-only solution, while in a solution with four LEOs LEO-specific characteristics are significantly reduced. We compare the derived geocenter coordinates w.r.t. LAGEOS results and external solutions based on GPS and SLR data. We found good agreement in the amplitudes of all components; however, the phases in x- and z-direction do not agree well.     

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

The biogeochemical sulfur cycle in the marine boundary layer over the Northeast Pacific Ocean

The major components of the marine boundary layer biogeochemical sulfur cycle were measured simultaneously onshore and off the coast of Washington State, U.S.A. during May 1987. Seawater dimethylsulfide (DMS) concentrations on the continental shelf were strongly influenced by coastal upwelling. Concentration further offshore were typical of summer values (2.2 nmol/L) at this latitude. Although seawater DMS concentrations were high on the biologically productive continental shelf (2–12 nmol/L), this region had no measurable effect on atmospheric DMS concentrations. Atmospheric DMS concentrations (0.1–12 nmol/m³), however, were extremely dependent upon wind speed and boundary layer height. Although there appeared to be an appreciable input of non-sea-salt sulfate to the marine boundary layer from the free troposphere, the local flux of DMS from the ocean to the atmosphere was sufficient to balance the remainder of the sulfur budget.     

Impacts of the leading modes of tropical Indian Ocean sea surface temperature anomaly on sub-seasonal evolution of the circulation and rainfall over East Asia during boreal spring and summer

The two leading modes of the interannual variability of the tropical Indian Ocean (TIO) sea surface temperature (SST) anomaly are the Indian Ocean basin mode (IOBM) and the Indian Ocean dipole mode (IODM) from March to August. In this paper, the relationship between the TIO SST anomaly and the sub-seasonal evolution of the circulation and rainfall over East Asia during boreal spring and summer is investigated by using correlation analysis and composite analysis based on multi-source observation data from 1979 to 2013, together with numerical simulations from an atmospheric general circulation model. The results indicate that the impacts of the IOBM on the circulation and rainfall over East Asia vary remarkably from spring to summer. The anomalous anticyclone over the tropical Northwest Pacific induced by the warm IOBM is closely linked with the Pacific–Japan or East Asia–Pacific teleconnection pattern, which persists from March to August. In the upper troposphere over East Asia, the warm phase of the IOBM generates a significant anticyclonic response from March to May. In June and July, however, the circulation response is characterized by enhanced subtropical westerly flow. A distinct anomalous cyclone is found in August. Overall, the IOBM can exert significant influence on the western North Pacific subtropical high, the South Asian high, and the East Asian jet, which collectively modulate the precipitation anomaly over East Asia. In contrast, the effects of the IODM on the climate anomaly over East Asia are relatively weak in boreal spring and summer. Therefore, studying the impacts of the TIO SST anomaly on the climate anomaly in East Asia should take full account of the different sub-seasonal response during boreal spring and summer.     

Nutrient pore water chemistry,Great Bay,New Hampshire: Benthic fluxes

Benthic fluxes of C, N, P and Si have been measured at two sites in Great Bay Estuary, New Hampshire. Higher fluxes of reactive phosphate, nitrate and reactive silicate were observed at the site where bioturbation is known to occur and the fluxes of NH₄ ⁺, PO₄ ^?3 and reactive silicate at this location were from 3 to 6 times higher than that calculated by simple pore water diffusion models. *** DIRECT SUPPORT *** A01BY019 00010     

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.