Insolation, erosion, and morphology of comet 67P/Churyumov-Gerasimenko

Context. The complex shape of comet 67P and its oblique rotation axis cause pronounced seasonal effects. Irradiation and hence activity vary strongly. Aims. We investigate the insolation of the cometary surface in order to predict the sublimation of water ice. The strongly varying erosion levels are correlated with the topography and morphology of the present cometary surface and its evolution. Methods. The insolation as a function of heliocentric distance and diurnal (spin dependent) variation is calculated using >10(5) facets of a detailed digital terrain model. Shading, but also illumination and thermal radiation by facets in the field of view of a specific facet are iteratively taken into account. We use a two-layer model of a thin porous dust cover above an icy surface to calculate the water sublimation, presuming steady state and a uniform surface. Our second model, which includes the history of warming and cooling due to thermal inertia, is restricted to a much simpler shape model but allows us to test various distributions of active areas. Results. Sublimation from a dirty ice surface yields maximum erosion. A thin dust cover of 50 pm yields similar rates at perihelion. Only about 6% of the surface needs to be active to match the observed water production rates at perihelion. A dust layer of 1 mm thickness suppresses the activity by a factor of 4 to 5. Erosion on the south side can reach more than 10 m per orbit at active spots. The energy input to the concave neck area (Hapi) during northern summer is enhanced by about 50% owing to self-illumination. Here surface temperatures reach maximum values along the foot of the Hathor wall. Integrated over the whole orbit this area receives the least energy input. Based on the detailed shape model, the simulations identify "hot spots" in depressions and larger pits in good correlation with observed dust activity. Three-quarters of the total sublimation is produced while the sub-solar latitude is south, resulting in a distinct dichotomy in activity and morphology. Conclusions. The northern areas display a much rougher morphology than what is seen on Imhotep, an area at the equator that will be fully illuminated when 67P is closer to the Sun. Self-illumination in concave regions enhance the energy input and hence erosion. This explains the early activity observed at Hapi. Cliffs are more prone to erosion than horizontal, often dust covered, areas, which leads to surface planation. Local activity can only persist if the forming cliff walls are eroding. Comet 67P has two lobes and also two distinct sides. Transport of material from the south to the north is probable. The morphology of the Imhotep plain should be typical for the terrains of the yet unseen southern hemisphere.

Fast Algorithms Using Minimal Data Structures for Common Topological Relationships in Large, Irregularly-spaced Topographic Data Sets

Digital terrain models (DTM) typically contain large numbers of postings, from hundreds of thousands to billions. Many algorithms that run on DTMs require topological knowledge of the postings, such as finding nearest neighbors, finding the posting closest to a chosen location, etc. If the postings are arranged irregu- larly, topological information is costly to compute and to store. This paper offers a practical approach to organizing and searching irregularly-space data sets by presenting a collection of efficient algorithms (O(N),O(lgN)) that compute important topological relationships with only a simple supporting data structure. These relationships include finding the postings within a window, locating the posting nearest a point of interest, finding the neighborhood of postings nearest a point of interest, and ordering the neighborhood counter-clockwise. These algorithms depend only on two sorted arrays of two-element tuples, holding a planimetric coordinate and an integer identification number indicating which posting the coordinate belongs to. There is one array for each planimetric coordinate (eastings and northings). These two arrays cost minimal overhead to create and store but permit the data to remain arranged irregularly.

Satellite image map of King George Island, Antarctica

A new topographic database for King George Island, one of the most visited areas in Antarctica, is presented. Data from differential GPS surveys, gained during the summers 1997/98 and 1999/2000, were combined with up to date coastlines from a SPOT satellite image mosaic, and topographic information from maps as well as from the Antarctic Digital Database. A digital terrain model (DTM) was generated using ARC/INFO GIS. From contour lines derived from the DTM and the satellite image mosaic a satellite image map was assembled. Extensive information on data accuracy, the database as well as on the criteria applied to select place names is given in the multilingual map. A lack of accurate topographic information in the eastern part of the island was identified. It was concluded that additional topographic surveying or radar interferometry should be conducted to improve the data quality in this area. In three case studies, the potential applications of the improved topographic database are demonstrated. The first two examples comprise the verification of glacier velocities and the study of glacier retreat from the various input data-sets as well as the use of the DTM for climatological modelling. The last case study focuses on the use of the new digital database as a basic GIS (Geographic Information System) layer for environmental monitoring and management on King George Island.

