- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Subscribe to alerts
- Peer review
- For authors
- For reviewers
- EGU publications

Journal cover
Journal topic
**Atmospheric Measurement Techniques**
An interactive open-access journal of the European Geosciences Union

Journal topic

- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Subscribe to alerts
- Peer review
- For authors
- For reviewers
- EGU publications

- Abstract
- Introduction
- Modification of instrument
- ROLO irradiance model
- Data
- Calibration of POM-02 using MLO data
- Results of comparison
- Summary and conclusion
- Data availability
- Appendix A
- Author contributions
- Competing interests
- Special issue statement
- Acknowledgements
- Financial support
- Review statement
- References
- Supplement

AMT | Articles | Volume 12, issue 12

Atmos. Meas. Tech., 12, 6465–6488, 2019

https://doi.org/10.5194/amt-12-6465-2019

© Author(s) 2019. This work is distributed under

the Creative Commons Attribution 4.0 License.

https://doi.org/10.5194/amt-12-6465-2019

© Author(s) 2019. This work is distributed under

the Creative Commons Attribution 4.0 License.

Special issue: SKYNET – the international network for aerosol, clouds,...

**Research article**
09 Dec 2019

**Research article** | 09 Dec 2019

Nocturnal aerosol optical depth measurements with modified sky radiometer POM-02 using the moon as a light source

^{1}Center for Global Environmental Research, National Institute for Environmental Studies, Tsukuba, Ibaraki 305-8506, Japan^{2}National Institute of Polar Research, Midoricho, Tachikawa, Tokyo 190-8518, Japan^{3}University of Yamanashi, Takeda, Kofu, Yamanashi 400-8510, Japan^{4}Meteorological Research Institute, Japan Meteorological Agency, Tsukuba, Ibaraki 305-0052, Japan^{5}Prede Co., Ltd., Kamidaira, Fussa, Tokyo 197-0012, Japan

^{1}Center for Global Environmental Research, National Institute for Environmental Studies, Tsukuba, Ibaraki 305-8506, Japan^{2}National Institute of Polar Research, Midoricho, Tachikawa, Tokyo 190-8518, Japan^{3}University of Yamanashi, Takeda, Kofu, Yamanashi 400-8510, Japan^{4}Meteorological Research Institute, Japan Meteorological Agency, Tsukuba, Ibaraki 305-0052, Japan^{5}Prede Co., Ltd., Kamidaira, Fussa, Tokyo 197-0012, Japan

**Correspondence**: Uchiyama Akihiro (uchiyama.akihiro@nies.go.jp)

**Correspondence**: Uchiyama Akihiro (uchiyama.akihiro@nies.go.jp)

Abstract

Back to toptop
The majority of aerosol data are obtained from daytime measurements, and there are few datasets available for studying nighttime aerosol characteristics. In order to estimate the aerosol optical depth (AOD) and the precipitable water vapor (PWV) during the nighttime using the moon as a light source, a sky radiometer (POM-02, Prede Ltd., Japan) was modified. The amplifier was adjusted so that POM-02 could measure lower levels of input irradiance. In order to track the moon based on the calculated values, a simplified formula was incorporated into the firmware. A new position sensor with a four-quadrant detector to adjust the tracking of the Sun and moon was also developed.

The calibration constant, which is the sensor output for the extraterrestrial solar and lunar irradiance at the mean Earth–Sun distance, was determined by using the Langley method. The measurements for the Langley calibration were conducted at the National Oceanic and Atmospheric Administration/Mauna Loa Observatory (NOAA/MLO) from 28 September 2017 to 7 November 2017. By assuming that the correct reflectance is proportional to the reflectance estimated by the Robotic Lunar Observatory (ROLO) irradiance model, the calibration constant for the lunar direct irradiance was successfully determined using the Langley method. The ratio of the calibration constant for the moon to that of the Sun was often greater than 1; the value of the ratio was 0.95 to 1.18 in the visible and near-infrared wavelength regions. This indicates that the ROLO model often underestimates the reflectance. In addition, this ratio depended on the phase angle. In this study, this ratio was approximated by a quadratic equation of the phase angle. By using this approximation, the reflectance of the moon can be calculated to within an accuracy of 1 % or less.

In order to validate the estimates of the AOD and PWV, continuous measurements with POM-02 were conducted at the Japan Meteorological Agency/Meteorological Research Institute (JMA/MRI) from January 2018 to May 2018, and the AOD and PWV were estimated. The results were compared with the AOD and PWV obtained by independent methods. The AOD was compared with that estimated by the National Institute for Environmental Studies (NIES) High Spectral Resolution Lidar measurements (wavelength: 532 nm), and the PWV was compared with the PWV obtained from a radiosonde and the Global Positioning System. In addition, the continuity of the AOD (PWV) before and after sunrise and sunset in Tsukuba was examined, and the AOD (PWV) of AERONET and that of POM-02 at MLO were compared. In the results, the daytime and nighttime AOD (PWV) measurements are shown to be statistically almost equivalent. The AODs (PWVs) during the daytime and nighttime for POM-02 are presumed to have the same degree of precision and accuracy within the measurement uncertainty.

How to cite

Back to top
top
How to cite.

Uchiyama, A., Shiobara, M., Kobayashi, H., Matsunaga, T., Yamazaki, A., Inei, K., Kawai, K., and Watanabe, Y.: Nocturnal aerosol optical depth measurements with modified sky radiometer POM-02 using the moon as a light source, Atmos. Meas. Tech., 12, 6465–6488, https://doi.org/10.5194/amt-12-6465-2019, 2019.

1 Introduction

Back to toptop
Atmospheric aerosols are an important constituent of the atmosphere. Aerosols change the radiation budget directly by absorbing and scattering solar radiation and indirectly through their role as cloud condensation nuclei (CCN), thereby increasing cloud reflectivity and lifetime (e.g., Ramanathan et al., 2001; Lohmann and Feichter, 2005). Aerosols also affect human health as one of the main components of air pollution (Dockery et al., 1993; WHO, 2006, 2013).

Atmospheric aerosols have a large variability in time and space. Therefore, measurement networks covering an extensive area on the ground and from space have been developed and established to determine the spatiotemporal distribution of aerosols. Well-known ground-based networks include AERONET (AErosol RObotic NETwork) (Holben et al., 1998), SKYNET (Takamura et al., 2004), and PFR-GAW (Precision Filter Radiometer – Global Atmosphere Watch) (Wehrli, 2005). These observation networks use passive radiometers which measure sunlight in the region from the ultraviolet to shortwave infrared wavelengths and the column-average effective aerosol characteristics such as aerosol optical depth (AOD) are retrieved.

Using lidar, which is an active remote sensing instrument, several networks have also been constructed: for example, the Micropulse Lidar Network (MPLNET) by NASA (National Aeronautics and Space Administration) (Welton et al., 2001; Levis et al., 2016), the European Aerosol Research Lidar Network (EARLINET) (Pappalardo et al., 2014) in Europe, the Asian Dust and aerosol lidar observation network (AD-Net) (Shimizu et al., 2016) in east Asia, and the Latin American Lidar Network (LALINET) (Guerrero-Rascado et al., 2016) in South America.

Several satellite programs provide aerosol optical depth data on a global scale: for example, the Moderate Resolution Imaging Spectroradiometer (MODIS) (Remer et al., 2005), Multiangle Imaging Spectroradiometer (MISR) (Kahn et al., 2005), Geostationary Operational Environmental Satellite (GOES) Aerosol/Smoke Product (GASP) (Prados et al., 2007), Sea-viewing Wide Field-of-view Sensor (SeaWiFS) (Wang et al., 2000), Advanced Himawari Imager (AHI) (Yoshida et al., 2018; Kikuchi et al., 2018), and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) (Winker et al., 2007).

With the exception of active sensor measurements such as lidar systems, to estimate aerosol characteristics, direct solar irradiance and scattered solar radiance measured with a passive sensor are required. Therefore, the majority of aerosol property data are obtained by daytime measurements, and there are few datasets of nighttime aerosol characteristics available.

To advance the understanding of the diurnal behavior of aerosols, and nocturnal mixing layer dynamics, nighttime continuous AOD measurements are necessary. In particular, in high-latitude regions during the winter polar night, aerosol properties cannot be measured using sunlight, and this results in gaps in the long-term aerosol data. Such nocturnal aerosol data would also contribute to the understanding of aerosol transport to polar regions, the influence of aerosol on cloud formation, and the cloud effect on the radiation budget.

Lidar instruments can be used to obtain aerosol data during the night. However, in many cases, lidar data retrieval requires some physical or mathematical constraints in inversion algorithms to allow the quantitative interpretation of the lidar backscatter signal (Fernald, 1984; Klett, 1985). In order to improve the accuracy of the analysis, constraining AOD is necessary.

In order to measure the optical depth of aerosol at night, research has been conducted using the moon and stars as light sources (Herber et al., 2002; Esposito et al., 1998, 2003; Pérez-Ramírez et al., 2008). Since the reflectivity of the moon changes depending on the observation angle, the determination of the calibration coefficient is an important obstacle to overcome (Herber et al., 2002). Instruments for observing stars are large, expensive, and complicated to use due to the low level of incoming energy from stars. Therefore, stellar measurements are limited in use, and no large-scale observation network has been established.

The moon is a bright light source at night, and the reflectance properties of
the moon's surface are virtually invariant (< 10^{−8} yr^{−1};
Kieffer, 1997). However, since the surface of the moon is not spatially
uniform and has non-Lambertian reflectance, the brightness of the moon as
seen by an observer on the Earth varies depending on the relationship
between the moon, the Sun, and the observer, that is, the phase and the
lunar libration. Therefore, it is difficult to use the moon as a light
source.

However, starting from the 2000s, the quality of reflectance data for the moon has improved. The empirical model known as ROLO (Robotic Lunar Observatory) was developed by the United States Geological Survey (USGS) (Kieffer and Stone, 2005). ROLO is a NASA-funded program aimed at using the moon for on-orbit calibration of Earth Observing System (EOS) satellite instruments. Furthermore, the Spectral Profiler (SP) aboard the Japanese Selenological and Engineering Explorer (SELENE, nicknamed Kaguya) measures lunar photometric properties in the region of visible, near-infrared, and shortwave infrared wavelengths (Yokota et al., 2011). These data made it possible to estimate the reflectance of the moon, and thus the moon can be used as a light source for aerosol optical depth estimation.

The Cimel Sun photometer used in AERONET has been modified for lunar observation and the aerosol optical depth at night can be estimated (Berkoff et al., 2011; Barreto et al., 2013, 2016, 2017). In addition, a lunar photometer – the Moon Precision Filter Radiometer, LunarPFR (Kouremeti et al., 2016) – has been developed by the Physical Meteorological Observatory in Davos (PMOD), which serves as the World Radiation Center (WRC), based on the Sun-PFR experience. Using these instruments and stellar photometers, a multi-instrument nocturnal intercomparison campaign was conducted to evaluate nighttime aerosol measurements and lunar irradiance models (Barreto et al., 2019).

In SKYNET, the POM-01 and POM-02 radiometers, manufactured by Prede Co. Ltd., Japan, are used. These radiometers are called “sky radiometers” and measure both the solar direct irradiance and sky radiances (Takamura et al., 2004). The sky radiometers (POM-01 and POM-02) can measure solar direct irradiance and sky radiances during the daytime, and the measured data are used for estimating aerosol characteristics during the daytime (Takamura et al., 2004). In this study, we will aim to measure the optical depth of aerosol using the moon as a light source by modifying POM-02.

