Millimeter and sub-millimeter atmospheric performance at Dome C combining radiosoundings and ATM synthetic spectra

The reliability of astronomical observations at millimeter and sub-millimeter wavelengths closely depends on a low vertical content of water vapor as well as on high atmospheric emission stability. Although Concordia station at Dome C (Antarctica) enjoys good observing conditions in this atmospheric spectral windows, as shown by preliminary site-testing campaigns at different bands and in, not always, time overlapped periods, a dedicated instrument able to continuously determine atmospheric performance for a wide spectral range is not yet planned. In the absence of such measurements, in this paper we suggest a semi-empirical approach to perform an analysis of atmospheric transmission and emission at Dome C to compare the performance for 7 photometric bands ranging from 100 GHz to 2 THz. Radiosoundings data provided by the Routine Meteorological Observations (RMO) Research Project at Concordia station are corrected by temperature and humidity errors and dry biases and then employed to feed ATM (Atmospheric Transmission at Microwaves) code to generate synthetic spectra in the wide spectral range from 100 GHz to 2 THz. To quantify the atmospheric contribution in millimeter and sub-millimeter observations we are considering several photometric bands in which atmospheric quantities are integrated. The observational capabilities of this site at all the selected spectral bands are analyzed considering monthly averaged transmissions joined to the corresponding fluctuations. Transmission and pwv statistics at Dome C derived by our semi-empirical approach are consistent with previous works. It is evident the decreasing of the performance at high frequencies. We propose to introduce a new parameter to compare the quality of a site at different spectral bands, in terms of high transmission and emission stability, the Site Photometric Quality Factor.Comment: accepted to MNRAS with minor revision

Terra and Aqua MODIS TEB Inter-Comparison Using Himawari-8/AHI as Reference

Intercomparison between the two MODIS instruments is very useful for both the instrument calibration and its uncertainty assessment. Terra and Aqua MODIS have almost identical relative spectral response, spatial resolution, and dynamic range for each band, so the site-dependent effect from spectral mismatch for their comparison is negligible. Major challenges in cross-sensor comparison of instruments on different satellites include differences in observation time and view angle over selected pseudoinvariant sites. The simultaneous nadir overpasses (SNO) between the two satellites are mostly applied for comparison and the scene under SNO varies. However, there is a dearth of SNO between the Terra and Aqua. This work focuses on an intercomparison method for MODIS thermal emissive bands using Himawari-8 Advanced Himawari Imager (AHI) as a reference. Eleven thermal emissive bands on MODIS are at least to some degree spectrally matched to the AHI bands. The sites selected for the comparison are an ocean area around the Himawari-8 suborbital point and the Strzelecki Desert located south of the Himawari-8 suborbital point. The time difference between the measurements from AHI and MODIS is <5 min. The comparison is performed using 2017 collection 6.1 L1B data for MODIS. The MODISAHI difference is corrected to remove the view angle dependence. The TerraAqua MODIS difference for the selected TEB is up to 0.6 K with the exception of band 30. Band 30 has the largest difference, which is site dependent, most likely due to a crosstalk effect. Over the ocean, the band 30 difference between the two MODIS instruments is around 1.75 K, while over the desert; the difference is around 0.68 K. The MODIS precision is also compared from the Gaussian regression of the double difference. Terra bands 27 to 30 have significant extra noise due to crosstalk effects on these bands. These TerraAqua comparison results are used for MODIS calibration assessments and are beneficial for future calibration algorithm improvement. The impact of daytime measurements and the scene dependence are also discussed

Early Analysis of Landsat-8 Thermal Infrared Sensor Imagery of Volcanic Activity

