The influence of island generated eddies on the carbon dioxide system, south of the Canary Islands

Measurements of surface partial pressure of CO₂ and water column alkalinity, pH_T, nutrients, oxygen, fluorescence and hydrography were carried out, south of the Canary Islands during September 1998. Cyclonic and anticyclonic eddies were alternatively observed from the northwestern area to the central area of the Canary Islands. Nutrient pumping and vertical uplifting of the deep chlorophyll maximum by cyclonic eddies were also ascertained by upward displacement of dissolved inorganic carbon. A model was applied to determine the net inorganic carbon balance in the cyclonic eddy. The fluxes were determined considering both the diffusive and convective contributions from the upward pumping and the corresponding horizontal transport of water outside the area. An increase in the total inorganic carbon concentration in the upper layers inside the eddy field of 133 mmol C m^− 2 d^− 1 was determined. The upward flux of inorganic carbon decreased the effect of the increased primary production on the carbon dioxide chemistry. The reduced fCO₂ inside the cyclonic eddy, 15 μatm lower than that observed in non-affected surface water, was explained by thermodynamic aspects, biological activity, eddy upward pumping and diffusion and air–sea water exchange effects.     

Applications of AnnAGNPS model for soil loss estimation and nutrient loading for Malaysian conditions

The study evaluated the performance and suitability of AnnAGNPS model in assessing runoff, sediment loading and nutrient loading under Malaysian conditions. The watershed of River Kuala Tasik in Malaysia, a combination of two sub-watersheds, was selected as the area of study. The data for the year 2004 was used to calibrate the model and the data for the year 2005 was used for validation purposes. Several input parameters were computed using methods suggested by other researchers and studies carried out in Malaysia. The study shows that runoff was predicted well with an overall R² value of 0.90 and E value of 0.70. Sediment loading was able to produce a moderate result of R² = 0.66 and E = 0.49, nitrogen loading predictions were slightly better with R² = 0.68 and E = 0.53, and phosphorus loading performance was slightly poor with an R² = 0.63 and E = 0.33. The erosion map developed was in agreement with the erosion risk map produced by the Department of Agriculture, Malaysia. Rubber estates and urban areas were found to be the main contributors to soil erosion. The simulation results showed that AnnAGNPS has the potential to be used as a valuable tool for planning and management of watersheds under Malaysian conditions.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

Study on Recovering the Earth＇s Potential Field Based on GOCE Gradiometry

Given the second radial derivative Vrr（P） ｜δs of the Earth＇s gravitational potential V（P） on the surface δS corresponding to the satellite altitude, by using the fictitious compress recovery method, a fictitious regular harmonic field rrVrr（P）^＊ and a fictitious second radial gradient field V：（P） in the domain outside an inner sphere Ki can be determined, which coincides with the real field V（P） in the domain outside the Earth. Vrr^＊（P）could be further expressed as a uniformly convergent expansion series in the domain outside the inner sphere, because rrV（P）^＊ could be expressed as a uniformly convergent spherical harmonic expansion series due to its regularity and harmony in that domain. In another aspect, the fictitious field V^＊（P） defined in the domain outside the inner sphere, which coincides with the real field V（P） in the domain outside the Earth, could be also expressed as a spherical harmonic expansion series. Then, the harmonic coefficients contained in the series expressing V^＊（P） can be determined, and consequently the real field V（P） is recovered. Preliminary simulation calculations show that the second radial gradient field Vrr（P） could be recovered based only on the second radial derivative V（P）｜δs given on the satellite boundary. Concerning the final recovery of the potential field V（P） based only on the boundary value Vrr （P）｜δs, the simulation tests are still in process.     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

大洋粘土制备白炭黑后的铝回收实验研究

为了回收大洋粘土制备白炭黑后残液中的Al3 ,并获得氧化铝超细粉体,实验选用NH3.H2O为沉淀剂,在常温下采用沉淀法使残液中Al3 以氢氧化铝形式沉淀并去除其它杂质离子,然后将得到的沉淀在1200℃灼烧,保温1h。利用IR和XRD测试了产物的物相组成和结构,证实实验可以获得纯度较高的α-Al2O3;而SEM和粒度分析结果显示产物粒度细小,部分颗粒可以达到纳米级;BET结果显示Al2O3的比表面积可达32m2/g。根据氧化铝结构和性能等优选的实验条件为大洋粘土的进一步开发利用提供了参考依据。     

Application of membrane technology for water treatment： Ultrasound-associated cleaning of fouled microfiltration membranes

The use of membranes is a widely employed, versatile, and effective separation process. One of the limiting aspects in applying microfiltration （MF） for wastewater treatment is that of problems with membrane fouling and consequent flux reduction. Membrane fouling occurs by the irreversible deposition of retained particles, colloids, macromolecules, salts, etc. at the membrane surface and/or inside the membrane. The predominant fouling mechanisms observed with MF membranes are classified as three categories： the build-up of a cake layer on the membrane surface, blocking of membrane pores, and adsorption of fouling material on the membrane surface or in the pore walls. Although many techniques have been developed to overcome fouling, studies on membrane cleaning still seem to be insufficient for practical membrane filtration systems. Current membrane cleaning technologies include hydraulic, chemical, and mechanical methods. Ultrasound （US） has been widely used as a method of cleaning materials because of the cavitation phenomenon. In this study, US cleaning technique was applied to removing the fouling of polyvinglidenefluoride （PVDF） MF membrane, which was used to treat yeast cell and isolated soybean protein （ISP） solution, respectively. The US employed 40 kHz frequency and the output power of 1.43-2.85 W/cm^2. The evolution of the cleaning effect is followed by the measurements of the flux recovery rate （FRstat） and the cleaning time cycle. Results showed that the membrane property, which was fouled by yeast cell solution, could be recovered by water cleaning with US irradiation. And the cleaning time with the same FRstat decreased with the increase of US intensity.     

Exploring the spatial and temporal dynamics of land use in Xizhuang watershed of Yunnan, southwest China

Over the last two decades, China has introduced a series of agricultural and forestland use reforms, aiming to feed the largest population in the world and maintain ecological services locally and nationally. This paper studies the impacts of local government-driven reforestation on land use and land cover change, as well as its further impacts on livelihoods of upland farmers in Xizhuang watershed. An analysis of aerial photographs and ASTER satellite imagery from 1987 to 2002, respectively, showed that the forest has significantly increased at the expense of decreasing farmland. However, the monoculture reforestation of pine has caused both biophysical and socio-economic consequences. This case study also shows forestry decentralization in China remains incomplete. Land use and land cover change is also a political economic issue. Some of the reforms designed to protect forest resources have had a negative impact on rural livelihoods.