Search NASA⌕ Search

SEARCH · Search NASA

Results for “solar zenith angle”

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 145 records · Page 8

A Reevaluation of Airborne HO(x) Observations from NASA Field Campaigns

In-situ observations of tropospheric HO(x) (OH and HO2) obtained during four NASA airborne campaigns (SUCCESS, SONEX, PEM-Tropics B and TRACE-P) are reevaluated using the NASA Langley time-dependent photochemical box model. Special attention is given to previously diagnosed discrepancies between observed and predicted HO2 which increase with higher NO(x) levels and at high solar zenith angles. This analysis shows that much of the model discrepancy at high NO(x) during SUCCESS can be attributed to modeling observations at time-scales too long to capture the nonlinearity of HO(x) chemistry under highly variable conditions for NO(x). Discrepancies at high NO(x) during SONEX can be moderated to a large extent by complete use of all available precursor observations. Differences in kinetic rate coefficients and photolysis frequencies available for previous studies versus current recommendations also explain some of the disparity. Each of these causes is shown to exert greater influence with increasing NO(x) due to both the chemical nonlinearity between HO(x) and NO(x) and the increased sensitivity of HO(x) to changes in sources at high NO(x). In contrast, discrepancies at high solar zenith angles will persist until an adequate nighttime source of HO(x) can be identified. It is important to note that this analysis falls short of fully eliminating the issue of discrepancies between observed and predicted HO(x) for high NO(x) environments. These discrepancies are not resolved with the above causes in other data sets from ground-based field studies. Nevertheless, these results highlight important considerations in the application of box models to observationally based predictions of HO(x) radicals.

Olson, Jennifer↗

3D Orbit Visualization for Earth-Observing Missions

This software visualizes orbit paths for the Orbiting Carbon Observatory (OCO), but was designed to be general and applicable to any Earth-observing mission. The software uses the Google Earth user interface to provide a visual mechanism to explore spacecraft orbit paths, ground footprint locations, and local cloud cover conditions. In addition, a drill-down capability allows for users to point and click on a particular observation frame to pop up ancillary information such as data product filenames and directory paths, latitude, longitude, time stamp, column-average dry air mole fraction of carbon dioxide, and solar zenith angle. This software can be integrated with the ground data system for any Earth-observing mission to automatically generate daily orbit path data products in Google Earth KML format. These KML data products can be directly loaded into the Google Earth application for interactive 3D visualization of the orbit paths for each mission day. Each time the application runs, the daily orbit paths are encapsulated in a KML file for each mission day since the last time the application ran. Alternatively, the daily KML for a specified mission day may be generated. The application automatically extracts the spacecraft position and ground footprint geometry as a function of time from a daily Level 1B data product created and archived by the mission s ground data system software. In addition, ancillary data, such as the column-averaged dry air mole fraction of carbon dioxide and solar zenith angle, are automatically extracted from a Level 2 mission data product. Zoom, pan, and rotate capability are provided through the standard Google Earth interface. Cloud cover is indicated with an image layer from the MODIS (Moderate Resolution Imaging Spectroradiometer) aboard the Aqua satellite, which is automatically retrieved from JPL s OnEarth Web service.

Jacob, Joseph C.↗

Reflected fluxes for broken clouds over a Lambertian surface

Reflected fluxes are calculated for broken cloudiness (i.e., nonplane parallel) as a function of cloud cover, cloud optical depth, solar zenith angle and surface albedo. These calculations extend previous results for broken cloud reflected fluxes over a black surface. The present study demonstrates that not only radiances but also radiative fluxes over high albedo surfaces may be decreased by the presence of broken cloudiness. Conventional wisdom states that cloud radiances (brightnesses) are always greater than the background. While most cloud retrieval schemes are built around this assumption, it is incorrect for clouds over high albedo surfaces such as found in polar regions. However, the most startling and counterintuitive conclusion of this study is that nonabsorbing finite clouds over a highly reflecting surface will decrease the system albedo. As a result, surface absorption is increased, the result of multiple scattering between surface and cloud layer, controlled by cloud morphology and cloud optical thickness. A simple parameterization of the effects of cloud contamination upon retrieved albedo is given in terms of solar zenith angle, cloud optical depth, surface albedo, cloud cover, and plane-parallel cloud albedo. In this way, the effects of broken cloudiness are modeled in terms of easily computed plane-parallel values.

