Deforestation and fragmentation of seasonal tropical forests in the southern Yucatán,Mexico (1990–2006)

The southern Yucatán (SY) has been recognized as a hotspot of biodiversity with great risk of deforestation. Land change analysis, based on classified Landsat TM and ETM?+?satellite imagery (1990, 2000 and 2006), was used to estimate the annual deforestation rates of 141 land management units of the SY, and spatial patterns of forest fragmentation around and within the Calakmul Biosphere Reserve (CBR), which comprises approximately one-third of the region. Results indicate a decrease in annual deforestation rates over 1990–2006, from 0.15% year^?1 to 0.06% year^?1, but with significant sub-regional variations in the quantity and rate of forest loss. Despite a decline in deforestation during this period, there was considerable fragmentation both inside and outside the CBR. While population pressures and the expansion of pasture have caused deforestation across the region, agricultural intensification, diversified income strategies and reserve conservation may have contributed to reduced forest loss during the study period.     

Spherical integral transforms of second-order gravitational tensor components onto third-order gravitational tensor components

New spherical integral formulas between components of the second- and third-order gravitational tensors are formulated in this article. First, we review the nomenclature and basic properties of the second- and third-order gravitational tensors. Initial points of mathematical derivations, i.e., the second- and third-order differential operators defined in the spherical local North-oriented reference frame and the analytical solutions of the gradiometric boundary-value problem, are also summarized. Secondly, we apply the third-order differential operators to the analytical solutions of the gradiometric boundary-value problem which gives 30 new integral formulas transforming (1) vertical-vertical, (2) vertical-horizontal and (3) horizontal-horizontal second-order gravitational tensor components onto their third-order counterparts. Using spherical polar coordinates related sub-integral kernels can efficiently be decomposed into azimuthal and isotropic parts. Both spectral and closed forms of the isotropic kernels are provided and their limits are investigated. Thirdly, numerical experiments are performed to test the consistency of the new integral transforms and to investigate properties of the sub-integral kernels. The new mathematical apparatus is valid for any harmonic potential field and may be exploited, e.g., when gravitational/magnetic second- and third-order tensor components become available in the future. The new integral formulas also extend the well-known Meissl diagram and enrich the theoretical apparatus of geodesy.     

Spin rate estimation of sounding rockets using GPS wind-up

Carrier phase wind-up is a well-known effect that arises from the relative rotation between a transmitting and receiving antenna. In GPS measurements at L1 frequency, this effect translates into an error of 19.029 cm per full relative rotation of antennas. Since this effect is independent of the satellite elevation for pure rotation about the antenna boresight axis, it is usually absorbed by the clock estimation in navigation algorithms. Therefore, the impact of wind-up is usually neglected for applications that do not require accuracies to the cm level like RTK. However, in receiving platforms with high rotation rate, the accumulated wind-up value can be important and actually be larger than receiver noise or even ionospheric variations. Therefore, in such scenarios, the wind-up contribution can be isolated and used as a source of information to compute the spin rate of such platforms using an appropriate combination of GPS observables. This work shows some results of a coarse, yet simple, approach to monitor the rotation angle and spin-rate of spin stabilized sounding rockets flown by DLR.     

The ionosphere: effects,GPS modeling and the benefits for space geodetic techniques

The main goal of this paper is to provide a summary of our current knowledge of the ionosphere as it relates to space geodetic techniques, especially the most informative technology, global navigation satellite systems (GNSS), specifically the fully deployed and operational global positioning system (GPS). As such, the main relevant modeling points are discussed, and the corresponding results of ionospheric monitoring are related, which were mostly computed using GPS data and based on the direct experience of the authors. We address various phenomena such as horizontal and vertical ionospheric morphology in quiet conditions, traveling ionospheric disturbances, solar flares, ionospheric storms and scintillation. Finally, we also tackle the question of how improved knowledge of ionospheric conditions, especially in terms of an accurate understanding of the distribution of free electrons, can improve space geodetic techniques at different levels, such as higher-order ionospheric effects, precise GNSS navigation, single-antenna GNSS orientation and real-time GNSS meteorology.     

Synergetic efficiency of Lidar and WorldView-2 for 3D urban cartography in Northeast Mexico

Three-dimensional urban cartography is needed for city changes’ assessment. The variety of studies using 3D calculations of urban elements grows each year. Building and vegetation volumes are necessary to assess and understand spatio-temporal urban changeable environments. However, there are technical questions as to which method can improve 3D urban cartographic accuracy. The innovative part of this current study is the creation of a six-band hybrid obtained from LIDAR and WorldView2 synergy. Two different enhancement algorithms demonstrated the most important spectral features for the urban development and vegetation classes. Results indicated an improvement in accuracy by up to 21.3%, according to the Kappa coefficient. Both infra-red band and intensity band were the most significant, according to the principal components analysis. The synergy delimited classes and polygons, as well as the direct display of information regarding heights of elements and improving the extraction of roads, buildings and vegetation classes.     

A calibration system for superconducting gravimeters

