scholarly journals Penetration of interferometric radar signals in Antarctic snow

2021 ◽  
Vol 15 (9) ◽  
pp. 4399-4419
Author(s):  
Helmut Rott ◽  
Stefan Scheiblauer ◽  
Jan Wuite ◽  
Lukas Krieger ◽  
Dana Floricioiu ◽  
...  

Abstract. Synthetic aperture radar interferometry (InSAR) is an efficient technique for mapping the surface elevation and its temporal change over glaciers and ice sheets. However, due to the penetration of the SAR signal into snow and ice, the apparent elevation in uncorrected InSAR digital elevation models (DEMs) is displaced versus the actual surface. We studied relations between interferometric radar signals and physical snow properties and tested procedures for correcting the elevation bias. The work is based on satellite and in situ data over Union Glacier in the Ellsworth Mountains, West Antarctica, including interferometric data of the TanDEM-X mission, topographic data from optical satellite sensors and field measurements on snow structure, and stratigraphy undertaken in December 2016. The study area comprises ice-free surfaces, bare ice, dry snow and firn with a variety of structural features related to local differences in wind exposure and snow accumulation. Time series of laser measurements of NASA's Ice, Cloud and land Elevation Satellite (ICESat) and ICESat-2 show steady-state surface topography. For area-wide elevation reference we use the Reference Elevation Model of Antarctica (REMA). The different elevation data are vertically co-registered on a blue ice area that is not affected by radar signal penetration. Backscatter simulations with a multilayer radiative transfer model show large variations for scattering of individual snow layers, but the vertical backscatter distribution can be approximated by an exponential function representing uniform absorption and scattering properties. We obtain estimates of the elevation bias by inverting the interferometric volume correlation coefficient (coherence), applying a uniform volume model for describing the vertical loss function. Whereas the mean values of the computed elevation bias and the elevation difference between the TanDEM-X DEMs and the REMA show good agreement, a trend towards overestimation of penetration is evident for heavily wind-exposed areas with low accumulation and towards underestimation for areas with higher accumulation rates. In both cases deviations from the uniform volume structure are the main reason. In the first case the dense sequence of horizontal structures related to internal wind crust, ice layers and density stratification causes increased scattering in near-surface layers. In the second case the small grain size of the top snow layers causes a downward shift in the scattering phase centre.

2021 ◽  
Author(s):  
Helmut Rott ◽  
Stefan Scheiblauer ◽  
Jan Wuite ◽  
Lukas Krieger ◽  
Dana Floricioiu ◽  
...  

Abstract. Synthetic aperture radar interferometry (InSAR) is an efficient technique for mapping the surface elevation and its temporal change over glaciers and ice sheets. However, due to the penetration of the SAR signal into snow and ice the apparent elevation in uncorrected InSAR digital elevation models (DEMs) is displaced versus the actual surface. We studied relations between interferometric radar signals and physical snow properties and tested procedures for correcting the elevation bias. The work is based on satellite and in-situ data over Union Glacier in the Ellsworth Mountains, West Antarctica, including interferometric data of the TanDEM-X mission, topographic data from optical satellite sensors and field measurements on snow structure and stratigraphy undertaken in December 2016. The study area comprises ice-free surfaces, bare ice, dry snow and firn with a variety of structural features related to local differences in wind exposure and snow accumulation. Time series of laser measurements of NASA’s Ice, Cloud and land Elevation Satellite (ICESat) and ICESat-2 show steady state surface topography. For area-wide elevation reference we use the Reference Elevation Model of Antarctica (REMA). The different elevation data are vertically co-registered on a blue ice area and an ice-free slope, surfaces not affected by radar signal penetration. The backscatter simulations with a multi-layer radiative transfer model show large variations for scattering of individual snow layers due to different size and structure of the scattering elements. The average depth-dependent backscatter contributions can be approximated by an exponential function. We obtain estimates of the elevation bias by inverting the interferometric volume correlation coefficient (coherence) applying a uniform volume model for describing the vertical loss function. Whereas the mean values of the computed elevation bias and the elevation difference between the TDM DEMs and the REMA show good agreement, a trend towards overestimation of penetration is evident for heavily wind-exposed areas and towards underestimation for areas with higher accumulation rates. The angular gradients of the backscatter intensity show also distinct differences between these two domains. This behaviour can be attributed to the anisotropy of the snow/firn volume structure showing differences in the size and shape of the scattering elements and in stratification related to snow accumulation and wind-driven erosion and deposition.