Welch, Ronald M.↗

Observations of "Hystereris" in Backscattered Ultraviolet Ozone Data

The backscattered ultraviolet (BUV) technique has been used for almost 3 decades to monitor global total ozone and the distribution of ozone in the stratosphere. Satellite BUV measurements in the 250-340 mn wavelength region are technically challenging because the Earth's radiance varies by approximately 4 orders of magnitude during a single scan. Further, the observed signal increases by over three orders of magnitude in about 2 minutes as the satellite emerges into daylight. The gain of the instrument's photomultiplier tube (PMT) detector is low when the spacecraft first emerges into the sunlit portion of the orbit relative to the gain observed after the PMT has experienced moderately high current levels. This "hysteresis" effect was first observed on the Nimbus-7 SBUV and TOMS instruments. The effect is difficult to characterize prelaunch because of the high signal levels and rapid variations required. We have recently observed and quantified the hysteresis effect for the NOAA-9 SBUV/2 instrument, which collected ozone data from February 1985 to February 1998. The instrument gain is observed to be up to 3% low at high solar zenith angles [Chi = 85-90 degrees] in the emergent hemisphere (i.e. Southern Hemisphere at launch). The gain error decreases as the SZA decreases and average PMT current increases, and is generally negligible for Chi < 65 degrees. The magnitude of the hysteresis effect varies with season, and exhibits long-term changes as the NOAA-9 sun-synchronous orbit drifts. In the latter portion of the record, when the spacecraft emerged from the dark in the Northern Hemisphere, hysteresis effects were then observed in the North. NOAA-9 total ozone errors due to the hysteresis effect are typically on the order of 2%, but can reach 5% in extreme cases. We have developed a quantitative correction for the hysteresis effect that incorporates both seasonal and long-term variations in magnitude. Results of similar analyses for the NOAA-11 and NOAA-14 SBUV/2 instruments will also be discussed. The characterization of the hysteresis effect in high solar zenith angle SBUV/2 ozone data represents a significant step towards reconciling polar ozone measurements from different satellite instruments.

DeLand, Matthew T.↗

Surface Net Solar Radiation Estimated from Satellite Measurements: Comparisons with Tower Observations

A parameterization that relates the reflected solar flux at the top of the atmosphere to the net solar flux at the surface in terms of only the column water vapor amount and the solar zenith angle was tested against surface observations. Net surface fluxes deduced from coincidental collocated satellite-measured radiances and from measurements from towers in Boulder during summer and near Saskatoon in winter have mean differences of about 2 W/sq m, regardless of whether the sky is clear or cloudy. Furthermore, comparisons between the net fluxes deduced from the parameterization and from surface measurements showed equally good agreement when the data were partitioned into morning and afternoon observations. This is in contrast to results from an empirical clear-sky algorithm that is unable to account adequately for the effects of clouds and that shows, at Boulder, a distinct morning to afternoon variation, which is presumably due to the predominance of different cloud types throughout the day. It is also demonstrated that the parameterization may be applied to irradiances at the top of the atmosphere that have been temporally averaged by using the temporally averaged column water vapor amount and the temporally averaged cosine of the solar zenith angle. The good agreement between the results of the parameterization and surface measurements suggests that the algorithm is a useful tool for a variety of climate studies.

Li, Zhanqing↗

Influence of the Surface and Cloud Nonuniformities in the Solar Energy Fluxes in the Arctic

Solar energy fluxes reaching the surface and absorbed by it are basic components of the energy balance of the Arctic. They depend mainly on the solar zenith angle, a state of the atmosphere, especially the cloudiness, and the surface albedo. However, they can also be modified by variabilities in the surface albedo and cloud optical thickness. The surface of the Arctic can be highly nonuniform. The surface of the Arctic Ocean, which covers the huge part of the Arctic can be view as a mosaic of sea water, sea ice, snow and, in the melting period, melting ponds. In our paper, results are presented of Monte Carlo simulations of the expected influence of nonuniform cloud structure and nonuniform surface albedo on radiative fluxes at the Arctic surface. In particular, the plane parallel biases in the surface absorptance and atmospheric transmittance are studied. The bias is defined as the difference between the real absorptance or transmittance (i.e. nonuniform conditions) averaged over a given area, and the uniform or plane parallel case with the same mean cloud optical thickness and the same mean surface albedo. The dependence of the biases is analysed with respect to the following: domain averaged values of the cloud optical thickness and surface albedo, scales of their spatial variabilities, correlation between cloud optical thickness and cloud albedo variabilities, cloud height, and the solar zenith angle. Ranges of means and standard deviations of the input parameters typical of Arctic conditions are obtained from the SHEBA experiment.

