Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface fluxes”

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 523 records · Page 29

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

Christensen, Matthew [Pacific Northwest National L↗

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

54 ENVIRONMENTAL SCIENCES↗

A Multiyear Gridded Data Ensemble of Surface Biogenic Carbon Fluxes for North America: Evaluation and Analysis of Results

Accurate and fine-scale estimates of biogenic carbon fluxes are critical for measuring and monitoring the biosphere's responses and feedback to the climate system. Currently available data products from flux towers and model-intercomparison projects struggle to adequately represent spatiotemporal dynamics of surface biogenic carbon fluxes, and to quantify their uncertainties, which also are crucial to atmospheric inversion systems. To address these gaps, we introduce a new perturbed-parameter model ensemble with the CASA model to estimate surface biogenic carbon fluxes at monthly and 3-hourly scales for North America at ~500-m and 5-km resolutions. We first use the Extended Fourier Amplitude Sensitivity Testing to choose the three most sensitive parameters to be perturbed, maximum light use efficiency (E max ), optimal temperature of photosynthesis (T opt ), and temperature response of respiration (Q 10 ). The initial range for each parameter is broadly sampled for the L1 ensemble, but then we pruned E max with site-level primary productivity to derive an L2 ensemble with narrower uncertainty ranges. Ensembles are strongly correlated with site-level results at both monthly and 3-hourly scales, and the spread across L1/L2 ensemble members encompasses the range of AmeriFlux observations. Monthly variability in the L2 ensemble mean is 85% of the observed variability. The L2 ensemble outperforms diverse data products with the highest Taylor skill scores at diurnal to annual scales. Finally, the ensemble's seasonality agrees well with other models for most biome types and in high and middle latitudes, but inconsistencies are found in subtropical and tropical ecoregions and for annual totals over North America.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Contributions of Land Carbon Fluxes Variability and Atmospheric Transport Variability to Atmospheric CO2 Variability

Much of the uncertainty in predicting variations in greenhouse gases originates in the complex dynamics of the land ecosystem and the atmosphere. To reduce the uncertainties, it is useful to decouple the contributions of land carbon fluxes and atmospheric transport to atmospheric carbon variability. Here we isolate these contributions using a version of the NASA GEOS model that couples carbon, energy and water cycles between the land and the atmosphere. Our current study is a follow-on to a preliminary analysis that suggested that an extreme event (e.g., imposed drought) in a free running AGCM simulation affects atmospheric CO2 more through its impact on atmospheric transport than through a modification of land carbon fluxes. In order to more carefully isolate the effects of the land carbon variability and atmospheric transport variability on atmospheric CO2 variability, we conducted two coupled-AGCM simulations in replay mode, a mode that forces the model's evolution of weather to match that of the MERRA-2 reanalysis. In the control simulation, the land carbon fluxes and the atmospheric CO2 concentrations, as well as the meteorology, are simulated over 2001-2015. We compute the climatological seasonal cycles of net land carbon production from this control simulation and then prescribe, in a second simulation, these climatological carbon fluxes to the atmosphere in the same replay mode. By comparing the atmospheric CO2 variability produced in the two simulations, we fully isolate the part of this variability associated with land surface fluxes. Relative contributions of land flux variability and atmospheric transport variability to CO2 variability are quantified both on a regional basis and as a function of height into the atmosphere to support interpretation of both ground-based and satellite observations.

Lee, Eunjee↗

Vacuum vessel design with lofted toroidal surfaces for a QHS configuration

A vacuum vessel in between the hot plasma and the magnetic coils is crucial for a fusion reactor. In this study, a tool is presented to perform the initial design of such a vessel automatically. Starting from the last closed flux surface, an algorithm is developed which lofts this manifold out as much as possible. At the same time, the algorithm makes sure that no intersection with the coils takes place. In the presented work, this algorithm is applied to the HSX stellarator and it has been verified that the resulting vessel wall can be built from an engineering perspective. Using field line tracing, the strikepoints on the vessel wall indicate the initial locations for a divertor.

HSX↗

TPSAS-NF1676L-30053-DND

The 2016-17 Arctic sea ice growth season (October-March) exhibited one of the lowest end-of-season sea ice volume and extent of any year since 1979. An analysis of MERRA2 atmospheric reanalysis data and CERES radiative flux data reveals that a record warm and moist Arctic atmosphere supported the reduced sea ice growth through two pathways. First, numerous regional episodes of increased atmospheric temperature and moisture, transported from lower latitudes, increased the cumulative energy input from downwelling longwave surface fluxes. Second, in those same episodes, the efficiency that the atmosphere cooled radiatively to space was reduced, increasing the amount of energy retained in the Arctic atmosphere and reradiated back toward the surface. Overall, the Arctic radiative cooling efficiency shows a decreasing trend since 2000. The results presented highlight the increasing importance of atmospheric forcing on sea ice variability demonstrating that episodic Arctic atmospheric rivers, regions of elevated poleward water vapor transport, and the subsequent surface energy budget response is a critical mechanism actively contributing to the evolution of Arctic sea ice.

Bradley Hegyi↗

Marine Boundary Layer Decoupling and the Stable Isotopic Composition of Water Vapor

Decoupling of the subcloud layer in stratocumulus-topped marine boundary layers (STMBL) influences low-cloud cover by limiting the supply of water vapor from the surface. However, the relative importance of mixing between surface fluxes and other reservoirs of water vapor as a function of the degree of decoupling is poorly understood. Water vapor transport within the STMBL and its response to decoupling is explored using surface measurements of water vapor isotopic composition that were obtained from Graciosa Island, Azores, during the summer and fall of 2018. The data show an inverse relationship between the decoupling metric, Δq, and the degree of drying from an isotopically depleted water vapor source. The isotopic data require some degree of drying with an isotopically depleted source for coupled conditions, and likely require a small amount of mixing even for strongly decoupled conditions. The data are consistent with mixing between purely local reservoirs of water vapor, subjected to a small amount of condensation and fractionation, and do not require the invocation of large-scale transport of water vapor aloft or additional cloud-formation effects aloft.

54 ENVIRONMENTAL SCIENCES↗

How similar are temperature and humidity profiles in the unstable boundary layer?

A comparison was made between temperature and humidity profiles obtained from radiosonde measurements under unstable conditions over hilly terrain in northeastern Kansas during the summer and fall of 1987. The summer profiles were found to be similar in a layer over the height range given by z = 65(plus or minus 45) to 124(plus or minus 54)m, where z is the height above the ground surface reference level. Evidence for this similarity was the closeness of the scaled ratio of the temperature and humidity gradients at any given time to the mean Bowen ratio determined at ground-based surface flux stations. Above the layer, in the outer region of the atmospheric boundary layer, it was found that the profiles of humidity and potential temperature are completely dissimilar in terms of the surface Bowen ratio. The profiles of the fall campaign when soil moisture content was low and vegetation senescent displayed a less pronounced similarity.

Sugita, Michiaki↗

Improving a Multilevel Turbulence Closure Model for a Shallow Lake in Comparison With Other 1-D Models

Lakes differ from lands in water availability, heat capacity, albedo, and roughness, which affect local surface-atmospheric interactions. This study modified a multilevel upper ocean model (UOM) for lake applications and evaluated its performance in Lake Taihu (China) with comprehensive measurements against three popular one-dimensional (1-D) lake models. These models were based on different concepts, including the self-similarity (FLake), the wind-driven eddy diffusion (LISSS), the k-ε turbulence closure (SIMSTRAT), and a simplified turbulence closure (UOM). The surface flux scheme in these models was unified to exclude the discrepancies in representing air-lake exchanges. All models in their default formulations presented obvious cold water temperature biases and largely underestimated the lake surface temperature (LST) diurnal range. For each model, these deficiencies were significantly reduced by incorporating new physics schemes or calibrated tunable parameters based on systematic sensitivity tests. The primary modifications for UOM included (1) a new scheme of decreased surface roughness lengths to better characterize the shallow lake, (2) a solar radiation penetration scheme with increased light extinction coefficient and surface absorption fraction to account for the high water turbidity, and (3) turbulent Prandtl number increased by a factor of 20 to reduce the turbulent vertical mixing. All other models were improved in these three aspects (roughness, extinction, and mixing) within their original formulations. Given these improvements, UOM showed superior performance to other models in capturing LST diurnal cycle and daily to seasonal variations, as well as summer-autumn vertical stratification changes. The new UOM is well suited for application in shallow lakes.

54 ENVIRONMENTAL SCIENCES↗

Dynamics of Downdrafts Around a Growing Convective Cloud: A Numerical Study

We examine the dynamics of cloud-edge downdrafts over the growth phase of isolated cumuli, combining Eulerian and Lagrangian analyses. As in previous studies, our results show that growing cumuli are surrounded by downdrafts linked to cloud-scale quasi-toroidal circulations at all times at middle and upper cloud levels consistent with the thermal chain description of convective clouds. These toroidal circulations are responsible for the most intense cloud-edge downdrafts in our simulations. In the upper cloud half, roughly 30%–50% of the upward mass flux is typically compensated within a radius of about twice the updraft radius in quasi-laminar simulations forced by a warm bubble in an initially quiescent flow. In a turbulent cloud forced by surface fluxes, this compensation fraction is around 10%–30% over the same region. In contrast to the buoyancy-centered view of subsiding shells, Eulerian and Lagrangian vertical momentum budget analyses show that the most intense cloud-edge downdrafts in the turbulent setup, and after spin-up of the toroidal circulation in the quasi-laminar experiments, are predominantly mechanically forced (i.e., driven by dynamic pressure accelerations). This is consistent throughout the entire growth phase of the cumulus clouds and across tests with varying assumptions, including drier and moister environments. Despite dynamic pressure perturbations being the main driver of toroidal downdrafts, the downdraft speed (relative to the corresponding updraft velocity) exceeds the prediction of the non-buoyant Hill's spherical vortex—a simple model frequently used for cloud-scale circulations—by more than 30%.

Pardo, Lianet Hernández [Goethe Univ., Frankfurt (↗

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Initial Results from High-Field-Side Transient CHI Start-Up on QUEST

Transient coaxial helicity injection (t-CHI) current start-up using a new design simple electrode configuration has been implemented on the QUEST. Discharges injected from the low field side (LFS) and from the high field side (HFS) were examined. Compared to the LFS injection case, the HFS injection has the advantages of providing access to a higher toroidal field and better controlling the location of the injector flux footprint location. Although the present PF coils on QUEST are not well positioned to form the injector flux on the HFS injector region and there has been a frequent occurrence of the spurious arcs, known as absorber arcs, HFS injection has shown flux evolution in a shape that is suitable for the formation of closed flux surfaces. The discharges were improved by installing an in-vessel-coil and adding a new cylindrical electrode to the existing CHI electrode. Overall, the results show that the new cylindrical electrode allowed the flux to evolve stably while allowing both the inner and the outer injector flux footprint to remain in the vicinity of the cylindrical electrode. This configuration which inherently generates a narrow injector flux footprint width resulted in discharges that strongly suggested the persistence of the CHI generated plasma after the injector current was reduced to zero. These studies have informed us of the need to improve the CHI gas injection system so that the absorber arcs could be better controlled in the HFS injection configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Io: Mountains and crustal extension

It is argued that there is good reason to conclude that mountains on Io, like those on Earth, are subject to growth and decay. The decay of mountains will be assisted by the ability of SO sub 2 to rot silicate rock and by explosive escape of sub-surface SO sub 2 from aquifers (Haemus Mons is seen to be covered by bright material, presumably fallout from a SO sub 2 rich plume which had been active on the mountain flanks). On the west side of the massif at 10 degrees S, 270 degrees W a rugged surface consists of long ridges running perpendicular to the downslope direction, suggesting tectonic denudation with crustal blocks sliding down the mountain flank. Tectonic denudation may be assisted, as in the case of the Bearpaw Mountains, Montana by overloading mountain flanks with volcanic products. The surfaces of some massifs exhibit a well developed, enigmatic corrugated terrain, consisting of complex ridge systems. Ridges may bifurcate, anastomose to form closed depressions and form concentric loops. Taken together, observations of morphology, heat flux, surface deposits and styles of volcanism may point to the existence of lithosphere domains with distinct compositions and tectonic regimes.

Heath, M. J.↗

Model Stellar Atmospheres and Real Stellar Atmospheres and Status of the ATLAS12 Opacity Sampling Program and of New Programs for Rosseland and for Distribution Function Opacity

I discuss errors in theory and in interpreting observations that are produced by the failure to consider resolution in space, time, and energy. I discuss convection in stellar model atmospheres and in stars. Large errors in abundances are possible such as the factor of ten error in the Li abundance for extreme Population II stars. Finally I discuss the variation of microturbulent velocity with depth, effective temperature, gravity, and abundance. These variations must be dealt with in computing models and grids and in any type of photometric calibration. I have also developed a new opacity-sampling version of my model atmosphere program called ATLAS12. It recognizes more than 1000 atomic and molecular species, each in up to 10 isotopic forms. It can treat all ions of the elements up through Zn and the first 5 ions of heavier elements up through Es. The elemental and isotopic abundances are treated as variables with depth. The fluxes predicted by ATLAS12 are not accurate in intermediate or narrow bandpass intervals because the sample size is too small. A special stripped version of the spectrum synthesis program SYNTHE is used to generate the surface flux for the converged model using the line data on CD-ROMs 1 and 15. ATLAS12 can be used to produce improved models for Am and Ap stars. It should be very useful for investigating diffusion effects in atmospheres. It can be used to model exciting stars for H II regions with abundances consistent with those of the H II region. These programs and line files will be distributed on CD-ROMs.

Kurucz, Robert L.↗

What causes the spread of model projections of ocean dynamic sea-level change in response to greenhouse gas forcing?

Abstract Sea levels of different atmosphere–ocean general circulation models (AOGCMs) respond to climate change forcing in different ways, representing a crucial uncertainty in climate change research. We isolate the role of the ocean dynamics in setting the spatial pattern of dynamic sea-level ( ζ ) change by forcing several AOGCMs with prescribed identical heat, momentum (wind) and freshwater flux perturbations. This method produces a ζ projection spread comparable in magnitude to the spread that results from greenhouse gas forcing, indicating that the differences in ocean model formulation are the cause, rather than diversity in surface flux change. The heat flux change drives most of the global pattern of ζ change, while the momentum and water flux changes cause locally confined features. North Atlantic heat uptake causes large temperature and salinity driven density changes, altering local ocean transport and ζ . The spread between AOGCMs here is caused largely by differences in their regional transport adjustment, which redistributes heat that was already in the ocean prior to perturbation. The geographic details of the ζ change in the North Atlantic are diverse across models, but the underlying dynamic change is similar. In contrast, the heat absorbed by the Southern Ocean does not strongly alter the vertically coherent circulation. The Arctic ζ change is dissimilar across models, owing to differences in passive heat uptake and circulation change. Only the Arctic is strongly affected by nonlinear interactions between the three air-sea flux changes, and these are model specific.

54 ENVIRONMENTAL SCIENCES↗

Diagnosing Tropical Divergent Circulations From Satellite-Derived Diabatic Forcing and Variational Integral Constraints

Large-scale divergent circulations are part of the atmospheric dynamic response to diabatic heating from condensation, radiative processes, and surface energy fluxes. Vertical motion and the associated divergent wind is thus intimately tied to the hydrologic cycle and the global beat balance. Despite its importance, the divergent circulation is too small in comparison to the rotational flow to measure directly with any accuracy. Vertical motions are recovered diagnostically from reanalyses and, as such, are subject to shortcomings in model physics, numerics, and data availability. While reanalysis estimates of tropical divergent circulations are much improved over those from the Global Weather Experiment, there are still substantial differences between products from operational centers. This is because these circulations are still forced largely by model physics and only secondarily by observations. In order to produce a refined estimate of tropical divergence and its interannual variability we have used a number of remotely-sensed data sets along with variational constraints to improve upon reanalysis estimates. Among these are: precipitation from SSM/I and GPCP, TOVS Path-A vertical cloud distributions; ISCCP radiative cooling rates; TOA radiative fluxes from ERBS, surface radiative fluxes from the SRB project, and surface latent and sensible flux estimates from SSM/I. The TOVS Path-A data constrain the divergent outflow in precipitating regions to have the same vertical structure as observed cloudiness. Using integral constraints for moisture, heat, and mass balance, we retrieve consistent divergent wind flows. We examine the ability of this type analysis to capture regional details of ENSO related perturbations to the divergent wind and associated tropical energy balance. Precipitation from satellite is found to be the major constraint in supplying this horizontal structure. We also consider the ability of this analysis to quantify integrated (land vs. ocean) fluctuations in the global monsoonal flow and energy balance. This depends, of course, on the error characteristics of the data constraints. Some estimate of these properties and their relative importance is provided.

Robertson, Franklin R.↗

Satellite Estimates of Surface Short-wave Fluxes: Issues of Implementation

Surface solar radiation reaching the Earth's surface is the primary forcing function of the land surface energy and water cycle. Therefore, there is a need for information on this parameter, preferably, at global scale. Satellite based estimates are now available at accuracies that meet the demands of many scientific objectives. Selection of an approach to estimate such fluxes requires consideration of trade-offs between the use of multi-spectral observations of cloud optical properties that are more difficult to implement at large scales, and methods that are simplified but easier to implement. In this study, an evaluation of such trade-offs will be performed. The University of Maryland Surface Radiation Model (UMD/SRB) has been used to reprocess five years of GOES-8 satellite observations over the United States to ensure updated calibration and improved cloud detection over snow. The UMD/SRB model was subsequently modified to allow input of information on aerosol and cloud optical depth with information from independent satellite sources. Specifically, the cloud properties from the Atmospheric Radiation Measurement (ARM) Satellite Data Analysis Program (Minnis et al., 1995) are used to drive the modified version of the model to estimate surface short-wave fluxes over the Southern Great Plain ARM sites for a twelve month period. The auxiliary data needed as model inputs such as aerosol optical depth, spectral surface albedo, water vapor and total column ozone amount were kept the same for both versions of the model. The estimated shortwave fluxes are evaluated against ground observations at the ARM Central Facility and four satellite ARM sites. During summer, the estimated fluxes based on cloud properties derived from the multi-spectral approach were in better agreement with ground measurements than those derived from the UMD/SRB model. However, in winter, the fluxes derived with the UMD/SRB model were in better agreement with ground observations than those estimated from cloud properties provided by the ARM Satellite Data Analysis Program. During the transition periods, the results were comparable.

Wang, H.↗

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↗