A survey of debris trails from short-period comets

We observed 34 comets using the 24 micron camera on the Spitzer Space Telescope. Each image contains the nucleus and covers at least 10^6 km of each comet's orbit. Debris trails due to mm-sized or larger particles were found along the orbits of 27 comets; 4 comets had small-particle dust tails and a viewing geometry that made debris trails impossible to distinguish; and only 3 had no debris trail despite favorable observing conditions. There are now 30 Jupiter-family comets with known debris trails, of which 22 are reported in this paper for the first time. The detection rate is >80%, indicating that debris trails are a generic feature of short-period comets. By comparison to orbital calculations for particles of a range of sizes ejected over 2 yr prior to observation, we find that particles comprising 4 debris trails are typically mm-sized while the remainder of the debris trails require particles larger than this. The lower-limit masses of the debris trails are typically 10^11 g, and the median mass loss rate is 2 kg/s. The mass-loss rate in trail particles is comparable to that inferred from OH production rates and larger than that inferred from visible-light scattering in comae.Comment: accepted by Icarus; figures compressed for astro-p

Volcanic Hot-Spot Detection Using SENTINEL-2: A Comparison with MODIS−MIROVA Thermal Data Series

In the satellite thermal remote sensing, the new generation of sensors with high-spatial resolution SWIR data open the door to an improved constraining of thermal phenomena related to volcanic processes, with strong implications for monitoring applications. In this paper, we describe a new hot-spot detection algorithm developed for SENTINEL-2/MSI data that combines spectral indices on the SWIR bands 8a-11-12 (with a 20-meter resolution) with a spatial and statistical analysis on clusters of alerted pixels. The algorithm is able to detect hot-spot-contaminated pixels (S2Pix) in a wide range of environments and for several types of volcanic activities, showing high accuracy performances of about 1% and 94% in averaged omission and commission rates, respectively, underlining a strong reliability on a global scale. The S2-derived thermal trends, retrieved at eight key-case volcanoes, are then compared with the Volcanic Radiative Power (VRP) derived from MODIS (Moderate Resolution Imaging Spectroradiometer) and processed by the MIROVA (Middle InfraRed Observation of Volcanic Activity) system during an almost four-year-long period, January 2016 to October 2019. The presented data indicate an overall excellent correlation between the two thermal signals, enhancing the higher sensitivity of SENTINEL-2 to detect subtle, low-temperature thermal signals. Moreover, for each case we explore the specific relationship between S2Pix and VRP showing how different volcanic processes (i.e., lava flows, domes, lakes and open-vent activity) produce a distinct pattern in terms of size and intensity of the thermal anomaly. These promising results indicate how the algorithm here presented could be applicable for volcanic monitoring purposes and integrated into operational systems. Moreover, the combination of high-resolution (S2/MSI) and moderate-resolution (MODIS) thermal timeseries constitutes a breakthrough for future multi-sensor hot-spot detection systems, with increased monitoring capabilities that are useful for communities which interact with active volcanoes

Long time series (1984–2020) of albedo variations on the Greenland ice sheet from harmonized Landsat and Sentinel 2 imagery

Albedo is a key factor in modulating the absorption of solar radiation on ice surfaces. Satellite measurements have shown a general reduction in albedo across the Greenland ice sheet over the past few decades, particularly along the western margin of the ice sheet, a region known as the Dark Zone (albedo < 0.45). Here we chose a combination of Landsat 4–8 and Sentinel 2 imagery to enable us to derive the longest record of albedo variations in the Dark Zone, running from 1984 to 2020. We developed a simple, pragmatic and efficient sensor transformation to provide a long time series of consistent, harmonized satellite imagery. Narrow to broadband conversion algorithms were developed from regression models of harmonized satellite data and in situ albedo from the Program for Monitoring of the Greenland Ice Sheet (PROMICE) automatic weather stations. The albedo derived from the harmonized Landsat and Sentinel 2 data shows that the maximum extent of the Dark Zone expanded rapidly between 2005 and 2007, increasing to ~280% of the average annual maximum extent of 2900 km2 to ~8000 km2 since. The Dark Zone is continuing to darken slowly, with the average annual minimum albedo decreasing at a rate of (p = 0.16, 2001–2020)