Sistemas multiespectrales e hiperespectrales para la observación del territorio : análisis y aplicación a la prospección de hidrocarburos.

La observación de la Tierra es una herramienta de gran utilidad en la actualidad para el estudio de los fenómenos que se dan en la misma. La observación se puede realizar a distintas escalas y por distintos métodos dependiendo del propósito. El actual Trabajo Final de Grado persigue exponer la observación del territorio mediante técnicas de Teledetección, o Detección Remota, y su aplicación en la exploración de hidrocarburos. Desde la Segunda Guerra Mundial el capturar imágenes aéreas de regiones de la Tierra estaba restringido a usos cartográficos en el sentido estricto. Desde aquellos tiempos, hasta ahora, ha acontecido una serie de avances científicos que permiten deducir características intrínsecas de la Tierra mediante mecanismos complejos que no apreciamos a simple vista, pero que, están configurados mediante determinados parámetros geométricos y electrónicos, que permiten generar series temporales de fenómenos físicos que se dan en la Tierra. Hoy en día se puede afirmar que el aprovechamiento del espectro electromagnético está en un punto máximo. Se ha pasado del análisis de la región del espectro visible al análisis del espectro en su totalidad. Esto supone el desarrollo de nuevos algoritmos, técnicas y procesos para extraer la mayor cantidad de información acerca de la interacción de la materia con la radiación electromagnética. La información que generan los sistemas de captura va a servir para la aplicación directa e indirecta de métodos de prospección de hidrocarburos. Las técnicas utilizadas en detección por sensores remotos, aplicadas en campañas geofísicas, son utilizadas para minimizar costes y maximizar resultados en investigaciones de campo. La predicción de anomalías en la zona de estudio depende del analista, quien diseña, calcula y evalúa las variaciones de la energía electromagnética reflejada o emitida por la superficie terrestre. Para dicha predicción se revisarán distintos programas espaciales, se evaluará la bondad de registro y diferenciación espectral mediante el uso de distintas clasificaciones (supervisadas y no supervisadas). Por su influencia directa sobre las observaciones realizadas, se realiza un estudio de la corrección atmosférica; se programan distintos modelos de corrección atmosférica para imágenes multiespectrales y se evalúan los métodos de corrección atmosférica en datos hiperespectrales. Se obtendrá temperatura de la zona de interés utilizando los sensores TM-4, ASTER y OLI, así como un Modelo Digital del Terreno generado por el par estereoscópico capturado por el sensor ASTER. Una vez aplicados estos procedimientos se aplicarán los métodos directos e indirectos, para la localización de zonas probablemente afectadas por la influencia de hidrocarburos y localización directa de hidrocarburos mediante teledetección hiperespectral. Para el método indirecto se utilizan imágenes capturadas por los sensores ETM+ y ASTER. Para el método directo se usan las imágenes capturadas por el sensor Hyperion. ABSTRACT The observation of the Earth is a wonderful tool for studying the different kind of phenomena that occur on its surface. The observation could be done by different scales and by different techniques depending on the information of interest. This Graduate Thesis is intended to expose the territory observation by remote sensing acquiring data systems and the analysis that can be developed to get information of interest. Since Second World War taking aerials photographs of scene was restricted only to a cartographic sense. From these days to nowadays, it have been developed many scientific advances that make capable the interpretation of the surface behavior trough complex systems that are configure by specific geometric and electronic parameters that make possible acquiring time series of the phenomena that manifest on the earth’s surface. Today it is possible to affirm that the exploitation of the electromagnetic spectrum is on a maxim value. In the past, analysis of the electromagnetic spectrum was carry in a narrow part of it, today it is possible to study entire. This implicates the development of new algorithms, process and techniques for the extraction of information about the interaction of matter with electromagnetic radiation. The information that has been acquired by remote sensing sensors is going to be a helpful tool for the exploration of hydrocarbon through direct and vicarious methods. The techniques applied in remote sensing, especially in geophysical campaigns, are employed to minimize costs and maximize results of ground-based geologic investigations. Forecasting of anomalies in the region of interest depends directly on the expertise data analyst who designs, computes and evaluates variations in the electromagnetic energy reflected or emanated from the earth’s surface. For an optimal prediction a review of the capture system take place; assess of the goodness in data acquisition and spectral separability, is carried out by mean of supervised and unsupervised classifications. Due to the direct influence of the atmosphere in the register data, a study of the minimization of its influence has been done; a script has been programed for the atmospheric correction in multispectral data; also, a review of hyperspectral atmospheric correction is conducted. Temperature of the region of interest is computed using the images captured by TM-4, ASTER and OLI, in addition to a Digital Terrain Model generated by a pair of stereo images taken by ASTER sensor. Once these procedures have finished, direct and vicarious methods are applied in order to find altered zones influenced by hydrocarbons, as well as pinpoint directly hydrocarbon presence by mean of hyperspectral remote sensing. For this purpose ETM+ and ASTER sensors are used to apply the vicarious method and Hyperion images are used to apply the direct method.

