The core quantity at stake is aerosol optical depth, or AOD, a measure of how much sunlight or moonlight is attenuated by particles suspended in the atmosphere. AOD is a fundamental climate and air quality variable, and AERONET’s fleet of Cimel sun-sky radiometers has long been the gold standard for measuring it during daylight hours. But the atmosphere does not stop changing at sunset. Nighttime AOD can differ significantly from daytime values because the boundary layer shallows, relative humidity rises and drives hygroscopic particle growth, emission patterns shift, and chemical evolution proceeds without photochemistry. Capturing this 24-hour cycle matters because aerosols that cool the surface by day can actually warm it at night, and nocturnal radiative effects can even promote radiation fog that further processes aerosols.

Measuring AOD by moonlight, however, is a fundamentally harder problem than measuring it by sunlight. The Sun is an extraordinarily stable calibration source whose irradiance varies by only about 3.4 percent over a year, an effect that is almost entirely geometric and can be corrected with precision. The Moon is another story entirely. Lunar brightness at full moon is ten times or more that at quarter moon, and near full moon the irradiance can change by a percentage comparable to the annual solar variation in just two to three hours. This is partly due to the opposition effect, in which the porous lunar regolith produces constructive backscatter, and shadow hiding, a geometric effect that boosts lunar brightness by roughly 40 percent during the full moon window. Earth-moon distance adds another 12 percent swing within two weeks from apogee to perigee.

The network’s solution builds on the ROLO model, the Robotic Lunar Observatory model developed by the United States Geological Survey to predict lunar spectral irradiance. ROLO alone was never accurate enough for absolute calibration, and early attempts that relied on it produced unphysical results, including large negative AOD values during clean atmospheric conditions. The AERONET team instead developed an empirical correction. At the high-altitude calibration facilities at Mauna Loa in Hawaii and Izaña in the Canary Islands, reference instruments undergo both solar and lunar Langley calibrations, and the ratio of the resulting calibration coefficients is compared against the nominal sun-to-moon gain factor of 4096. Deviations from that ideal ratio, binned as a function of lunar phase angle, yield a correction factor for the ROLO model at each of six wavelengths from 440 to 1640 nanometers.

The updated assessment relies on a vastly larger set of calibration pairs than the original provisional lunar product released in 2016, with automated Langley pair counts ranging from 634 to 1537 depending on the channel. Each pair was retained only if the associated AOD fell below strict thresholds and the correlation coefficient of the fit was high enough. The corrections were then fitted with a nonparametric local regression, or LOESS, model to generate lookup tables at 0.8-degree resolution in lunar phase angle. Notably, the team chose to include data from earlier versions of the Cimel CE318-T Sun-Sky-Moon photometer despite their greater scatter, because restricting the dataset to newer instruments would have cut the number of Langley pairs by a factor of two to three, and the additional noise proved largely random.

To validate the corrected product, the researchers exploited the brief transition windows between day and night, comparing solar and lunar AOD measured by the same instrument within three-hour intervals under low and stable aerosol conditions. The analysis drew on 372 sites and 454 different instruments, yielding more than 280,000 comparison events. The result is striking: mean AOD differences at 440 nanometers were within plus or minus 0.005 for nearly all lunar phase angle bins, and at 675 nanometers and longer wavelengths the differences were within plus or minus 0.002. At low AOD, the mean solar-lunar difference of roughly 0.002 sits well inside the AERONET instrumental uncertainty of 0.01. The transition statistics were comparable to or smaller than lagged comparisons made entirely within day or entirely within night, confirming the absence of systematic calibration bias.

The assessment also probed what happens at high aerosol loading, where real atmospheric change over the transition interval dominates the comparison. Using both standard normal and bootstrap statistical methods, the team found a small but statistically significant tendency for daytime AOD to exceed nighttime AOD by about 0.013 to 0.016 in the 0.6 to 0.8 AOD bin. Rather than a calibration artifact, the researchers attribute this to genuine diurnal aerosol production, including daytime photochemistry and greater emissions from human activity and biomass burning. Scatter plots spanning all AOD values up to 1.0 showed nearly identical slopes of about 0.95 and concordance correlation coefficients of about 0.96 for both day-to-night and night-to-day transitions.

One subtlety of lunar photometry is that the Moon’s dimmer signal imposes a lower ceiling on measurable AOD. Under the network’s minimum signal threshold, solar photometry can theoretically capture total optical depth up to about 7.31 at an airmass of one, and observed solar maxima in the dataset, including extreme biomass burning events in Indonesia in 2015, came close to that limit. Lunar measurements under bright moon conditions reach roughly two to three times lower maximum AOD, and under minimum illumination the ceiling drops to about half the bright-moon value. The consequence is a disproportionate exclusion of high-AOD observations at shorter wavelengths: on average 3.2 percent of Level 2 lunar AOD at 440 nanometers is removed, three times the solar rate, with some dusty and smoky sites excluding more than 8 percent. At 675 and 870 nanometers the removal rates are 0.3 percent or less.

