Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition   总被引：3，自引：0，他引：3  

Oliver Baur  Gerrit Austen  Jürgen Kusche 《Journal of Geodesy》2008,82(4-5):207-221

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     

Wenbin Shen  Jin Li  Jiancheng Li  Zhengtao Wang  Jinsheng Ning  Dingbo Chao 《地球空间信息科学学报》2008,11(4):273-278

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     

M. Gerritsen  J. V. Lambers 《Computational Geosciences》2008,12(2):193-208

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.  相似文献   

Environmental geochemistry of the River Gomti: A tributary of the Ganges River     

L. P. Gupta  V. Subramanian 《Environmental Geology》1994,24(4):235-243

Water and sediment samples collected from the Gomti River, a tributary of the Ganges River system, during the postmonsoon season have been analyzed to estimate major elemental chemistry. Water chemistry of the River Gomti shows almost monotonous spatial distribution of various chemical species, especially because of uniform presence of alluvium Dun gravels throughout the basin. The river annually transports 0.34×10⁶ tonnes of total suspended material (TSM) and 3.0×10⁶ tonnes of total dissolved solids (TDS), 69 percent of which is accounted for by bicarbonate ions only. Samples collected downstream of the city of Lucknow show the influence of anthropogenic loadings for a considerable distance in the river water. Na⁺, Cl^–, and SO₄ ^2– concentrations build up downstream. The bed sediment chemistry is dominated by Si (36 percent), reflecting a high percentage of detrital quartz, which makes up about 74 percent of the mineralogy of the bed sediments in the River Gomti. The average Kjeldahl nitrogen concentration (234 g/g) indicates indirectly the amount of organic matter in the sediments. The Hg concentration in sediments has been found to be higher (average 904 ppb) than the background value. The suspended sediments are well sorted, very finely skewed, and extremely leptokurtic, indicating a low energy condition of flow in the Gomti River. The influence of chemical loads in the Gomti has been found to be small or nonexistent on the Ganges River, perhaps because the water discharge of the Gomti (1.57 percent) to the Ganges is quite low.  相似文献   

Sediment budget-based estimates of trace metal inputs to a Chesapeake estuary     

W. A. Marcus  C. C. Nielsen  J. C. Cornwell 《Environmental Geology》1993,22(1):1-9

This article evaluates whether a sediment budget for the South River, Maryland, can be coupled with metals data from sediment cores to identify and quantify sources of historic metal inputs to marsh and subtidal sediments along the estuary. Metal inputs to estuarine marsh sediments come from fluvial runoff and atmospheric deposition. Metal inputs to subtidal sediments come from atmospheric deposition, fluvial runoff, coastal erosion, and estuarine waters. The metals budget for the estuary indicates that metal inputs from coastal erosion have remained relatively constant since 1840. Historical variations in metal contents of marsh sediments have probably resulted primarily from increasing atmospheric deposition in this century, but prior to 1900 may reflect changing fluvial sources, atmospheric inputs, or factors not quantified by the budget. Residual Pb, Cu, and Zn in the marsh sediments not accounted for by fluvial inputs was low to moderate in 1840, decreased to near zero circa 1910, and by 1987 had increased to levels that were one to ten times greater than those of 1840. Sources of variability in subtidal cores could not be clearly discerned because of geochemical fluxes, turbulent mixing, and bioturbation within the cores. The sediment-metal budgeting approach appears to be a viable method for delineating metal sources in small, relatively simple estuarine systems like the South River and in systems where recent deposition (for example, prograding marshes) prevents use of deep core analysis to identify background levels of metal. In larger systems or systems with more variable sources of sediment and metal input, however, assumptions and measurement errors in the metal budgeting approach suggest that deep core analysis and normalization techniques are probably preferable for identifying anthropogenic impacts.Field and laboratory research conducted at the Department of Geography, University of Maryland, College Park, Maryland, 20742, USAField and laboratory research conducted at the Marine and Estuarine Environmental Science Program, University of Maryland, College Park, Maryland, 20742, USA  相似文献   

Hydrometallurgical processing and recovery of molybdenum trioxide from spent catalyst   总被引：5，自引：0，他引：5  

Kyung Ho Park  B. Ramachandra Reddy  D. Mohapatra  Chul-Woo Nam 《International Journal of Mineral Processing》2006,80(2-4):261-265

Hydrometallurgical processing of spent hydrodesulphurisation (HDS) catalyst for the recovery of molybdenum using sodium carbonate and hydrogen peroxide mixtures was investigated. The results indicated that the recovery of molybdenum was largely dependent on the concentrations of Na₂CO₃ and H₂O₂ in the reaction medium, which controls the acidity of the leach liquor and carry over of impurities such as Al, Ni, P, Si and V. Leaching process was exothermic and leaching efficiency of molybdenum decreased with increasing solid to liquid ratio. Large scale leaching of spent catalyst, under optimum conditions: 20% pulp density, 85 g/L Na₂CO₃, 10 vol.% H₂O₂ and 1 h reaction, resulted a leaching efficiency of 84% Mo. The obtained leach liquor contained (g/L): Mo — 22.0, Ni — 0.015 and Al — 0.82, P — 1.1, Si — 0.094 and minor quantities of V — 8 mg/L, As and Co — < 1 mg/L. Recovery of Mo from leach solution as MoO₃ of 97.30% purity was achieved by ammonium molybdate precipitation method.  相似文献   

The research of phosphorus of Xiangxi River nearby the Three Gorges, China     