Two-dimensional direct rainfall hydraulic models applied to direct calculation of hydrologic events in dams to prevent overtopping

One of the main concerns when conducting a dam test is the acute determination of the hydrograph for a specific flood event. The use of 2D direct rainfall hydraulic mathematical models on a finite elements mesh, combined with the efficiency of vector calculus that provides CUDA (Compute Unified Device Architecture) technology, enables nowadays the simulation of complex hydrological models without the need for terrain subbasin and transit splitting (as in HEC-HMS). Both the Spanish PNOA (National Plan of Aereal Orthophotography) Digital Terrain Model GRID with a 5 x 5 m accuracy and the CORINE GIS Land Cover (Coordination of INformation of the Environment) that allows assessment of the ground roughness, provide enough data to easily build these kind of models

Boston, Massachusetts Region LIDAR First Return Elevation Data, 2002

LIDAR (LIght Detection And Ranging) first return elevation data of the Boston, Massachusetts region from MassGIS at 1-meter resolution. This LIDAR data was captured in Spring 2002. LIDAR first return data (which shows the highest ground features, e.g. tree canopy, buildings etc.) can be used to produce a digital terrain model of the Earth's surface. This dataset consists of 74 First Return DEM tiles. The tiles are 4km by 4km areas corresponding with the MassGIS orthoimage index. This data set was collected using 3Di's Digital Airborne Topographic Imaging System II (DATIS II). The area of coverage corresponds to the following MassGIS orthophoto quads covering the Boston region (MassGIS orthophoto quad ID: 229890, 229894, 229898, 229902, 233886, 233890, 233894, 233898, 233902, 233906, 233910, 237890, 237894, 237898, 237902, 237906, 237910, 241890, 241894, 241898, 241902, 245898, 245902). The geographic extent of this dataset is the same as that of the MassGIS dataset: Boston, Massachusetts Region 1:5,000 Color Ortho Imagery (1/2-meter Resolution), 2001 and was used to produce the MassGIS dataset: Boston, Massachusetts, 2-Dimensional Building Footprints with Roof Height Data (from LIDAR data), 2002 [see cross references].

AWI Bathymetric Chart of the Fram Strait (BCFS) (Scale 1:100,000)

Based on data from R/V Polarstern multibeam sonar surveys between 1984 and 1997 a high resolution bathymetry has been generated for the central Fram Strait. The area ensonified covers approx. 36,500 sqkm between 78°N - 80°N and 0°E - 7.5°E. Basic outcome of the investigation is a Digital Terrain Model (DTM) with 100 m grid spacing which was utilized for contouring and generation of a new series of bathymetric charts at a scale of 1:100,000, the AWI Bathymetric Chart of the Fram Strait (AWI BCFS).

Evolução espaço-temporal da ocupação urbana sobre áreas naturais no município de Natal-RN