Cloud screening posed another question, since lunar measurements cannot perform the cross-scan curvature check that helps identify thin cirrus during solar operation. The team estimated the impact by comparing solar data processed with and without the curvature test across 259 sites. The differences proved negligible in most cases: aggregated multi-year monthly AOD differences ranged from only 0.0004 to 0.0014, and of nearly 15,000 individual months, just nine showed a mean AOD difference greater than 0.02. Even at Singapore, where the curvature check removes 5.1 percent more data points, the effect on monthly means was small. The researchers note that future versions may narrow remaining screening disparities, potentially with the help of machine learning algorithms.

The broader significance is considerable. Nighttime AOD from AERONET is already being used to validate satellite retrievals from the VIIRS Day-Night Band, which achieved correlations of roughly 0.8 over the United States, the Middle East, and the Indian Subcontinent. Lunar photometry is also the only option for monitoring Arctic haze during the polar night, which can last nearly four months at high-latitude sites such as Ny-Ålesund in Svalbard. With the new assessment demonstrating that the empirically corrected lunar AOD product agrees with solar measurements across diverse aerosol conditions worldwide, scientists now have a validated, continuous day-night record of atmospheric aerosols, opening the door to a truly complete picture of the 24-hour aerosol cycle and its consequences for air quality and climate.

Subject of Research: An assessment of lunar photometry in AERONET

Article Title: An assessment of lunar photometry in AERONET

Article References: Schafer, J., Eck, T. F., Slutsker, I., Gupta, P., Sorokin, M., Sinyuk, A., & Smirnov, A. (2026). An assessment of lunar photometry in AERONET. Atmospheric Measurement Techniques, 19(20), 6487-6502. https://doi.org/10.5194/amt-19-6487-2026

Image Credits: AI Generated

DOI: 10.5194/amt-19-6487-2026

Keywords: assessment, lunar, photometry, AERONET, scientific research, peer-reviewed research

What makes MISR extraordinary among Earth-observing instruments is its geometry. Flying on Terra, the flagship spacecraft of NASA’s Earth Observing System, the push-broom imager has been collecting radiance data since late February 2000, viewing the planet at nine angles ranging from straight down to 70.5 degrees fore and aft along the flight path, in four spectral bands centered at 446, 558, 672, and 866 nanometers. Because aerosol particles scatter and absorb light differently depending on the angle of observation, this multi-angle approach can constrain particle size, shape, and light-absorbing properties in ways that single-view imagers simply cannot. The instrument covers the globe roughly once a week, and even though Terra began constellation exit procedures in October 2022, MISR leaves behind a legacy of more than 24 years of multi-angle, multi-spectral observations of the entire planet.

The new MRPA is the culmination of a long-running parallel effort. NASA’s operational MISR Standard Algorithm, currently at version 23, was designed from the start to run automatically on the entire dataset, producing aerosol retrievals at 4.4 kilometer resolution. Alongside it, the research team built the MISR Research Algorithm, a more experimental code designed to push the sensitivity of the retrieval further, exploring innovations such as corrections for radiometric calibration errors like flat-field anomalies and ghosting in high-contrast scenes, and improved performance over shallow, turbid, and biologically productive waters. The catch was that the research algorithm demanded considerable hands-on user involvement and could only be run case by case, on limited subsets of data. The MRPA changes that equation by streamlining and automating the research approach so that any user can generate research-quality results across most of the MISR data record, spanning February 2000 through October 2022.

At the heart of the algorithm lies a clever dual strategy for dealing with the hardest problem in satellite aerosol retrieval: separating the signal of particles in the air from the light reflected by the ground below. The Retrieved-Surface Algorithm solves for surface reflectance self-consistently along with the aerosol properties, using a shape-similarity assumption in which the angular shape of surface reflectance is treated as independent of wavelength. In the updated MRPA, the shape-similarity coefficient is now calculated exclusively from the near-infrared band, because near-infrared light penetrates the atmosphere more reliably than the shorter blue and green wavelengths, giving better constraint on the surface. The Prescribed Surface Algorithm, by contrast, takes surface brightness from an external dataset, the Multi-Angle Implementation of Atmospheric Correction product derived from NASA’s MODIS instrument, and then solves for aerosol quantities directly using a bounded-variable least squares approach. This prescribed-surface route excels precisely where the retrieved-surface method struggles, namely in high-aerosol-optical-depth cases where the thick haze itself makes it hard to pin down what the surface is doing.

