Recent Retreat of Wilkins Ice Shelf Reveals New Insights in Ice Shelf Breakup Mechanisms

The disintegration of various ice shelves on the Antarctic Peninsula has demonstrated their vulnerability and impacts on tributary glaciers. A satellite image of Wilkins Ice Shelf (WIS) from July of 2007 reveals the formation of a large new double fracture, accompanied by numerous small fractures. We show that bending stresses induced by buoyancy forces were responsible for fracture formation. On February 28-29, 2008, an area of about 425 km² broke up at a narrow connection of the WIS to one of its confining islands. In contrast to Larsen B Ice Shelf, melt ponds that drain into crevasses played no role in this breakup process. A further breakup of 160 km² in the same area occurred on May 30-31, 2008 and documented that breakup can occur during austral winter. Radar images reveal a frozen surface, which demonstrates that in this breakup, surface melt water did not play a role. We conclude that ice shelves with strong thickness contrasts carry potential for disintegration. The fact that the WIS experienced two breakup events under two widely contrasting surface conditions (one during the melt season and one during winter) reveals that there may be several reasons for the disintegration of ice shelves that operate under differing circumstances.     

Characterization of diurnal cycles in ZTD from a decade of global GPS observations

The diurnal cycle of the tropospheric zenith total delay (ZTD) is one of the most obvious signals for the various physical processes relating to climate change on a short time scale. However, the observation of such ZTD oscillations on a global scale with traditional techniques (e.g. radiosondes) is restricted due to limitations in spatial and temporal resolution. Nowadays, the International GNSS Service (IGS) provides an important data source for investigating the diurnal and semidiurnal cycles of ZTD and related climatic signals. In this paper, 10 years of ZTD data from 1997 to 2007 with a 2-hour temporal resolution are derived from global positioning system (GPS) observations taken at 151 globally distributed IGS reference stations. These time series are used to investigate diurnal and semidiurnal oscillations. Significant diurnal and semidiurnal oscillations of ZTD are found for all GPS stations used in this study. The diurnal cycles (24 hours period) have amplitudes between 0.2 and 10.9 mm with an uncertainty of about 0.5 mm and the semidiurnal cycles (12 h period) have amplitudes between 0.1 and 4.3 mm with an uncertainty of about 0.2 mm. The larger amplitudes of the diurnal and semidiurnal ZTD cycles are observed in the low-latitude equatorial areas. The peak times of the diurnal cycles spread over the whole day, while the peak value of the semidiurnal cycles occurs typically about local noon. These GPS-derived diurnal and semidiurnal ZTD signals are similar with the surface pressure tides derived from surface synoptic pressure observations, indicating that atmospheric tides are the main driver of the diurnal and semidiurnal ZTD variations.     

On the use of iterative re-weighting least-squares and outlier detection for empirically modelling rates of vertical displacement

The proper identification and removal of outliers in the combination of rates of vertical displacements derived from GPS, tide gauges/satellite altimetry, and GRACE observations is presented. Outlier detection is a necessary pre-screening procedure in order to ensure reliable estimates of stochastic properties of the observations in the combined least-squares adjustment (via rescaling of covariance matrices) and to ensure that the final vertical motion model is not corrupted and/or distorted by erroneous data. Results from this study indicate that typical data snooping methods are inadequate in dealing with these heterogeneous data sets and their stochastic properties. Using simulated vertical displacement rates, it is demonstrated that a large variety of outliers (random scattered and adjacent, as well as jointly influential) can be dealt with if an iterative re-weighting least-squares adjustment is combined with a robust median estimator. Moreover, robust estimators are efficient in areas weakly constrained by the data, where even high quality observations may appear to be erroneous if their estimates are largely influenced by outliers. Four combined models for the vertical motion in the region of the Great Lakes are presented. The computed vertical displacements vary between  − 2 mm/year (subsidence) along the southern shores and 3 mm/year (uplift) along the northern shores. The derived models provide reliable empirical constraints and error bounds for postglacial rebound models in the region.     

Land suitability assessment for selected crops in Vellore district based on agro-ecological characterisation

Planning based on agro-ecological zoning aims at scientific management of regional resources to meet the food, fibre, fodder and fuel wood requirements without adversely affecting the status of natural resources and environment. An attempt has been made to map the agro-ecological units for Vellore district of Tamil Nadu and derive the crop-zone map for the four major crops namely, paddy, sugarcane, groundnut and millets. The basic theory of FAO framework for Land Evaluation was adopted to define the suitability of crops. Land quality details necessary for evaluating the agro-land suitability of crops and for delineating the agroecological units include the terrain, soil and climatic characteristics. Agro-ecological units map was generated by overlaying the agro-edaphic and agroclimatic map layers in GIS. The agro-land suitability map was generated by matching the crop requirement details with the land qualities. The results of the suitability evaluation, when compared with the current land use statistics of these crops showed that area cultivated is less than the area suitable for these crops.     