Chang-Ying Fu  Tao Fang  Nan Sheng Deng 《Environmental Geology》2006,49(6):923-928

The Xiangxi River is the first middling tributary of the Changjiang River near the Three Gorges Dam. The River is subject to phosphorus pollution mainly from industrial wastewater. As the water quality of the Xiangxi River could directly influence the water quality of the Three Gorges Reservoir, the research on phosphorus levels and its change in the sediment profile of the Xiangxi River could provide useful information in the dynamic changes in the system, thereby offering options for mitigative measures. Water and sediment samples from lower reaches of Xiangxi River were collected and the different forms of phosphorus in sediments of the Xiangxi River were studied. The concentrations of total phosphorus in sediment ranged from 757.67 to 1438.54 mg/kg. Inorganic phosphorus concentrations ranged from 684.63 to 1055.58 mg/kg. Phosphorus contamination was serious in some parts of the Xiangxi River. With an average concentration of 635.17 mg/kg, calcium-bound phosphorus is the main form among different inorganic phosphorus forms. Labile phosphorus and iron/aluminum-bound phosphorus measured 3.40, 0.05and 35.28 mg/kg, respectively. The mobilization potential of phosphorus of sediments was studied through adsorption and release experiments. The equilibrium concentration of phosphorus adsorption and release was around 0.1 mg/L. The initial concentrations of phosphorus in the overlying water and the sediments have obvious effect on phosphorus mobilization potential. In addition, the release rate of phosphorus in sediment increased with water depth.  相似文献   

Factors influencing the presence or absence of tributary-junction fans in the Iberian Range, Spain     

Amelia Gmez-Villar  Javier lvarez-Martínez  Jos M. García-Ruiz 《Geomorphology》2006,81(3-4):252-264

This paper analyses the factors which influence the presence or absence of tributary-junction fans in the Iberian Range, northern Spain. Two valleys were selected, both characterised by wide variations in lithology, altitude, land use and plant cover. Two groups of factors were studied: those related to the internal characteristics of the drainage basins, which particularly control sediment generation; and those related to the characteristics of the depositional area which control accommodation space and main river power. Among the internal factors, the development of alluvial fans was related to: (i) the capacity of the basin to yield large volumes of sediment, (ii) the occurrence of intense human pressure until recent times, a good indicator of sediment yield, and (iii) the capacity of the basin to quickly increase discharge during rainstorms (discharge density and torrentiality). It is suggested that the areas that were intensively cultivated in the past, and have therefore been affected by intense erosion, have played a decisive role on the development of alluvial fans. This would imply that many of these alluvial fans have a relatively recent origin, perhaps related to the beginning of a widespread deforestation. The basins without alluvial fans are characterised by relatively steep hillslope gradients (that is, slopes that never were subjected to historical cultivation), low drainage densities and dense forest and shrub cover, mostly coinciding with high altitude basins composed of quartzite and shale bedrocks. Regarding the external factors, the shape, size and longitudinal gradient of the main river to which the fans are tributary are the most relevant conditioning factors determining the development of alluvial fans.  相似文献   

The Record of an Extreme Flood in the Sediments of Montane Lillooet Lake, British Columbia: Implications for Paleoenvironmental Assessment     

Robert Gilbert  Sarah Crookshanks  Kyle R. Hodder  John Spagnol  Roland B. Stull 《Journal of Paleolimnology》2006,35(4):737-745

Severe rainfall in mid October, 2003 produced the largest floods in almost a century of record on rivers in the Cordillera of southwestern British Columbia. Sediment deposited in Lillooet Lake as a result of this event is clearly distinguished by stratigraphy, colour, texture, magnetic properties, and organic content. Each of these physical properties is related to the lacustrine processes, especially turbid underflow, that distributed the sediment through the lake. The flood, which lasted less than a week, delivered 8–12 times the amount of sediment that accumulates in most entire years in the deepest, central parts of the lake. Recognition of events of this type in the stratigraphic record offers a means of assessing the changing nature of extreme hydroclimatic events, and their relation to more ubiquitous, lower-energy processes.  相似文献   

A Comparison of Sediment Varves (1950–2003 AD) in Two Adjacent Lakes in Northern Sweden     

V. Gälman  G. Petterson  I. Renberg 《Journal of Paleolimnology》2006,35(4):837-853

Koltjärnen and Nylandssjön are two closely situated lakes (<2 km apart) in northern Sweden. During the past century, distinct varved sediments have formed in these lakes. Nylandssjön has two varved, deep basins. Since lake and catchment characteristics superficially appear very similar for the two lakes and they are exposed to the same climate, one would expect the sediment varves to be similar. This investigation compares the varves in the two deep basins of Nylandssjön, and in the two lakes. The comparison of basins of Nylandssjön shows that varve thickness, water content and annual accumulation rates of organic matter and nitrogen are correlated for the period (1970–2003). The grey-scale curves are only clearly similar in about 50% of the varves. In the between-lake comparison varve thickness, water content and annual accumulation rates of organic matter and nitrogen are correlated for the period (1950–1996). However, the annual accumulation rates of dry mass, minerogenic matter and biogenic silica differ between the lakes, as well as within-varve successions in grey-scale. A general explanation to the differences is that the prerequisites for varve formation are not totally similar because of differences in catchment size, catchment- to-lake material fluxes, lake productivity and land-use influence. This study illustrates the complex relationships that exist between a lake, its catchment, in-lake productivity and formation of sediment varves. As a consequence, we must not apply a too simplistic view of the potential of varves as past climatic indicators, especially if the lakes are affected by land-use.  相似文献