The team blends the two approaches into a Combined Surface Algorithm, weighting the contributions of each according to the Normalized Difference Vegetation Index of the underlying terrain. Bright, barren landscapes such as deserts are where small errors in prescribed surface reflectance can balloon into large aerosol errors, so the new weighting scheme, tuned against an extensive validation dataset, leans more heavily on the retrieved-surface results in low-vegetation scenes. The payoff is striking: in scenes where the prescribed-surface NDVI falls below 0.1, the new combined algorithm shows a root-mean-square error 45 percent lower than the previous version, with a 43 percent lower bias. The MRPA also introduces new cloud masking interpolated to the MISR swath from the MAIAC cloud mask over land, and over deep ocean it compares 12-micrometer brightness temperatures from MODIS-Terra against MERRA-2 reanalysis sea surface temperatures, flagging pixels that run 7.5 kelvin warmer as cloudy. Column ozone for absorption corrections now comes from MERRA-2 as well, and ocean surface roughness is prescribed from Cross-Calibrated Multiplatform wind speeds.

Validation against ground truth is where the new product truly shines. The researchers compared the retrievals against the AErosol RObotic NETwork, or AERONET, a global federation of calibrated sun photometers whose direct-sun aerosol optical depth measurements carry an uncertainty of about 0.01. Using more than 75,000 over-land coincidences between February 2000 and October 2022, the combined algorithm achieved an AOD root-mean-square error of 0.088, compared with 0.154 for the MISR V23 standard product, roughly half the error, along with a correlation coefficient of 0.938 versus 0.81. The bias of 0.008 is somewhat higher than the standard product’s −0.002, one of the compromises of automation, but the overall statistical improvement is substantial. Over water, where more than 29,000 coincidences were available, the algorithm performed even better, with an RMSE of 0.069 and a correlation of 0.923, and the expected error envelope works out to roughly 17 percent of the retrieved AOD plus 0.01 over land.

Particle-type information, the real prize of multi-angle imaging, is harder to validate quantitatively because the ground-based inversion retrievals of size, sphericity, and single-scattering albedo carry their own assumptions and limitations, particularly at low aerosol load. Still, the Ångström exponent, a robust measure of particle size derived from how aerosol optical depth varies with wavelength, shows two clear clusters in the validation data, around 0.25 and 1.5, corresponding to coarse-mode dust and fine-mode smoke or pollution respectively, and the agreement with AERONET improves as aerosol loading increases, with RMSE dropping from 0.43 to 0.25 at higher AOD thresholds. When the team examined regions and seasons dominated by a single aerosol source, the retrievals behaved exactly as physics demands. Smoke-dominated scenes, once the mid-visible AOD exceeded roughly 0.2 over water and 0.7 over land, showed a preponderance of small-to-medium, spherical, light-absorbing particles with single-scattering albedo between about 0.9 and 0.975. Dust-dominated regions near the Sahara in late spring and early summer yielded distinctly larger, weakly absorbing, non-spherical particles, with the coarse-mode signal especially clear over the high-contrast backdrop of dark ocean water.

The research algorithm behind the MRPA has already proven its worth in a remarkable range of scientific applications, from characterizing the properties and evolution of wildfire smoke plumes across North America, Alaska, and Siberia, to tracking volcanic eruption plumes including the catastrophic 2022 Hunga Tonga-Hunga Ha’apai event, to air quality assessments and studies of how dust particles age and change during trans-Atlantic transport. With the algorithm now automated, publicly released, and archived for general use, researchers anywhere can run it on any portion of the MISR archive that interests them, opening the door to systematic studies of how particle properties vary and evolve across the entire 22-year record rather than only in hand-picked case studies. The code’s availability marks a shift from a specialist tool requiring expert knowledge to a community resource.

The broader significance of this work lies in what it preserves and what it promises. As single-view instruments like MODIS have built their own long aerosol records, MISR’s multi-angle perspective offers a complementary and uniquely information-rich dataset, one that can distinguish, for example, a plume of absorbing smoke from a layer of scattering desert dust even when both produce similar optical depths. The compromises required to automate the research algorithm do come at a cost, with somewhat diminished particle-type sensitivity compared to the fully manual version, and the MRPA is currently configured to process data only through October 2022 to maintain continuity with its calibration corrections, though new Level 1 radiance products accounting for Terra’s orbit drift are planned. Even so, with about half the retrieval error of the operational standard product, four times finer pixel resolution, and demonstrated ability to separate smoke from dust by size, shape, and absorption, the MRPA hands atmospheric scientists a powerful new lens on the invisible particles that shape our climate and our air, and it ensures that the full scientific value of one of NASA’s most inventive instruments will be realized for years to come.

Subject of Research: Automated multi-angle satellite retrieval of global aerosol optical depth and particle properties from MISR