In Sect. 2, we describe our modification of the instrument. In Sect. 3, the ROLO model is briefly explained. In Sect. 4, we briefly describe the data used in this study. In Sect. 5, the calibration method and corresponding results are described. In Sect. 6, we show the results of comparing the aerosol optical depth and precipitable water vapor obtained by continuous observation with those obtained by other independent instruments. We also show the results of comparing the aerosol optical depth and precipitable water vapor before and after sunrise and sunset using the continuous observation data. Furthermore, we show the results of a comparison between the AERONET and POM-02 data during the period of the MLO calibration measurement.

2 Modification of instrument

Back to toptop
In the modification of the POM-02 for solar observation, only the amplifier and the position sensor were changed. The other components, e.g., detectors, filters, and lenses, are not changed. Therefore, the magnitude of the solid view angle (field of view) for the new POM-02 is the same as in the non-modified POM-02. Measurements can still be obtained in the daytime using the modified POM-02.

The sky radiometer POM-02 is designed to measure the direct solar irradiance
and the scattered sky radiance with a single radiometer. An example of the
calibration constant, which is the sensor output for the extraterrestrial
solar irradiance at the mean Earth–Sun distance (one astronomical unit; AU)
at the reference temperature, is shown in Table 1. The calibration constant
is $\mathrm{1.8}\times {\mathrm{10}}^{-\mathrm{5}}$ to $\mathrm{3.4}\times {\mathrm{10}}^{-\mathrm{4}}$ A in the visible and
near-infrared region, and $\mathrm{7.9}\times {\mathrm{10}}^{-\mathrm{5}}$ to $\mathrm{1.3}\times {\mathrm{10}}^{-\mathrm{4}}$ A in
the short-wavelength infrared region. Figure 1 shows an example of
measurements of scattered radiances in the visible and near-infrared
wavelength regions. The output for the scattered radiance from the sky is
$\mathrm{1}\times {\mathrm{10}}^{-\mathrm{7}}$ to $\mathrm{1}\times {\mathrm{10}}^{-\mathrm{10}}$ A, and this value is
$\mathrm{1}\times {\mathrm{10}}^{-\mathrm{6}}$ smaller than the output for the direct solar
irradiance. The direct lunar irradiance is $\mathrm{1}\times {\mathrm{10}}^{-\mathrm{5}}$ as strong
as the direct solar irradiance during a full moon, and $\mathrm{1}\times {\mathrm{10}}^{-\mathrm{6}}$
during a half-moon (Berkoff et al., 2011). From Table 1, the calibration
constants at 340 and 380 nm are $\mathrm{1.8}\times {\mathrm{10}}^{-\mathrm{5}}$ and $\mathrm{1.9}\times {\mathrm{10}}^{-\mathrm{5}}$ (about $\mathrm{2}\times {\mathrm{10}}^{-\mathrm{5}})$, respectively. Therefore, the output
for the direct lunar irradiance during the half-moon is about $\mathrm{2}\times {\mathrm{10}}^{-\mathrm{5}}\times {\mathrm{10}}^{-\mathrm{6}}=\mathrm{2}\times {\mathrm{10}}^{-\mathrm{11}}$ in the 340 and 380 nm channels.
This is close to the detectable limits of the current POM-02.
Without modification, it is possible to measure the direct lunar irradiance
with the current POM-02 except for wavelengths between 340 and 380 nm, where
the sensitivity of the detector is low, and wavelengths of 1225, 1627, and
2200 nm with poor *S*∕*N*.

Table 2 shows the measurement ranges before and after modification of POM-02. POM-02 measures input energy in seven ranges according to the magnitude of the input energy, and the measured value is digitized with 15 bits. After modification, the measurement ranges are slightly expanded, and the measurement limit depends on the magnitude of the dark current and the magnitude of the noise. The sensor output takes into account the magnification of the amplifier, and the same amplifier was used for both the solar and lunar measurements.

The dark current of the detector in the visible and near-infrared region was about $\mathrm{5}\times {\mathrm{10}}^{-\mathrm{13}}$ A, and the root mean square (rms) of the random component of the noise was $\mathrm{4}\times {\mathrm{10}}^{-\mathrm{14}}$ A. In consideration of these values, the new POM-02 can use amplifiers for measurement ranges 1 to 7 and the minimum meaningful current is about $\mathrm{4}\times {\mathrm{10}}^{-\mathrm{13}}$ A (∼ rms ×10) in the visible and near-infrared region. This value is smaller than the output for the direct solar irradiance by a factor of $\mathrm{1}\times {\mathrm{10}}^{-\mathrm{8}}$ to $\mathrm{1}\times {\mathrm{10}}^{-\mathrm{9}}$.

The dark current of the detector in the shortwave infrared wavelength region was about $\mathrm{1.5}\times {\mathrm{10}}^{-\mathrm{8}}$ A, and the rms of the random component of the noise was $\mathrm{4}\times {\mathrm{10}}^{-\mathrm{11}}$ A. The measurement range is limited due to the large dark current. The new POM-02 can use amplifiers for measurement ranges 1 to 5 and the minimum meaningful current is about $\mathrm{4}\times {\mathrm{10}}^{-\mathrm{10}}$ A (∼ rms ×10). This value and the magnitude of the measured value of the direct lunar irradiance are comparable. Therefore, it is difficult to measure the direct lunar irradiance even with the new POM-02 in the shortwave infrared wavelength region.

The tracking of the Sun and the moon is based on the calculated position. The moon positions are calculated with the simplified formula in Nagasawa (1981). The necessary software is installed in the firmware of POM-02. Deviations may occur even if the instrument is pointed in the calculated direction due to errors in the moon position calculation, instrument installation errors, misalignment of the rotation axis, and so on. A position sensor is used to correct this deviation.

A position sensor with a four-quadrant detector is used to adjust the tracking of the Sun and the moon. In order to adjust the tracking of the moon, a position sensor incorporating a new electronic circuit to amplify the signal and new software to process the signal data were developed. The new position sensor can be used to track both the Sun and the moon.

When the input energy to the position sensor is small, it is difficult to
adjust the tracking with the position sensor. The magnitude of the input
energy to the position sensor varies depending on the lunar phase and the
aerosol optical depth. It was confirmed that the function of the moon
tracking adjustment works during the period of the full moon ± about
90^{∘} of the phase angle (half-moon).

Whether the position sensor can be used can be determined by a user-specified threshold value. That is, the position sensor can be deactivated when the input energy to the position sensor becomes less than the threshold value. For phase angles larger than the half-moon, the signal of the position sensor was small, and the position sensor was deactivated.

When the position sensor is not functioning, tracking is performed based on
the calculated values. When comparing the moon position calculated by this
simplified formula with that calculated using the NASA SPICE toolkit (Acton,
1996), the difference in the zenith angle is less than 0.01^{∘}, and the
difference in the azimuth angle is less than 0.04^{∘}. The center of the
field of view has a flat region of ±0.5^{∘}; the flat
region is ±0.25^{∘} in the solar disk scan. The apparent
diameters of the Sun and the moon are about 0.5^{∘}. Since the
calculation error of the moon position is less than 0.25^{∘}, if the
misalignment of the rotation axis is negligible and POM-02 is installed
correctly, it is possible to track the moon using only the calculated
positions. In fact, measurements could be made on the day of a full moon
±10 d (phase angle about 120^{∘}). Figure 2 shows an example
of the measurements on 14 October 2017 at NOAA/MLO. In this example, the phase
angle of the moon is 117.6 to 118.0^{∘}.

3 ROLO irradiance model

Back to toptop
In order to estimate the aerosol optical thickness using the moon as a light
source, measurement of the extraterrestrial irradiance of the moon is
necessary. In this study, a model known as the ROLO irradiance model
(Kieffer and Stone, 2005) was used. This model was developed at USGS and is based on an extensive database of radiance
images acquired by the ground-based ROLO over more than 8 years. ROLO is a
NASA-funded program designed to use the moon for on-orbit calibration of
Earth Observing System (EOS) satellite instruments. The empirical irradiance
model was developed for 32 wavelengths from 350 to 2450 nm and has the same
form for each wavelength. The average residual is less than 1 %. The
coefficients of the empirical formula were constrained and determined using
data with a phase angle between 1.55 and 97^{∘}. The empirically derived
analytic form based on the primary geometric variables is as follows:

$$\begin{array}{}\text{(1)}& \begin{array}{rl}\mathrm{ln}{A}_{k}& =\sum _{i=\mathrm{0}}^{\mathrm{3}}{a}_{ik}{g}^{i}+\sum _{j=\mathrm{1}}^{\mathrm{3}}{b}_{jk}{\mathrm{\Phi}}^{\mathrm{2}j-\mathrm{1}}+{c}_{\mathrm{1}}\mathit{\varphi}+{c}_{\mathrm{2}}\mathit{\theta}+{c}_{\mathrm{3}}\mathrm{\Phi}\mathit{\varphi}\\ & +{c}_{\mathrm{4}}\mathrm{\Phi}\mathit{\theta}+{d}_{\mathrm{1}k}{e}^{-g/{p}_{\mathrm{1}}}+{d}_{\mathrm{2}k}{e}^{-g/{p}_{\mathrm{2}}}\\ & +{d}_{\mathrm{3}k}\mathrm{cos}\left(\right(g-{p}_{\mathrm{3}})/{p}_{\mathrm{4}}),\end{array}\end{array}$$

where *A*_{k} is the disk-equivalent reflectance, *g* is the absolute phase
angle in radians, *θ* and *ϕ* are the selenographic latitude and
longitude of the observer in degrees, and Φ is the selenographic
longitude of the Sun in radians.

This formula must be used with caution. The equation in Kieffer and Stone (2005)
has well-known typographical errors. In Eq. (1), *θ* and *ϕ* in the original
expression by Kieffer and Stone (2005) are exchanged. In
addition, the units of the coefficients *p*_{1}, *p*_{2}, *p*_{3}, and
*p*_{4} are degrees. Therefore, in order to make the dimensions the same,
*g* in the exponent and the cosine terms must be converted into units of
degrees.

The astronomical parameter was calculated using our own software developed using the NASA SPICE toolkit; an observation geometry information system named SPICE is offered by NASA's Navigation and Ancillary Information Facility (NAIF) (Acton, 1996). SPICE is widely used in the NASA and international planetary exploration communities (for more information about SPICE, refer to the NAIF webpage at http://naif.jpl.nasa.gov, last access: 2 December 2019).

In this study, only the values of the reflectance are used, and it is assumed that there is an error in the ROLO reflectance and that the correct reflectance is proportional to the ROLO reflectance. This indicates that the relative variation in the ROLO reflectance is assumed to be correct. The reflectance values are not converted to irradiance values by assuming the extraterrestrial solar spectral irradiance. The wavelength of POM-02 used in this study does not necessarily match the wavelength of the ROLO model. Here, the reflectance at the wavelength of POM-02 was calculated by linearly interpolating from the reflectance of the ROLO model at two adjacent wavelengths. Information on the filters used in the ROLO measurement was not available. Here, the wavelength is represented by the center wavelength. In addition, the ROLO model does not have reflectance data for the 340 nm wavelength. The reflectance at the 340 nm wavelength was obtained by extrapolating linearly from the values at the two end wavelengths.

4 Data

Back to toptop
The aerosol optical thickness is estimated by measuring the attenuation of the direct solar or lunar irradiance. Therefore, in order to estimate the aerosol optical thickness, the output of the instrument for the input irradiance at the top of the atmosphere is necessary. The determination of this constant is referred to as calibration, and the output of the instrument for the extraterrestrial solar or lunar irradiance at the mean Earth–Sun distance (1 AU) at the reference temperature is called the calibration constant. In this study, the calibration constant was determined by the Langley method.