The Landsat-8 satellite of the Landsat Data Continuity Mission was launched by the National Aeronautics and Space Administration (NASA) in April 2013. Just weeks after it entered active service, its sensors observed activity at Paluweh Volcano, Indonesia. Given that the image acquired was in the daytime, its shortwave infrared observations were contaminated with reflected solar radiation; however, those of the satellite’s Thermal Infrared Sensor (TIRS) show thermal emission from the volcano’s summit and flanks. These emissions detected in sensor’s band 10 (10.60–11.19 µm) have here been quantified in terms of radiant power, to confirm reports of the actual volcanic processes operating at the time of image acquisition, and to form an initial assessment of the TIRS in its volcanic observation capabilities. Data from band 11 have been neglected as its data have been shown to be unreliable at the time of writing. At the instant of image acquisition, the thermal emission of the volcano was found to be 345 MW. This value is shown to be on the same order of magnitude as similarly timed NASA Earth Observing System (EOS) Moderate Resolution Imaging Spectroradiometer thermal observations. Given its unique characteristics, the TIRS shows much potential for providing useful, detailed and accurate volcanic observations in the future

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

Overview of Intercalibration of Satellite Instruments

Intercalibration of satellite instruments is critical for detection and quantification of changes in the Earth’s environment, weather forecasting, understanding climate processes, and monitoring climate and land cover change. These applications use data from many satellites; for the data to be interoperable, the instruments must be cross-calibrated. To meet the stringent needs of such applications, instruments must provide reliable, accurate, and consistent measurements over time. Robust techniques are required to ensure that observations from different instruments can be normalized to a common scale that the community agrees on. The long-term reliability of this process needs to be sustained in accordance with established reference standards and best practices. Furthermore, establishing physical meaning to the information through robust Système International d’unités traceable calibration and validation (Cal/Val) is essential to fully understand the parameters under observation. The processes of calibration, correction, stabilitymonitoring, and quality assurance need to be underpinned and evidenced by comparison with “peer instruments” and, ideally, highly calibrated in-orbit reference instruments. Intercalibration between instruments is a central pillar of the Cal/Val strategies of many national and international satellite remote sensing organizations. Intercalibration techniques as outlined in this paper not only provide a practical means of identifying and correcting relative biases in radiometric calibration between instruments but also enable potential data gaps between measurement records in a critical time series to be bridged. Use of a robust set of internationally agreed upon and coordinated intercalibration techniques will lead to significant improvement in the consistency between satellite instruments and facilitate accurate monitoring of the Earth’s climate at uncertainty levels needed to detect and attribute the mechanisms of change. This paper summarizes the state-of-the-art of postlaunch radiometric calibration of remote sensing satellite instruments through intercalibration

Observing CMB polarisation through ice

Ice crystal clouds in the upper troposphere can generate polarisation signals at the uK level. This signal can seriously affect very sensitive ground based searches for E- and B-mode of Cosmic Microwave Background polarisation. In this paper we estimate this effect within the ClOVER experiment observing bands (97, 150 and 220 GHz) for the selected observing site (Llano de Chajnantor, Atacama desert, Chile). The results show that the polarisation signal from the clouds can be of the order of or even bigger than the CMB expected polarisation. Climatological data suggest that this signal is fairly constant over the whole year in Antarctica. On the other hand the stronger seasonal variability in Atacama allows for a 50% of clean observations during the dry season.Comment: 7 Pages, 4 figure

OGLE-2014-BLG-0289: Precise Characterization of a Quintuple-peak Gravitational Microlensing Event

We present the analysis of the binary-microlensing event OGLE-2014-BLG-0289. The event light curve exhibits five very unusual peaks, four of which were produced by caustic crossings and the other by a cusp approach. It is found that the quintuple-peak features of the light curve provide tight constraints on the source trajectory, enabling us to precisely and accurately measure the microlensing parallax πE. Furthermore, the three resolved caustics allow us to measure the angular Einstein radius θE. From the combination of πE and θE, the physical lens parameters are uniquely determined. It is found that the lens is a binary composed of two M dwarfs with masses M1 = 0.52 ± 0.04 M⊙ and M2 = 0.42 ± 0.03 M⊙ separated in projection by a⊥ = 6.4 ± 0.5 au. The lens is located in the disk with a distance of DL = 3.3 ± 0.3 kpc. The reason for the absence of a lensing signal in the Spitzer data is that the time of observation corresponds to the flat region of the light curve