2019 ◽  
Vol 19 (18) ◽  
pp. 11651-11668 ◽  
Author(s):  
Francisco Navas-Guzmán ◽  
Giovanni Martucci ◽  
Martine Collaud Coen ◽  
María José Granados-Muñoz ◽  
Maxime Hervo ◽  
...  

Abstract. This study focuses on the analysis of aerosol hygroscopicity using remote sensing techniques. Continuous observations of aerosol backscatter coefficient (βaer), temperature (T) and water vapor mixing ratio (r) have been performed by means of a Raman lidar system at the aerological station of MeteoSwiss at Payerne (Switzerland) since 2008. These measurements allow us to monitor in a continuous way any change in aerosol properties as a function of the relative humidity (RH). These changes can be observed either in time at a constant altitude or in altitude at a constant time. The accuracy and precision of RH measurements from the lidar have been evaluated using the radiosonde (RS) technique as a reference. A total of 172 RS profiles were used in this intercomparison, which revealed a bias smaller than 4 % RH and a standard deviation smaller than 10 % RH between both techniques in the whole (in lower) troposphere at nighttime (at daytime), indicating the good performance of the lidar for characterizing RH. A methodology to identify situations favorable to studying aerosol hygroscopicity has been established, and the aerosol hygroscopicity has been characterized by means of the backscatter enhancement factor (fβ). Two case studies, corresponding to different types of aerosol, are used to illustrate the potential of this methodology. The first case corresponds to a mixture of rural aerosol and smoke particles (smoke mixture), which showed a higher hygroscopicity (fβ355=2.8 and fβ1064=1.8 in the RH range 73 %–97 %) than the second case, in which mineral dust was present (fβ355=1.2 and fβ1064=1.1 in the RH range 68 %–84 %). The higher sensitivity of the shortest wavelength to hygroscopic growth was qualitatively reproduced using Mie simulations. In addition, a good agreement was found between the hygroscopic analysis done in the vertical and in time for Case I, where the latter also allowed us to observe the hydration and dehydration of the smoke mixture. Finally, the impact of aerosol hygroscopicity on the Earth's radiative balance has been evaluated using the GAME (Global Atmospheric Model) radiative transfer model. The model showed an impact with an increase in absolute value of 2.4 W m−2 at the surface with respect to the dry conditions for the hygroscopic layer of Case I (smoke mixture).


2017 ◽  
Vol 34 (8) ◽  
pp. 1723-1741 ◽  
Author(s):  
Sebastian Schemm ◽  
Aleksi Nummelin ◽  
Nils Gunnar Kvamstø ◽  
Øyvind Breivik

AbstractThe Lagrangian Analysis Tool (LAGRANTO) is adopted and applied to ECMWF’s latest ocean reanalysis. The primary motivation behind this study is to introduce and document LAGRANTO Ocean (LAGRANTO.ocean) and explore its capabilities in combination with an eddy-permitting ocean reanalysis. The tool allows for flexibly defining starting points, within circles, cylinders, or any user-defined region or volume. LAGRANTO.ocean also offers a sophisticated way to refine a set of computed trajectories according to a wide range of mathematical operations that can be combined into a single refinement criterion. Tools for calculating—for example, along-trajectory cross sections or trajectory densities—are further provided. After introducing the tool, three case studies are presented, which were chosen to reflect a selection of phenomena on different spatial and temporal scales. The case studies also serve as hands-on examples. For the first case study, at the mesoscale, ocean trajectories are computed during the formation of a Gulf Stream cold-core ring to study vertical motion in the developing eddy. In the second example, source waters are traced to the East Greenland Spill Jet. This example highlights the usefulness of a Lagrangian method for identifying sources or sinks of buoyancy. The third example, on annual time scales, focuses on the temporal evolution of extreme potential temperature anomalies in the South Pacific and the memory of the involved water parcels. Near-surface trajectories reveal that it takes approximately 5 months after the highest temperature anomaly before the involved water parcels cool to their climatological mean values at their new positions. LAGRANTO.ocean will be made publicly available.


