Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “invariant”

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.

At least 1,063 records · Page 59

Generating Vegetation Leaf Area Index Earth System Data Record from Multiple Sensors: Theory - Part 1

The generation of multi-decade long Earth System Data Records (ESDRs) of Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation absorbed by vegetation (FPAR) from remote sensing measurements of multiple sensors is key to monitoring long-term changes in vegetation due to natural and anthropogenic influences. Challenges in developing such ESDRs include problems in remote sensing science (modeling of variability in global vegetation, scaling, atmospheric correction) and sensor hardware (differences in spatial resolution, spectral bands, calibration, and information content). In this paper, we develop a physically based approach for deriving LAI and FPAR products from the Advanced Very High Resolution Radiometer (AVHRR) data that are of comparable quality to the Moderate resolution Imaging Spectroradiometer (MODIS) LAI and FPAR products, thus realizing the objective of producing a long (multi-decadal) time series of these products. The approach is based on the radiative transfer theory of canopy spectral invariants which facilitates parameterization of the canopy spectral bidirectional reflectance factor (BRF). The methodology permits decoupling of the structural and radiometric components and obeys the energy conservation law. The approach is applicable to any optical sensor, however, it requires selection of sensor-specific values of configurable parameters, namely, the single scattering albedo and data uncertainty. According to the theory of spectral invariants, the single scattering albedo is a function of the spatial scale, and thus, accounts for the variation in BRF with sensor spatial resolution. Likewise, the single scattering albedo accounts for the variation in spectral BRF with sensor bandwidths. The second adjustable parameter is data uncertainty, which accounts for varying information content of the remote sensing measurements, i.e., Normalized Difference Vegetation Index (NDVI, low information content), vs. spectral BRF (higher information content). Implementation of this approach indicates good consistency in LAI values retrieved from NDVI (AVHRRmode) and spectral BRF (MODIS-mode). Specific details of the implementation and evaluation of the derived products are detailed in the second part of this two-paper series.

Ganguly, Sangram↗

Maximizing TDRS Command Load Lifetime

The GNC software onboard ISS utilizes TORS command loads, and a simplistic model of TORS orbital motion to generate onboard TORS state vectors. Each TORS command load contains five "invariant" orbital elements which serve as inputs to the onboard propagation algorithm. These elements include semi-major axis, inclination, time of last ascending node crossing, right ascension of ascending node, and mean motion. Running parallel to the onboard software is the TORS Command Builder Tool application, located in the JSC Mission Control Center. The TORS Command Builder Tool is responsible for building the TORS command loads using a ground TORS state vector, mirroring the onboard propagation algorithm, and assessing the fidelity of current TORS command loads onboard ISS. The tool works by extracting a ground state vector at a given time from a current TORS ephemeris, and then calculating the corresponding "onboard" TORS state vector at the same time using the current onboard TORS command load. The tool then performs a comparison between these two vectors and displays the relative differences in the command builder tool GUI. If the RSS position difference between these two vectors exceeds the tolerable lim its, a new command load is built using the ground state vector and uplinked to ISS. A command load's lifetime is therefore defined as the time from when a command load is built to the time the RSS position difference exceeds the tolerable limit. From the outset of TORS command load operations (STS-98), command load lifetime was limited to approximately one week due to the simplicity of both the onboard propagation algorithm, and the algorithm used by the command builder tool to generate the invariant orbital elements. It was soon desired to extend command load lifetime in order to minimize potential risk due to frequent ISS commanding. Initial studies indicated that command load lifetime was most sensitive to changes in mean motion. Finding a suitable value for mean motion was therefore the key to achieving this goal. This goal was eventually realized through development of an Excel spreadsheet tool called EMMIE (Excel Mean Motion Interactive Estimation). EMMIE utilizes ground ephemeris nodal data to perform a least-squares fit to inferred mean anomaly as a function of time, thus generating an initial estimate for mean motion. This mean motion in turn drives a plot of estimated downtrack position difference versus time. The user can then manually iterate the mean motion, and determine an optimal value that will maximize command load lifetime. Once this optimal value is determined, the mean motion initially calculated by the command builder tool is overwritten with the new optimal value, and the command load is built for uplink to ISS. EMMIE also provides the capability for command load lifetime to be tracked through multiple TORS ephemeris updates. Using EMMIE, TORS command load lifetimes of approximately 30 days have been achieved.