Article Title: The MISR Research Product Algorithm – producing global, self-consistent, pixel-level aerosol retrievals for the Multi-angle Imaging SpectroRadiometer

Article References: Anstett, M., Limbacher, J. A., & Kahn, R. A. (2026). The MISR Research Product Algorithm – producing global, self-consistent, pixel-level aerosol retrievals for the Multi-angle Imaging SpectroRadiometer. Atmospheric Measurement Techniques, 19(20), 6503-6516. https://doi.org/10.5194/amt-19-6503-2026

Image Credits: AI Generated

DOI: 10.5194/amt-19-6503-2026

Keywords: MISR, aerosol retrieval, Terra satellite, NASA, AERONET, aerosol optical depth, remote sensing, wildfire smoke, mineral dust, air quality, atmospheric measurement, multi-angle imaging

The measurement site, a small terminal depression called L’Bour near M’Hammid El Ghizlane, sits at the transition between the southern foothills of the Atlas Mountains and the open Saharan Desert. It is an ideal natural laboratory for dust emission studies. Mornings were typically calm, but as solar heating destabilized the near-surface air in the afternoons, winds strengthened and frequent local saltation events injected fresh dust into the lower atmosphere. On two days, 4 and 6 September 2019, the campaign captured intense cold-pool outflows known as haboobs, dramatic walls of dust that provided the researchers with abundant material for their instruments. This diurnal rhythm of quiet mornings and dusty afternoons allowed the team to build a rich dataset spanning a wide range of emission intensities.

The instrumentation combined a multi-angle integrating nephelometer, which measured light scattering at three wavelengths and eight angular positions, with a seven-wavelength dual-spot aethalometer that tracked absorption from the near-ultraviolet to the near-infrared. A Fine Dust Aerosol Spectrometer recorded particle size distributions in 63 logarithmically spaced bins covering optical diameters from roughly 0.2 to 19 micrometers. During the first part of the campaign, the optical instruments sampled through a PM2.5 inlet; later they switched to a PM10 inlet, a change that profoundly influenced what the team observed, because most of the mass of mineral dust resides in the coarse fraction. All measurements were averaged into 15-minute intervals to smooth out instrumental noise while preserving the turbulent structures that carry dust aloft.

The headline optical results are strikingly consistent with previous Saharan dust studies, which is itself reassuring for the field. During dust events, the single scattering albedo, the ratio of scattering to total extinction, was 0.95 for PM2.5 particles at 370 nanometers, rising to 0.97 at 660 nanometers, with slightly lower values of 0.94 and 0.96 for the PM10 fraction. These numbers tell us that Saharan dust is predominantly a scatterer of sunlight, but that its absorption is far from negligible, especially at short wavelengths. The absorption Ångström exponent, a measure of how strongly absorption varies with wavelength, reached values of up to 2.5, far above the values near 1 that characterize soot-like anthropogenic particles. This steep spectral signature is the fingerprint of iron oxide minerals such as hematite and goethite, which absorb much more strongly in the ultraviolet and blue than in the red and near-infrared.

Particle size emerged as a dominant control on nearly every optical property the team measured. The asymmetry parameter, which describes whether light is preferentially scattered forward or backward, was 0.7 at 520 nanometers for PM2.5 dust, and the backscatter fraction was correspondingly low at 0.11, both hallmarks of large particles that funnel scattered light strongly into the forward direction. The scattering Ångström exponent dropped to negative values when coarse dust dominated, a well-known indicator of coarse-mode particles. Mass absorption efficiency fell from about 0.30 to 0.15 square meters per gram at 370 nanometers as particle size increased, showing that larger grains absorb less light per unit mass. Meanwhile, the mass scattering efficiency shifted its peak toward longer wavelengths for larger particles, exactly as scattering theory predicts when the size parameter, the ratio of particle circumference to wavelength, grows beyond the regime of wavelength-independent scattering.

The most consequential part of the study, however, goes beyond describing the dust. The researchers used their combined measurements to retrieve the imaginary part of the complex refractive index, the quantity usually written as k that directly controls how much light a particle absorbs. Retrieving k from field data is an exercise in conditional reasoning: the answer depends on assumptions about particle density, the real part of the refractive index, and, critically, the shape of the particles. The team therefore ran their inversion repeatedly under a menu of prescribed shape assumptions, including perfect spheres, biaxial spheroids using aspect ratios measured on individual particles by electron microscopy during the campaign, the spheroid distribution commonly used in AERONET satellite retrievals, triaxial ellipsoids, and a set of six irregular particle geometries from the MOPSMAP modeling database.