Rozwadowska, A.↗

The Use of the BSRN Data as A Benchmark for the POWER Hourly DHI and DNI and In Validating Derived Hourly GTI

The satellite-based CERES SYN1deg hourly data is the source data of the POWER GIS solar data that covers 2001 to near present. The SYN1deg(Ed4.1) hourly GHI agrees well with the BSRN data, but the hourly DHI and DirHI (Direct Horizontal Irradiance) are positively and negatively, respectively, biased with appreciable magnitudes. The hourly DNI, derived by dividing the DirHI by cos(SZA), or the cosine of the solar zenith angle, is therefore negatively biased. Based on the statistics of comparisons with the BSRN data, we performed bias corrections on the hourly DHI and DNI. The corrections were executed in the 3-D phase space of latitude, cos(SZA), and cloud fraction (CLFR). The isotropic model is then used to derive the hourly global tilted irradiance (GTI). For validation purpose, we applied the isotropic model to the BSRN data at the original 1-, 2-, 3- or 5-minute interval. The satellite-based hourly GTI shows good agreement with their BSRN counterpart. We also examined two monthly-mean-based methods that empirically derive monthly mean GTI and DNI from monthly mean GHI and from both monthly mean GHI and DHI. The monthly-mean-based results compare favorably with the hourly-mean-based results. The GEWEX SRB (V4-IP) provides POWER with daily mean GHI for the years before the CERES era, and the data were corrected using quantile mapping by referencing the CERES SYN1deg data. We used the Kolmogorov -Smirnov test (K-S test) and Cramer-von Mises test to examine how well the results agree with the BSRN data. We found that if we set the lower limit for the daily mean GHI to 30 W m-2, the data can pass the K-S test at 0.01 significance level and the Cramer-von Mises test at 0.001 significance level. If no lower limit is set on the daily means, the data fail both tests. The satellite-based CERES SYN1deg hourly data is the source data of the POWER GIS solar data that covers 2001 to near present. The SYN1deg(Ed4.1) hourly GHI agrees well with the BSRN data, but the hourly DHI and DirHI (Direct Horizontal Irradiance) are positively and negatively, respectively, biased with appreciable magnitudes. The hourly DNI, derived by dividing the DirHI by cos(SZA), or the cosine of the solar zenith angle, is therefore negatively biased. Based on the statistics of comparisons with the BSRN data, we performed bias corrections on the hourly DHI and DNI. The corrections were executed in the 3-D phase space of latitude, cos(SZA), and cloud fraction (CLFR). The isotropic model is then used to derive the hourly global tilted irradiance (GTI). For validation purpose, we applied the isotropic model to the BSRN data at the original 1-, 2-, 3- or 5-minute interval. The satellite-based hourly GTI shows good agreement with their BSRN counterpart. We also examined two monthly-mean-based methods that empirically derive monthly mean GTI and DNI from monthly mean GHI and from both monthly mean GHI and DHI. The monthly-mean-based results compare favorably with the hourly-mean-based results. The GEWEX SRB (V4-IP) provides POWER with daily mean GHI for the years before the CERES era, and the data were corrected using quantile mapping by referencing the CERES SYN1deg data. We used the Kolmogorov -Smirnov test (K-S test) and Cramer-von Mises test to examine how well the results agree with the BSRN data. We found that if we set the lower limit for the daily mean GHI to 30 W m-2, the data can pass the K-S test at 0.01 significance level and the Cramer-von Mises test at 0.001 significance level. If no lower limit is set on the daily means, the data fail both tests.

Taiping Zhang↗

The Re-Analysis of Ozone Profile Data from a 41-Year Series of SBUV Instruments