To calibrate the POM-02 by the Langley method, measurements were conducted
at NOAA/MLO during the period from 28 September 2017
to 7 November 2017; the full moon was on 4 October and 3 November 2017. The MLO
(19.5362^{∘} N, 155.5763^{∘} W) is located at an elevation of
3397.0 m a.m.s.l. on the northern slope of Mauna Loa, island of Hawaii,
Hawaii, USA. The atmospheric pressure is about 680 hPa. The MLO is one of
the most suitable places to obtain data for a Langley plot for the solar
direct irradiance measurement (Shaw, 1983). Though the air at MLO is highly
transparent, it is affected in the late morning and afternoon hours by
marine aerosol that reaches the observatory during the marine inversion
boundary layer breakdown under solar heating (Shaw, 1983; Perry et al., 1999).
Therefore, using data taken in the morning is recommended (Shaw, 1982; Dutton
et al., 1994; Holben et al., 1998).

However, during the nighttime, the upslope winds change to downslope winds, which bring low moisture and aerosol-poor air above the marine boundary layer down to the observatory. As a result, daytime orographic clouds at the observatory disappear and the atmosphere stratification becomes stable. These atmospheric conditions are suitable for obtaining data for the Langley plot from the lunar direct irradiance measurement.

During the calibration period, the data obtained for the moon over 18 nights for the visible and near-infrared region, and 13 nights for the short wavelength infrared region and water vapor channel (940 nm) were used to determine the calibration constants. The data obtained for the Sun over 22 d for the visible, near-infrared, and short wavelength infrared regions, and 24 d for the water vapor channel (940 nm) were used to determine the calibration constants.

The measurements for the estimation of the aerosol optical depth and
precipitable water vapor were performed at 1 min intervals at the Japan
Meteorological Agency/Meteorological Research Institute (JMA/MRI)
(36.05^{∘} N, 140.13^{∘} E) in Tsukuba, which is located about
50 km northeast of Tokyo. The comparison was made using data obtained during
the period from 1 January to 31 May 2018. During this period, the AOD and the
precipitable water vapor (PWV) were estimated assuming the calibration
constant was unchanged.

The optical depth estimated from POM-02 was compared with the value of the
National Institute for Environmental Studies (NIES) High Spectral Resolution
Lidar (HSRL; wavelength: 532 nm). The NIES/HSRL is one of the lidars operated
by the lidar measurement group of the NIES (Shimizu et al., 2016). The NIES
and MRI observation sites are located about 800 m apart. Since POM-02
was not measured at the 532 nm wavelength, the AOD at 532 nm was
interpolated from the values of 500 and 675 nm by assuming that AOD is
proportional to *λ*^{−α}, where *λ* is the wavelength.
Furthermore, since the AOD of NIES/HSRL is the 15 min average, the value
of POM-02 was also averaged over 15 min.

The PWV estimated from POM-02 was compared with that obtained from the vertical profile of a radiosonde and that obtained from the Global Positioning System (GPS) receiver. The radiosonde observation is operated from the JMA Aerological Observatory, which is adjacent to JMA/MRI. The GPS receiver is installed at JMA/MRI, and GPS data were processed by one of the JMA/MRI researchers (Shoji et al., 2013). The comparison of the PWV was performed using the 30 min average values.

5 Calibration of POM-02 using MLO data

Back to toptop
In this study, the calibration constant was determined by the Langley method (Uchiyama et al., 2018). Here, we do not consider the temperature dependence of the sensor output for the POM-02. Under these observation conditions in Tsukuba, the temperature dependence of the sensor output can be ignored except for the 340, 380, and 2200 nm channels (Uchiyama et al., 2018).

The sensor output when measuring the direct solar irradiance can be written as follows:

$$\begin{array}{}\text{(2)}& V\left({\mathit{\lambda}}_{\mathrm{0}}\right)={\displaystyle \frac{{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)}{{R}_{S}^{\mathrm{2}}}}\mathrm{exp}(-m(\mathit{\theta}\left)\mathit{\tau}\right({\mathit{\lambda}}_{\mathrm{0}}\left)\right){\stackrel{\mathrm{\u203e}}{T}}_{\mathrm{gas}}({\mathit{\lambda}}_{\mathrm{0}},\mathit{\theta}),\end{array}$$

where *V*(*λ*_{0}) is the sensor output in the *λ*_{0}
wavelength channel, *R*_{S} is the Earth–Sun distance in AU, *m*(*θ*)
is the total air mass, *τ*(*λ*) is the total optical depth, *θ* is the solar zenith angle, and ${\stackrel{\mathrm{\u203e}}{T}}_{\mathrm{gas}}({\mathit{\lambda}}_{\mathrm{0}},\mathit{\theta})$
is the channel average transmittance of the gas line absorption.
Furthermore, *V*_{S0}(*λ*_{0}) is the sensor output for the
extraterrestrial solar irradiance at 1 AU and is called the calibration
constant. *τ*(*λ*) consists of the optical thickness for molecular
scattering (Rayleigh scattering), aerosol, and the continuous absorption of
gas. In this study, it is assumed that air mass *m*(*θ*) is the same for
all components. The air mass *m*(*θ*) for molecular scattering is used
(Schmid and Wehrli, 1995; Holben et al., 1998).

In the case of no “gas absorption”, the following equation is used:

$$\begin{array}{}\text{(3)}& V\left({\mathit{\lambda}}_{\mathrm{0}}\right)={\displaystyle \frac{{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)}{{R}_{S}^{\mathrm{2}}}}\mathrm{exp}(-m(\mathit{\theta}\left)\mathit{\tau}\right({\mathit{\lambda}}_{\mathrm{0}}\left)\right).\end{array}$$

Taking the logarithm of the equation leads to

$$\begin{array}{}\text{(4)}& \begin{array}{rl}\mathrm{ln}\left(V\left({\mathit{\lambda}}_{\mathrm{0}}\right){R}_{S}^{\mathrm{2}}\right)& =\mathrm{ln}{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)-m\left(\mathit{\theta}\right)\mathit{\tau}\left({\mathit{\lambda}}_{\mathrm{0}}\right)\\ & ={C}_{\mathrm{1}}m\left(\mathit{\theta}\right)+{C}_{\mathrm{2}}.\end{array}\end{array}$$

The parameters on the left-hand side are known: *V* is the measurement
value, and *R*_{S} and *m*(*θ*) can be calculated from the solar zenith
angle. For example, *R*_{S} can be calculated with the simplified formula
in Nagasawa (1981), and *m*(*θ*) can be calculated as in Kasten and
Young (1989). In the case of POM-02, the sensor output is the current, and
the unit of the measurements of *V* is the ampere (A). *C*_{2}=ln *V*_{S0} is
determined from the ordinate intercept of a least-squares fit when one plots
the left-hand side of the above equation versus air mass *m*(*θ*).

For the water vapor absorption band at a wavelength of 940 nm, the Beer–Lambert–Bouguer law is not valid. Calibration methods for the 940 nm channel, which is in the water vapor absorption band, have been considered extensively in previous studies (Reagan et al., 1987a, b, 1995; Bruegge et al., 1992; Thome et al., 1992, 1994; Michalsky et al., 1995, 2001; Schmid et al., 1996, 2001; Shiobara et al., 1996; Halthore et al., 1997; Cachorro et al., 1998; Plana-Fattori et al., 1998, 2004; Ingold et al., 2000; Kiedron et al., 2001, 2003; Uchiyama et al., 2014, Campanelli et al., 2014).

In this study, the modified Langley method is used (Reagan et al., 1987a; Bruegge et al., 1992; Schmid and Wehrli, 1995). In the modified Langley method, the transmittance is approximated by an empirical formula. The water vapor transmittance is approximated as follows:

$$\begin{array}{}\text{(5)}& Tr\left({\mathrm{H}}_{\mathrm{2}}\mathrm{O}\right)=\mathrm{exp}(-a(m\left(\mathit{\theta}\right)\cdot \mathrm{pwv}{)}^{b}),\end{array}$$

where *a* and *b* are fitting coefficients, and pwv is PWV.

Coefficients *a* and *b* were determined by computing the transmittance for
several atmospheric models (Uchiyama et al., 2014).

The output of the 940 nm channel can be written as follows:

$$\begin{array}{}\text{(6)}& \begin{array}{rl}V\left({\mathit{\lambda}}_{\mathrm{0}}\right)& ={\displaystyle \frac{{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)}{{R}_{S}^{\mathrm{2}}}}\mathrm{exp}(-m(\mathit{\theta}\left)\mathit{\tau}\right({\mathit{\lambda}}_{\mathrm{0}}\left)\right)Tr\left({\mathrm{H}}_{\mathrm{2}}\mathrm{O}\right)\\ & ={\displaystyle \frac{{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)}{{R}_{S}^{\mathrm{2}}}}\mathrm{exp}(-m(\mathit{\theta}\left)\mathit{\tau}\right({\mathit{\lambda}}_{\mathrm{0}}\left)\right)\\ & \mathrm{exp}(-a(m\left(\mathit{\theta}\right)\cdot \mathrm{pwv}{)}^{b}).\end{array}\end{array}$$

Taking the logarithm of the equation leads to

$$\begin{array}{}\text{(7)}& \begin{array}{rl}\mathrm{ln}V{R}_{s}^{\mathrm{2}}+m\left(\mathit{\theta}\right)({\mathit{\tau}}_{\mathrm{aer}}+{\mathit{\tau}}_{R})& =\mathrm{ln}{V}_{S\mathrm{0}}-a\left(\mathrm{pwv}{)}^{b}m\right(\mathit{\theta}{)}^{b}\\ & ={C}_{\mathrm{1}}^{\prime}m(\mathit{\theta}{)}^{b}+{C}_{\mathrm{2}}^{\prime}.\end{array}\end{array}$$

In the same way as the normal Langley method, the parameters on the
left-hand side are known: *V* is the measurement value, and *R*_{S} and
*m*(*θ*) can be calculated from the solar zenith angle. *τ*_{R} is
also estimated from the surface pressure; for example, *τ*_{R} can be
calculated as in Asano et al. (1983). In addition, *τ*_{aer} is the
aerosol optical depth at the 940 nm wavelength, which is interpolated from
the aerosol optical depth from the values at the 870 and 1020 nm
wavelengths.

If pwv is constant, then the right-hand side of the equation is a linear
function of *m*(*θ*)^{b}. Therefore, the values on the left-hand side
can be fitted by a linear function of *m*(*θ*)^{b}, and the intersection
of the *y* axis and the fitted line is ln *V*_{S0}.

The sensor output when measuring the direct lunar irradiance can be written as follows:

$$\begin{array}{}\text{(8)}& \begin{array}{rl}V\left({\mathit{\lambda}}_{\mathrm{0}}\right)& ={\displaystyle \frac{{\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}}{\mathit{\pi}}}{\mathrm{\Omega}}_{\mathrm{M}}{\displaystyle \frac{{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)}{{R}_{S}^{\mathrm{2}}}}\cdot {\displaystyle \frac{\mathrm{1}}{{R}_{\mathrm{m}}^{\mathrm{2}}}}\\ & \mathrm{exp}(-m(\mathit{\theta}\left)\mathit{\tau}\right({\mathit{\lambda}}_{\mathrm{0}}\left)\right){\stackrel{\mathrm{\u203e}}{T}}_{\mathrm{gas}}({\mathit{\lambda}}_{\mathrm{0}},\mathit{\theta}),\end{array}\end{array}$$