The outcome was unambiguous. For central assumptions about density and the real refractive index, the retrieved k at 520 nanometers ranged from approximately 0.00104 for spherical particles to approximately 0.00145 for triaxial ellipsoids, a difference of about 40 percent. Non-spherical shapes systematically demanded larger values of k to reproduce the same observed absorption-to-scattering ratio, and the more elongated or irregular the assumed particle, the higher the retrieved absorption. The team found a nearly linear relationship between k and a shape descriptor, the ratio of volume-equivalent to cross-section-equivalent particle radius, with the effect weakening at longer wavelengths. Particle density assumptions also mattered, because density controls how the measured size distribution is converted between optical and aerodynamic diameters, whereas the assumed real refractive index had only a weak influence across its plausible range.

To put their in-situ results in context, the researchers compared them with two decades of dust-filtered retrievals from the AERONET network, using stations across western and eastern North Africa, the Sahel, the Arabian Peninsula, central Asia, and India, including Ouarzazate, the station closest to their field site. The comparison showed qualitative agreement with the Ouarzazate climatology, especially when the more realistic non-spherical shape assumptions were used, with median values converging closely at 550 and 675 nanometers. But the analysis also exposed pronounced regional variability in dust absorption, with retrieved k increasing from western toward eastern North Africa and further still across Asia, likely reflecting regional differences in iron oxide content and mineralogy. Differences in particle size distributions between surface measurements and column-integrated retrievals, along with dust ageing and residual aerosol contamination, complicate any strict comparison, so the authors frame the AERONET analysis as context rather than validation.

Why does all this matter? Mineral dust is one of the most abundant aerosol types in the atmosphere by mass, and its direct radiative effect is estimated at a central value of minus 0.40 watts per square meter, but with an uncertainty of plus or minus 0.30 that leaves open the possibility of a substantially different sign or magnitude. Dust absorption is the least constrained ingredient in these calculations, and most existing measurements come either from laboratory chambers using resuspended soils or from remote sensing of dust that has already been transported, aged, and mixed with other aerosols. Direct measurements of dust at the point of emission remain scarce. By showing that the shape assumptions embedded in retrieval algorithms and climate models can shift inferred absorption by tens of percent, the study provides a concrete explanation for part of the persistent spread among dust absorption estimates, and it suggests a practical path forward: past datasets built on spherical assumptions could potentially be reinterpreted using the relationship between particle asphericity and retrieved k.

The authors are careful about the limitations of their work. Their retrieval is explicitly conditional, restricted to the PM2.5 dust-event period, and their observing system simulation experiments revealed that size-dependent distortions in the measured size distributions, arising from different sampling inlets and possible particle losses, cannot be fully eliminated even by basing the cost function on the ratio of truncated scattering to extinction. Future campaigns, they argue, should prioritize co-located and size-consistent measurements, well-characterized inlets, and independent constraints on particle density and morphology. Still, the message of this desert field campaign is clear and broadly resonant: the light absorbed by a billowing Saharan dust plume depends not only on what the dust is made of, but on what shape we assume its grains to be, and getting that shape right may be essential for pinning down one of climate science’s most stubborn uncertainties.

Subject of Research: Optical properties and shape-dependent refractive index retrieval of freshly emitted Saharan mineral dust

Article Title: Optical properties and shape-dependent retrievals of the imaginary part of the refractive index of freshly emitted Saharan dust

Article References: Yus-Díez, J., Escribano, J., Pandolfi, M., Alastuey, A., González-Flórez, C., Gonzalez-Romero, A., Gonçalves Ageitos, M., Ivančič, M., Klose, M., Kandler, K., Obiso, V., Panta, A., Reche, C., Rigler, M., Querol, X., & Perez García-Pando, C. (2026). Optical properties and shape-dependent retrievals of the imaginary part of the refractive index of freshly emitted Saharan dust. Atmospheric Chemistry and Physics, 26(19), 14531-14557. https://doi.org/10.5194/acp-26-14531-2026

Image Credits: AI Generated

DOI: 10.5194/acp-26-14531-2026

Keywords: Saharan dust, mineral dust aerosol, refractive index, light absorption, single scattering albedo, particle shape, iron oxides, AERONET, Morocco, dust emission, radiative forcing, aerosol optics

AERONET, the Aerosol Robotic Network, has been measuring aerosol optical depth (AOD) with ground-based sun photometers since 1993. Its data underpin estimates of how aerosols influence Earth’s energy balance and serve as the gold standard for checking satellite retrievals. But the entire enterprise rests on a fragile assumption: that the measurements describe aerosols, not clouds. A thin veil of cirrus can masquerade as an aerosol layer, inflating AOD values and distorting inferred particle properties. Cloud screening has therefore always been a critical component of the network’s data processing, and any undetected cloud contamination quietly corrupts downstream science.