In this study we present the validation of ozone profiles from a number of Solar Back Scattered Ultra Violet (SBUV) and SBUV/2 instruments that were recently reprocessed using an updated (Version 8.6) algorithm. The SBUV dataset provides the longest available record of global ozone profiles, spanning a 41-year period from 1970 to 2011 (except a 5-year gap in the 1970s) and includes ozone profile records obtained from the Nimbus-4 BUV and Nimbus-7 SBUV instruments, and a series of SBUV(/2) instruments launched on NOAA operational satellites (NOAA 09, 11, 14, 16, 17, 18, 19). Although modifications in instrument design were made in the evolution from the BUV instrument to the modern SBUV(/2) model, the basic principles of the measurement technique and retrieval algorithm remain the same. The long term SBUV data record allows us to create a consistent, calibrated dataset of ozone profiles that can be used for climate studies and trend analyses. In particular, we focus on estimating the various sources of error in the SBUV profile ozone retrievals using independent observations and analysis of the algorithm itself. For the first time we include in the metadata a quantitative estimate of the smoothing error, defined as the error due to profile variability that the SBUV observing system cannot inherently measure. The magnitude of the smoothing error varies with altitude, latitude, season and solar zenith angle. Between 10 and 1 hPa the smoothing errors for the SBUV monthly zonal mean retrievals are of the order of 1 %, but start to increase above and below this layer. The largest smoothing errors, as large as 15-20%, were detected in in the troposphere. The SBUV averaging kernels, provided with the ozone profiles in version 8.6, help to eliminate the smoothing effect when comparing the SBUV profiles with high vertical resolution measurements, and make it convenient to use the SBUV ozone profiles for data assimilation and model validation purposes. The smoothing error can also be minimized by combining layers of data, and we will discuss recommendations for this approach as well. The SBUV ozone profiles have been intensively validated against satellite profile measurements obtained from the Microwave Limb Sounders (MLS) (on board the UARS and AURA satellites), Stratospheric Aerosol and Gas Experiment (SAGE) and Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). Also, we compare coincident and collocated SBUV ozone retrievals with observations made by ground-based instruments, such as microwave spectrometers, lidars, Umkehr instruments and balloon-borne ozonosondes. Finally, we compare the SBUV ozone profiles with output from the NASA GSFC GEOS-CCM model. In the stratosphere between 25 and 1 hPa the mean biases and standard deviations are within 5% for monthly mean ozone profiles. Above and below this layer the vertical resolution of the SBUV algorithm decreases and the effects of vertical smoothing should be taken into account. Though the SBUV algorithm has a coarser vertical resolution in the lower stratosphere and troposphere, it is capable of precisely estimating the integrated ozone column between the surface and 25 hPa. The time series of the tropospheric - lower stratospheric ozone column derived from SBUV agrees within 5% with the corresponding values observed by an ensemble of ozone sonde stations in North Hemisphere. Drift of the ozone time series obtained from each SBUV(/2) instrument relative to ground based and satellite measurements are evaluated and some features of individual SBUV(l2) instruments are discussed. In addition to evaluating individual instruments against independent observations, we also focus on the instrument to instrument consistency in the series. Overall, Version 8.6 ozone profiles obtained from two different SBUV(l2) instruments compare within a couple of percent during overlap periods and are consistently varying in time, with some exceptions. Some of the noted discrepancies might bssociated with ozone diurnal variations, since the difference in the local time of the observations for a pair of SBUV(l2) instruments could be several hours. Other issues include the potential short-term drift in measurements as the instrument orbit drifts, and measurements are obtained at high solar zenith angles (>85 ). Based on the results of the validation, a consistent, calibrated dataset of SBUV ozone profiles has been created based on internal calibration only.

Kramarova, Natalya↗

MISR Geometric Parameters (MIB2GEOP_V2)