In this work the landscape morphodynamics was used to check the strength and importance of the changes carried out by man on the environment over time, in Natal-RN municipality. The occupation of partially preserved natural areas was analyzed, but environmentally fragile, such as riparian forests, vegetation on the banks of waterways, which play regulatory role of the water flow, and the dunes, which guarantee the rapid recharge of aquifers. The impacts of urban sprawl in Natal Southern and West zones Were identified and characterized, through a detailed mapping in the period between 1969 and 2013 the main Permanent Preservation Areas - PPA (banks of rivers and lagoons, and dunes remaining) and their temporal changes. For this were used aerial photographs and satellite imagery, altimetry data, and pre-existing information, which allowed the creation of a spatial database, and evolution of maps of impervious areas, evolution of the use and occupation and Digital Terrain Model (DTM) from contour lines with contour interval of 1 meter. Based on this study presents a diagnosis of the environmental situation and the state of conservation of natural areas, over the last 44 years, compared to human pressures. In general, it was found that the urban settlement has advanced about 60% of studied natural areas. This advance was growing by the year 2006, when there was a slowdown in the process, except for the Environmental Protection Zone (EPZ) 03, where the river Pitimbú and your PPA, which experienced a more significant loss area. The urban occupation affected the natural drainage and contributed to the contamination of groundwater Natal, due to increased sealed area, the release of liquid and solid waste, as well as the removal of riparian vegetation. Changed irreversibly the natural landscape, and reduced the quality and quantity of water resources necessary for the population. Thus, it is necessary to stimulate the adoption of use and protection of PPA planning measures, to the preservation of the San Valley Region inserted into the EPZ 01, and integrate more remaining dunes, in good condition, this EPZ, due to the importance of those remaining on the environment and the maintenance of quality of life. It is suggested, also, protection of catchment areas, such as PPA ponds and Pitimbú River. Finally, it is expected that this study can assist the managers in making decisions in urban and environmental planning of the municipality

Mapping of the Eltanin impact area in the South-East Pacific

Two Polarstern expeditions were conducted in 1995 (ANT-XII/4) and 2001 (ANT-XVIII/5a) to the Bellingshausen Sea and Amundsen Sea and the suspected Eltanin meteorite impact in the SE-Pacific. A survey of the sediment distribution and its acoustic structure along the cruise track was performed. The seafloor topography was sampled using the multibeam sonar system Hydrosweep DS2 which operates on a frequency of 15.5 kHz. The resulting AWI Bathymetric Chart of the Eltanin Meteorite Impact Area is based on a Digital Terrain Model of this area. The mapping was performed using ArcGIS. The Eltanin impact area which covers the 4.100 m high Freden Seamount is visualized by one overview sheet of the scale 1:200,000 and four 1:100,000 subsheets.

AWI Bathymetric Chart of the Gakkel Ridge (AWI BCGR) (Scale 1:150,000)

The 1 : 1,500,000 AWI Bathymetric Chart of the Gakkel Ridge (AWI BCGR) has been developed from multibeam data measured during the Arctic Mid-Ocean Ridge Expedition in 2001 (AMORE 2001, ARK-XVII/2). This expedition was conducted to investigate the Gakkel Ridge in the Arctic Ocean and was carried out by the icebreaking research vessels RV Polarstern and USCGC Healy. Polarstern is equipped with the multibeam sonar system Hydrosweep DS-2, whereas Healy carries Seabeam 2112. During the expedition an area of 8890 km length and 18 - 46 km width, situated between 82°N/8°W and 87°N/75°E, was surveyed simultaneously by both vessels. Water depths ranged from 566 to 5673 meters. Dense sea ice cover derogated the sonar measurements and decreased data quality. Data errors were corrected in an extensive post-processing. The data of two different sonar systems had to be consolidated in order to derive a high resolution bathymetry of the Gakkel Ridge. Final result was a digital terrain model (DTM) with a grid spacing of 100 meters, which was utilized for generating the map series AWI Bathymetric Chart of the Gakkel Ridge, consisting of ten map sheets.

Geomorphology and sedimentology of a modern isolated carbonate platform: The glorieuses archipelago, SW Indian Ocean

The study of modern carbonate systems is commonly helps in improving facies interpretation in fossil reefs and in providing analogues of sediment distribution depending on the specific platform configuration (i.e. rimmed shelves and isolated carbonate platforms). This paper deals with a geomorphological and sedimentological study of the Glorieuses Archipelago, an isolated carbonate platform located between the northern tip of Madagascar and Mayotte. The dataset consists of Digital Terrain Model, satellite imagery, and box-sediment samples. Analyses of grain-size and composition of carbonate grains are used to characterize the distribution and heterogeneity of sediment accumulated on the isolated platform. Main results show that the Glorieuses Archipelago is organized in distinctive morphological units, including a reef flat developed along the windward side, an apron, and a semi-enclosed (< 12 m water depth) to open lagoon (> 12 m and up to 15 m water depth). The lack of carbonate mud in sediments deposited on the archipelago can be explained by the direct connection between the lagoon and the open ocean. The main carbonate grains include Halimeda segments, coral fragments, large benthic foraminifers, red algae, and molluscs. According to the shape and the position of intertidal sandwaves, the current arrangement of moderately sorted fine to medium sands appears to be strongly influenced by tidal currents. The in-situ sediment production, accumulation and transport on the platform finally contribute to carbonate sand export to distinct deep marine areas depending on wind regimes and currents.