Automatic main road extraction from high resolution satellite imageries by means of particle swarm optimization applied to a fuzzy-based mean calculation approach

Manual extraction of road network by human operator is an expensive and time-consuming procedure. Alternatively, automation of the extraction process would be a great advancement. For this purpose, an automatic method is proposed to extract roads from high resolution satellite images. In this study, using few samples from road surface, a particle swarm optimization is applied to a fuzzy-based mean calculation system to obtain road mean values in each band of high resolution satellite colour images. Then, the images are segmented using the calculated mean values from the fuzzy system. Optimizing the fuzzy cost function by particle swarm optimization enables the fuzzy approach to be the best mean value of road with sub-grey level precision. Initially, this method was applied to simulated images where the calculated mean values are consistent with the hypothetic mean values. Application of the method to IKONOS satellite images has shown a prospective outcome for automatic road extraction. Mathematical morphology is subsequently used to extract an initial main road centreline from the segmented image. Then, small redundant segments are automatically removed. The quality of the extracted road centreline indicates the effectiveness of the proposed approach.     

Classical globally reflected gravity field determination in modern locally oriented multiscale framework

The purpose of this paper is the canonical connection of classical global gravity field determination following the concept of Stokes (Trans Camb Philos Soc 8:672–712, 1849), Bruns (Die Figur der Erde, Publikation Königl. Preussisch. Geodätisches Institut, P. Stankiewicz Buchdruckerei, Berlin, 1878), and Neumann (Vorlesungen über die Theorie des Potentials und der Kugelfunktionen. Teubner, Leipzig, pp 135–154, 1887) on the one hand and modern locally oriented multiscale computation by use of adaptive locally supported wavelets on the other hand. The essential tools are regularization methods of the Green, Neumann, and Stokes integral representations. The multiscale approximation is guaranteed simply as linear difference scheme by use of Green, Neumann, and Stokes wavelets. As an application, gravity anomalies caused by plumes are investigated for the Hawaiian and Iceland areas.     

Flood inundation modeling using MIKE FLOOD and remote sensing data

A coupled 1D-2D hydrodynamic model, MIKE FLOOD was used to simulate the flood inundation extent and flooding depth in the delta region of Mahanadi River basin in India. Initially, the 1D model MIKE 11 was calibrated using river water level and discharge data of various gauging sites for the monsoon period (June to October) of the year 2002. Subsequently, the calibrated set up was validated using both discharge and water level data for the same period of the year 2001. The performance of calibration and validation results of MIKE 11 were evaluated using different performance indices. A bathymetry of the study area with a spatial resolution of 90m was prepared from SRTM DEM and provided as an input to the 2D model, MIKE 21. MIKE 11 and MIKE 21 models were then coupled using lateral links to form the MIKE FLOOD model set up for simulating the two dimensional flood inundations in the study area. Flood inundation is simulated for the year 2001 and the maximum flood inundation extent simulated by the model was compared with the corresponding actual inundated area obtained from IRS-1D WiFS image.     

The IGS VTEC maps: a reliable source of ionospheric information since 1998

The International GNSS Service (IGS) Working Group on Ionosphere was created in 1998. Since then, the Scientific community behind IGS, in particular CODE, ESA, JPL and UPC, have been continuosly contributing to reliable IGS combined vertical total electron content (VTEC) maps in both rapid and final schedules. The details on how these products are being generated, performance numbers, proposed improvement as far as VTEC evolution trends during near one Solar Cycle, are summarized in this paper. The confirmation of (1) the good performance of the IGS combined VTEC maps, and (2) the characteristic VTEC variability periods, are two main results of this work.     

CPLat: first operational experimental processing center for SIRGAS in Argentina

The SIRGAS permanent GPS network which is in fact the IGS network densification for the American continent, consists today of more than 200 stations covering the continent and islands. It is currently processed by the IGS RNAAC SIR centre at Deutsches Geodätisches Forschungsinstitut producing weekly free solutions relying on IGS final orbits and EOP that contribute to the ITRF through IGS. By August 2006, the SIRGAS Working Group I had accepted five proposals for experimental processing centers within the region that would collaborate with IGS RNAAC SIR. One of them, Centro de Procesamiento La Plata (CPLat) in Argentina, began processing 60 stations on October 2006. By January 2007 CPLat reached operational capability, delivering weekly free solution SINEX files, with an internal consistency of 1.5 mm average for the horizontal components, and 3 mm in the vertical. Comparisons with IGS global and IGS RNAAC SIR weekly solutions were taken as external consistency indications, showing average RMS residuals of 1.8, 2.4 and 5 mm for the north, east, and vertical component, respectively. Analysis and comparison of adjusted solution time series from CPLat and other processing centers has proved to be highly valuable for solution QC, namely detection and identification of station anomalous behavior or modelling problems. These procedures will ensure the maintenance of the performance specifications for CPLat solutions. Action is being taken in order to guarantee the continuity of this effort beyond the experimental phase.