Ocean surface currents reconstruction from microwave radiometers measurements

Premi Extraordinari de Doctorat, promoció 2014-2015. Àmbit d'Enginyeria de les TICOcean currents are a key component to understanding many oceanic and climatic phenomena and knowledge of them is crucial for both navigation and operational applications. Therefore, a key problem in oceanography is the estimation of the synoptic velocity field. Currently, global ocean surface velocities are routinely estimated from Sea Surface Height (SSH) measurements provided by altimeters. However, the separation between passes, as well as and the limited number of available altimeters leads to errors in the accurate location of oceanic currents when these measurements are used exclusively. Contrarily, satellite images of Sea Surface Temperature (SST) provide a good qualitative view of the location of ocean patterns, which has encouraged the investigation of alternative methodologies to reconstruct the velocity field based on these observations. This Ph.D. thesis has assessed the capability of SST microwave radiometers observations to retrieve ocean surface currents. The reconstruction of the ocean surface currents from SST observations can be expressed in terms of a transfer function notation, that allows to convert SST maps into SSH, and thus into currents. Because under geostrophic balance, the slope of SSH is proportional to ocean surface currents. This transfer function can be theoretically derived using the Surface Quasi-Geostrophic equations (SQG). Two different approaches were analyzed at a global scale: on one side, the analysis of the validity of the SQG approach has been performed, and on the other, an approach based on the synergetic properties between simultaneous SST and SSH observations has been analyzed. Both approaches have been compared with ocean surface currents retrieved from merged altimetric observations. The study has been focused on the period from October 2002 to May 2005, since during that period there were available four different altimeters, and the quality of the merged altimetric observations was enhanced. The analysis of the validity of SQG at a global scale revealed that this dynamical model is valid near the major extratropic current system such us the Gulf Stream, the Antartic Circumpolar Current, Kuroshio currents. Besides, the potential of MW SST observations to reconstruct ocean surface currents was analyzed using a synergetic approach: the combination of the SST phase with the SSH spectra. Actually, we explored under which environmental conditions the phase of the MW SST is close to the SSH phase. Results showed that the phase of the MW SST can be used to retrieve ocean currents during winter, near the major extratropical current systems, which are characterized by an intense mesoscale activity and the presence of strong thermal gradients, and deep ML. Furthermore, the reconstruction of the velocity fields from an ideal transfer function built up from simultaneous SST and SSH observations revealed that the SQG approach can be enhanced. The spectral properties of this ideal transfer function derived from simultaneous SST and SSH observations were characterized at a global scale. The analysis of spectral properties of the transfer function between SST and SSH observations revealed that despite daily spectral can be flatter or steeper than the k^{-1} predicted by SQG theory, in mean eSQG is a good statistically approach to retrieve ocean currents, when no simultaneous observations of SSH and SST are available.Las corrientes oceánicas son clave en muchos procesos oceánicos y climáticos, y su conocimiento es crucial para aplicaciones operacionales y de navegación. Por lo tanto, un aspecto importante en oceanografía es la estimación de campos sinópticos del campo de velocidades superficiales del mar. Actualmente, las velocidades superficiales el mar se estiman rutinariamente a partir de medidas del nivel del mar proporcionadas por altimetros, denotadas a partir de ahora con sus siglas en inglés SSH. Sin embargo, la llocalización de las corrientes puede no ser la correcta si solo se utilizan este tipo de medidas para su estimación, debido a la separación entre trazas del satélite. Por contra, las imágenes de temperatura superficial del mar, SST, proporiconan una visión cualitativa de la localización de las estructruas oceánicas. Este hecho ha motivado la investigación de metodologías alternativas para reconstruir los campos de velocidades superficiales del mar basados en estas observaciones. Esta tesis doctoral ha investigado la capacidad de las observaciones de SST proporcionadas por radiometros de microondas para recuperar las corrientes oceánicas superficiales. La reconstrucción de estas velocidades a partir de observaciones de SST se puede expresar en términos de una función de transferencia que relacione las observaciones de SST con las observaciones de SSH. Con lo que la estimación del campo de velocidades es directa, puesto que bajo la condición de equilibrio geostrófico la pendiente de la SSH es proporcional a las corrientes oceánicas. Esta función de transferencia se puede derivar teóricamente mediante las equaciones superficiales cuasi-geotróficas, denotadas con sus siglas en inglés SQG a partir de ahora. Una pregunta clave, es si las ecuaciones de este modelo dinámico son válidas. En esta tesis, se han llevado a cabo dos aproximaciones diferentes para la reconstrucción del campo de velocidades superficiales del mar: por un lado, el análisis de la validez de las ecuaciones SQG, y por otro, una aproximación basada en las propiedades espectrales de medidas simultáneas de SST y SSH. El estudio se ha centrado en el período comprendido entre Octubre del 2002 y Mayo del 2005, puesto que durante este período había disponibles hasta cuatro altmímetros, y consecuentemente la calidad de las observaciones es mayor. El análisis de la validez de SQG a escala global reveló que este modelo dinámico es válido en las regiones cerca de los sistemas de corrientes extratropicales, como la corriente del Golfo, la Corriente Circumpolar Antártica (ACC), o la Kuroshio. Además, el potencial de las observaciones de SST en el rango de las microondas para la recuperación del campo de velocidades superficiales del mar, ha sido analizado utilizando un método que combina la fase de la SST con el espectro de SSH. En realidad, se ha investigado bajo que condiciones la SST y SSH están en fase. Los resultados mostraron que la fase de la SST de microondas puede utilizarse para para la reconstrucción en invierno, cerca de los sistemas de corrientes extratropicales, caracterizados por una intensa actividad de mesoscala y la presencia de fuertes gradientes termales, así como de capas de mezcla profundas. Asimismo, la reconstrucción del campo de velocidades a partir de una función de transferencia ideal, construida a partir de imágenes simultaneas de SST y SSH, reveló que la aproximación SQG puede ser mejorada. Las propiedades espectrales de esta función de tranferencia ideal han sido estudiadas., así como su variabilidad temporal. Este análisis desveló que para escalas pequeñas y zonas enegéticas, la aproximación SQG es una buena aproximación, al menos, desde un punto de vista estádistico.Award-winningPostprint (published version