Brown, Aaron J.↗

The Characterization of Deep Convective Cloud Albedo as a Calibration Target Using MODIS Reflectances

There are over 25 years of historical satellite data available to climate analysis. The historical satellite data needs to be well calibrated, especially in the visible, where there is no onboard calibration on operational satellites. The key to the vicarious calibration of historical satellites relies on invariant targets, such as the moon, Dome C, and deserts. Deep convective clouds (DCC) also show promise of being a stable invariant or predictable target viewable by all satellites, since they behave as solar diffusers. However DCC have not been well characterized for calibration. Ten years of well-calibrated MODIS is now available. DCC can easily be identified using IR thresholds, where the IR calibration can be traced to the onboard black-bodies. The natural variability of DCC albedo will be analyzed geographically and seasonally, especially difference of convection initiated over land or ocean. Functionality between particle size and ozone absorption with DCC albedo will be examined. Although DCC clouds are nearly Lambertion, the angular distribution of reflectances will be sampled and compared with theoretical models. Both Aqua and Terra MODIS DCC angular models will be compared for consistency. Normalizing angular geostationary DCC reflectances, which were calibrated against MODIS, with SCIAMACHY spectral reflectances and comparing them to MODIS DCC reflectances will inspect the usage of DCC albedos as an absolute calibration target.

Doelling, David R.↗

Statistical and Scaling Properties of Remotely-Sensed Soil Moisture in Two Contrasting Domains in the North American Monsoon Region

Characterizing soil moisture (theta) variability is important for inferring high-resolution information from coarse estimates provided by remote sensors. In this study, we analyze the spatial variability and scale invariance of high-resolution theta estimates collected in two contrasting semiarid areas, Arizona (AZ) and Sonora (SON), during the Soil Moisture Experiment - North American Monsoon in 2004 (SMEX04- NAME). Results reveal that as the mean theta condition ( ) becomes drier, the spatial standard deviation becomes smaller in both domains. The coefficient of variation of theta decreases with in SON, but does not display a clear tendency with in AZ. We also found the presence of scale invariance and multifractality in the range of support scales from 51.2 km to 0.8 km for all soil moisture fields in the two regions. The multifractal properties of theta are clearly linked to in SON, while the relation is affected by more dispersion in AZ. We argue this is due to differences in the dynamic (rainfall) and static (vegetation) controls on theta in the two domains.

Mascaro, Giuseppe↗

Image Segmentation, Registration, Compression, and Matching

A novel computational framework was developed of a 2D affine invariant matching exploiting a parameter space. Named as affine invariant parameter space (AIPS), the technique can be applied to many image-processing and computer-vision problems, including image registration, template matching, and object tracking from image sequence. The AIPS is formed by the parameters in an affine combination of a set of feature points in the image plane. In cases where the entire image can be assumed to have undergone a single affine transformation, the new AIPS match metric and matching framework becomes very effective (compared with the state-of-the-art methods at the time of this reporting). No knowledge about scaling or any other transformation parameters need to be known a priori to apply the AIPS framework. An automated suite of software tools has been created to provide accurate image segmentation (for data cleaning) and high-quality 2D image and 3D surface registration (for fusing multi-resolution terrain, image, and map data). These tools are capable of supporting existing GIS toolkits already in the marketplace, and will also be usable in a stand-alone fashion. The toolkit applies novel algorithmic approaches for image segmentation, feature extraction, and registration of 2D imagery and 3D surface data, which supports first-pass, batched, fully automatic feature extraction (for segmentation), and registration. A hierarchical and adaptive approach is taken for achieving automatic feature extraction, segmentation, and registration. Surface registration is the process of aligning two (or more) data sets to a common coordinate system, during which the transformation between their different coordinate systems is determined. Also developed here are a novel, volumetric surface modeling and compression technique that provide both quality-guaranteed mesh surface approximations and compaction of the model sizes by efficiently coding the geometry and connectivity/topology components of the generated models. The highly efficient triangular mesh compression compacts the connectivity information at the rate of 1.5-4 bits per vertex (on average for triangle meshes), while reducing the 3D geometry by 40-50 percent. Finally, taking into consideration the characteristics of 3D terrain data, and using the innovative, regularized binary decomposition mesh modeling, a multistage, pattern-drive modeling, and compression technique has been developed to provide an effective framework for compressing digital elevation model (DEM) surfaces, high-resolution aerial imagery, and other types of NASA data.