The MISR instrument consists of nine pushbroom cameras which measure radiance in four spectral bands. Global coverage is achieved in nine days. The cameras are arranged with one camera pointing toward the nadir, four cameras pointing forward and four cameras pointing aftward. It takes 7 minutes for all nine cameras to view the same surface location. The view angles relative to the surface reference ellipsoid, are 0, 26.1, 45.6, 60.0, and 70.5 degrees. The spectral band shapes are nominally gaussian, centered at 443, 555, 670, and 865 nm. The MISR Level 1B2 Geometric Parameters data product is part of the georectified radiance product. It contains the geometric parameters which measure the sun and view angles at the reference ellipsoid. [Location=GLOBAL] [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=17.6 km; Longitude_Resolution=17.6 km; Horizontal_Resolution_Range=10 km - < 50 km or approximately .09 degree - < .5 degree; Temporal_Resolution=about 15 orbits/day; Temporal_Resolution_Range=Daily - < Weekly].

AM-1↗

Parameterization of albedo, thermal inertia, and surface roughness of desert scrub/sandy soil surface

Spectral albedo, A sub n, for the direct solar beam is defined as A sub n (r sub i,s, theta sub 0) = r sub i exp(-s tan theta sub 0)1-I(s) where I(s) is the integral over all reflection angles describing the interception by the absorbing plants of the flux reflected from the soil, r sub i soil reflectance, assumed Lambertian, S the projection on a vertical plane of plants per unit surface area, and theta sub 0 is the solar zenith angle. Hemispheric reflectance for the direct solar beam equals 1-I(s) times the reflectance to the zenith. The values of s of 0.1, 0.2, and 0.3 respectively quantify sparse, moderately dense, and very dense desert scrub. Thin plants are assumed to be of negligible thermal inertia, and thus directly yield the absorbed insolation to the atmosphere. Surface thermal inertia is therefore effectively reduced. The ratio of surface roughness height to plant height is parameterized for sparse, moderately dense, and very dense desert-scrub as a function of s based on data expressing the dependence of this ratio on plant silhouette.

Otterman, J.↗

Evaluation of SLAR and thematic mapper MSS data for forest cover mapping using computer-aided analysis techniques

The column normalizing technique was used to adjust the data for variations in the amplitude of the signal due to look angle effects with respect to solar zenith angle along the scan lines (i.e., across columns). Evaluation of the data set containing the geometric and radiometric adjustments, indicates that the data set should be satisfactory for further processing and analysis. Software was developed for degrading the spatial resolution of the aircraft data to produce a total of four data sets for further analysis. The quality of LANDSAT 2 CCT data for the test site is good for channels four, five, and six. Channel seven was not present on the tape. The data received were reformatted and analysis of the test site area was initiated.

Hoffer, R. M.↗

Off Median Phenomena and International Reference Ionosphere; COSPAR International Scientific Symposium, Trieste, Italy, Oct. 19-22, 1993

The main subject of the symposium was 'off-median phenomena'. The title denotes a range of problems that are rarely considered in 'pure science' studies of the ionosphere. The appearance of regular ionospheric variations is well known. Most of these depend on evident solar-geophysical influences like day and night, solar zenith angle, the seasons, geomagnetic control, solar activity, etc. Applicants and theoreticians as well used to work with monthly medians so that the in fact existing and quite important day-to-day variability is systematically overlooked. It is evident that a descriptive model like the International Reference Ionosphere (IRI) would be inadequate if this variability were denied. Interesting contributions from the symposium on 'off-median phenomena' and the IRI are presented.

Rawer, K.↗

Scene identification probabilities for evaluating earth radiation measurements

Alongtrack data from the Earth Radiation Budget Experiment (ERBE) satellite are used here to derive scene identification probabilities empirically. These error probabilities are shown to be functions of view zenith angle and scene type, but not of other parameters such as solar zenith angle and relative azimuth angle, except in the sunglint region. The standard deviations of the resulting errors in radiative fluxes are also presented.

Manalo, N.↗

Magnetization of the ionospheres of Venus and Mars - Results from radio occultation measurements

Remote sensing radio occultation measurements are used here to study magnetization of the ionospheres of Venus and Mars. For Venus, the measurements yield results on frequency of occurrence of magnetization during solar maximum that are similar to those obtained from Pioneer Venus in situ magnetic field measurements. During solar minimum, magnetization of the Venus ionosphere is more pervasive than at solar maximum. Magnetization extends to higher solar zenith angles and appears stronger than at solar maximum. These results confirm that during solar minimum the high solar wind dynamic pressure state is more prevalent at Venus because the ionospheric plasma pressure is weaker than at solar maximum. Comparison of a large number of electron density profiles of Mars with those of Venus shows an absence of the ledge and disturbed topside plasma observed in the Venus profiles. These results do not constitute evidence against magnetization of the ionosphere of Mars.