2010 ◽  
Vol 10 (5) ◽  
pp. 13373-13405 ◽  
Author(s):  
B. Mayer ◽  
S. W. Hoch ◽  
C. D. Whiteman

Abstract. The MYSTIC three-dimensional Monte-Carlo radiative transfer model has been extended to simulate solar and thermal irradiances with a rigorous consideration of topography. Forward as well as backward Monte Carlo simulations are possible for arbitrarily oriented surfaces and we demonstrate that the backward Monte Carlo technique is superior to the forward method for applications involving topography, by greatly reducing the computational demands. MYSTIC is used to simulate the short- and longwave radiation fields during a clear day and night in and around Arizona's Meteor Crater, a bowl-shaped, 165-m-deep basin with a diameter of 1200 m. The simulations are made over a 4 by 4 km domain using a 10-m horizontal resolution digital elevation model and meteorological input data collected during the METCRAX (Meteor Crater Experiment) field experiment in 2006. Irradiance (or radiative flux) measurements at multiple locations inside the crater are then used to evaluate the simulations. MYSTIC is shown to realistically model the complex interactions between topography and the radiative field, resolving the effects of terrain shading, terrain exposure, and longwave surface emissions. The effects of surface temperature variations and of temperature stratification within the crater atmosphere on the near-surface longwave irradiance are then evaluated with additional simulations.


2020 ◽  
Author(s):  
Yang Wang ◽  
Arnoud Apituley ◽  
Alkiviadis Bais ◽  
Steffen Beirle ◽  
Nuria Benavent ◽  
...  

Abstract. We present the inter-comparison of delta slant column densities (SCDs) and vertical profiles of nitrous acid (HONO) derived from measurements of different MAX-DOAS instruments and using different inversion algorithms during the Second Cabauw Inter-comparison campaign for Nitrogen Dioxide measuring Instruments (CINDI-2), in September 2016, at Cabauw, The Netherlands (51.97° N, 4.93° E). Systematic discrepancies of HONO delta SCDs are observed in the range of ±0.3 × 1015 molecules cm−2, which is half of the typical random discrepancy of 0.6 × 1015 molecules cm−2. For a typical high HONO delta SCD of 2 × 1015 molecules cm−2, the relative systematic and random discrepancies are about 15 % and 30 %, respectively. The inter-comparison of HONO profiles shows that both systematic and random discrepancies of HONO VCDs and near-surface volume mixing ratios (VMRs) are mostly in the range of ~ ±0.5 × 1015 molecules cm−2 and ~ ±0.1 ppb (typically ~ 20 %). Further we find that the discrepancies of the retrieved HONO profiles are dominated by discrepancies of the HONO delta SCDs. The profile retrievals only contribute to the discrepancies of the HONO profiles by ~ 5 %. However, some data sets with substantial larger discrepancies than the typical values indicate that inappropriate implementations of profile inversion algorithms and configurations of radiative transfer models in the profile retrievals can also be an important uncertainty source. In addition, estimations of measurement uncertainties of HONO dSCDs, which can significantly impact profile retrievals using the optimal estimation method, need to consider not only DOAS fit errors, but also atmospheric variability, especially for an instrument with a DOAS fit error lower than ~ 3 × 1015 molecules cm−2. The MAX-DOAS results during the CINDI-2 campaign indicate that the peak HONO levels (e.g. near-surface VMRs of ~ 0.4 ppb) often appeared in the early morning and below 0.2 km. The near-surface VMRs retrieved from the MAX-DOAS observations are compared with those measured using a co-located long-path DOAS instrument. The systematic differences are smaller than 0.15 ppb and 0.07 ppb during early morning and around noon, respectively. Since true HONO values at high altitudes are not known in the absence of real measurements, in order to evaluate the abilities of profile inversion algorithms to respond to different HONO profile shapes, we performed sensitivity studies using synthetic HONO delta SCDs simulated by a radiative transfer model with assumed HONO profiles. The tests indicate that the profile inversion algorithms based on the optimal estimation method with proper configurations can well reproduce the different HONO profile shapes. Therefore we conclude that the feature of HONO accumulated near the surface derived from MAX-DOAS measurements are expected to well represent the ambient HONO profiles.