The current Version 3 screening algorithm became fully automatic by exploiting the solar aureole, the bright region of sky surrounding the sun. The angular shape of sky radiance in the aureole depends on particle size, so the algorithm parameterizes that shape using the curvature of the radiance curve and compares it against thresholds calibrated with lidar observations. Modern Cimel instruments support this with a curvature scan (CCS) that measures sky radiance from 3 to 7.5 degrees scattering angle in 0.3-degree steps after every AOD measurement. Yet the method is fundamentally threshold-based, and when cirrus parameters fall outside the pre-set limits, the clouds go undetected. The problem is most acute at high latitudes, where in situ aircraft studies have shown that cirrus ice crystals can shrink to sizes comparable to aerosol particles, with integral crystal radii as low as roughly 3 micrometers in the Arctic.

The failure mode is vividly illustrated at the OPAL site in Eureka, Canada, at about 80 degrees north latitude. There, the researchers found that Level 2 AOD time series, already passed through Version 3 screening, still displayed sharp upward spikes in AOD paired with sharp downward spikes in the Angstrom exponent, a size-sensitive parameter. Both signatures point to large ice crystals contaminating the record. The curvature parameters showed telltale spikes too, but most fell short of the operational thresholds, so the clouds slipped through. Similar failures appeared across northern Canada, Greenland, and the Arctic, where cold temperatures favor small ice crystals that scatter light like aerosols.

The team’s insight came from a different way of looking at the same data. In cloud-free conditions, the diurnal variation of sky radiance measured at the 3.3-degree scattering angle, the second smallest angle of the curvature scan, should be smooth, because it is driven almost entirely by the slow, predictable change in solar zenith angle. Radiative transfer calculations confirm this smooth behavior for an aerosol-laden atmosphere. The researchers showed that cloud-free days at sites as different as Ilorin in Nigeria, dominated by heavy dust, and the DEWA Research Center in the United Arab Emirates yielded radiance time series that fit polynomials with high correlation coefficients. On partially cloudy days, by contrast, the smooth aerosol signature is interrupted by irregular spikes, and those spikes become the fingerprint of cloud contamination.

From this observation the team built what they call the FD algorithm. It takes the full day of CCS(3.3) radiance measurements and computes first differences between values at neighboring time stamps. The standard deviation of these differences is compared against a threshold, ultimately set at 3.0, estimated from cloud-free days. If the standard deviation exceeds the threshold, the algorithm flags the largest value in each offending pair as cloud contaminated and removes the associated AOD measurement, then recalculates and iterates. The process stops when the standard deviation drops below the threshold or begins to increase, the latter signaling that the algorithm has started eating into genuine aerosol variability or that the day is completely overcast with no smooth aerosol signature to preserve. Crucially, no polynomial detrending is applied beforehand; the raw measurements speak for themselves.

Applying the algorithm across all AERONET sites for 2022 revealed both its power and its limits. At fine-mode-dominated sites such as Thule in Greenland and Hokkaido University in Japan, the Angstrom exponent rose systematically after screening, while the number of AOD measurements fell, with the two changes strongly anti-correlated, exactly what one expects when large cloud particles are removed. But at Capo Verde, where transported desert dust dominates, the algorithm stripped out large numbers of observations with almost no change in the Angstrom exponent, indicating it was deleting cloud-free data. The team’s remedy was an elegant second criterion: track how the Angstrom exponent changes with each algorithm iteration. Cloudy cases show a steep, highly correlated linear increase, while cloud-free cases show a nearly flat, uncorrelated drift. The slope of a linear regression between iteration number and Angstrom exponent, with a threshold of 0.01, now vetoes the FD results whenever the exponent barely moves, restoring wrongly eliminated data. Satellite imagery from Terra and Aqua MODIS confirmed the logic on test days at Capo Verde, where the slope criterion saved roughly 25 percent of mostly cloud-free observations.

Validation came from an entirely independent source: the Micro-Pulse Lidar Network, or MPLNET, whose vertically resolved backscatter profiles are particularly sensitive to cirrus. The team collocated AERONET and MPLNET data in space and time, averaged over one-hour periods, and drew on eight stations from GSFC in the United States to Kanpur in India and Barcelona in Spain. Comparing the frequency distributions of FD standard deviations for 90 percent clear and cloudy conditions fixed the threshold at 3.0, consistent with the case-study estimate. Performance was scored with the Matthews Correlation Coefficient, a stringent metric that rewards accuracy across true positives, true negatives, false positives, and false negatives simultaneously. On average, the FD algorithm outperformed Version 3 Level 1.5 screening by about 0.02 in MCC, a modest but consistent gain that held across solar zenith angle thresholds of 30, 60, and 90 degrees.