Earth Resources: A continuing bibliography with indexes (Issue 37)

This bibliography lists 512 reports, articles, and other documents introduced into the NASA scientific and technical information system between January 1 and March 31, 1983. Emphasis is placed on the use of remote sensing and geophysical instrumentation in spacecraft and aircraft to survey and inventory natural resources and urban areas. Subject matter is grouped according to agriculture and forestry, environmental changes and cultural resources, geodesy and cartography, geology and mineral resources, hydrology and water management, data processing and distribution systems, instrumentation and sensors, and economic analysis

The Burst and Transient Source Experiment Earth Occultation Technique

An Earth orbiting detector sensitive to gamma ray photons will see step-like occultation features in its counting rate when a gamma ray point source crosses the Earth's limb. This is due to the change in atmospheric attenuation of the gamma rays along the line of sight. In an uncollimated detector, these occultation features can be used to locate and monitor astrophysical sources provided their signals can be individually separated from the detector background. We show that the Earth occultation technique applied to the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory (CGRO) is a viable and flexible all-sky monitor in the low energy gamma ray and hard X-ray energy range (20 keV - 1 MeV). The method is an alternative to more sophisticated photon imaging devices for astronomy, and can serve well as a cost-effective science capability for monitoring the high energy sky. Here we describe the Earth occultation technique for locating new sources and for measuring source intensity and spectra without the use of complex background models. Examples of transform imaging, step searches, spectra, and light curves are presented. Systematic uncertainties due to source confusion, detector response, and contamination from rapid background fluctuations are discussed and analyzed for their effect on intensity measurements. A sky location-dependent average systematic error is derived as a function of galactic coordinates. The sensitivity of the technique is derived as a function of incident photon energy and also as a function of angle between the source and the normal to the detector entrance window. Occultations of the Crab Nebula by the Moon are used to calibrate Earth occultation flux measurements independent of possible atmospheric scattering effects.Comment: 39 pages, 24 figures. Accepted for publication in the Astrophysical Journal Supplement