2010 ◽  
Vol 10 (2) ◽  
pp. 331-340 ◽  
Author(s):  
I. Ialongo ◽  
V. Buchard ◽  
C. Brogniez ◽  
G. R. Casale ◽  
A. M. Siani

Abstract. The aerosol Single Scattering Albedo (SSA) and Absorbing Aerosol Optical Depth (AAOD) at 320.1 nm are derived at Rome site by the comparison between Brewer and modelled spectra. The UVSPEC radiative transfer model is used to calculate the UV irradiances for different SSA values, taking into account as input data total ozone and Aerosol Optical Depth (AOD) obtained from Brewer spectral measurements. The accuracy in determining SSA depends on the aerosol amount and on Solar Zenith Angle (SZA) value: SSA uncertainty increases when AOD and SZA decrease. The monthly mean values of SSA and AAOD during the period January 2005–June 2008 are analysed, showing a monthly and seasonal variability. It is found that the SSA and AAOD averages are 0.80±0.08 and 0.056±0.028, respectively. AAOD retrievals are also used to quantify the error in the Ozone Monitoring Instrument (OMI) surface UV products due to absorbing aerosols, not included in the current OMI UV algorithm. OMI and Brewer UV irradiances at 324.1 nm and Erythemal Dose Rates (EDRs) under clear sky conditions, are compared as a function of AAOD. Three methods are considered to investigate on the applicability of an absorbing aerosol correction on OMI UV data at Rome site. Depending on the correction methodology, the bias value decreases from 18% to 2% for spectral irradiance at 324.1 nm and from 25% to 8% for EDR.


2013 ◽  
Vol 6 (6) ◽  
pp. 1521-1532 ◽  
Author(s):  
R. Sinreich ◽  
A. Merten ◽  
L. Molina ◽  
R. Volkamer

Abstract. We present a novel parameterization method to convert multi-axis differential optical absorption spectroscopy (MAX-DOAS) differential slant column densities (dSCDs) into near-surface box-averaged volume mixing ratios. The approach is applicable inside the planetary boundary layer under conditions with significant aerosol load, and builds on the increased sensitivity of MAX-DOAS near the instrument altitude. It parameterizes radiative transfer model calculations and significantly reduces the computational effort, while retrieving ~ 1 degree of freedom. The biggest benefit of this method is that the retrieval of an aerosol profile, which usually is necessary for deriving a trace gas concentration from MAX-DOAS dSCDs, is not needed. The method is applied to NO2 MAX-DOAS dSCDs recorded during the Mexico City Metropolitan Area 2006 (MCMA-2006) measurement campaign. The retrieved volume mixing ratios of two elevation angles (1° and 3°) are compared to volume mixing ratios measured by two long-path (LP)-DOAS instruments located at the same site. Measurements are found to agree well during times when vertical mixing is expected to be strong. However, inhomogeneities in the air mass above Mexico City can be detected by exploiting the different horizontal and vertical dimensions probed by the MAX-DOAS and LP-DOAS instruments. In particular, a vertical gradient in NO2 close to the ground can be observed in the afternoon, and is attributed to reduced mixing coupled with near-surface emission inside street canyons. The existence of a vertical gradient in the lower 250 m during parts of the day shows the general challenge of sampling the boundary layer in a representative way, and emphasizes the need of vertically resolved measurements.