Overview of Intercalibration of Satellite Instruments

Intercalibration of satellite instruments is critical for detection and quantification of changes in the Earth’s environment, weather forecasting, understanding climate processes, and monitoring climate and land cover change. These applications use data from many satellites; for the data to be interoperable, the instruments must be cross-calibrated. To meet the stringent needs of such applications, instruments must provide reliable, accurate, and consistent measurements over time. Robust techniques are required to ensure that observations from different instruments can be normalized to a common scale that the community agrees on. The long-term reliability of this process needs to be sustained in accordance with established reference standards and best practices. Furthermore, establishing physical meaning to the information through robust Système International d’unités traceable calibration and validation (Cal/Val) is essential to fully understand the parameters under observation. The processes of calibration, correction, stabilitymonitoring, and quality assurance need to be underpinned and evidenced by comparison with “peer instruments” and, ideally, highly calibrated in-orbit reference instruments. Intercalibration between instruments is a central pillar of the Cal/Val strategies of many national and international satellite remote sensing organizations. Intercalibration techniques as outlined in this paper not only provide a practical means of identifying and correcting relative biases in radiometric calibration between instruments but also enable potential data gaps between measurement records in a critical time series to be bridged. Use of a robust set of internationally agreed upon and coordinated intercalibration techniques will lead to significant improvement in the consistency between satellite instruments and facilitate accurate monitoring of the Earth’s climate at uncertainty levels needed to detect and attribute the mechanisms of change. This paper summarizes the state-of-the-art of postlaunch radiometric calibration of remote sensing satellite instruments through intercalibration

GNSS-R as a source of opportunity for remote sensing of the cryosphere