Spitzer Observations of Comet 67P/Churyumov-Gerasimenko at 5.5-4.3 AU From the Sun

We report Spitzer Space Telescope observations of comet 67P/Churyumov-Gerasimenko at 5.5 and 4.3 AU from the Sun, post-aphelion. Comet 67P is the primary target of the European Space Agency's Rosetta mission. The Rosetta spacecraft will rendezvous with the nucleus at heliocentric distances similar to our observations. Rotationally resolved observations at 8 and 24 microns (at a heliocentric distance, rh, of 4.8 AU) that sample the size and color-temperature of the nucleus are combined with aphelion R-band light curves observed at the Very Large Telescope (VLT) and yield a mean effective radius of 2.04 +/- 0.11 km, and an R-band geometric albedo of 0.054 +/- 0.006. The amplitudes of the R-band and mid-infrared light curves agree, which suggests that the variability is dominated by the shape of the nucleus. We also detect the dust trail of the comet at 4.8 and 5.5 AU, constrain the grain sizes to be less than or similar to 6 mm, and estimate the impact hazard to Rosetta. We find no evidence for recently ejected dust in our images. If the activity of 67P is consistent from orbit to orbit, then we may expect the Rosetta spacecraft will return images of an inactive or weakly active nucleus as it rendezvous with the comet at rh = 4 AU in 2014.Comment: 19 pages, 2 tables, 10 figures. Accepted for publication in the Astronomical Journa

Observational and Dynamical Characterization of Main-Belt Comet P/2010 R2 (La Sagra)

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out.Comment: 23 pages, 13 figures, accepted for publication in A

Exploring the Venus global super-rotation using a comprehensive General Circulation Model

The atmospheric circulation in Venus is well known to exhibit strong super-rotation. However, the atmospheric mechanisms responsible for the formation of this super-rotation are still not fully understood. In this work, we developed a new Venus general circulation model to study the most likely mechanisms driving the atmosphere to the current observed circulation. Our model includes a new radiative transfer, convection and suitably adapted boundary layer schemes and a dynamical core that takes into account the dependence of the heat capacity at constant pressure with temperature. The new Venus model is able to simulate a super-rotation phenomenon in the cloud region quantitatively similar to the one observed. The mechanisms maintaining the strong winds in the cloud region were found in the model results to be a combination of zonal mean circulation, thermal tides and transient waves. In this process, the semi-diurnal tide excited in the upper clouds has a key contribution in transporting axial angular momentum mainly from the upper atmosphere towards the cloud region. The magnitude of the super-rotation in the cloud region is sensitive to various radiative parameters such as the amount of solar radiative energy absorbed by the surface, which controls the static stability near the surface. In this work, we also discuss the main difficulties in representing the flow below the cloud base in Venus atmospheric models. Our new radiative scheme is more suitable for 3D Venus climate models than those used in previous work due to its easy adaptability to different atmospheric conditions. This flexibility of the model was crucial to explore the uncertainties in the lower atmospheric conditions and may also be used in the future to explore, for example, dynamical-radiative-microphysical feedbacks.Comment: Accepted for publication in Planet. Space Sc