The detailed metrics tell a nuanced story. The FD algorithm identified a higher percentage of true positives and fewer false negatives than Version 3 alone, meaning it caught more cloud-contaminated measurements, at the cost of only slightly more false positives. Both algorithms maintained high true negative rates, confirming reliable clear-sky identification. Both showed declining MCC at larger solar zenith angles, largely because the lidar points straight up while the sun photometer looks toward the horizon, so the two instruments can sample different patches of a variable sky. Globally, the FD test removed an additional average of about 3.1 percent of data points from Version 3 Level 1.5, a small fraction that nonetheless matters enormously for field campaigns and high-cirrus regions where every clean observation counts.

The new technique is slated for inclusion in AERONET Version 4, expected in 2027, where it will run as an augmentation to the existing curvature algorithm rather than a replacement. Its beauty lies in parsimony: no new instrument, no new thresholds on angular shape, just the recognition that the atmosphere’s aerosol burden changes slowly and gracefully while clouds punctuate the record with spikes. For a network whose data feed climate assessments, satellite calibration, and aerosol retrievals on every continent, teaching an old scan a new trick may be one of the quiet but consequential upgrades in modern Earth observation.

Subject of Research: Cloud screening of aerosol optical depth measurements using time series smoothness of solar aureole sky radiances in the AERONET network

Article Title: Employing smoothness of the time series of sky radiances measured in the solar aureole for cloud screening

Article References: Sinyuk, A., Eck, T. F., Slutsker, I., Lewis, J., Grigorov, P., Smirnov, A., Schafer, J. S., Sorokin, M., Lind, E., & Gupta, P. (2026). Employing smoothness of the time series of sky radiances measured in the solar aureole for cloud screening. Atmospheric Measurement Techniques, 19(18), 6125-6144. https://doi.org/10.5194/amt-19-6125-2026

Image Credits: AI Generated

DOI: 10.5194/amt-19-6125-2026

Keywords: AERONET, aerosol optical depth, cloud screening, cirrus clouds, solar aureole, sky radiance, curvature scan, MPLNET, Angstrom exponent, Matthews correlation coefficient, sun photometer, atmospheric measurement

At the heart of the study is a measurement called aerosol optical depth, or AOD, which quantifies how much sunlight is removed from a beam passing through the atmosphere by particles suspended in the air. An AOD of 0.1 describes relatively clean skies, while values above 1.0 indicate conditions so hazy that the sun itself can appear dimmed and diffuse. The researchers worked with daily retrievals from the MODIS instrument aboard NASA’s Terra satellite, which has been scanning the Earth since early 2000. Raw satellite AOD products, however, carry well-documented biases, particularly over the bright surfaces and complex aerosol mixtures of West Africa. To correct these errors, the team built a machine learning framework known as Gradient Boosting Quantile Regression, trained against ground truth from fourteen AERONET sun photometer stations spread across nine West African countries.

The correction proved dramatic. At the Nigeria-Ilorin AERONET site, the machine learning adjustment shrank the raw mean bias of MODIS retrievals from minus 0.252 to just minus 0.015, a near-perfect agreement on average. More importantly, the fraction of satellite observations falling within the expected error envelope of ground measurements jumped from a dismal 23.0 percent to 95 percent. When the team tested the model’s ability to generalize to locations it had never seen, using a leave-one-site-out validation strategy, the expected error agreement remained a respectable 74.5 percent. This bias-corrected, boundary-clipped daily record gives Nigeria something it has never had before: a trustworthy, continuous, quarter-century chronicle of its atmospheric particle loading, suitable for everything from climate model evaluation to public health epidemiology.

So what does 25 years of data actually show? Perhaps surprisingly, the national annual trend is essentially flat. The linear slope across the full record is minus 0.0020 per year, with a p-value of 0.226, meaning the change cannot be distinguished from random year-to-year variability. But seasonal decomposition told a subtler story. The dry season, when Saharan dust and harmattan winds dominate, shows a suggestive decline of minus 0.0043 per year, equivalent to roughly a 16 percent reduction over the study period. This trend sits just at the edge of conventional statistical significance, with a p-value of 0.059 after accounting for autocorrelation through lag-1 pre-whitening, so the authors are careful to describe it as indicative rather than definitive. The wet season shows only a marginal, non-significant decrease of minus 0.0012 per year. In other words, Nigeria’s aerosol environment has been remarkably stable, with at most a gentle cleaning of its dustiest months.

The seasonal rhythm itself, however, is unmistakable. The monthly climatology is predominantly unimodal, peaking in February at a mean AOD of 0.865, when dry northeasterly harmattan winds carry desert dust deep into the country, and bottoming out in November at 0.370, as the monsoon transition washes the atmosphere clean. That February peak is more than double the November minimum, a swing that shapes everything from solar power output to aviation visibility to respiratory disease patterns across the region. Previous research has linked these dust episodes to meningitis outbreaks in the Sahel, and the harmattan season has been implicated in altered transmission dynamics of respiratory infections, making the timing and intensity of this annual maximum a matter of direct public health consequence.