Utilização de restituição aerofotogramétrica para análise da influência das alterações da rede de drenagem e desenvolvimento de voçorocas devido a urbanização no Município de Jacareí (SP)

Dissertação (mestrado)—Universidade de Brasília, Instituto de Ciências Humanas, Departamento de Geografia, 2015.

Spatial Modeling of a Soil Fertility Index using Visible-Near-Infrared Spectra and Terrain Attributes

Our objective was to develop a methodology to predict soil fertility using visible near-infrared (vis-NIR) diffuse reflectance spectra and terrain attributes derived from a digital elevation model (DEM). Specifically, our aims were to: (i) assemble a minimum data set to develop a soil fertility index for sugarcane (Sarcharum officinarum L.) (SFI-SC) for biofuel production in tropical soils; (ii) construct a model to predict the SFI-SC using soil vis-NIR spectra and terrain attributes; and (iii) produce a soil fertility map for our study area and assess it by comparing it with a green vegetation index (GVI). The study area was 185 ha located in sao Paulo State, Brazil. In total, 184 soil samples were collected and analyzed for a range of soil chemical and physical properties. Their vis-NIR spectra were collected from 400 to 2500 nm. The Shuttle Radar Topographic Mission 3-arcsec (90-m resolution) DEM of the area was used to derive 17 terrain attributes. A minimum data set of soil properties was selected to develop the SFI-SC. The SFI-SC consisted of three classes: Class 1, the highly fertile soils; Class 2, the fertile soils; and Class 3, the least fertile soils. It was derived heuristically with conditionals and using expert knowledge. The index was modeled with the spectra and terrain data using cross-validated decision trees. The cross-validation of the model correctly predicted Class 1 in 75% of cases, Class 2 in 61%, and Class 3 in 65%. A fertility map was derived for the study area and compared with a map of the GVI. Our approach offers a methodology that incorporates expert knowledge to derive the SFI-SC and uses a versatile spectro-spatial methodology that may be implemented for rapid and accurate determination of soil fertility and better exploration of areas suitable for production.

Spatial prediction of monthly wind speeds in complex terrain with adaptive general regression neural networks

This paper presents the general regression neural networks (GRNN) as a nonlinear regression method for the interpolation of monthly wind speeds in complex Alpine orography. GRNN is trained using data coming from Swiss meteorological networks to learn the statistical relationship between topographic features and wind speed. The terrain convexity, slope and exposure are considered by extracting features from the digital elevation model at different spatial scales using specialised convolution filters. A database of gridded monthly wind speeds is then constructed by applying GRNN in prediction mode during the period 1968-2008. This study demonstrates that using topographic features as inputs in GRNN significantly reduces cross-validation errors with respect to low-dimensional models integrating only geographical coordinates and terrain height for the interpolation of wind speed. The spatial predictability of wind speed is found to be lower in summer than in winter due to more complex and weaker wind-topography relationships. The relevance of these relationships is studied using an adaptive version of the GRNN algorithm which allows to select the useful terrain features by eliminating the noisy ones. This research provides a framework for extending the low-dimensional interpolation models to high-dimensional spaces by integrating additional features accounting for the topographic conditions at multiple spatial scales. Copyright (c) 2012 Royal Meteorological Society.

Linear and nonlinear terrain deformation maps from a reduced set of interferometric SAR images

In this paper, an advanced technique for the generation of deformation maps using synthetic aperture radar (SAR) data is presented. The algorithm estimates the linear and nonlinear components of the displacement, the error of the digital elevation model (DEM) used to cancel the topographic terms, and the atmospheric artifacts from a reduced set of low spatial resolution interferograms. The pixel candidates are selected from those presenting a good coherence level in the whole set of interferograms and the resulting nonuniform mesh tessellated with the Delauney triangulation to establish connections among them. The linear component of movement and DEM error are estimated adjusting a linear model to the data only on the connections. Later on, this information, once unwrapped to retrieve the absolute values, is used to calculate the nonlinear component of movement and atmospheric artifacts with alternate filtering techniques in both the temporal and spatial domains. The method presents high flexibility with respect to the required number of images and the baselines length. However, better results are obtained with large datasets of short baseline interferograms. The technique has been tested with European Remote Sensing SAR data from an area of Catalonia (Spain) and validated with on-field precise leveling measurements.