Polarimetric properties of the lunar surface and its interpretation. VI - Albedo determinations from polarimetric characteristics
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
A method for unsupervised segmentation of polarimetric synthetic aperture radar (SAR) data into classes of homogeneous microwave polarimetric backscatter characteristics is presented. Classes of polarimetric backscatter are selected on the basis of a multidimensional fuzzy clustering of the logarithm of the parameters composing the polarimetric covariance matrix. The clustering procedure uses both polarimetric amplitude and phase information, is adapted to the presence of image speckle, and does not require an arbitrary weighting of the different polarimetric channels; it also provides a partitioning of each data sample used for clustering into multiple clusters. Given the classes of polarimetric backscatter, the entire image is classified using a maximum a posteriori polarimetric classifier. Four-look polarimetric SAR complex data of lava flows and of sea ice acquired by the NASA/JPL airborne polarimetric radar (AIRSAR) are segmented using this technique. The results are discussed and compared with those obtained using supervised techniques.
The use of synthetic aperture radar (SAR) polarimetric processing for enhancing the detection, classification, and/or identification of scene scattering features is described. This paper describes and compares the results obtained from applying the polarimetric techniques that have been developed and/or advanced by researchers at MIT Lincoln Laboratory, and basic polarimetric match filter (PMF) techniques that have been modified to treat distributed scatterers. The paper describes the individual polarimetric processing techniques and the formulation used fro their application to polarimetric imagery obtained from the NAWC P-3 SAR. In all cases, the polarimetric techniques have been applied to enhance the detection of distributed scatterers in clutter The emphasis arises from considering most complex scatterers of interest (either man made or natural) as distributed scatterers consisting of multiple scattering centers, and many modern polarimetric radar systems have the resolution performance to resolve the multiple scattering centers. Furthermore, the individual scattering centers of a distributed target can often exhibit different polarimetric scattering characteristics, and consequently do not respond favorably to polarimetric processing techniques derived for individual point scatterers. The treatment of distributed scatterers also includes concepts for tuning the polarimetric ensemble response of individual scattering centers.
The Global Precipitation Measurement (GPM) Mission Iowa Flood Studies (IFloodS) campaign was conducted in central and northeastern Iowa during the months of April-June, 2013. Specific science objectives for IFloodS included quantification of uncertainties in satellite and ground-based estimates of precipitation, 4-D characterization of precipitation physical processes and associated parameters (e.g., size distributions, water contents, types, structure etc.), assessment of the impact of precipitation estimation uncertainty and physical processes on hydrologic predictive skill, and refinement of field observations and data analysis approaches as they pertain to future GPM integrated hydrologic validation and related field studies. In addition to field campaign archival of raw and processed satellite data (including precipitation products), key ground-based platforms such as the NASA NPOL S-band and D3R Ka/Ku-band dual-polarimetric radars, University of Iowa X-band dual-polarimetric radars, a large network of paired rain gauge platforms, and a large network of 2D Video and Parsivel disdrometers were deployed. In something of a canonical approach, the radar (NPOL in particular), gauge and disdrometer observational assets were deployed to create a consistent high-quality distributed (time and space sampling) radar-based ground "reference" rainfall dataset, with known uncertainties, that could be used for assessing the satellite-based precipitation products at a range of space/time scales. Subsequently, the impact of uncertainties in the satellite products could be evaluated relative to the ground-benchmark in coupled weather, land-surface and distributed hydrologic modeling frameworks as related to flood prediction. Relative to establishing the ground-based "benchmark", numerous avenues were pursued in the making and verification of IFloodS "reference" dual-polarimetric radar-based rain maps, and this study documents the process and results as they pertain specifically to efforts using the NPOL radar dataset. The initial portions of the "process" involved dual-polarimetric quality control procedures which employed standard phase and correlation-based approaches to removal of clutter and non-meteorological echo. Calculation of a scale-adaptive KDP was accomplished using the method of Wang and Chandrasekar (2009; J. Atmos. Oceanic Tech.). A dual-polarimetric blockage algorithm based on Lang et al. (2009; J. Atmos. Oceanic Tech.) was then implemented to correct radar reflectivity and differential reflectivity at low elevation angles. Next, hydrometeor identification algorithms were run to identify liquid and ice hydrometeors. After the quality control and data preparation steps were completed several different dual-polarimetric rain estimation algorithms were employed to estimate rainfall rates using rainfall scans collected approximately every two to three minutes throughout the campaign. These algorithms included a polarimetrically-tuned Z-R algorithm that adjusts for drop oscillations (via Bringi et al., 2004, J. Atmos. Oceanic Tech.), and several different hybrid polarimetric variable approaches, including one that made use of parameters tuned to IFloodS 2D Video Disdrometer measurements. Finally, a hybrid scan algorithm was designed to merge the rain rate estimates from multiple low level elevation angle scans (where blockages could not be appropriately corrected) in order to create individual low-level rain maps. Individual rain maps at each time step were subsequently accumulated over multiple time scales for comparison to gauge network data. The comparison results and overall error character depended strongly on rain event type, polarimetric estimator applied, and range from the radar. We will present the outcome of these comparisons and their impact on constructing composited "reference" rainfall maps at select time and space scales.
In the past, wind measurements from space using fan-beam antennas, such as Seasat Scatterometer (SASS-1), ERS-1 &2, and NASA scatterometer (NSCAT), required up to six large stick-like antennas and suffered a nadir gap of up to 400 km. In the near future, a spinning pencil-beam scatterometer system is to be used for the SeaWinds scatterometer on QuikSCAT (QSCAT) and on ADEOS-2 (SeaWinds). This scatterometer, though offering wind measurements in the nadir region, still suffers from degraded performance in the nadir and outer swath. The purpose of this paper is to present an advanced polarimetric spinning pencil-beam scatterometer system, which can significantly improve the wind performance across the entire swath. The polarimetric scatterometer simultaneously measures co-polarized backscatter and the polarimetric correlation of co- and cross-polarized radar returns from the ocean surface. The advantage over the conventional scatterometer system is that, while the co-polarization radar returns are even function of the wind direction, the polarimetric correlation is an odd function of wind direction due to the reflection symmetry of the wind roughened surface. Therefore, this polarimetric scatterometer system can provide additional, equivalent measurements at azimuth angle 45degree away from the corresponding co-polarization measurements. The combined co-polarization and correlation measurements enable good wind performance across the whole swath to be obtained. In this paper, we will first present the theoretical formulation of all of the key components required for designing a polarimetric scatterometer. Then, we show that good wind performance can be achieved by a slight improvement in the signal-to-noise ratio of the current QSCAT/SeaWinds design. We then present the predicated wind performance using computer simulation based on a model function for the co-polarized backscatter obtained from actual spaceborne scatterometer data and an estimated model function for the polarimetric correlation based on the asymmetry observed in backscatter data. Finally, we will show that, aside from ocean applications, this polarimetric scatterometer can also be used for ice and land applications.
An unsupervised selection of polarimetric features useful for the segmentation and analysis of polarimetric synthetic aperture radar (SAR) data is presented. The technique is based on multidimensional clustering of the parameters composing the polarimetric covariance matrix of the data. Clustering is performed on the logarithm of these quantities. Once the polarimetric cluster centers have been determined, segmentation of the polarimetric data into regions is performed using a maximum likelihood polarimetric classifier. Segmentation maps are further improved using a Markov random field to describe the statistics of the regions and computing the maximum of the product of the local conditional densities. Examples with real polarimetric SAR imagery are given to illustrate the potential of this method.
Many passive remote-sensing techniques have been developed to retrieve cloud microphysical properties from satellite-based sensors, with the most common approaches being the bispectral and polarimetric techniques. These two vastly different retrieval techniques have been implemented for a variety of polar-orbiting and geostationary satellite platforms, providing global climatological data sets. Prior instrument comparison studies have shown that there are systematic differences between the droplet size retrieval products (effective radius) of bispectral (e.g., MODIS, Moderate Resolution Imaging Spectroradiometer) and polarimetric (e.g., POLDER, Polarization and Directionality of Earth's Reflectances) instruments. However, intercomparisons of airborne bispectral and polarimetric instruments have yielded results that do not appear to be systematically biased relative to one another. Diagnosing this discrepancy is complicated, because it is often difficult for instrument intercomparison studies to isolate differences between retrieval technique sensitivities and specific instrumental differences such as calibration and atmospheric correction. In addition to these technical differences the polarimetric retrieval is also sensitive to the dispersion of the droplet size distribution (effective variance), which could influence the interpretation of the droplet size retrieval. To avoid these instrument-dependent complications, this study makes use of a cloud remote-sensing retrieval simulator. Created by coupling a large-eddy simulation (LES) cloud model with a 1-D radiative transfer model, the simulator serves as a test bed for understanding differences between bispectral and polarimetric retrievals. With the help of this simulator we can not only compare the two techniques to one another (retrieval intercomparison) but also validate retrievals directly against the LES cloud properties. Using the satellite retrieval simulator, we are able to verify that at high spatial resolution (50 m) the bispectral and polarimetric retrievals are highly correlated with one another within expected observational uncertainties. The relatively small systematic biases at high spatial resolution can be attributed to different sensitivity limitations of the two retrievals. In contrast, a systematic difference between the two retrievals emerges at coarser resolution. This bias largely stems from differences related to sensitivity of the two retrievals to unresolved inhomogeneities in effective variance and optical thickness. The influence of coarse angular resolution is found to increase uncertainty in the polarimetric retrieval but generally maintains a constant mean value.
Many passive remote-sensing techniques have been developed to retrieve cloud microphysical properties from satellite-based sensors, with the most common approaches being the bispectral and polarimetric techniques. These two vastly different retrieval techniques have been implemented for a variety of polar-orbiting and geostationary satellite platforms, providing global climatological data sets. Prior instrument comparison studies have shown that there are systematic differences between the droplet size retrieval products (effective radius) of bispectral (e.g., MODIS, Moderate Resolution Imaging Spectroradiometer) and polarimetric (e.g., POLDER, Polarization and Directionality of Earth's Reflectances) instruments. However, intercomparisons of airborne bispectral and polarimetric instruments have yielded results that do not appear to be systematically biased relative to one another. Diagnosing this discrepancy is complicated, because it is often difficult for instrument intercomparison studies to isolate differences between retrieval technique sensitivities and specific instrumental differences such as calibration and atmospheric correction. In addition to these technical differences the polarimetric retrieval is also sensitive to the dispersion of the droplet size distribution (effective variance), which could influence the interpretation of the droplet size retrieval. To avoid these instrument-dependent complications, this study makes use of a cloud remote-sensing retrieval simulator. Created by coupling a large-eddy simulation (LES) cloud model with a 1-D radiative transfer model, the simulator serves as a test bed for understanding differences between bispectral and polarimetric retrievals. With the help of this simulator we can not only compare the two techniques to one another (retrieval intercomparison) but also validate retrievals directly against the LES cloud properties. Using the satellite retrieval simulator, we are able to verify that at high spatial resolution (50m) the bispectral and polarimetric retrievals are highly correlated with one another within expected observational uncertainties. The relatively small systematic biases at high spatial resolution can be attributed to different sensitivity limitations of the two retrievals. In contrast, a systematic difference between the two retrievals emerges at coarser resolution. This bias largely stems from differences related to sensitivity of the two retrievals to unresolved inhomogeneities in effective variance and optical thickness. The influence of coarse angular resolution is found to increase uncertainty in the polarimetric retrieval but generally maintains a constant mean value.
The normalized polarimetric classifier (NPC) which uses only the relative magnitudes and phases of the polarimetric data is proposed for discrimination of terrain elements. The probability density functions (PDFs) of polarimetric data are assumed to have a complex Gaussian distribution, and the marginal PDF of the normalized polarimetric data is derived by adopting the Euclidean norm as the normalization function. The general form of the distance measure for the NPC is also obtained. It is demonstrated that for polarimetric data with an arbitrary PDF, the distance measure of NPC will be independent of the normalization function selected even when the classifier is mistrained. A complex Gaussian distribution is assumed for the polarimetric data consisting of grass and tree regions. The probability of error for the NPC is compared with those of several other single-feature classifiers. The classification error of NPCs is shown to be independent of the normalization function.
Based on symmetry properties, polarimetric remote sensing of geophysical media is studied. From the viewpoint of symmetry groups, media with reflection, rotation, azimuthal, and centrical symmetries are considered. The symmetries impose relations among polarimetric scattering coefficients, which are valid to all scattering mechanisms in the symmetrical configurations. Various orientation distributions of non-spherical scatterers can be identified from the scattering coefficients by a comparison with the symmetry calculations. Experimental observations are then analyzed for many geophysical scenes acquired with the Jet Propulsion Laboratory (JPL) airborne polarimetric SAR at microwave frequencies over sea ice and vegetation. Polarimetric characteristics of different ice types are compared with symmetry behaviors. The polarimetric response of a tropical rain forest reveals characteristics close to the centrical symmetry properties, which can be used as a distributed target to relatively calibrate polarimetric radars without any deployment of manmade calibration targets.
Relationships among polarimetric backscattering coefficients are derived from the viewpoint of symmetry groups. For both reciprocal and non-reciprocal media, symmetry encountered in remote sensing due to reflection, rotation, azimuthal, and centrical symmetry groups is considered. The derived properties are general and valid to all scattering mechanisms, including volume and surface scatterings and their interactions, in a given symmetrical configuration. The scattering coefficients calculated from theoretical models for layer random media and rough surfaces are shown to obey the symmetry relations. Use of symmetry properties in remote sensing of structural and environmental responses of scattering media is also discussed. Orientations of spheroidal scatterers described by spherical, uniform, planophile, plagiothile, erectophile, and extremophile distributions are considered to derive their polarimetric backscattering characteristics. These distributions can be identified from the observed scattering coefficients by comparison with theoretical symmetry calculations. A new parameter is then defined to study scattering structures in geophysical media. Observations from polarimetric data acquired by the Jet Propulsion Laboratory airborne synthetic aperture radar over forests, sea ice, and sea surface are presented. Experimental evidences of the symmetry relationships are shown and their use in polarimetric remote sensing is illustrated. For forests, the coniferous forest in Mt. Shasta area (California) and mixed forest near Presque Isle (Maine) exhibit characteristics of the centrical symmetry at C-band. For sea ice in the Beaufort Sea, multi-year sea ice has a cross-polarized ratio e close to e(sub 0), calculated from symmetry, due to the randomness in the scattering structure. First-year sea ice has e much smaller than e(sub 0) due to the preferential alignment of the columnar structure of the ice. From polarimetric data of a sea surface in the Bering Sea, it is observed that e and e(sub 0) are increasing with incident angle and e is greater than e(sub 0) at L-band because of the directional feature of sea surface waves. Symmetry properties of geophysical media can also be used to calibrate polarimetric radars.
NASA Marshall Space Flight Center and the University of Alabama Huntsville are collaborating with the 45th Weather Squadron (45WS) to develop improved lightning prediction capabilities for the new C-band dual-polarimetric weather radar being acquired for use by 45WS and launch weather forecasters at Cape Canaveral Air Force Station (CCAFS). In particular, these algorithms will focus on lightning onset, cessation and combined lightning-radar applications for convective winds assessment. Research using radar reflectivity (Z) data for prediction of lightning onset has been extensively discussed in the literature and subsequently applied by launch weather forecasters as it pertains to lightning nowcasting. Currently the forecasters apply a relatively straight forward but effective temperature-Z threshold algorithm for assessing the likelihood of lightning onset in a given storm. In addition, a layered VIL above the freezing level product is used as automated guidance for the onset of lightning. Only limited research and field work has been conducted on lightning cessation using Z and vertically-integrated Z for determining cessation. Though not used operationally vertically-integrated Z (basis for VIL) has recently shown promise as a tool for use in nowcasting lightning cessation. The work discussed herein leverages and expands upon these and similar reflectivity-threshold approaches via the application/addition of over two decades of polarimetric radar research focused on distinct multi-parameter radar signatures of ice/mixed-phase initiation and ice-crystal orientation in highly electrified convective clouds. Specifically, our approach is based on numerous previous studies that have observed repeatable patterns in the behavior of the vertical hydrometeor column as it relates to the temporal evolution of differential reflectivity and depolarization (manifested in either LDR or p(sub hv)), development of in-situ mixed and ice phase microphysics, electric fields, and ensuing lightning in the sub-tropical/tropical convection typical of the southeastern U.S., Maritime Continent, and southwestern Amazon. The polarimetric signatures detected in this setting provide a basis for automated 3-D detection of hydrometeor types in fuzzy logic hydrometeor identification algorithms (HID). Our working hypothesis is that improvement in lightning onset warning lead time and specificity for a given storm, relative to application of a Z-threshold algorithm, should arise as a consequence of the ability of dual-polarimetric radar to unambiguously detect and identify (through HID algorithms) the updraft elevation of rain-water cores above the freezing level and subsequent onset of drop freezing, riming, and robust mixed phase processes leading to significant charge separation and lightning. This type of algorithm, though dependent on the quality of the polarimetric data should be less susceptible to variable Z-calibration that can impact a given Z-threshold approach. To facilitate development of the algorithm while the 45WS dual-pol radar is in its current test stages and to evaluate the impact of polarimetric data quality (e.g., modified scan parameters and sampling) on the ensuing algorithms, we are using the ARMOR C-band dual-pol radar in Huntsville combined with N. Alabama LMA data and ARMOR HID algorithms [NCAR algorithm modified for application at C-band] in a testbed fashion. For lightning cessation we are revisiting the application of differential propagation phase variables for the monitoring of ice crystal alignment driven by in-cloud electric fields combined with metrics of ice water path (i.e., vertically integrated reflectivity). Importantly it should be noted that this approach is still very much a research topic and as such, we will explore operational applications that involve radar frequencies other than C-Band by using the UAH MAX X-band dual-pol radar in slow staring modes.
The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made
Estimation of rainfall integral parameters, including radar observables, and empirical relations between them are sensitive to the truncation of the drop size distribution (DSD), particularly at the large drop end. The sensitivity of rainfall integral parameters to the maximum drop diameter (D(sub max)) is exacerbated at C‐band since resonance effects are pronounced for large drops in excess of 5 mm diameter (D). Due to sampling limitations, it is often difficult to reliably estimate D(sub max) with disdrometers. The resulting uncertainties in D(sub max0 potentially increase errors in radar retrieval methods, particularly at C‐band, that rely on disdrometer observations for DSD input to radar models. In fact, D(sub max) is typically an assumed DSD parameter in the development of radar retrieval methods. Because of these very uncertainties, it is difficult to independently confirm disdrometer estimates of D(sub max) with polarimetric radar observations. A couple of approaches can be taken to reduce uncertainty in large drop measurement. Longer integration times can be used for the collection of larger disdrometer samples. However, integration periods must be consistent with a radar resolution volume (RRV) and the temporal and spatial scales of the physical processes affecting the DSD therein. Multiple co‐located disdrometers can be combined into a network to increase the sample size within a RRV. However, over a reasonable integration period, a single disdrometer sample volume is many orders of magnitudes less than a RRV so it is not practical to devise a network of disdrometers that has an equivalent volume to a typical RRV. Since knowledge of DSD heterogeneity and large drop occurrence in time and space is lacking, the specific accuracy or even general representativeness of disdrometer based D(sub max) and large drop concentration estimates within a RRV are currently unknown. To address this complex issue, we begin with a simpler question. Is the frequency of occurrence of large rain drops (D > 5 mm) in disdrometer observations, either stand alone or networked, generally representative and consistent with polarimetric radar observations? We first show from simulations that the concentration of large (D > 5 mm) rain drops (N(sub T5)) can be estimated from polarimetric observations of specific differential phase (K(sub dp)) and differential reflectivity (Z(sub dr)), N(sub T5)=F(K(sub dp),Z(sub dr)), or horizontal reflectivity (Z(sub h)) and Z(sub dr), N(sub T5)=(Z(sub h),Z(sub dr)). We assess the error associated with polarimetric retrieval of N(sub T5), including sensitivity to D(sub max) parameterization assumptions and measurement error in the radar simulations. Polarimetric measurements at S‐band and C‐band will then be used to retrieve estimates of N(sub T5) and compared to disdrometer estimates of N(sub T5). After careful consideration of retrieval error, we will check consistency between disdrometer and polarimetric radar estimates of N(sub T5) and the frequency of occurrence of large rain drops in a variety of precipitating regimes using data from NASA's Global Precipitation Measurement (GPM) Ground Validation (GV) program, including field campaigns such as MC3E (Oklahoma) and IFloodS (Iowa) and extended measurements over Huntsville, Alabama and NASA Wallops Flight Facility in coastal Virginia.
An experiment was carried out in January 1994 at the Geophysical Research Facility in the Cold Regions Research and Engineering Laboratory. To investigate effects on polarimetric scattering signatures of sea ice growth under diurnal temperature variations, an ice sheet was grown for 2.5 days for the thickness of 10 cm and a polarimetric radar operating at C-band was used to obtain backscattering data in conjunction with ice-characterization measurements. The ice sheet was grown in the late morning of January 19, 1994. The initial growth rate was slow due to high insolation and temperature. As the air temperature dropped during the night, the growth rate increased significantly. The air temperature changed drastically from about -10(deg)C to -35(deg)C between day and night. The temperature cycle was repeated during the next day and the growth rate varied in the same manner. The surface of the ice was partially covered by frost flowers and the areal coverage increased as the ice became thicker. Throughout the ice growth duration of 2.5 days, polarimetric backscatter data were collected at roughly every centimeter of ice growth. For each set of radar measurements of saline ice, a set of calibration measurements was carried out with trihedrial corner reflectors and a metallic sphere. Measured polarimetric backscattering coefficients of the ice sheet reveal a strong correlation between radar data and temperature variations. As the temperature increased (decreased), the backscatter increased (decreased) correspondingly. From the ice-characterization data, temperatures of the air, at the ice-air interface, and in the ice layer had the same variation trend. Another interesting experimental observation is that the salinity measured as a function of ice depth from a sample of 10-cm thich ice indicated that the salinity variations had a similar cycle as the temperature; i.e., the salinity profile recorded the history of the temperature variations. Characterization data of the ice sheet are used in a theoretical model for scattering from saline ice with frost cover to explain the observed polarimetric signatures.
A polarimetric scattering model is developed to study effects of snow cover and frost flowers with brine infiltration on thin sea ice. Leads containing thin sea ice in the Artic icepack are important to heat exchange with the atmosphere and salt flux into the upper ocean. Surface characteristics of thin sea ice in leads are dominated by the formation of frost flowers with high salinity. In many cases, the thin sea ice layer is covered by snow, which wicks up brine from sea ice due to capillary force. Snow and frost flowers have a significant impact on polarimetric signatures of thin ice, which needs to be studied for accessing the retrieval of geophysical parameters such as ice thickness. Frost flowers or snow layer is modeled with a heterogeneous mixture consisting of randomly oriented ellipsoids and brine infiltration in an air background. Ice crystals are characterized with three different axial lengths to depict the nonspherical shape. Under the covering multispecies medium, the columinar sea-ice layer is an inhomogeneous anisotropic medium composed of ellipsoidal brine inclusions preferentially oriented in the vertical direction in an ice background. The underlying medium is homogeneous sea water. This configuration is described with layered inhomogeneous media containing multiple species of scatterers. The species are allowed to have different size, shape, and permittivity. The strong permittivity fluctuation theory is extended to account for the multispecies in the derivation of effective permittivities with distributions of scatterer orientations characterized by Eulerian rotation angles. Polarimetric backscattering coefficients are obtained consistently with the same physical description used in the effective permittivity calculation. The mulitspecies model allows the inclusion of high-permittivity species to study effects of brine infiltrated snow cover and frost flowers on thin ice. The results suggest that the frost cover with a rough interface significantly increases the backscatter from thin saline ice and the polarimetric signature becomes closer to the isotropic characteristics. The snow cover also modifies polarimetric signatures of thin sea ice depending on the snow mixture and the interface condition.
This paper introduces the concept and potential value of polarimetric maps and the techniques for generating these maps in operational remote sensing. The application-oriented polarimetric signature analyses in the literature are compiled, and several optical models are illustrated to bring out requirements of a sensor system for polarimetric mapping. By use of the concepts of Stokes parameters the descriptive specification of one sensor system is refined. The descriptive specification for a multichannel digital photometric-polarimetric mapper is based upon our experience with the present single channel device which includes the generation of polarimetric maps and pictures. High photometric accuracy and stability coupled with fast, accurate digital output has enabled us to overcome the handicap of taking sequential data from the same terrain.
A new technique is presented for calibrating polarimetric synthetic aperture radar (SAR) images using only the responses from natural distributed targets. The model for polarimetric radars is assumed to be X = cRST where X is the measured scattering matrix corresponding to the target scattering matrix S distorted by the system matrices T and R (in general T does not equal R(sup t)). To allow for the polarimetric calibration using only distributed targets and corner reflectors, van Zyl assumed a reciprocal polarimetric radar model with T = R(sup t); when applied for JPL SAR data, a heuristic symmetrization procedure is used by POLCAL to compensate the phase difference between the measured HV and VH responses and then take the average of both. This heuristic approach causes some non-removable cross-polarization responses for corner reflectors, which can be avoided by a rigorous symmetrization method based on reciprocity. After the radar is made reciprocal, a new algorithm based on the responses from distributed targets with reflection symmetry is developed to estimate the cross-talk parameters. The new algorithm never experiences problems in convergence and is also found to converge faster than the existing routines implemented for POLCAL. When the new technique is implemented for the JPL polarimetric data, symmetrization and cross-talk removal are performed on a line-by-line (azimuth) basis. After the cross-talks are removed from the entire image, phase and amplitude calibrations are carried out by selecting distributed targets either with azimuthal symmetry along the looking direction or with some well-known volume and surface scattering mechanisms to estimate the relative phases and amplitude responses of the horizontal and vertical channels.