This work evaluates the potential use of signals from the Global Navigation Satellite Systems (GNSS) that scatter off the Earth surface for the retrieval of geophysical information from the cryosphere. For this purpose, the present study is based on data collected with a dedicated reflectometry GNSS receiver during two field campaigns, which were focused on two types of characteristic surfaces of the cryosphere: thin sea ice covers and thick dry snow accumulations. During the first experiment, the complete process of formation, evolution and melting of sea ice was monitorized for more than seven months in a bay located in Greenland. This type of ice is typically characterized by its thickness, concentration and roughness. Different observables from GNSS reflections are analyzed to try to infer these properties. The ice thickness is linked to the free-board level, defined as the height of the sea ice surface. Accurate phase altimetry is achieved, showing good agreement with an Arctic tide model. In addition, the long term results of ellipsoidal height retrievals are consistent with the evolution of the ice surface temperature product given by MODIS, which is a key parameter in the rate of growth of sea ice. On the other hand, the presence of salinity in the sea ice modifies its dielectric properties, resulting in different amplitude and phase for the co- and cross-polar components of the complex Fresnel coefficients. The polarimetric measurements obtained show good agreement with visual inspections of ice concentration from an Arctic weather station. Finally, the shape of the reflected signals and its phase dispersion are tested as potential signatures of surface roughness. For comparison, ice charts of the experimental area are employed. In particular, maximums in roughness given by the GNSS observables coincide with fast ice events. Fast ice is defined as ice anchored to the coast, where the tidal movements contribute to the development of strange patterns, cracks, and fissures on its surface, thus consistent with the GNSS-R roughness retrievals. The second experiment took place on Antarctica, monitoring a pristine snow area which is well-known for the calibration of remote sensing instruments. Due to the relative stability of the snow layers, the data acquisition was limited to ten continuous days. Interferometric beats were found after a first analysis of the amplitude from the collected signals, which were consistent with a multipath model where the reflector lies below the surface level. Motivated by these results, a forward model is developed that reconstructs the complex received signal as a sum of a finite number of reflections, coming from different snow layers (a snow density profile obtained from in-situ measurements). The interferometric information is then retrieved from the spectral analysis applied to time series from both real and modeled signals (lag-holograms). We find that the frequency bands predicted by the model are in general consistent with the data and the lag-holograms show repeatability for different days. Then, we attempt a proper inversion of the collected data to determine the dominant layers of the dry snow profile that contribute to L-band reflections, which are related to significant gradients of snow density/permittivity.Aquest treball avalua el possible ús dels senyals dels sistemes mundials de navegació per satèl lit (GNSS) que es reflecteixen a la superfície terrestre, per a l’extracció de la informació geofísica de la criosfera. Amb aquest propòsit, el present estudi es basa en dades recollides amb un reflectòmetre GNSS durant dues campanyes experimentals, centrades en dos tipus de superfícies característiques de la criosfera: cobertes de gel marí i gruixudes acumulacions de neu seca. En el primer experiment, el procés complet de formació, evolució i fusió del gel marí va ser monitoritzat durant més de set mesos a una badia situada a Groenlàndia. Aquest tipus de gel es caracteritza típicament amb el seu gruix, concentració i rugositat. Diferents observables de les reflexions GNSS són analitzats per tractar de fer una estimació d’aquestes propietats. El gruix de gel està relacionat amb el nivell de francbord, que a la seva vegada està relacionat amb l’alçada de la superfície de gel marí. S’ha aconseguit altimetria de fase precisa, que mostra correlació amb un model de marea de l’Àrtic. A més, els resultats a llarg termini de l’alçada elipsoidal segueixen l’evolució de les mesures de temperatura de superfície de gel donades per MODIS. La temperatura és un paràmetre clau en el ritme de creixement del gel marí. Per altra banda, la presència de sal a aquest tipus de gel modifica les seves propietats dielèctriques, el que implica variacions d’amplitud i fase per als coeficients de Fresnel complexos amb polaritzacions oposades. Les mesures polarimètriques obtingudes mostren concordança amb els valors de concentració de gel obtinguts des d’una estació meteorològica propera. Finalment, la forma de la senyal reflectida i la dispersió de la seva fase s’evaluen com a potencials indicadors de la rugositat de superfície. Per a la seva comparació, es fan servir mapes del gel de la zona experimental. En concret, els valors màxims a la rugositat estimada a partir pels observables GNSS coincideixen amb el gel fixe, que es refereix a gel ancorat a la costa, on els moviments de les marees contribueixen al desenvolupament de patrons estranys, esquerdes i fissures en la seva superfície. El segon experiment es va dur a terme a l’Antàrtida, monitoritzant una àrea de neu pristina que és ben coneguda per al calibratge d’instruments de teledetecció. A causa de la relativa estabilitat de les capes de neu, l’adquisició de dades es va limitar a deu dies consecutius. Es van trobar pulsacions interferomètriques a partir d’un primer anàlisi de l’amplitud de les senyals recollides, les quals eren compatibles amb un model de propagació multicamí a on el reflector es troba per sota del nivell de superfície. Com a conseqüència d’aquests resultats, s’ha desenvolupat un model que reconstrueix el senyal complexe rebut com la suma d’un nombre finit de reflexions, procedents de diferents capes de neu (determinat per mesures locals). La informació interferomètrica es recupera després de l’anàlisi espectral aplicat a les sèries temporals tant de les senyals reals, com de les modelades (lag-hologrames). Trobem que les bandes de freqüències predites pel model són en general consistents amb les dades i que els lag-hologrames mostren repetibilitat per dies diferents. Posteriorment, es realitza un anàlisi de les dades recollides per determinar les capes dominants del perfil de neu seca que contribueixen a les reflexions en banda L, i que a la seva vegada, estan relacionades amb gradents significatius de densitat/permitivitat.Este trabajo evalúa el posible uso de las señales de los sistemas globales de navegación por satélite (GNSS) que se reflejan en la superficie terrestre para la extracción de información geofísica de la criosfera. Con este propósito, el presente estudio se basa en datos recogidos con un reflectómetro GNSS durante dos campañas experimentales, centradas en dos tipos de superficies características de la criosfera: capas de hielo marino y gruesas acumulaciones de nieve seca. Durante el primer experimento, el proceso completo de formación, evolución y fusión del hielo marino fue monitorizado durante más de siete meses en una bahía ubicada en Groenlandia. Este tipo de hielo se caracteriza típicamente por su grosor, concentración y rugosidad. Diferentes observables de las reflexiones GNSS son analizados para tratar de estimar dichas propiedades. El espesor de hielo está relacionado con el nivel de francobordo o borda libre, que a su vez está relacionado con la altura de la superficie de hielo marino. Se ha logrado altimetría de fase precisa, mostrando correlación con un modelo de marea del Ártico. Además, los resultados a largo plazo de la altura elipsoidal siguen la evolución de las mediciones de temperatura de superficie de hielo proporcionadas por MODIS. La temperatura es un parámetro clave en el ritmo de crecimiento del hielo marino. Por otro lado, la presencia de sal en este tipo de hielo modifica sus propiedades dieléctricas, lo que implica variaciones en las amplitudes y fases de los coeficientes complejos de Fresnel con polarizaciones opuestas. Los resultados polarimétricos concuerdan con los valores de concentración de hielo obtenidos mediante inspección visual desde una estación meteorológica cercana. Por último, la forma de la señal reflejada y la dispersión de su fase son evaluadas como potenciales indicadores de la rugosidad de superficie. Para su comparación, se emplean mapas del hielo de la zona experimental. En particular, valores máximos de rugosidad estimada por los observables GNSS coinciden con hielo fijo, que se refiere al hielo anclado a la costa, donde los movimientos de las mareas contribuyen al desarrollo de patrones extraños, grietas y fisuras en su superficie. El segundo experimento se llevó a cabo en la Antártida, monitorizando una área de nieve pristina que es bien conocida para la calibración de instrumentos de teledetección. Debido a la relativa estabilidad de las capas de nieve, la adquisición de datos se limitó a diez días consecutivos. Se encontraron pulsaciones interferométricas a partir de un primer análisis de la amplitud de las señales recibidas, las cuales eran compatibles con un modelo de propagación multicamino donde el reflector se encuentra por debajo del nivel de la superficie. Como consecuencia de estos resultados, se ha desarrollado un modelo que reconstruye la señal recibida como la suma de un número finito de reflexiones, procedentes de diferentes capas de nieve (caracterizados por mediciones locales). La información interferométrica se recupera después del análisis espectral aplicado a las series temporales tanto de las señales reales, como de las modeladas (lag-hologramas). Encontramos que las bandas de frecuencias predichas por el modelo son en general consistentes con los datos y que los lag-hologramas muestran repetibilidad para días diferentes. Posteriormente, se realiza un análisis de los datos recogidos para determinar las capas dominantes del perfil de nieve seca que contribuyen a las reflexiones en banda L, y que a su vez, están relacionadas con gradientes significativos de densidad/permitivida

Astroclimatic Characterization of Vallecitos: A candidate site for the Cherenkov Telescope Array at San Pedro Martir

We conducted an 18 month long study of the weather conditions of the Vallecitos, a proposed site in Mexico to harbor the northern array of the Cherenkov Telescope Array (CTA). It is located in Sierra de San Pedro Martir (SPM) a few kilometers away from Observatorio Astron\'omico Nacional. The study is based on data collected by the ATMOSCOPE, a multi-sensor instrument measuring the weather and sky conditions, which was commissioned and built by the CTA Consortium. Additionally, we compare the weather conditions of the optical observatory at SPM to the Vallecitos regarding temperature, humidity, and wind distributions. It appears that the excellent conditions at the optical observatory benefit from the presence of microclimate established in the Vallecitos.Comment: 16 pages, 16 figures, Publication of the Astronomical Society of the Pacific, accepte