Application of remote sensing data in characterization and mapping of soil resources for watershed planning in arid Western Rajasthan

Visual interpretation of IRS-L1SS-II (January, 1995) FCC (1:50,000 scale) of spectral bands 2, 3 and 4 was carried out for the identification and mapping of major physiographic units in an arid watershed of Jodhpur district (Rajasthan). Based on image characteristics and field traverses, seven major physiographic units identified are (1) hills (2) pediments, flat to undulating (3) buried pediments, moderately deep to deep, coarse textured (4) buried pediment, shallow to moderately deep and deep, medium to fine textured, saline (5) older alluvial plains, deep and very deep, coarse textured (6) younger alluvial plains, deep to very deep, very coarse textured and (7) dune complexes. Based on physiographicvariatton and soil or site characteristics such as texture, depth, slope, erosion and underneath substrata, 41 soil mapping units were identified and mapped. Final physiography, soil, slope, drainage and landuse maps were prepared on 1:5,000 scale. Taxonomically, the soils of the watershed were classified as Para-Lithic Torriorthents, coarse-loamy, Lithic/Typic Haplocambids, fine-loamy, Lithic/Typic Haplosalids and Typic Torrifluvents and Typic Torripsamments. Land suitability for various mapping units in the watershed have been assessed on the basis of soil physico-chemical characteristics.     

Non-stationary Approaches for Mapping Terrain and Assessing Prediction Uncertainty

It is well known that terrain may vary markedly over small areas and that statistics used to characterise spatial variation in terrain may be valid only over small areas. In geostatistical terminology, a non-stationary approach may be considered more appropriate than a stationary approach. In many applications, local variation is not accounted for sufficiently. This paper assesses potential benefits in using non-stationary geostatistical approaches for interpolation and for the assessment of uncertainty in predictions with implications for sampling design. Two main non-stationary approaches are employed in this paper dealing with (1) change in the mean and (2) change in the variogram across the region of interest. The relevant approaches are (1) kriging with a trend model (KT) using the variogram of residuals from local drift and (2) locally-adaptive variogram KT, both applied to a sampled photogrammetrically derived digital terrain model (DTM). The fractal dimension estimated locally from the double-log variogram is also mapped to illustrate how spatial variation changes across the data set. It is demonstrated that estimation of the variogram of residuals from local drift is worthwhile in this case for the characterisation of spatial variation. In addition, KT is shown to be useful for the assessment of uncertainty in predictions. This is shown to be true even when the sample grid is dense as is usually the case for remotely-sensed data. In addition, both ordinary kriging (OK) and KT are shown to provide more accurate predictions than inverse distance weighted (IDW) interpolation, used for comparative purposes.     