Yadegar, Jacob↗

Wavespace-Based Coherent Deconvolution

Array deconvolution is commonly used in aeroacoustic analysis to remove the influence of a microphone array's point spread function from a conventional beamforming map. Unfortunately, the majority of deconvolution algorithms assume that the acoustic sources in a measurement are incoherent, which can be problematic for some aeroacoustic phenomena with coherent, spatially-distributed characteristics. While several algorithms have been proposed to handle coherent sources, some are computationally intractable for many problems while others require restrictive assumptions about the source field. Newer generalized inverse techniques hold promise, but are still under investigation for general use. An alternate coherent deconvolution method is proposed based on a wavespace transformation of the array data. Wavespace analysis offers advantages over curved-wave array processing, such as providing an explicit shift-invariance in the convolution of the array sampling function with the acoustic wave field. However, usage of the wavespace transformation assumes the acoustic wave field is accurately approximated as a superposition of plane wave fields, regardless of true wavefront curvature. The wavespace technique leverages Fourier transforms to quickly evaluate a shift-invariant convolution. The method is derived for and applied to ideal incoherent and coherent plane wave fields to demonstrate its ability to determine magnitude and relative phase of multiple coherent sources. Multi-scale processing is explored as a means of accelerating solution convergence. A case with a spherical wave front is evaluated. Finally, a trailing edge noise experiment case is considered. Results show the method successfully deconvolves incoherent, partially-coherent, and coherent plane wave fields to a degree necessary for quantitative evaluation. Curved wave front cases warrant further investigation. A potential extension to nearfield beamforming is proposed.

Bahr, Christopher J.↗

Dipole Alignment in Rotating MHD Turbulence

We present numerical results from long-term CPU and GPU simulations of rotating, homogeneous, magnetohydrodynamic (MHD) turbulence, and discuss their connection to the spherically bounded case. We compare our numerical results with a statistical theory of geodynamo action that has evolved from the absolute equilibrium ensemble theory of ideal MHD turbulence, which is based on the ideal MHD invariants are energy, cross helicity and magnetic helicity. However, for rotating MHD turbulence, the cross helicity is no longer an exact invariant, although rms cross helicity becomes quasistationary during an ideal MHD simulation. This and the anisotropy imposed by rotation suggests an ansatz in which an effective, nonzero value of cross helicity is assigned to axisymmetric modes and zero cross helicity to non-axisymmetric modes. This hybrid statistics predicts a large-scale quasistationary magnetic field due to broken ergodicity , as well as dipole vector alignment with the rotation axis, both of which are observed numerically. We find that only a relatively small value of effective cross helicity leads to the prediction of a dipole moment vector that is closely aligned (less than 10 degrees) with the rotation axis. We also discuss the effect of initial conditions, dissipation and grid size on the numerical simulations and statistical theory.

Shebalin, John V.↗

A New Formulation for Hybrid LES-RANS Computations