A new method for the calibration of a superconducting gravity meter is described, in which a 273 Kg annular mass is placed around the meter and is moved up and down. The geometry of the apparatus is easy to model and the accuracy in the computation of the gravity variation induced by the mass, 6.7µgal, is limited only by the accuracy in the knowledge of value of the gravitational constant. Measurements done in 91 and 92 for the calibration of the instrument GWR-T015 are described. The calibration factor has been determined with a precision of about 0.3%.     

Crustal structure of Iraq from receiver functions and surface wave dispersion: implications for understanding the deformation history of the Arabian–Eurasian collision

We report the crustal structure for two locations in Iraq estimated by joint inversion of P -wave receiver functions (RFs) and surface (Rayleigh) wave group velocity dispersion. RFs were computed from teleseismic recordings at two temporary broad-band seismic stations located in Mosul (MSL) in the Zagros Fold Belt and Baghdad (BHD) in the Mesopotamian Foredeep. Group velocity dispersion curves at the sites were derived from continental-scale tomography. The inversion results show that the crustal thicknesses are 39 km at MSL and 43 km at BHD. We observe a strong Ps _Moho at BHD consistent with a sharp Moho discontinuity. However, at MSL we observe a weak Ps _Moho suggesting a transitional Moho where crustal thickening is likely to be occurring in the deep crust. Both sites reveal low velocity surface layers consistent with sedimentary thickness of about 3 km at station MSL and 7 km at BHD and agreeing well with the previous reports. Ignoring the sediments, the crystalline crustal velocities and thicknesses are remarkably similar at both stations. The similarity of crustal structure suggests that the crust of the northeastern proto-Arabian Platform was uniform before subsidence and deposition of the sediments in the Cenozoic. If crystalline crustal structure is uniform across the northern Arabian Platform then crustal thickness variations in the Zagros Fold Belt and Thrust Zone should reveal the history of deformation and crustal shortening in the Arabian–Eurasian collision zone and not reflect pre-existing crustal thickness variations in the Arabian Plate.     

Spatial and temporal patterns of dry spell lengths in the Iberian Peninsula for the second half of the twentieth century

Summary Daily pluviometric records of 43 meteorological stations across the Iberian Peninsula have permitted a detailed analysis of dry spell patterns for the period 1951–2000 by distinguishing daily amount thresholds of 0.1, 1.0, 5.0 and 10.0 mm/day. The analyses are based on three annual series, namely the number of dry spells, N, the average dry spell length, L, and the extreme dry spell length, L _max. First, the statistical significance of local trends for the annual series of N, L and L _max has been investigated by means of the Mann-Kendall test and significant field trends have been established by means of Monte Carlo simulations. Clear signs of negative field trends are detected for N (1.0 and 10.0 mm/day) and L (0.1 mm/day). Second, the Weibull model fits well the empirical distributions of dry spell lengths for all the rain gauges, whatever the daily amount threshold, with a well ranged spatial distribution of their parameters u and k. On the basis of the Weibull distribution, return period maps for 2, 5, 10, 25 and 50 years have been obtained for dry spell lengths with respect to the four daily threshold levels. While for 0.1 and 1.0 mm/day the longest dry spells are expected at the south of the Iberian Peninsula, for 5.0 and 10.0 mm/day they are mostly detected at the southeast. Finally, the elapsed time between consecutive dry spells has been analysed by considering the same rain amount thresholds and different dry spell lengths at increasing intervals of 10 days. This analysis makes evident a significant negative field trend of the elapsed time between consecutive dry spells of lengths ranging from 10 to 20 days for daily amount thresholds of 1.0, 5.0 and 10.0 mm/day. Authors’ addresses: X. Lana, C. Serra, Departament de Física i Enginyeria Nuclear, ETSEIB, Universitat Politècnica de Catalunya, Av. Diagonal 647 planta 11, 08028 Barcelona, Spain; M. D. Marínez, Departament de Física Aplicada, Universitat Politècnica de Catalunya, 08028 Barcelona, Spain; A. Burgue?o, Departament de Meteorologia i Astronomia, Universitat de Barcelona, 08028 Barcelona, Spain; J. Martín-Vide, L. Gómez, Grup de Climatologia, Universitat de Barcelona, 08028 Barcelona, Spain.     

Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography

The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth’s topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth’s surface. It is based on an evolution of a surface, which approximates the Earth’s topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth’s surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth’s topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model.     

An open source GIS‐based Planning Support System: Application to the land use plan of La Troncal,Ecuador

Planning Support Systems (PSS) comprise a wide variety of geo‐technological tools related to GIS and spatial modeling aimed at addressing land planning processes. This article describes the OpenRules system, a PSS based on a previous system called RULES. Among OpenRules new features are its architecture, based exclusively on free and open source software, and its applicability to all land use types, including rural and urban uses. In addition, OpenRules incorporates an unlimited number of land evaluation factors and a new objective in land use spatial allocation. OpenRules has been programmed in Java and implemented as a module of the free GIS software gvSIG, with full integration between the GIS and the decision support tools. Decision support tools include multicriteria evaluation, multiobjective linear programming and heuristic techniques, which support three basic stages of land use planning processes, namely land suitability evaluation, land use area optimization and land use spatial allocation. The application of OpenRules to the region of La Troncal, Ecuador, demonstrates its capability to generate alternative and coherent solutions through a scientific and justified procedure at low cost in terms of time and resources.