Somigliana–Pizzetti gravity: the international gravity formula accurate to the sub-nanoGal level

 The Somigliana–Pizzetti gravity field (the International gravity formula), namely the gravity field of the level ellipsoid (the International Reference Ellipsoid), is derived to the sub-nanoGal accuracy level in order to fulfil the demands of modern gravimetry (absolute gravimeters, super conducting gravimeters, atomic gravimeters). Equations (53), (54) and (59) summarise Somigliana–Pizzetti gravity Γ(φ,u) as a function of Jacobi spheroidal latitude φ and height u to the order ?(10⁻¹⁰ Gal), and Γ(B,H) as a function of Gauss (surface normal) ellipsoidal latitude B and height H to the order ?(10⁻¹⁰ Gal) as determined by GPS (`global problem solver'). Within the test area of the state of Baden-Württemberg, Somigliana–Pizzetti gravity disturbances of an average of 25.452 mGal were produced. Computer programs for an operational application of the new international gravity formula with (L,B,H) or (λ,φ,u) coordinate inputs to a sub-nanoGal level of accuracy are available on the Internet. Received: 23 June 2000 / Accepted: 2 January 2001     

Local geoid determination combining gravity disturbances and GPS/levelling: a case study in the Lake Nasser area, Aswan, Egypt

 The use of GPS for height control in an area with existing levelling data requires the determination of a local geoid and the bias between the local levelling datum and the one implicitly defined when computing the local geoid. If only scarse gravity data are available, the heights of new data may be collected rapidly by determining the ellipsoidal height by GPS and not using orthometric heights. Hence the geoid determination has to be based on gravity disturbances contingently combined with gravity anomalies. Furthermore, existing GPS/levelling data may also be used in the geoid determination if a suitable general gravity field modelling method (such as least-squares collocation, LSC) is applied. A comparison has been made in the Aswan Dam area between geoids determined using fast Fourier transform (FFT) with gravity disturbances exclusively and LSC using only the gravity disturbances and the disturbances combined with GPS/levelling data. The EGM96 spherical harmonic model was in all cases used in a remove–restore mode. A total of 198 gravity disturbances spaced approximately 3 km apart were used, as well as 35 GPS/levelling points in the vicinity and on the Aswan Dam. No data on the Nasser Lake were available. This gave difficulties when using FFT, which requires the use of gridded data. When using exclusively the gravity disturbances, the agreement between the GPS/levelling data were 0.71 ± 0.17 m for FFT and 0.63 ± 0.15 for LSC. When combining gravity disturbances and GPS/levelling, the LSC error estimate was ±0.10 m. In the latter case two bias parameters had to be introduced to account for a possible levelling datum difference between the levelling on the dam and that on the adjacent roads. Received: 14 August 2000 / Accepted: 28 February 2001     

New analytical and numerical approaches for geopotential modeling

 The standard analytical approach which is applied for constructing geopotential models OSU86 and earlier ones, is based on reducing the boundary value equation to a sphere enveloping the Earth and then solving it directly with respect to the potential coefficients _n,m . In an alternative procedure, developed by Jekeli and used for constructing the models OSU91 and EGM96, at first an ellipsoidal harmonic series is developed for the geopotential and then its coefficients _n,m ^e are transformed to the unknown _n,m . The second solution is more exact, but much more complicated. The standard procedure is modified and a new simple integral formula is derived for evaluating the potential coefficients. The efficiency of the standard and new procedures is studied numerically. In these solutions the same input data are used as for constructing high-degree parts of the EGM96 models. From two sets of _n,m (n≤360,|m|≤n), derived by the standard and new approaches, different spectral characteristics of the gravity anomaly and the geoid undulation are estimated and then compared with similar characteristics evaluated by Jekeli's approach (`etalon' solution). The new solution appears to be very close to Jekeli's, as opposed to the standard solution. The discrepancies between all the characteristics of the new and `etalon' solutions are smaller than the corresponding discrepancies between two versions of the final geopotential model EGM96, one of them (HDM190) constructed by the block-diagonal least squares (LS) adjustment and the other one (V068) by using Jekeli's approach. On the basis of the derived analytical solution a new simple mathematical model is developed to apply the LS technique for evaluating geopotential coefficients. Received: 12 December 2000 / Accepted: 21 June 2001     

Primary production and fluxes of organic carbon to the seabed in the Eurasian arctic seas, 2003–2012

The primary production and fluxes of organic matter to the seabed and their variations were estimated in the Greenland, Barents, Kara, Laptev, East Siberian, and Chukchi seas in 2003–2012 on the basis of satellite and in situ data. When counting the open water area with the assumptions made for the assessment of the primary production in the sea areas hidden under clouds, an increase in primary production was recorded in all these seas, the total production (phytoplankton and ice algae) has grown from 250 × 10⁶ to 400 × 10⁶ t of C per year over the last ten years. The calculation of the OM flux to the seabed showed growth for certain seas from 4 to 12% per year.     

Sediment transfer in an extremely low-gradient,low-relief and highly buffered system: Darwin Harbour catchment,northern Australia

Abstract

Sediment yields from and sediment transfer within catchments of very low relief and gradient, which make up about 50% of Earth’s surface, are poorly documented and their internal sediment dynamics are poorly known. Sediment sources, their proportionate contributions to valley floors and sediment yield, and storage are estimated using fallout radionuclides ²¹⁰Pb_(ex) and ¹³⁷Cs in the catchments that drain into Darwin Harbour, northern Australia, an example of this understudied catchment type that appears to be globally at the extreme end of this category of catchments. Unchannelled grassy valley floors (dambos, or seasonal wetlands) trap ～90% of the sediment delivered from hillslopes by sheet and rill erosion. Further down valley, small channels transport ～10% of the sediment that escapes from the dambos, and the remaining sediment comes from erosion of the channels. In this case, the fractional sediment storage is very high as a result of the existence of dambos, a landform that depends for its existence on low gradients.