Ideally, a hybrid LES-RANS computation would employ LES only where necessary to make up for the failure of the RANS model to provide sufficient accuracy or to provide time-dependent information. Current approaches are fairly restrictive in the placement of LES and RANS regions; an LES-RANS transition in a boundary layer, for example, yields an unphysical log-layer shift. A hybrid computation is formulated here to allow greater control over the placement of LES and RANS regions and the transitions between them. The concept of model invariance is introduced, which provides a basis for interpreting hybrid results within an LES-RANS transition zone. Consequences of imposing model invariance include the addition of terms to the governing equations that compensate for unphysical gradients created as the model changes between RANS and LES. Computational results illustrate the increased accuracy of the approach and its insensitivity to the location of the transition and to the blending function employed.

Woodruff, Stephen L.↗

Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes a hybrid approach for determining physical displacements by augmenting the original quadratically constrained least squares (LSQI) algorithm with Direct Shape Mapping (DSM) and modifying the energy constraints. The approach presented is applicable to simulation of the elastic behavior of launch vehicles and other structures that utilize discrete LTI finite element model (FEM) derived mode sets (eigenvalues and eigenvectors) that are propagated throughout time. The time invariant nature of the elastic data presents a problem of how to properly transition elastic states from the prior to the new model while preserving motion across the transition and ensuring there is no truncation or excitation of the system. A previous approach utilizes a LSQI algorithm with an energy constraint to effect smooth transitions between eigenvector sets with no requirement that the models be of similar dimension or have any correlation. This approach assumes energy is conserved across the transition, which results in significant non-physical transients due to changing quasi-steady state energy between mode sets, a phenomenon seen when utilizing a truncated mode set. The computational burden of simulating a full mode set is significant so a subset of modes is often selected to reduce run time. As a result of this truncation, energy between mode sets may not be constant and solutions across transitions could produce non-physical transients. In an effort to abate these transients an improved methodology was developed based on the aforementioned approach, but this new approach can handle significant changes in energy across mode set transitions. It is proposed that physical velocities due to elastic behavior be solved for using the LSQI algorithm, but solve for displacements using a two-step process that independently addresses the quasi-steady-state and non-steady-state contributions to the elastic displacement. For structures subject to large external forces, such as thrust or atmospheric drag, it is imperative to capture these forces when solving for elastic displacement. To simplify the mathematical formulation, assumptions are made regarding mass matrix normalization, constant external forcing, and constant viscous damping. These simplifications allow for direct solutions to the quasi-steady-state displacements through a process titled Direct Shape Mapping. DSM solves for the displacements using the eigenvalues of the elastic modes and the external forcing and returns a set of elastic displacements dictated by the eigenvectors of the post-transition mode set. For the non-steady-state contributions to displacement we formulate a LSQI problem that is constrained by energy of the non-steady state terms. The contributions from the quasi-steady-state and non-steady state solutions are then combined to obtain the physical displacements associated with the new set of eigenvectors. Results for the LSQI-DSM approach show significant reduction/complete removal of transients across mode set transitions while maintaining elastic motion from the prior state. For time propagation applications employing discrete elastic models that need to be transitioned in time and where running with full a full mode set is not feasible, the method developed offers a practical solution to simulating vehicle elasticity.

Hannan, Mike R.↗

Elastic Model Transitions: a Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes both quadratically constrained least squares (LSQI) and Direct Shape Mapping (DSM) algorithms to determine physical displacements. This approach is applicable to the simulation of the elastic behavior of launch vehicles and other structures that utilize multiple LTI finite element model (FEM) derived mode sets that are propagated throughout time. The time invariant nature of the elastic data for discrete segments of the launch vehicle trajectory presents a problem of how to properly transition between models while preserving motion across the transition. In addition, energy may vary between flex models when using a truncated mode set. The LSQI-DSM algorithm can accommodate significant changes in energy between FEM models and carries elastic motion across FEM model transitions. Compared with previous approaches, the LSQI-DSM algorithm shows improvements ranging from a significant reduction to a complete removal of transients across FEM model transitions as well as maintaining elastic motion from the prior state.