To bring order to the daily chaos of aerosol variability, the researchers applied an unsupervised K-means clustering algorithm to the full record, letting the data sort itself into distinct atmospheric states. Three regimes emerged. Clean days, with a mean AOD of 0.399, account for 57.7 percent of the record. Transitional days, averaging 0.647, make up 31.6 percent. And then there are the Extreme days, averaging a staggering AOD of 1.163, which occupy 10.7 percent of all days in the record. These extreme events are overwhelmingly a dry season phenomenon: they occur on 18.1 percent of dry season days but only 3.2 percent of wet season days. That asymmetry confirms quantitatively what residents of Kano or Maiduguri know viscerally: the harmattan is when the air becomes a hazard.

One of the most practically useful contributions of the study is its characterization of how often extreme aerosol events recur. Using empirical return period analysis, the team calculated that days with AOD exceeding 1.0 arrive on average every 12.5 days, based on 711 exceedances across the record, with a 95 percent confidence interval of 11.6 to 13.4 days. Thresholds above 1.5 are crossed roughly every 61 days, with 145 exceedances recorded. And truly exceptional events, with AOD above 2.0, occur about once every 6.06 years, though with only four such events in 25 years the confidence interval is wide, spanning 3.03 to 24.25 years. These recurrence benchmarks give atmospheric hazard planners, health officials, and solar energy operators a climatological yardstick for what counts as a normal bad day versus a once-in-a-decade crisis.

The spatial structure of Nigeria’s aerosol burden is equally striking. A global spatial autocorrelation analysis yielded a Moran’s I of 0.949, with a p-value of 0.001 and a z-score of 163.9, confirming that aerosol loading is not randomly distributed but intensely clustered in space. Local Indicators of Spatial Association, corrected for the false discovery rate, flagged 40.1 percent of grid cells as statistically significant hotspots or coldspots. The pattern reveals a stable High-High hotspot belt across Southern and Southwestern Nigeria, where dense populations, urban pollution, and biomass burning combine, and persistent Low-Low spots in the Northeast. This geographic fingerprint has remained stable across the record, suggesting that the underlying drivers, from dust transport corridors to emission sources, are entrenched features of the regional environment rather than transient anomalies.

The study also probed the atmosphere’s memory, asking whether today’s aerosol conditions predict tomorrow’s. The answer is a strong yes, at least over short horizons. The lag-1 daily autocorrelation is 0.75, decaying to 0.21 by lag 30, indicating that aerosol states persist for days to weeks before the signal fades. An ARIMA (2,0,2) time series model successfully captured these synoptic dynamics, achieving a pseudo R-squared of 0.58. This short-range predictability has real-world value: it means that once a major dust outbreak begins, forecasters can anticipate with reasonable confidence that hazy conditions will persist for several days, giving hospitals, airports, and power grid operators time to prepare.

Taken together, these findings establish a rigorous long-term atmospheric baseline for a country whose air quality has historically been monitored far less intensively than its population warrants. The stability of the national trend, punctuated by a possible dry season improvement, offers a measure of reassurance, but the sheer frequency of extreme events, one in every nine days nationally and nearly one in five during the dry season, underscores the scale of the ongoing exposure. As West Africa’s combustion emissions continue to grow and climate change reshapes dust mobilization across the Sahara, records like this one will serve as the indispensable reference point against which future change is measured. For Nigeria’s policymakers, the message is clear: the dustiest skies are predictable, clustered, and recurring on a known clock, and that predictability is itself a tool for protecting public health.

Subject of Research: Long-term satellite analysis of aerosol optical depth dynamics, seasonal regimes, and extreme dust events over Nigeria from 2000 to 2024

Article Title: A 25-year analysis of aerosol optical depth dynamics over Nigeria (2000–2024): spatial-temporal dynamics, seasonal regimes, and extreme event characterization

Article References: Ekpe, E. E., Shi, C., Yin, S., Wang, H., Zhao, Y., Letu, H., & Raihan, M. Z. U. (2026). A 25-year analysis of aerosol optical depth dynamics over Nigeria (2000–2024): spatial-temporal dynamics, seasonal regimes, and extreme event characterization. Air Quality, Atmosphere & Health, 19(10), Article 224. https://doi.org/10.1007/s11869-026-02111-4

Image Credits: AI Generated

DOI: 10.1007/s11869-026-02111-4

Keywords: aerosol optical depth, Nigeria, MODIS Terra, machine learning bias correction, AERONET, Saharan dust, harmattan, air quality, spatial autocorrelation, extreme events, seasonal climatology, West Africa