2005 ◽  
Vol 62 (4) ◽  
pp. 1053-1071 ◽  
Author(s):  
Zhonghai Jin ◽  
Thomas P. Charlock ◽  
Ken Rutledge ◽  
Glenn Cota ◽  
Ralph Kahn ◽  
...  

Abstract Spectral and broadband radiances and irradiances (fluxes) were measured from surface, airborne, and spaceborne platforms in the Chesapeake Lighthouse and Aircraft Measurements for Satellites (CLAMS) campaign. The radiation data obtained on the 4 clear days over ocean during CLAMS are analyzed here with the Coupled Ocean–Atmosphere Radiative Transfer (COART) model. The model is successively compared with observations of broadband fluxes and albedos near the ocean surface from the Clouds and the Earth's Radiant Energy System (CERES) Ocean Validation Experiment (COVE) sea platform and a low-level OV-10 aircraft, of near-surface spectral albedos from COVE and OV-10, of broadband radiances at multiple angles and inferred top-of-atmosphere (TOA) fluxes from CERES, and of spectral radiances at multiple angles from Airborne Multiangle Imaging Spectroradiometer (MISR), or “AirMISR,” at 20-km altidude. The radiation measurements from different platforms are shown to be consistent with each other and with model results. The discrepancies between the model and observations at the surface are less than 10 W m−2 for downwelling and 2 W m−2 for upwelling fluxes. The model–observation discrepancies for shortwave ocean albedo are less than 8%; some discrepancies in spectral albedo are larger but less than 20%. The discrepancies between low-altitude aircraft and surface measurements are somewhat larger than those between the model and the surface measurements; the former are due to the effects of differences in height, aircraft pitch and roll, and the noise of spatial and temporal variations of atmospheric and oceanic properties. The discrepancy between the model and the CERES observations for the upwelling radiance is 5.9% for all angles; this is reduced to 4.9% if observations within 15° of the sun-glint angle are excluded. The measurements and model agree on the principal impacts that ocean optical properties have on upwelling radiation at low levels in the atmosphere. Wind-driven surface roughness significantly affects the upwelling radiances measured by aircraft and satellites at small sun-glint angles, especially in the near-infrared channel of MISR. Intercomparisons of various measurements and the model show that most of the radiation observations in CLAMS are robust, and that the coupled radiative transfer model used here accurately treats scattering and absorption processes in both the air and the water.


2021 ◽  
Author(s):  
Boris Barja González ◽  
Jorge Rosas Santana ◽  
Victoria Eugenia Cachorro ◽  
Carlos Toledano ◽  
Juan Carlos Antuña-Marrero ◽  
...  

The effects of cumulus clouds (Cu) and the combination of stratocumulus-cumulus clouds (Sc-Cu) on the solar radiation at the earth’s surface were evaluated at Camagüey (Cuba) for a 6-year period (from June 2010 to May 2016). Two methods to calculate the cloud radiative effect (CRE) were employed. The first method (CREm) uses solar irradiances in cloudy conditions from actinometric observations, where cloud information was also reported by visual observation. In the second method (CRE0) surface solar irradiances were estimated for both cloudy and clear sky conditions using a 1-D radiative transfer model, and cloud optical depth (COD) retrieved from an AERONET sun-photometer as the main input. A temporal correspondence criterion between COD retrievals and actinometric observations was performed in order to classify COD of each cloud type. After the application of this criterion, the COD belonging to the optically thin clouds was removed. Finally, 255 and 732 COD observations for Cu and Sc-Cu respectively, were found. Results show statistically significant difference at the 95% confidence level between CRE calculated for Sc-Cu and Cu, using the both methods. Mean values of CREm and CRE0 for Cu (Sc-Cu) were −442 W/m2 (−390 W/m2) and −460 W/m2 (−417 W/m2), respectively. CRE0 shows a linear relation with ln(COD), with stronger correlation at lower solar zenith angle. The shortwave cloud effect efficiency (CEE) for the two cloud types sharply decreases with the increase of the COD up to the value of 20. For larger COD, the CEE is less sensitive to the increase of COD.


Sign in / Sign up

Export Citation Format

Share Document