Jurenko, Robert J.↗

Coordinates for Mapping the Distribution of Magnetically Trapped Particles

Dipole representations of the earth's magnetic field have been found to have insufficient accuracy for the study of magnetically trapped particles. A coordinate system consisting of the magnitude of the magnetic field B, and the integral invariant I has been found to adequately organize, measurements made at different geographic locations. It is shown in the present paper, that a parameter L = f(B,I) can be defined which retains most of the desirable properties of I and which has the additional property of organizing measurements along lines of force. Since the parameter L is the analog of a physical distance in a dipole field (the equatorial radius of a magnetic shell), it is usually found to present fewer conceptual difficulties than the integral invariant I.

McIlwain, Carl E.↗

BRDF Characterization and Calibration Inter-Comparison between Terra MODIS, Aqua MODIS, and S-NPP VIIRS

The inter-comparison of reflective solar bands (RSB) between Terra MODIS, Aqua MODIS, and SNPP VIIRS is very important for assessment of each instrument's calibration and to identify calibration improvements. One of the limitations of using their ground observations for the assessment is a lack of the simultaneous nadir overpasses (SNOs) over selected pseudo-invariant targets. In addition, their measurements over a selected Earth view target have significant difference in solar and view angles, and these differences magnify the effects of Bidirectional Reflectance Distribution Function (BRDF). In this work, an inter-comparison technique using a semi-empirical BRDF model is developed for reflectance correction. BRDF characterization requires a broad coverage of solar and view angles in the measurements over selected pseudo-invariant targets. Reflectance measurements over Libya 1, 2, and 4 desert sites from both the Aqua and Terra MODIS are regressed to a BRDF model with an adjustable coefficient accounting for the calibration difference between the two instruments. The BRDF coefficients for three desert sites for MODIS bands 1 to 9 are derived and the wavelength dependencies are presented. The analysis and inter-comparison are for MODIS bands 1 to 9 and VIIRS moderate resolution radiometric bands (M bands) M1, M2, M4, M5, M7, M8, M10 and imaging bands (I bands) I1-I3. Results show that the ratios from different sites are in good agreement. The ratios between Terra and Aqua MODIS from year 2003 to 2014 are presented. The inter-comparison between MODIS and VIIRS are analyzed for year 2014.

Chang, Tiejun↗

A Consistent EPIC Visible Channel Calibration using VIIRS and MODIS as a Reference

The Earth Polychromatic Imaging Camera (EPIC) aboard the Deep Space Climate Observatory (DSCOVR) satellite constantly images the sunlit disk of Earth from the Lagrange-1 (L1) point in 10 spectral channels spanning the UV, VIS, and NIR spectrums. Recently, the DSCOVR EPIC team has publicly released version 2 dataset, which has implemented improved navigation, stray-light correction, and flat-fielding of the CCD array. The EPIC 2-year data record must be well-calibrated for consistent cloud, aerosol, trace gas, land use and other retrievals. Because EPIC lacks onboard calibrators, the observations made by EPIC channels must be calibrated vicariously using the coincident measurements from radiometrically stable instruments that have onboard calibration systems. MODIS and VIIRS are best-suited instruments for this task as they contain similar spectral bands that are well-calibrated onboard using solar diffusers and lunar tracking. We have previously calibrated the EPIC version 1 dataset by using EPIC and VIIRS angularly matched radiance pairs over both all-sky ocean and deep convective clouds (DCC). We noted that the EPIC image required navigations adjustments, and that the EPIC stray-light correction provided an offset term closer to zero based on the linear regression of the EPIC and VIIRS ray-matched radiance pairs. We will evaluate the EPIC version 2 navigation and stray-light improvements using the same techniques. In addition, we will monitor the EPIC channel calibration over the two years for any temporal degradation or anomalous behavior. These two calibration methods will be further validated using desert and DCC invariant Earth targets. The radiometric characterization of the selected invariant targets is performed using multiple years of MODIS and VIIRS measurements. Results of these studies will be shown at the conference.