Woo, Richard↗

Extraction of soil and vegetation parameters from high resolution bi-directional reflectance spectra

High spectral resolution reflectance spectra were collected over a semi-desert grassland at both dry and wet season periods. Spectral reflectance measurements were made from several viewing angles at both low and high solar zenith angles. A mixture model was used to separate and extract green vegetation from dry/dead vegetation and soil. The extracted vegetation signal varied greatly with view and sun angle variations such that off-nadir viewing and illuminating angles resulted in the highest vegetation loadings. These variations were normalized with cosine functions for both sun and view angle. These results offer a methodology for standardizing multi-temporal and multi-angular satellite measurements of vegetation activity.

Huete, A. R.↗

A Chapman-Layers Ionspheric Model for Mars

A numerical model (CLIMM) is developed that adopts functions of two Chapman layers to compute Mars ionospheric electron densities at given local solar zenith angle and height. Electron density profiles derived from Mars Global Survey (MGS)-to-Earth radio occultation measurements collected during 1998 through 2005 are used to fit the model. The present model does not include variations with solar extreme ultraviolet (EUV) radiation cycles and seasons, and may have increased errors at lower latitudes. A more sophisticated model taking into account these variations is being developed and will be available in the future.

solar extreme ultraviolet (EUV)↗

Nimbus-7 total ozone algorithm

The total ozone algorithm used for the processing of Nimbus-4 BUV measurements is modified to account for off-nadir viewing of the Total Ozone Mapping Spectrometer onboard Nimbus-7 satellite. TOMS is a cross-track scanning instrument with a 50 x 50-km field of view, increasing to 150 x 150 km at the extreme off-nadir angle of 51 deg and six spectral bands in the 310-380-nm region. Consideration is given to the determination of the pressure at the reflecting surface and profile selection scheme at very high solar zenith angles. The algorithm is also applicable to the Solar Backscattered UV Experiment aboard the same satellite. The modifications provide for a better cloud height, snow correction, and improved retrievals near the terminator.

Bhartia, P. K.↗

Monthly Mean DNI and GTI Derived from Monthly Mean GHI and DHI Using Two Methods: Comparisons with the BSRN Data and the Results Derived from the CERES Hourly Data

Monthly mean Global Horizontal Irradiances (GHI) and Diffuse Horizontal Irradiances (DHI) are more widely available than monthly mean Direct Normal Irradiances (DNI) and Global Tilted Irradiances (GTI). Empirical methods have been developed to derive monthly mean DNI and GTI from monthly mean GHI or from GHI and DHI. In this paper, we evaluate two such methods. The first one was the Whitlock Method developed by Charles H. Whitlock (2005) for the NASA POWER project by means of regression of the BSRN data. The method expresses the monthly mean DHI-to-GHI ratio as polynomial functions of monthly mean clearness index, sunset hour angle and noon solar elevation angle on the monthly-average-day. The monthly mean DNI is calculated by dividing the monthly mean GHI-DHI difference, or DirHI, by the cosine of the solar zenith angle at the mid-time between sunrise and solar noon on the monthly-average-day. The second method is the LJCR Method developed by Liu and Jorden (1960) and Collares-Pereira and Rabl (1979), and this method empirically splits monthly mean GHI and DHI into hourly means on the monthly-average-day, and the resulting hourly mean GHI and DHI and their difference, DirHI, can then be used to compute the monthly mean DNI, GTI and the global solar tracker irradiance (GTrI). This method is also used by RETScreen. We recently produced a set of hourly DNI and DHI by bias-correcting the CERES hourly DNI and DHI, and computed hourly GTI and GTrI as well. The data span twenty plus years from March 2000 to near present on a 1 by 1 grid system. The monthly mean CERES GHI and the corrected DHI are used as inputs to the above two methods to compute monthly mean DNI, GTI and GTrI. Through comparisons with the BSRN data, it is found that the Whitlock Method, with slight modification, and the LJCR Method can produce results that are nearly as good as the results derived from the CERES hourly data.

Taiping Zhang↗