where Ω_{M} is the solid angle of the moon, *R*_{S} is the
distance between the moon and the Sun in AU, and *R*_{m} is the distance
between the moon and the observer normalized by 384 400 km (the mean radius
of the moon's orbit around the Earth). ${\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}$ is the smoothed
ROLO reflectance adjusted to the laboratory reflectance spectra of the
Apollo 16 samples. ${\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}$ is calculated using the lunar
reflectance *A*_{ROLO} with the ROLO irradiance model by the method shown
in Kieffer and Stone (2005) (see Appendix A).

Let ${\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}={F}_{C}{A}_{\mathrm{ROLO}}$, where *F*_{C} is a constant for
smoothing (see Appendix A). Using this equation, Eq. (8) becomes

$$\begin{array}{}\text{(9)}& \begin{array}{rl}V\left({\mathit{\lambda}}_{\mathrm{0}}\right)& ={\displaystyle \frac{{F}_{C}{A}_{\mathrm{ROLO}}}{\mathit{\pi}}}{\mathrm{\Omega}}_{\mathrm{M}}{\displaystyle \frac{{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)}{{R}_{S}^{\mathrm{2}}}}\cdot {\displaystyle \frac{\mathrm{1}}{{R}_{\mathrm{m}}^{\mathrm{2}}}}\\ & \mathrm{exp}(-m(\mathit{\theta}\left)\mathit{\tau}\right({\mathit{\lambda}}_{\mathrm{0}}\left)\right){\stackrel{\mathrm{\u203e}}{T}}_{\mathrm{gas}}({\mathit{\lambda}}_{\mathrm{0}},\mathit{\theta}).\end{array}\end{array}$$

It is known that the aerosol optical depth retrieved using the ROLO
reflectance contains an error, which is dependent on the phase angle
(Barreto et al., 2016, 2017, 2019; Juryšek and Prouza, 2017). We assume
that there is an error in the ROLO reflectance and that the correct lunar
reflectance is proportional to the ROLO reflectance. This indicates that the
relative variation in the ROLO model reflectance is assumed to be correct.
Let the proportional constant be denoted *C*^{′}, and *A*_{ROLO} in Eq. (9)
be replaced with ${C}^{\prime}\cdot {A}_{\mathrm{ROLO}}$. Equation (9) then becomes

$$\begin{array}{}\text{(10)}& \begin{array}{rl}V\left({\mathit{\lambda}}_{\mathrm{0}}\right)& ={\displaystyle \frac{{F}_{C}{C}^{\prime}{A}_{\mathrm{ROLO}}}{\mathit{\pi}}}{\mathrm{\Omega}}_{\mathrm{M}}{\displaystyle \frac{{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)}{{R}_{S}^{\mathrm{2}}}}\cdot {\displaystyle \frac{\mathrm{1}}{{R}_{\mathrm{m}}^{\mathrm{2}}}}\\ & \mathrm{exp}(-m(\mathit{\theta}\left)\mathit{\tau}\right({\mathit{\lambda}}_{\mathrm{0}}\left)\right){\stackrel{\mathrm{\u203e}}{T}}_{\mathrm{gas}}({\mathit{\lambda}}_{\mathrm{0}},\mathit{\theta})\\ & ={\displaystyle \frac{C{A}_{\mathrm{ROLO}}}{\mathit{\pi}}}{\mathrm{\Omega}}_{\mathrm{M}}{\displaystyle \frac{{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)}{{R}_{S}^{\mathrm{2}}}}\cdot {\displaystyle \frac{\mathrm{1}}{{R}_{\mathrm{m}}^{\mathrm{2}}}}\\ & \mathrm{exp}(-m(\mathit{\theta}\left)\mathit{\tau}\right({\mathit{\lambda}}_{\mathrm{0}}\left)\right){\stackrel{\mathrm{\u203e}}{T}}_{\mathrm{gas}}({\mathit{\lambda}}_{\mathrm{0}},\mathit{\theta}),\end{array}\end{array}$$

where *F*_{C}*C*^{′} is substituted with *C*.

In the case of no “gas absorption”, taking the logarithm of the equation leads to

$$\begin{array}{}\text{(11)}& \begin{array}{rl}\mathrm{ln}\left({\displaystyle \frac{\mathit{\pi}V\left({\mathit{\lambda}}_{\mathrm{0}}\right)}{{A}_{\mathrm{ROLO}}{\mathrm{\Omega}}_{\mathrm{M}}}}{R}_{S}^{\mathrm{2}}{R}_{\mathrm{m}}^{\mathrm{2}}\right)& =\mathrm{ln}C{V}_{S\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)-m\left(\mathit{\theta}\right)\mathit{\tau}\left({\mathit{\lambda}}_{\mathrm{0}}\right)\\ & =\mathrm{ln}{V}_{m\mathrm{0}}\left({\mathit{\lambda}}_{\mathrm{0}}\right)-m\left(\mathit{\theta}\right)\mathit{\tau}\left({\mathit{\lambda}}_{\mathrm{0}}\right)\\ & ={C}_{\mathrm{1}}^{\prime \prime}m\left(\mathit{\theta}\right)+{C}_{\mathrm{2}}^{\prime \prime},\end{array}\end{array}$$

where *V*_{m0}(*λ*_{0})=*C**V*_{S0}(*λ*_{0}). ${C}_{\mathrm{2}}^{\prime \prime}=\mathrm{ln}{V}_{m\mathrm{0}}$ is determined from the ordinate intercept of a
least-squares fit when one plots the left-hand side of the above equation
versus air mass *m*(*θ*).

*V*_{S0} can be determined by applying the Langley method to data taken
during the daytime. If *V*_{S0} is determined, the coefficient *C* can be
determined by taking the ratio of *V*_{m0} and *V*_{S0}. If the
coefficient *C* is 1, the reflectance of the ROLO model will be correct. If
the coefficient *C* is greater than 1 (less than 1), the reflectance in the
ROLO model is underestimated (overestimated).

Examples of Langley plots in the visible and near-infrared wavelengths are shown in Fig. 3. In these examples, the regression lines can be well determined for any wavelength. ${C}_{\mathrm{2}}^{\prime \prime}=\mathrm{ln}{V}_{m\mathrm{0}}$ is determined from the ordinate intercept of the regression line (see Eq. 11). At the 340 nm wavelength, the regression line tends to deviate from the measured values in the region of air masses larger than 6. It is presumed that the detector output at the 340 nm wavelength is small and hence may be nonlinear. The output at the time of observation was about $\mathrm{1}\times {\mathrm{10}}^{-\mathrm{12}}$ A. When using output values less than this, the user needs to treat their results with caution. At the 940 nm wavelength, the modified Langley method was applied. In this example, the regression line provides a good fit.

In Fig. 4, examples of the Langley plot in the shortwave infrared region (1225, 1627, and 2200 nm) are shown. The detector output of these channels ranges from $\mathrm{2}\times {\mathrm{10}}^{-\mathrm{10}}$ to $\mathrm{5}\times {\mathrm{10}}^{-\mathrm{10}}$ A, and the root mean square error of the random noise is $\mathrm{4}\times {\mathrm{10}}^{-\mathrm{11}}$ A. The ratio of noise to detector output is large, and it is difficult to use these channels for estimating the aerosol optical depth.

In Fig. 5, the relationship between the coefficient $C(={V}_{m\mathrm{0}}/{V}_{S\mathrm{0}})$
and the phase angle in the visible and near-infrared wavelength region (from
340 to 1020 nm) is shown. As shown in the previous section, the corrected
lunar reflectance is assumed to be proportional to the ROLO reflectance, and
the proportional coefficient *C* is the ratio of the calibration constant
for the moon and the Sun. That is, the coefficient *C* indicates the error
of the ROLO reflectance, and thus more accurate reflectance can be obtained
by multiplying the ROLO reflectance by the coefficient *C*. As can be seen
from this figure, the coefficient *C* is often greater than 1 and depends on
the phase angle. At most wavelengths, the coefficient *C* is small when the
absolute value of the phase angle is small (near the full moon) and
increases as the absolute value of the phase angle increases. The range of
*C* is 0.95 to 1.18. The absorption band of water vapor is at the 940 nm
wavelength. Water vapor in the atmosphere tends to fluctuate. Therefore, it
is difficult to make accurate Langley plots, and the accuracy of both
*V*_{S0} and *V*_{m0} is poor. Therefore, no clear relationship between
*C* and the phase angle is found, but the coefficient *C* is about 1.16. The
fact that *C* is larger than 1 means that the reflectance of the ROLO
irradiance model is underestimated.

In Fig. 6, the relationship between the coefficient $C(={V}_{m\mathrm{0}}/{V}_{S\mathrm{0}})$
and the phase angle in the shortwave infrared wavelength region (1225, 1627, and
2200 nm) is shown. In these channels, the error for *C* is large, but the
coefficient *C* depends on the phase angle as in the visible and
near-infrared wavelength region; *C* is small when the phase angle is near
zero and increases as the absolute value of the phase angle increases.

In this study, the phase angle dependence of the coefficient *C* is
approximated by a quadratic equation of the absolute value of the phase
angle:

$$\begin{array}{}\text{(12)}& C={A}_{c}\cdot {g}^{\mathrm{2}}+{B}_{c},\end{array}$$

where *g* is the phase angle.

That is,

$$\begin{array}{}\text{(13)}& {V}_{m\mathrm{0}}={V}_{S\mathrm{0}}\cdot ({A}_{c}\cdot {g}^{\mathrm{2}}+{B}_{c}).\end{array}$$

The coefficients *A*_{c} and *B*_{c} are shown in Table 3. The regression
line was plotted in Figs. 5 and 6. By using this approximation, the
reflectance of the ROLO model can be estimated to within 1 % in most
channels. By using this approximation, the data processing to estimate the
aerosol optical depth from the measured value becomes straightforward. The
coefficients, *F*_{C}, for smoothing the ROLO reflectance are also shown in
Table 3. The coefficients ${A}_{c}^{\prime}$ and ${B}_{c}^{\prime}$ of the
regression equation when using the smoothed ROLO reflectance are also given.

The size of the error in the reflectance in the ROLO irradiance model is
dependent on the phase angle. The ROLO reflectance was obtained by dividing
the lunar irradiance measured by Kieffer and Stone (2005) by the solar
spectral irradiance of the 1985 Wehrli Standard Extraterrestrial Solar
Irradiance Spectrum (Wehrli, 1985: Neckel and Labs, 1981). The solar spectral
irradiances are dependent on the solar spectral models. Therefore, the ROLO
reflectance includes an error due to the error in the solar spectral
irradiance of 1985 Wehrli. Instrument performance, data processing, and so
on are also sources of error. In this study, *C* is approximated as a
symmetric quadratic equation of the phase angle, but the phase angle
dependence of *C* is asymmetric (see Figs. 5 and 6). The applicable range of
the ROLO reflectance model is a phase angle of about 95^{∘} or less. In
order to improve the accuracy of the ROLO reflectance model and expand its
application range, it is necessary to further accumulate the reflectance
data of the moon.

6 Results of comparison

Back to toptop
In order to validate the estimations of AOD and PWV, we compared them with the AOD and PWV obtained by independent methods. We investigated whether there is a difference between daytime and nighttime measurements, and compared the measurements for the daytime and nighttime with measurement data which were recorded independently of POM-02 and have the same accuracy and precision in the daytime and nighttime.

Furthermore, the continuity of the AOD and PWV before and after sunrise and sunset was investigated, and the AOD and PWV of AERONET and POM-02 at MLO were also compared.

The AOD estimated from POM-02 was compared with the value of the NIES/HSRL (wavelength: 532 nm).

Figure 7a and b show the scatter plot of the aerosol optical depth during the daytime and nighttime, respectively. In Fig. 7c, the scatter plot during the nighttime is shown together with that during the daytime. Table 4 shows the results of the comparison between NIES/HSRL and POM-02 AOD: the statistics of the difference between the two AODs, the coefficients of the linear regression equation of NIES/HSRL and POM-02 AOD (${\mathit{\tau}}_{\mathrm{HSRL}}={C}_{\mathrm{1}}\cdot {\mathit{\tau}}_{\mathrm{POM}\mathrm{02}}+{C}_{\mathrm{2}})$, the root mean square error (RMSE) of the residual, the 95 % confidence interval of the coefficients, and the number of observations.

RMSE: root mean square error. CR: correlation coefficient. *C*_{1}
and *C*_{2}: coefficients of regression line
(${\mathit{\tau}}_{\mathrm{HSRL}}={C}_{\mathrm{1}}\cdot {\mathit{\tau}}_{POM-\mathrm{02}}+{C}_{\mathrm{2}}$). CI of *C*_{1} (95 %):
95 % confidential interval of *C*_{1}. CI of *C*_{2} (95 %): 95 %
confidential interval of *C*_{2}. RMSE of reg.: RMSE of regression line.

The difference in the slope value of the regression coefficients is 0.1600 ($=\mathrm{1.0477}-\mathrm{0.8877}$). The 95 % confidence interval of the coefficient is about ±0.04 during both the daytime and the nighttime. It cannot be said that the slopes of the two regression lines are equal based on their 95 % confidence intervals. However, the correlation between NIES/HSRL and POM-02 AOD is high, and the differences between them and their RMSEs are similar. Furthermore, as shown in Fig. 7c, the scatter diagrams for the daytime and nighttime are almost overlapping, and it seems that the two sets of measurements obtained similar results.

Examples of time series of the AOD from NIES/HSRL and POM-02 are shown in Fig. 8. As can be seen from these figures, the AODs of the daytime and nighttime estimated from POM-02 constitute a continuous series. The AOD from NIES/HSRL and that from POM-02 have qualitatively similar time variations. However, in these limited examples, while there are periods when the values are consistent, there are periods when there are systematic differences.

In the NIES/HSRL data processing, the AOD below an altitude of 500 m is calculated by using the value of the extinction coefficient for an altitude of 500 m. Since the height of the atmospheric boundary layer is typically 1500 to 2000 m, a large amount of aerosols exist at altitudes below 500 m. If the actual distribution deviates from the assumed distribution, the estimated AOD is shifted systematically.

In Fig. 8, only limited examples were shown, but in the Supplement, the time series of the AOD at 500 nm at Tsukuba for 5 months is shown in Fig. S1. In addition, the time series of the comparison between the HSRL and POM-02 AOD for 5 months is shown in Fig. S3.

The PWV estimated from POM-02 was compared with that obtained from the vertical profile of the radiosonde and that obtained from the GPS receiver. The PWV estimated from the radiosonde data has a frequency of two values per day, whereas the PWV obtained from GPS is continuous.

The PWV from a radiosonde is often used as a reference for the PWV measurement value. The PWV from the radiosonde and PWV from POM-02 are first compared. Figure 9 shows a scatter plot of the PWV from the radiosonde and from POM-02. The red symbol denotes 00:00 UTC (09:00 LST), and the blue symbol is 12:00 UTC (21:00 LST). Table 5 shows the results of the comparison between the radiosonde and POM-02 precipitable water vapor (Table 5 is the same as Table 4 except for radiosonde and POM-02 precipitable water vapor).

PWV, bias, RMSE, RMSE of reg.: g cm^{−2}. *C*_{1} and *C*_{2}:
coefficients of regression line (${\mathrm{PWV}}_{\mathrm{POM}-\mathrm{02}}={C}_{\mathrm{1}}\cdot {\mathrm{PWV}}_{\mathrm{Sonde}}+{C}_{\mathrm{2}}$).

The ratio of PWV estimated from POM-02 and the radiosonde in both daytime and nighttime is almost constant: the slope of the regression line is 0.80 in the daytime and 0.78 in the nighttime.

The empirical formula of the transmittance is expressed as Eq. (5). The
ratio of the two PWVs is almost constant. In addition, as shown in Fig. 5,
the modified Langley plot provides a good fit for the data. From these
facts, it seems that the value of the coefficient *b* in Eq. (5) is
appropriate but the value of the coefficient *a* in Eq. (5) was
inappropriate. It is possible that the filter characteristics of the 940 nm
channel have changed from the nominal characteristics due to degradation.

Let $\mathrm{pwv}=c\cdot {\mathrm{pwv}}^{\prime}$ and rewrite Eq. (5) as follows:

$$\begin{array}{}\text{(14)}& \begin{array}{rl}Tr\left({\mathrm{H}}_{\mathrm{2}}\mathrm{O}\right)& =\mathrm{exp}(-a(m\left(\mathit{\theta}\right)\cdot (c\cdot {\mathrm{pwv}}^{\prime}){)}^{b})\\ & =\mathrm{exp}(-a{c}^{b}(m\left(\mathit{\theta}\right)\cdot {\mathrm{pwv}}^{\prime}{)}^{b}).\end{array}\end{array}$$

Then, the PWV can be corrected by replacing *a* with *a**c*^{b}.

Figure 10 shows a scatter plot of the PWV from the radiosonde and the
corrected PWV from POM-02. For the correction coefficient *c*, the average
value of the coefficients *C*_{1} of the daytime and nighttime regression
equations was used. Table 6 shows the results of the comparison between the
radiosonde and corrected POM-02 precipitable water vapor (Table 6 is the
same as Table 4 except for radiosonde and corrected POM-02 precipitable
water vapor).

*C*_{1} and *C*_{2}: coefficients of regression line (${\mathrm{PWV}}_{\mathrm{POM}-\mathrm{02}}\left(\mathrm{corrected}\right)={C}_{\mathrm{1}}\cdot {\mathrm{PWV}}_{\mathrm{Sonde}}+{C}_{\mathrm{2}}$).

The slope *C*_{1} of the regression line during the daytime and nighttime
is 1.0160 and 0.9869, respectively, and the difference between them is
0.0291 ($=\mathrm{1.0160}-\mathrm{0.9869}$). The 95 % confidence intervals of the
slopes during the daytime and nighttime are ±0.0206 and ±0.0271, respectively. The difference between them is 0.0291, which is larger
than the respective 95 % confidence intervals. Therefore, the two slopes
are not equivalent based on the 95 % confidence intervals.

However, since the slope of the regression line determined using all of the data is 1.0042 and the 95 % confidence interval is ±0.0173, the three slopes of the regression lines can be regarded as equivalent at the 95 % confidence level. Furthermore, there are no large differences in the bias, RMSE, and correlation coefficient between PWV from the radiosonde and POM-02. Therefore, the PWVs of daytime and nighttime for POM-02 are statistically equivalent. That is, both PWVs are presumed to have the same degree of precision and accuracy within the measurement uncertainty.

Next, the result of the comparison between the PWV obtained from POM-02 and GPS is shown. Before that, the result of the comparison between the PWV obtained from GPS and the radiosonde is shown in Fig. 11. Table 7 shows the results of the comparison between GPS and radiosonde precipitable water vapor (Table 7 is the same as Table 4 except for GPS and radiosonde precipitable water vapor).

*C*_{1} and *C*_{2}: coefficients of regression line (${\mathrm{PWV}}_{\mathrm{Sonde}}={C}_{\mathrm{1}}\cdot {\mathrm{PWV}}_{\mathrm{GPS}}+{C}_{\mathrm{2}}$).

The slope of the regression line in Fig. 11 is about 0.94. In the region of
the PWV of less than 2 g cm^{−2}, the PWV from GPS tends to be smaller than
the PWV from the radiosonde. In the region of PWV of more than 3 g cm^{−2},
the difference between PWV from GPS and the radiosonde is more scattered.
Therefore, the slope of the regression line became smaller than 1. In a
previous comparison conducted by the authors, the slope of the regression
line was almost 1 (Uchiyama et al., 2014). There is a possibility that the
PWV from GPS used in this study has a larger error than the PWV used
previously.

Figure 12 shows a scatter diagram of the PWV from GPS and the corrected PWV from POM-02. Table 8 shows the results of the comparison between PWV from GPS and corrected PWV from POM-02 (Table 8 is the same as Table 4 except for GPS and corrected POM-02 precipitable water vapor).

*C*_{1} and *C*_{2}: coefficients of regression line (${\mathrm{PWV}}_{\mathrm{POM}-\mathrm{02}}\left(\mathrm{corrected}\right)={C}_{\mathrm{1}}\cdot {\mathrm{PWV}}_{\mathrm{GPS}}+{C}_{\mathrm{2}}$).

The slope of the regression line is about 0.91 for both the daytime and
nighttime. Similar to the results of the comparison between the PWV from the
radiosonde and GPS, in the region of PWV from GPS less than 2 g cm^{−2},
the PWV from GPS tends to be somewhat smaller than the PWV from POM-02
during both the daytime and nighttime. In the region of PWV greater than
3 g cm^{−2}, the difference between the PWV from GPS and the radiosonde is
more scattered.

The difference between the slopes of the regression lines is 0.0076 ($=\mathrm{0.9132}-\mathrm{0.9056}$) and the 95 % confidence intervals during the daytime and nighttime are ±0.0097 and ±0.0221, respectively. Therefore, the confidence intervals of the two slopes are overlapping, and the values of slopes can be regarded as equivalent at the 95 % confidence level.

In Fig. 12c, the scatter plot obtained using nighttime data is shown together with that obtained using daytime data. The data obtained during the daytime and nighttime overlap, and it seems that the PWV from POM-02 during the daytime and nighttime are estimated with the same degree of precision and accuracy.

Examples of time series of PWV from GPS and POM-02 are shown in Fig. 13. The PWV from GPS and that from POM-02 have qualitatively similar time variations. In these limited examples, although there are some systematic differences in Fig. 13b, the PWV from GPS and the PWV from POM-02 almost overlap in Figs. 13a and c. In addition, the PWV during the daytime and nighttime estimated from POM-02 are continuously connected.

In Fig. 13, only limited examples were shown, but in the Supplement, the time series of the PWV at Tsukuba for 5 months is shown in Fig. S2. In addition, the time series of the comparison between GPS and POM-02 PWV for 5 months is shown in Fig. S4.

The comparison of the AOD (PWV) before and after sunrise and sunset is used to evaluate the moon photometry (Berkoff et al., 2011; Barreto et al., 2013, 2016, 2017, 2019).

Before and after sunrise (sunset), the AOD before sunrise (after sunset) is
the average of the data with a solar altitude angle between −10 and −15^{∘},
with a lunar phase angle less than 100^{∘}, and with a lunar
altitude angle of more than 10^{∘}. The AOD after sunrise (before
sunset) is the average of the data with a solar altitude angle between 10
and 15^{∘}. Since this comparison is effective when the atmosphere is
stable, only data with small variations were selected; standard deviation/average
value is less than 0.1 or standard deviation is less than 0.02.

Figure 14 shows a scatter plot of the AOD at the wavelengths of 340, 380, 400, 500, 675, 870, and 1020 nm, and the PWV from the 940 nm channel. Table 9 shows the results of the comparison between the AOD (PWV) from the Sun and from the moon (the contents of Table 9 are the same as Table 4 except for the AOD (PWV) from the Sun and the moon).

*C*_{1} and *C*_{2}: coefficients of the regression line (${\mathit{\tau}}_{\mathrm{Moon}}={C}_{\mathrm{1}}\cdot {\mathit{\tau}}_{\mathrm{Sun}}+{C}_{\mathrm{2}}$, ${\mathrm{PWV}}_{\mathrm{Moon}}={C}_{\mathrm{1}}\cdot {\mathrm{PWV}}_{\mathrm{Sun}}+{C}_{\mathrm{2}}$).

The biases at wavelengths of 340 and 380 nm are relatively large, 0.05 and 0.03, respectively, but the biases at other wavelengths are 0.007 or less. The bias and RMSE of the PWV are 0.02 and 0.14, respectively, which are comparable to those from the comparison with POM-02 and the radiosonde or GPS. The correlation coefficient is high for all wavelengths; 0.65 at a wavelength of 340 nm, and 0.97 or higher at other wavelengths. Furthermore, the 95 % confidence interval of the slope value of the regression line includes 1, and the 95 % confidence interval of the intercept value includes 0. That is, the regression line is not different from a straight line with a slope of 1 and zero intercept at the 95 % confidence level. From these facts, the AOD and PWV retrieved using the moon as the light source are considered to be the same as those retrieved using the Sun as the light source at the 95 % confidence level.

There is an AERONET observation site at MLO. In the nighttime, the AODs at wavelengths of 500, 675, 870, and 1020 nm, and the PWV can be compared. In addition to these channels, the AOD at wavelengths of 340, 380, and 1627 nm can be compared in the daytime. The AERONET data used here are level 2.0 in the daytime and level 1.5 in the nighttime. There were no level 2.0 nighttime data. AERONET level 1.5 are cloud-screened data but may not have had the final calibration applied. Thus, these data are not quality assured. AERONET level 2.0 has pre- and post-field calibration applied, cloud-screened, and quality-assured data (see the AERONET homepage; https://aeronet.gsfc.nasa.gov/, last access: 2 December 2019). The nighttime comparison in this paper uses the AERONET data without quality assurance.

Figure 15 shows a scatter plot of the AERONET and POM-02 AOD (PWV). The blue (red) symbols show the daytime (nighttime) data. Both the daytime and the nighttime data are overlaid: the AOD at wavelengths of 500, 675, 870, 1020 nm, and the PWV from the 940 nm channel. The plotted data are the 15 min averages. The number of measurements for POM-02 in a 15 min interval is 10 to 16, and that for AERONET is 1 to 6. Only POM-02 data showing small variations were selected; standard deviation/average is less than 0.1 or standard deviation is less than 0.02.

Table 10 shows the results of the comparison between the AERONET and POM-02 aerosol optical depth (precipitable water vapor); Table 10 is the same as Table 4 except for AERONET and POM-02 aerosol optical depth (precipitable water vapor). The values at 940 nm are the precipitable water vapor.

In the daytime, from Fig. 15, it can be seen that the differences between AERONET and POM-02 AOD (PWV) are small. The 95 % confidence interval for the slope of the regression line does not necessarily include 1, but the slope value is nearly 1: between 0.97 and 1.11. The 95 % confidence interval for the intercept of the regression line does not necessarily include 0, but the magnitude of the intercept is 0.01 or less except for the 380 nm channel (0.015). The same can be said for the PWV of the 940 nm channel. In addition, the bias and RMSE are less than 0.01 except for the 380 nm channel (0.015), and those for the PWV at 940 nm are 0.018 and 0.022, respectively. Considering that the accuracy of the calibration constant is 0.5 % to 1 %, these values seem reasonable. Therefore, it can be inferred that in the daytime, POM-02 can measure the AOD (PWV) with the same level of accuracy as AERONET.

In the nighttime, the atmosphere observed at MLO was pristine, and most of the AODs at 500, 675, 870, and 1020 nm were below 0.02. Considering that the accuracy of the calibration constant is 0.5 % to 1 %, it is difficult to compare the AODs of AERONET and POM-02. In the nighttime, the slope of the regression line deviates from 1 at several wavelengths, but the bias and the RMSE are less than about 0.01. Therefore, the difference between AERONET and POM-02 is small. The slopes of the regression line for the PWV of 940 nm channel in the daytime and the nighttime are 1.07 and 1.16, respectively. Thus, the daytime and nighttime values differ. In the results of Sect. 6.3, there is almost no difference between the daytime and nighttime values. Therefore, this difference may be due to the lack of quality control in the nighttime data.

7 Summary and conclusion

Back to toptop
Aerosol data are often estimated using the solar direct irradiance and the solar scattered radiance. Therefore, the majority of data on aerosol properties are obtained using daytime measurements, and there are few data available on aerosol characteristics at night. In order to estimate the aerosol optical depth (AOD) and the precipitable water vapor (PWV) during the nighttime using the moon as a light source, POM-02 (Prede Ltd., Japan), which is used to estimate aerosol characteristics during the daytime, was modified.

The current version of POM-02 has the ability to measure the direct irradiance from the moon for some channels in the visible and near-infrared wavelength region without requiring modification. Several modifications were made to also be able to measure the AOD during the nighttime and expand the measurement ranges.

The amplifier was adjusted so that POM-02 could measure up to about $\mathrm{5}\times {\mathrm{10}}^{-\mathrm{13}}$ A, allowing the lunar direct irradiance to be measured in the wavelength range of 340 to 1020 nm.

In order to track the moon based on the calculated value, the simplified formula by Nagasawa (1981) was incorporated into the firmware.

A position sensor with a four-quadrant detector is used to adjust the tracking of the Sun and the moon. In order to adjust the tracking of the moon, a position sensor incorporating a new electronic circuit to amplify the signal and new software to process the signal data were developed. The new position sensor can be used to track both the Sun and the moon.

The calibration constant was determined by using the Langley method. The measurements of the solar and lunar direct irradiance were conducted at the NOAA/MLO during the period from 28 September to 7 November 2017. Assuming that the correct lunar reflectance is proportional to the ROLO reflectance, the calibration constant for the lunar direct irradiance was determined by using the Langley method. The calibration by the Langley method was successfully performed.

The ratio of the calibration constant for the moon to that of the Sun was often greater than 1, where the ratio is a coefficient for correcting the ROLO reflectance and includes a smoothing factor. This ratio shows the error of the ROLO irradiance model. The value of the ratio was 0.95 to 1.18 in the visible and near-infrared wavelength region. This means that the ROLO model often underestimates the reflectance. In addition, this ratio depended on the phase angle: when the phase angle was small (near the full moon), the ratio was small, and as the phase angle became larger, the ratio increased. In this study, this ratio was approximated by the quadratic equation of the phase angle. By using this approximation, the reflectance of the moon can be calculated to within an accuracy of 1 % or less.

The continuous measurement of POM-02 was conducted at JMA/MRI from January 2018 to May 2018 and the AOD and PWV were estimated. In order to validate the estimates of the AOD and PWV, we compared them with the AOD and PWV obtained by independent methods. The AOD was compared with the AOD (532 nm) estimated from NIES/HSRL, and the PWV was compared with the PWV from a radiosonde and GPS. In addition, the continuity of the AOD (PWV) before and after sunrise and sunset at Tsukuba was examined, and the AOD (PWV) of AERONET and that of POM-02 at MLO were compared.

Concerning the AOD, there were sometimes systematic differences between NIES/HSRL and POM-02. The cause of the systematic differences seems to be that NIES/HSRL assumes a constant extinction coefficient at altitudes of less than 500 m. The slopes of the linear regression lines during the daytime and nighttime could not be said to be equivalent at the 95 % confidence level, but the scatter diagrams of the daytime and nighttime were almost overlapping.

Concerning the PWV, the slopes of the linear regression lines during the daytime and nighttime were equivalent at the 95 % confidence level in the comparisons between the PWV from POM-02 and the radiosonde and in the comparison between the PWV from POM-02 and GPS. Furthermore, the scatter diagrams of the daytime and the nighttime data were almost overlapping.

In addition, the comparison of the AOD (PWV) before and after sunrise and sunset showed that the AOD and PWV retrieved using the moon as the light source are the same as those retrieved using the Sun as the light source at the 95 % confidence level.

The comparison of the AOD (PWV) between AERONET and POM-02 was performed using the data taken during the calibration measurements. The comparison in the daytime showed that POM-02 can measure AOD (PWV) with the same accuracy as AERONET. The comparison in the nighttime showed that the difference in the AOD between AERONET and POM-02 was small. However, since there were a lot of optically thin data and AERONET data are not quality assured, we cannot make a definite conclusion.

From these facts, the daytime and nighttime AOD (PWV) measurements are statistically almost equivalent. The AODs (PWVs) during the daytime and nighttime for POM-02 are presumed to have the same degree of precision and accuracy within the measurement uncertainty.

The accuracy of the nighttime calibration constant is lower than that for
the daytime. The measurement *S*∕*N* in the nighttime is also worse than that in
daytime. Considering these facts, even if there is no statistically
significant difference, the magnitude of the error in the AOD (PWV) during
the nighttime is not always the same as that during the daytime.

In this study, the calibration was performed using about 40 d of data
including two full-moon days. As a result, it was found that there was an
error in the reflectance of the ROLO irradiance model. In the future, it is
necessary to accumulate more data for calibration and to reduce the error of
the ROLO irradiance model. It is said that the ROLO model can be applied
over a phase angle range of about 90^{∘}. POM-02 has the ability to
measure the direct lunar irradiance up to a phase angle range of about 120^{∘}.
It is necessary to expand the ROLO irradiance model so that it can
be applied to larger phase angles.

It is now possible to estimate the aerosol optical depth during the nighttime. It is necessary to promote the adoption of this system in the existing observation network. After that, the data obtained by using this instrument can be used to better understand nighttime aerosol behavior, for the validation of aerosol transport models, and as input data in assimilation systems.

Data availability

Back to toptop
Data availability.

The data used in this study are available from the corresponding author.

Appendix A

Back to toptop
The smoothed ROLO reflectance ${\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}$ can be obtained by the procedure described in Kieffer and Stone (2005).

The calculated reflectance ${A}_{\mathrm{ROLO}}(g,\mathrm{\Phi},\mathit{\theta},\mathit{\varphi})$ at the 32 ROLO
wavelengths for a specific geometric configuration $(g=\mathrm{7}{}^{\circ},\mathrm{\Phi}=\mathrm{7}{}^{\circ},\mathit{\theta}=\mathrm{0},\mathit{\varphi}=\mathrm{0})$ is fitted to a composite spectrum of the samples
obtained by the Apollo 16 mission with a linear equation of wavelength
*λ*.

$$\begin{array}{}\text{(A1)}& {A}_{\mathrm{Apollo}}=(a+b\mathit{\lambda}){A}_{\mathrm{ROLO}}(\mathrm{7},\mathrm{7},\mathrm{0},\mathrm{0}),\end{array}$$

where *A*_{Apollo} is the composite laboratory reflectance spectrum for the
Apollo samples of soil (95 %) (Apollo 16 sample 62231; Pieters, 1999) and
breccia (5 %) (Apollo 16 sample 67455; Pieters and Mustard, 1988).

The Apollo sample 62231 spectrum is available at http://www.planetary.brown.edu/pds/AP62231.html (last access: 2 December 2019). The Apollo sample 67455 spectrum is shown in Fig. 8 in the paper of Pieters and Mustard (1988).

The values of the coefficients *a* and *b* are not shown in Kieffer and Stone (2005)
but were determined here with the least-squares method as follows:

$$\begin{array}{ll}{\displaystyle}a& {\displaystyle}=\mathrm{1.640875}\\ {\displaystyle}b& {\displaystyle}=-\mathrm{1.192034}\times {\mathrm{10}}^{-\mathrm{4}}\\ {\displaystyle},\end{array}$$

where the unit of the wavelengths is nanometers.

By dividing *A*_{Apollo} by *a*+*b**λ*, the smoothed ROLO reflectance
for a specific geometric configuration ${\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}(\mathrm{7},\mathrm{7},\mathrm{0},\mathrm{0})$ can be
obtained.

$$\begin{array}{}\text{(A2)}& {\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}(\mathrm{7},\mathrm{7},\mathrm{0},\mathrm{0})={A}_{\mathrm{Apollo}}/(a+b\mathit{\lambda})\end{array}$$

The smoothed ROLO reflectance ${\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}(g,\mathrm{\Phi},\mathit{\theta},\mathit{\varphi})$ for any viewing geometry is given by the following equation:

$$\begin{array}{}\text{(A3)}& \begin{array}{rl}{\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}& ={\displaystyle \frac{{\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}(\mathrm{7},\mathrm{7},\mathrm{0},\mathrm{0})}{{A}_{\mathrm{ROLO}}(\mathrm{7},\mathrm{7},\mathrm{0},\mathrm{0})}}{A}_{\mathrm{ROLO}}(g,\mathrm{\Phi},\mathit{\theta},\mathit{\varphi})\\ & ={F}_{C}{A}_{\mathrm{ROLO}}(g,\mathrm{\Phi},\mathit{\theta},\mathit{\varphi}),\end{array}\end{array}$$

where ${F}_{C}={\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}(\mathrm{7},\mathrm{7},\mathrm{0},\mathrm{0})/{A}_{\mathrm{ROLO}}(\mathrm{7},\mathrm{7},\mathrm{0},\mathrm{0})$.

The values of *F*_{C} are dependent on the interpolation method of the
reflectance table and the accuracy of the values read from the figure. The
smoothed and adjusted spectrum ${\stackrel{\mathrm{\u0303}}{A}}_{\mathrm{ROLO}}(\mathrm{7},\mathrm{7},\mathrm{0},\mathrm{0})$ is shown in
Fig. A1. The values of *F*_{C} determined by the authors are shown in Table A1 and Fig. A2.

Supplement

Back to toptop
Supplement.

The supplement related to this article is available online at: https://doi.org/10.5194/amt-12-6465-2019-supplement.

Author contributions

Back to toptop
Author contributions.

This study was designed by AU, MS, HK, and TM. The measurements for the sky radiometer were conducted by AU, AK, KI, and YW. The adjustment of the amplifier and the development of the position sensor were performed by MS, HK, KI, KK, and YW. The development of the related software and the data analyses were performed by AU. The manuscript was written by AU, and all authors contributed to editing and revision.

Competing interests

Back to toptop
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement

Back to toptop
Special issue statement.

This article is part of the special issue “SKYNET – the international network for aerosol, clouds, and solar radiation studies and their applications (AMT/ACP inter-journal SI)”. It is not associated with a conference.

Acknowledgements

Back to toptop
Acknowledgements.

We would like to thank Yoshitake Jin and Tomoaki Nishizawa of NIES for providing the NIES/HSRL data for the comparison of the aerosol optical depth. We also would like to thank Yoshinori Shoji of JMA/MRI for providing the GPS data for the comparison of precipitable water vapor. We thank Brent N. Holben and his staff for their effort in establishing and maintaining the AERONET Mauna Loa site. We would like to thank Thomas Stone and two anonymous reviewers for their useful comments.

Financial support

Back to toptop
Financial support.

This work was supported by the NIES GOSAT-2 project, Japan. This work was also supported by JSPS KAKENHI grant no. 17K00531.

Review statement

Back to toptop
Review statement.

This paper was edited by Monica Campanelli and reviewed by two anonymous referees.

References

Back to toptop
Acton Jr, C. H.: Ancillary data services of NASA's Navigation and Ancillary Information Facility, Planet. Space Sci., 44, 65–70, 1996.

Asano, S., Murai, K., and Yamauchi, T.: An improvement of the computation method of the atmospheric turbidity factors, J. Meteorol. Res., 35, 135–144, 1983. (in Japanese)

Barreto, A., Cuevas, E., Damiri, B., Guirado, C., Berkoff, T., Berjón, A. J., Hernández, Y., Almansa, F., and Gil, M.: A new method for nocturnal aerosol measurements with a lunar photometer prototype, Atmos. Meas. Tech., 6, 585–598, https://doi.org/10.5194/amt-6-585-2013, 2013.

Barreto, Á., Cuevas, E., Granados-Muñoz, M.-J., Alados-Arboledas, L., Romero, P. M., Gröbner, J., Kouremeti, N., Almansa, A. F., Stone, T., Toledano, C., Román, R., Sorokin, M., Holben, B., Canini, M., and Yela, M.: The new sun-sky-lunar Cimel CE318-T multiband photometer – a comprehensive performance evaluation, Atmos. Meas. Tech., 9, 631–654, https://doi.org/10.5194/amt-9-631-2016, 2016.

Barreto, Á., Román, R., Cuevas, E., Berjón, A. J., Almansa, A. F., Toledano, C., González, R., Hernández, Y., Blarel, L., Goloub, P., Guirado, C., and Yela, M.: Assessment of nocturnal aerosol optical depth from lunar photometry at the Izaña high mountain observatory, Atmos. Meas. Tech., 10, 3007–3019, https://doi.org/10.5194/amt-10-3007-2017, 2017.

Barreto, A., Román, R., Cuevas, E., Pérez-Ramírez, D., Berjón, A. J., Kouremeti, N., Kazadzis, S., Gröbner, J., Mazzola, M., Toledano, C., Benavent-Oltra, J. A., Doppler, L., Juryšek, J., Almansa, A. F., Victori, S., Maupin, F., Guirado-Fuentes, C., González, R., Vitale, V., Goloub, P., Blarel, L., Alados-Arboledas, L., Woolliams, E., Taylor, S., Antuña, J. C., and Yela, M.: Evaluation of night-time aerosols measurements and lunar irradiance models in the frame of the first multi-instrument nocturnal intercomparison campaign, Atmos. Environ., 202, 190–211, https://doi.org/10.1016/j.atmosenv.2019.01.006, 2019.

Berkoff, T. A., Sorokin, M., Stone, T., Eck, T. F., Hoff, R., Welton, E., and Holben, B.: Nocturnal aerosol optical depth measurements with a small-aperture automated photometer using the moon as a light source, J. Atmos. Ocean. Tech., 28, 1297–1306, https://doi.org/10.1175/JTECH-D-10-05036.1, 2011.

Bruegge, C. J., Conel, J. E., Green, R. O., Margolis, J. S., Holm, R. G., and Toon, G.: Water vapor column abundance retrievals during FIFE, J. Geophys. Res., 97, 18759–18768, 1992.

Cachorro, V. E., Utrillas, P., Vergaz, R., Duran, P., de Frutos, A. M., and Martinez-Lozano, J. A.: Determination of the atmospheric-water-vapor content in the 940 nm absorption band by use of moderate spectral-resolution measurements of direct solar irradiance, Appl. Opt., 37, 4678–4689, 1998.

Campanelli, M., Nakajima, T., Khatri, P., Takamura, T., Uchiyama, A., Estelles, V., Liberti, G. L., and Malvestuto, V.: Retrieval of characteristic parameters for water vapour transmittance in the development of ground-based sun–sky radiometric measurements of columnar water vapour, Atmos. Meas. Tech., 7, 1075–1087, https://doi.org/10.5194/amt-7-1075-2014, 2014.

Dockery, D. W., Pope, C. A., Xu, X., Spengler, J. D., Ware, J. H., Fay, M. E., Ferris, Jr., B. G., and Speizer, F. E.: An Association between Air Pollution and Mortality in Six U.S. Cities, New Engl. J. Med., 329, 1753–1759, 1993.

Dutton, E.G., Reddy, P., Ryan, S., and DeLuisi, J.: Features and Effects of Aerosol Optical Depth Observed at Mauna Loa, Hawaii: 1982–1992, J. Geophys. Res., 99, 8295–8306, 1994.

Esposito, F., Serio, C., Pavese, G., Auriemma, G., and Satriano, C.: Measurements of nighttime atmospheric optical depth preliminary data from a mountain site in southern Italy, J. Aerosol Sci., 29, 1213–1218, 1998.

Esposito, F., Mari, S., Pavese, G., and Serio, C.: Diurnal and Nocturnal Measurements of Aerosol Optical Depth at a Desert Site in Namibia, Aerosol Sci. Technol., 37, 392–400, https://doi.org/10.1080/02786820300972, 2003.

Fernald, F. G.: Analysis of atmospheric lidar observations: some comments, Appl. Opt., 5, 652–653, https://doi.org/10.1364/AO.23.000652, 1984.

Guerrero-Rascado, J. L., Landulfo, E., Antuña, J. C., Barbosa, H., de M. J., Barja, B., Bastidas, Á. E., Bedoya, A. E., da Costa, R. F., Estevan, R., Forno, R., Gouveia, D. A., Jiménez, C., Larroza E. G., da Silva Lopes, F. J., Montilla-Rosero, E., Moreira, G. de A., Nakaema, W. M., Nisperuza, D., Alegria, D., Múnera, M., Otero, L., Papandrea, S., Pallota, J. V., Pawelko, E. Quel, E. J., Ristori, P., Rodrigues, P. F., Salvador, J., Sánchez, M. F., and Silva, A.: Latin American Lidar Network (LALINET) for aerosol research: diagnosis on network instrumentation, J. Atmos. Sol. Terr. Phys., 138–139, 112–120, 2016.

Halthore, R. N., Eck, T. F., Holben, B. N., and Markham, B. L.: Sun photometric measurements of atmospheric water vapor column abundance in the 940-nm band, J. Geophys. Res., 102, 4343–4352, https://doi.org/10.1029/96JD03247, 1997.

Herber, A., Thomason, L. W., Gernandt, H., Leiterer, U., Nagel, D., Schulz, K-H., Kaptur, J., Albrecht, T., and Notholt, J.: Continuous day and night aerosol optical depth observations in the Arctic between 1991 and 1999, J. Geophys. Res., 107, 4097, https://doi.org/10.1029/2001JD000536, 2002.

Holben, B. N., Eck, T. F., Slutsker, I., Tanré, D., Buis, J. P., Setzer, A., Vermote, E., Reagan, J. A., Kaufman, Y. J., Nakajima, T., Lavenu, F., Jankowiak, I., and Smirnov, A.: AERONET-A federated instrument network and data archive for aerosol characterization, Remote Sens. Environ., 66, 1–16, 1998.

Ingold, T., Schmid, B., Matzler, C., Demoulin, P., and Kampfer, N.: Modeled and empirical approaches for retrieving columnar water vapor from solar transmittance measurements in the 0.72, 0.82, and 0.94 µm absorption bands, J. Geophys. Res., 105, 24327–24343, 2000.

Juryšek, J. and Prouza, M.: Sun/Moon photometer for the Cherenkov Telescope Array – first results, Proceedings of Science, 35th International Cosmic Ray Conference-ICRC2017, 10–20 July, 2017 Bexco, Busan, Korea, 2017.

Kahn, R. A., Gaitley, B. J., Martonchik, J. V., Diner, D. J., Crean, K. A., and Holben, B.: Multiangle Imaging Spectroradiometer (MISR) global aerosol optical depth validation based on 2 years of coincident Aerosol Robotic Network (AERONET) observations, J. Geophys. Res., 110, D10S04, https://doi.org/10.1029/2004JD004706, 2005.

Kasten, F. and Young, A. T.: Revised optical air mass tables and approximation formula, Appl. Opt., 28, 4735–4738, 1989.

Kiedron, P., Michalsky, J., Schmid, B., Slater, D., Berndt, J., Harrison, L., Racette, P., Westwater, E., and Han, Y.: A robust retrieval of water vapor column in dry Arctic conditions using the rotating shadowband spectroradiometer, J. Geophys. Res., 106, 24007–24016, 2001.

Kiedron, P., Berndt, J., Michalsky, J., and Harrison, L.: Column water vapor from diffuse irradiance, Geophys. Res. Lett., 30, 1565, https://doi.org/10.1029/2003GL016874, 2003.

Kieffer, H. H.: Photometric stability of the lunar surface, Icarus, 130, 323–327, 1997.

Kieffer, H. H. and Stone, T. C.: The spectral irradiance of the moon, Astron. J., 129, 2887–2901, 2005.

Kikuchi, M., Murakami, H., Suzuki, K., Nagao, T. M., and Higurashi, A.: Improved Hourly Estimates of Aerosol Optical Thickness Using Spatiotemporal Variability Derived From Himawari-8 Geostationary Satellite, IEEE T. Geosci. Remote, 56, 3442–3455, https://doi.org/10.1109/TGRS.2018.2800060, 2018.

Klett, J. D.: Lidar inversion with variable backscatter/extinction ratios, Appl. Opt., 11, 1638–1643, https://doi.org/10.1364/AO.24.001638, 1985.

Kouremeti, N., Gröbner, J., Kazadzis, S., Pfiffner, D., and Soder, R.: Development of a Lunar PFR, 21 pp., https://www.pmodwrc.ch/wp-content/uploads/2017/09/2015_Annual_Report.pdf (last access: 2 December 2019), 2016.

Levis, J. R., Campbell, J. R., Welton, E. J., Stewart, S. A., Phillip C., and Haftings, P. C.: Overview of MPLNET version 3 cloud detection, J. Atmos. Ocean. Technol., 33, 2113–2134, https://doi.org/10.1175/JTECH-D-15-0190.1, 2016.

Lohmann, U. and Feichter, J.: Global indirect aerosol effects: a review, Atmos. Chem. Phys., 5, 715–737, https://doi.org/10.5194/acp-5-715-2005, 2005.

Michalsky, J. J., Liljegren, J. C., and Harrison, L. C.: A comparison of sun photometer derivations of total column water vapor and ozone to standard measures of same at the Southern Great Plains atmospheric radiation measurement site, J. Geophys. Res., 100, 25995–26003, 1995.

Michalsky, J. J., Min, Q., Kiedron, P. W., Slater, D. W., and Barnard, J. C.: A differential technique to retrieve column water vapor using sun radiometry, J. Geophys. Res., 106, 17433–17442, 2001.

Nagasawa, K.: Tentai no ichi keisan (Position calculation of celestial bodies), Chjin Shokan, p. 239, 1981. (in Japanese)

Neckel, H. and Labs, D.: Improved Data of Solar Spectral Irradiance from 0.33 to 1.25 µm, Sol. Phy., 74, 231–249, https://doi.org/10.1007/BF00151293, 1981.

Pappalardo, G., Amodeo, A., Apituley, A., Comeron, A., Freudenthaler, V., Linné, H., Ansmann, A., Bösenberg, J., D'Amico, G., Mattis, I., Mona, L., Wandinger, U., Amiridis, V., Alados-Arboledas, L., Nicolae, D., and Wiegner, M.: EARLINET: towards an advanced sustainable European aerosol lidar network, Atmos. Meas. Tech., 7, 2389–2409, https://doi.org/10.5194/amt-7-2389-2014, 2014.

Pérez-Ramírez, D., Aceituno, J., Ruiza, B., Olmo, F. J., and Alados-Arboledas, L.: Development and calibration of a star photometer to measure the aerosol optical depth: Smoke observations at a high mountain site, Atmos. Environ., 42, 2733–2738, 2008.

Perry, K. D., Cahill, T. A., Schnell, R. C., and Harris, J. M.: Long-range transport of anthropogenic aerosols to the National Oceanic and Atmospheric Administration baseline station at Mauna Loa Observatory, Hawaii, J. Geophys. Res., 104, 18521–18533, 1999.

Pieters, C. M.: The Moon as a Calibration Standard Enabled by Lunar Samples, in: New Views of the Moon II: Understanding the Moon through the Integration of Diverse Datasets, edited by: Gaddis, L. and Shearer, C. K., 47, 1999.

Pieters, C. M. and Mustard, J. F.: Exploration of crustal/mantle material for the earth and moon using reflectance spectroscopy, Remote Sens. Environ., 24, 151–178, 1988.

Plana-Fattori, A., Legrand, M., Tanre, D., Devaux, C., Vermeulen, A., and Dubuisson, P.: Estimating the atmospheric water vapor content from sun photometer measurements, J. Appl. Meteorol., 37, 790–804, 1998.

Plana-Fattori, A., Dubuisson, P., Fomin, B. A., and de Paula Correa, M.: Estimating the atmospheric water vapor content from multi-filter rotating shadow-band radiometry at Sao Paulo, Brazil, Atmos. Res., 71, 171–192, 2004.

Prados, A. I., Kondragunta, S., Ciren, P., and Knapp, K. R.: GOES Aerosol/Smoke Product (GASP) over North America: Comparisons to AERONET and MODIS observations, J. Geophys. Res., 112, D15201, https://doi.org/10.1029/2006JD007968, 2007.

Ramanathan, V., Crutzen, P. J., Kiehl, J. T., and Rosenfeld, D.: Aerosols, Climate, and the Hydrological Cycle, Science, 294, 2119–2124, 2001.

Reagan, J. A., Thome, K., Herman, B., and Gall, R.: Water vapor measurements in the 0.94 micron absorption band-Calibration, measurements and data applications, in: Proc. Int. Geoscience and Remote Sensing '87 Symposium, Ann Arbor, Michigan, 63–67, 1987a.

Reagan, J. A., Pilewskie, P. A., Herman, B. M., and Ben-David, A.: Extrapolation of Earth-based solar irradiance measurements to exoatmospheric levels for broad-band and selected absorption-band observations, IEEE T. Geosci. Remote, 25, 647–653, 1987b.

Reagan, J., Thome, K., Herman, B., Stone, R., DeLuisi, J., and Snider, J.: A comparison of columnar water vapor retrievals obtained with near-IR solar radiometer and microwave radiometer measurements, J. Appl. Meteorol., 34, 1384–1391, 1995.

Remer, L. A., Kaufman, Y. J., Tanré, D., Mattoo, S., Chu, D. A., Martins, J. V., Li, R.-R., Ichoku, C., Levy, R. C., Kleidman, R. G., Eck, T. F., Vermote, E., and Holben, B. N.: The MODIS aerosol algorithm, products and validation, J. Atmos. Sci., 62, 947–973, https://doi.org/10.1175/JAS3385.1, 2005.

Schmid, B. and Wehrli, C.: Comparison of Sun photometer calibration by use of the Langley technique and the standard lamp, Appl. Opt., 34, 4500–4512, 1995.

Schmid, B., Thome, K. J., Demoulin, P., Peter, R., Matzler, C., and Sekler,
J.: Comparison of modeled and empirical approaches for retrieving columnar
water vapor from solar transmittance measurements in the 0.94-*μ*m
region, J. Geophys. Res., 101, 9345–9358, 1996.

Schmid, B., Michalsky, J. J., Slater, D. W., Barnard, J. C., Halthore, R. N., Liljegren, J. C., Holben, B. N., Eck, T. F., Livingston, J. M., Russell, P. B., Ingold, T., and Slutsker, I.: Comparison of columnar water-vapor measurements from solar transmittance methods, Appl. Opt., 40, 1886–1896, 2001.

Shaw, G. E.: Solar spectral irradiance and atmospheric transmission at Mauna Loa Observatory, Appl. Opt., 21, 2006–2011, 1982.

Shaw, G. E.: Sun photometry, B. Am. Meteor. Soc., 64, 4–11, 1983.

Shimizu, A., Nishizawa, T., Jin, Y., Kim, S.-W., Wang, Z., Batdorj, D., and Sugimoto, N.: Evolution of a lidar network for tropospheric aerosol detection in East Asia, Opt. Eng., 56, 031219, https://doi.org/10.1117/1.OE.56.3.031219, 2016.

Shiobara, M., Spinhirne, J. D., Uchiyama, A., and Asano, S.: Optical depth measurements of aerosol, cloud, and water vapor using sun photometers during FIRE Cirrus IFO II, J. Appl. Meteorol., 35, 36–46, 1996.

Shoji, Y.: Retrieval of Water Vapor Anisotropy using the Japanese Nationwide GPS Array and its Potential for Prediction of Convective Precipitation, J. Meteor. Soc. Japan, 91, 43–62, 2013.

Takamura, T., Nakajima, T., and SKYNET community group: Overview of SKYNET and its Activities. Proceedings of AERONET workshop, El Arenosillo. Optica Pura y Aplicada, 37, 3303–3308, 2004.

Thome, K., Herman, B. M., and Reagan, J. A.: Determination of precipitable water from solar transmission, J. Appl. Meteorol., 31, 157–165, 1992.

Thome, K. J., Smith, M. W., Palmer, J. M., and Reagan, J. A.: Three-channel solar radiometer for the determination of atmospheric columnar water vapor, Appl. Opt., 33, 5811–5819, 1994.

Uchiyama, A., Yamazaki, A., and Kudo, R.: Column Water Vapor Retrievals from Sky-radiometer (POM-02) 940 nm Data, J. Meteorol. Soc. Japan, 92A, 195–203, https://doi.org/10.2151/jmsj.2014-A13, 2014.

Uchiyama, A., Matsunaga, T., and Yamazaki, A.: The instrument constant of sky radiometers (POM-02) – Part 1: Calibration constant, Atmos. Meas. Tech., 11, 5363–5388, https://doi.org/10.5194/amt-11-5363-2018, 2018.

Wang, M., Bailey, S., and McClain, C. R.: SeaWiFS provides unique global aerosol optical property data, Eos, Trans. Amer. Geophys. Union, 81, 197–202, 2000.

Wehrli, C.: Extraterrestrial Solar Spectrum, Publication no. 615, Physikalisch-Meteorologisches Observatorium + World Radiation Center (PMO/WRC) Davos Dorf, Switzerland, July, 1985.

Wehrli, C.: GAW–PFR: A network of Aerosol Optical Depth observations with Precision Filter Radiometers, in: WMO/GAW Experts workshop on a global surface based network for long term observations of column aerosol optical properties, Tech. rep., GAW Report No. 162, WMO TD No. 1287, https://community.wmo.int/gaw-reports, (last accsses: 4 December 2019), 2005.

Welton, E. J., Campbell, J. R., Spinhirne, J. D., and Scott, V. S.: Global monitoring of clouds and aerosols using a network of micro-pulse lidar systems, in: Lidar Remote Sensing for Industry and Environmental Monitoring, edited by: Singh, U. N., Itabe, T., and Sugimoto, N., International Society for Optical Engineering, SPIE Proceedings, Vol. 4153, 151–158, 2001.

Winker, D. M., Hunt, W. H., and McGill, M. J.: Initial performance assessment of CALIOP, Geophys. Res. Lett., 34, L19803, https://doi.org/10.1029/2007GL030135, 2007.

World Health Organisation: WHO Air quality guidelines for particulate matter, ozone, nitrogen dioxide and sulfur dioxide; Global update 2005; Summary of risk assessment, World Health Organization, Geneva, Switzerland, p. 22, 2006.

World Health Organization Regional Office for Europe: Review of Evidence on Health Aspects of Air Pollution–REVIHAAP Project, Technical Report WHO, Copenhagen, p. 298, 2013.

Yokota, Y., Matsunaga, T., Ohtake, M., Haruyama, J., Nakamura, R., Yamamoto,
S., Ogawa, Y., Morota, T., Honda, C., Saiki, K., Nagasawa, K., Kitazato, K.,
Sasaki, S., Iwasaki, A., Demura, H., Hirata, N., Hiroi, T., Honda, R.,
Iijima, Y., and Mizutani, H.: Lunar photometric properties at wavelengths
0.5–1.6 *μ*m acquired by SELENE Spectral Profiler and their dependency
on local albedo and latitudinal zones, Icarus, 215, 639–660, 2011.

Yoshida, M., Kikuchi, M., Nagao, T. M., Murakami, H., Nomaki, T., and Higurashi, A.: Common Retrieval of Aerosol Properties for Imaging Satellite Sensors, J. Meteorol. Soc. Jpn., 96B, 193–209, https://doi.org/10.2151/jmsj.2018-039, 2018.

Short summary

The majority of aerosol data are obtained from daytime measurements using the Sun as a light source, and there are few datasets available for studying nighttime aerosol characteristics. To estimate the aerosol optical depth (AOD) during the nighttime using the moon as a light source, a radiometer for the daytime was modified, and a new calibration method was developed. As a result, the estimations of the nighttime AOD were made with the same degree of precision and accuracy during the daytime.

The majority of aerosol data are obtained from daytime measurements using the Sun as a light...

Atmospheric Measurement Techniques

An interactive open-access journal of the European Geosciences Union