Haney, C. O.↗

Terra and Aqua MODIS Intercomparison Using LEO-GEO Double Difference Method

The Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Terra and Aqua satellites havesuccessfully operated since their launch in 1999 and in 2002, providing more than 18 and 16 years ofcontinuous global observations, respectively. The inter-comparison between the two MODIS instruments can bevery supportive for the instrument calibration and uncertainty assessment. Aqua and Terra MODIS have almostidentical relative spectral response, spatial resolution, and dynamic range for each band. Therefore, a sitedependent correction for a sensor spectral band pair is not necessary for their comparison. However, Terra is inthe morning orbit with an equator crossing time of 10:30 am, and Aqua is in the afternoon orbit with equatorcrossing time of 1:30 pm. Consequently, there is a dearth of simultaneous nadir overpasses (SNOs)between the two satellites. Major challenges in cross-sensor comparison of instruments on different satellitesinclude differences in observation time, solar angle, and view angle over selected pseudo-invariant sites.In this work, the inter-comparisons of thermal emissive bands are performed over a pseudo-invariant target,using the observations from a sensor onboard a geostationary satellite as a bridge. Himawari8 was launched onOctober 7, 2014. The Advanced Himawari Imager (AHI) onboard Himawari8 can be used as a reference tobridge the comparison between Terra and Aqua MODIS. AHI has 16 channels; with spatial resolutions from 0.5km to 2 km at nadir and produces a full disk observations every 10 minutes. The band spectral coveragematchup, comparable spatial resolution and near-simultaneous observation between MODIS and AHI providefeasibility to implement a double difference method. This comparison method minimizes the impact of thedifference in observation time and solar angle. The comparison results will be used as an assessment for MODISinstrument calibration and will be helpful for future enhancement of the L1B product.

Himawari8↗

Cloud Edge Properties Measured by the ARM Shortwave Spectrometer over Ocean and Land

We use the spectrally-invariant method to study the variability of cloud optical thickness () and droplet effective radius (reff) in transition zones between the cloudy and clear sky columns observed by Shortwave Array Spectroradiometer-Zenith (SASZe) at the Southern Great Plains Central Facility site (SGP C1) and during the Marine ARM GPCI Investigation of Clouds (MAGIC) field campaign. The spectrally-invariant method approximates the spectra in the transition zone as a linear combination of definitely clear and definitely cloudy spectra. The slope and intercept of the linear relations characterize and reff in the transition region, respectively. The radiative transfer model simulations show that that (i) the slope of the visible (VIS) band is positively correlated with while (ii) the intercept of the near-infrared (NIR) band has a high negative correlation with reff. We have analyzed 22 cloud edge cases from the SGP and MAGIC and found that from cloud to clear in the transitions (a) the slopes of the VIS band decrease, indicating the decrease of towards cloud edges; (b) the intercepts of the NIR band show a much more significant increase at the SGP than from the MAGIC. The results from observed cases suggest that while decreases for all cases, the decrease in reff is much more significant for cloud over land at the SGP site compared to the ocean counterpart during the MAGIC campaign.

Yang, Weidong↗

TPSAS-NF1676L-11363-DND

The key to the vicarious calibration of historical satellites relies on invariant targets, such as the moon, Dome C, and deserts. Deep convective clouds (DCC) also show promise of being a stable invariant or predictable target viewable by all satellites, since they behave as solar diffusers. However DCC have not been well characterized for calibration. Ten years of well-calibrated MODIS is now available. DCC can easily be identified using IR thresholds, where the IR calibration can be traced to the onboard blackbodies. The natural variability of DCC albedo will be analyzed geographically and seasonally, especially difference of convection initiated over land or ocean.

D R Doelling↗

Radiometric Assessment of the First Three Years of the NOAA-20 VIIRS Reflective Solar Bands Calibration

NASA’s Clouds and the Earth’s Radiant Energy System (CERES) SYN1deg Ed4.1 product utilizes geostationary (GEO) satellite measured radiances and retrieved cloud properties to account for the regional diurnal fluctuations in the Earth’s radiant broadband fluxes for times between the CERES measurements gathered from the Aqua (1:30 PM) and Terra (10:30 AM) sun-synchronous satellites. In CERES Edition 4 products, a global uniformity in the cloud properties and computed surface fluxes across the GEO satellite domains is maintained by scaling the radiance observations from more than twenty GEO visible imagers in the CERES record to a common radiometric reference scale, i.e., Aqua-MODIS. With the new-generation GEO (Himawari-8/9 and GOES-16/17) imagers having multiple reflective solar bands (RSB) that are spectrally similar to those of VIIRS, the CERES Imager and Geostationary Calibration Group (IGCG) is preparing to use NOAA-20 VIIRS as the reference imager for maintaining the radiometric uniformity across the GEO imager constellation. Given the recent Aqua satellite anomaly, this transition may happen sooner than the projected date of Aqua’s de-orbit. This paper presents an independent performance evaluation of the first three years of the NOAA-20 VIIRS RSB calibration in the NASA VIIRS Land Science Investigator-led Processing System (Land SIPS) L1b Collection 2 dataset. The temporal radiometric stability is assessed using multiple invariant Earth targets, including tropical deep convective clouds and the Saharan desert. The invariant target anisotropic reflectance at the top of atmosphere was modeled using five years of stable satellite observations acquired from the previous VIIRS instrument onboard the Suomi National Polar-orbiting Partnership (S-NPP) satellite. The anisotropic corrections are essential for detecting temporal trends with a high statistical confidence. In addition, the radiometric consistency between the RSB of the two VIIRS instruments will be evaluated.

Rajendra Bhatt↗

Radiometric Assessment of the First Three Years of the NOAA-20 VIIRS Reflective Solar Bands Calibration

NASA’s Clouds and the Earth’s Radiant Energy System (CERES) SYN1deg Ed4.1 product utilizes geostationary (GEO) satellite measured radiances and retrieved cloud properties to account for the regional diurnal fluctuations in the Earth’s radiant broadband fluxes for times between the CERES measurements gathered from the Aqua (1:30 PM) and Terra (10:30 AM) sun-synchronous satellites. In CERES Edition 4 products, a global uniformity in the cloud properties and computed surface fluxes across the GEO satellite domains is maintained by scaling the radiance observations from more than twenty GEO visible imagers in the CERES record to a common radiometric reference scale, i.e., Aqua-MODIS. With the new-generation GEO (Himawari-8/9 and GOES-16/17) imagers having multiple reflective solar bands (RSB) that are spectrally similar to those of VIIRS, the CERES Imager and Geostationary Calibration Group (IGCG) is preparing to use NOAA-20 VIIRS as the reference imager for maintaining the radiometric uniformity across the GEO imager constellation. Given the recent Aqua satellite anomaly, this transition may happen sooner than the projected date of Aqua’s de-orbit. This paper presents an independent performance evaluation of the first three years of the NOAA-20 VIIRS RSB calibration in the NASA VIIRS Land Science Investigator-led Processing System (Land SIPS) L1b Collection 2 dataset. The temporal radiometric stability is assessed using multiple invariant Earth targets, including tropical deep convective clouds and the Saharan desert. The invariant target anisotropic reflectance at the top of the atmosphere was modeled using five years of stable satellite observations acquired from the previous VIIRS instrument onboard the Suomi National Polar-orbiting Partnership (S-NPP) satellite. The anisotropic corrections are essential for detecting temporal trends with high statistical confidence. In addition, the radiometric consistency between the RSB of the two VIIRS instruments will be evaluated.

Rajendra Bhatt↗