the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
The impact and estimation of uncertainty correlation for multiangle polarimetric remote sensing of aerosols and ocean color
Kirk Knobelspiesse
Bryan A. Franz
PengWang Zhai
Brian Cairns
Xiaoguang Xu
J. Vanderlei Martins
Multiangle polarimetric (MAP) measurements contain rich information for characterization of aerosol microphysical and optical properties that can be used to improve atmospheric correction in ocean color remote sensing. Advanced retrieval algorithms have been developed to obtain multiple geophysical parameters in the atmosphere–ocean system, although uncertainty correlation among measurements is generally ignored due to lack of knowledge on its strength and characterization. In this work, we provide a practical framework to evaluate the impact of the angular uncertainty correlation from retrieval results and a method to estimate correlation strength from retrieval fitting residuals. The Fast MultiAngular Polarimetric Ocean coLor (FastMAPOL) retrieval algorithm, based on neuralnetwork forward models, is used to conduct the retrievals and uncertainty quantification. In addition, we also discuss a flexible approach to include a correlated uncertainty model in the retrieval algorithm. The impact of angular correlation on retrieval uncertainties is discussed based on synthetic Airborne HyperAngular Rainbow Polarimeter (AirHARP) and HyperAngular Rainbow Polarimeter 2 (HARP2) measurements using a Monte Carlo uncertainty estimation method. Correlation properties are estimated using autocorrelation functions based on the fitting residuals from both synthetic AirHARP and HARP2 data and real AirHARP measurement, with the resulting angular correlation parameters found to be larger than 0.9 and 0.8 for reflectance and degree of linear polarization (DoLP), respectively, which correspond to correlation angles of 10 and 5^{∘}. Although this study focuses on angular correlation from HARP instruments, the methodology to study and quantify uncertainty correlation is also applicable to other instruments with angular, spectral, or spatial correlations and can help inform laboratory calibration and characterization of the instrument uncertainty structure.
Satellite remote sensing is important for the study of the earth system at a global scale (National Academies of Sciences, Engineering, and Medicine, 2018). Remote sensing instruments are evolving rapidly, with increasing accuracy and spatial, spectral, and angular resolutions (Kokhanovsky et al., 2015; Dubovik et al., 2019). Multiangle polarimeters (MAPs), measuring polarization states at multiple spectral bands and viewing angles, contain high information content for the study of aerosol and cloud optical and microphysical properties (Mishchenko and Travis, 1997; Chowdhary et al., 2001; Hasekamp and Landgraf, 2007; Knobelspiesse et al., 2012). The aerosol properties derived from MAP instruments can be used to assist atmospheric correction for ocean color remote sensing (Frouin et al., 2019; Gao et al., 2020; Hannadige et al., 2021).
Uncertainty quantification from MAP retrievals provides information on the quality of the data products and improves our understanding of retrieval sensitivities. These uncertainties depend on the retrieval algorithm as well as the instrument characterization, including the spectral bands, viewing angles, and polarization capability and the measurement accuracy. As shown in Fig. 1 and Table 1, MAP instruments collect a large number of highquality measurements with differing numbers of spectral bands and viewing angles (Gao et al., 2021b). The number of spectral bands are mostly within 4–14 for MAP instruments. An exception is the SPEX sensors (SPEX airborne and SPEXone, Smit et al., 2019; Hasekamp et al., 2019), which acquire up to 400 spectral bands. The number of viewing angles for many instruments vary between 5 and 16 (e.g., 5 for SPEXone and 9 for SPEX airborne and MISR) but can be on the order of 100 for several hyperangular instruments (90 for HARP2, 120 for AirHARP, 152 for RSP, and 250 for APS). More viewing angles are preferred for the observation of clouds (Waquet et al., 2009; McBride et al., 2020), and they can also be used to conduct multiangle cloud masking and data screening to increase aerosol retrieval accuracy and coverage (Gao et al., 2021b).
Deschamps et al. (1994)Tanré et al. (2011)Fougnie et al. (2018)Diner et al. (1998)Diner et al. (2013)Diner et al. (2018)Cairns et al. (1999)Mishchenko et al. (2007)Martins et al. (2018)Martins et al. (2018)Smit et al. (2019)Hasekamp et al. (2019)To understand the retrieval uncertainties, an uncertainty model is required to describe the combined uncertainties from the MAP measurements, forward model, and a priori assumptions. These combined uncertainty sources are often assumed to be independent and without correlations; however, measurements with high angular or spectral resolution are likely to have correlated uncertainty, depending on instrument design. For example, a sensor may use the same detector to scan through all measurement view angles (e.g. the Research Scanning Polarimeter, RSP, Cairns et al., 1999), and thus the systematic errors due to calibration will be correlated for different measurement angles. The correlation property should be part of the MAP uncertainty model, but often it has not been sufficiently characterized. The characterization of measurement uncertainty correlation is affected by instrument calibration and data processing protocols, which is challenging to quantify. A sensitivity study considering the angular uncertainty correlation of the RSP data was conducted by Knobelspiesse et al. (2012), which showed that information content is affected by correlation strength. However, the actual impacts of uncertainty correlation in a retrieval algorithm have not been well explored, as this requires better understanding of the correlation characteristics and efficient implementation in a retrieval algorithm.
Retrieval algorithms that exploit correlation information in retrieval parameters and measurement uncertainties have shown benefits in improving remote sensing capabilities. The Generalized Retrieval of Aerosol and Surface Properties (GRASP) algorithm retrieves multiple pixels simultaneously while considering the spatial correlation of the retrieval parameters (Dubovik et al., 2014, 2021). Xu et al. (2019) developed a correlated multipixel inversion approach (CIMAP), which further considers the correlation between different retrieval parameters. Theys et al. (2021) developed a CovarianceBased Retrieval Algorithm (COBRA) based on an error covariance matrix estimated from measurements with spectral correlation, applied their approach to sulfur dioxide (SO2) retrievals from the TROPOspheric Monitoring Instrument (TROPOMI) data, and demonstrated improved retrieval performance. To accurately evaluate pixellevel uncertainty in ocean color retrievals, spectral correlations associated with the uncertainties in topofatmosphere reflectance are also accounted for OLCI (Lamquin et al., 2013) and MODIS (Zhang et al., 2022) in the uncertainty propagation.
In this study, we provide a practical framework to understand the measurement uncertainty structure, study the impact of correlation in MAP retrievals, and demonstrate the potential for improvement in geophysical retrieval performance when proper correlation information is incorporated into the retrieval algorithm. Angular uncertainty correlations in measurements from the AirHARP and HARP2 instruments are studied as examples. Both instruments measure 60 angles at 670 nm. AirHARP measures 20 angles for the 440, 550, and 870 nm bands, while HARP2 measures at 10 angles for these bands. Angular correlation within each band is considered and modeled separately. Two methods are used to evaluate the retrieval uncertainties under different correlation strengths: (1) the error propagation method is used to evaluate the optimal retrieval uncertainties, by mapping the input uncertainty model describing the total uncertainty of the measurement and forward model to the retrieval parameter domain, and (2) comparative analyses are performed between the retrieval results from synthetic MAP measurements and the “truth data” that were assumed in the generation of that synthetic MAP data. The Monte Carlo error propagation method (MCEP) is adopted to compare the retrieval uncertainties from these two methods. To efficiently conduct retrieval and uncertainty analysis, the FastMAPOL retrieval algorithm is employed in this study, which uses neuralnetwork forward models for coupled atmosphere and ocean systems (Gao et al., 2021a). Analytical Jacobian matrices are derived based on the neural network and used to improve the efficiency of the retrieval (Gao et al., 2021b) and uncertainty quantification (Gao et al., 2022). To accurately evaluate the retrieval uncertainties of real measurements, an adaptive data screening approach is employed (Gao et al., 2021b). This ensures that only those measurements that can be sufficiently described by the forward model are used in this study, by avoiding uncharacterized uncertainty contributions due to contamination by cirrus clouds and other anomalies.
Furthermore, we study the angular uncertainty correlation in the measurements and demonstrate that the correlation property can be derived using the autocorrelation function from the retrieval fitting residuals. Studies on both synthetic data with various correlation strengths are conducted with results applied to the real measurement retrievals from AirHARP over multiple ocean scenes. Useful tools are provided to understand and analyze the angular correlated uncertainty structure and models. Note that autocorrelation analysis based on fitting residuals has been found useful in analyzing performance of machine learning algorithms such as using the Durbin–Watson test (Chatterjee and Simonoff, 2012).
In the following sections, we will discuss how to conveniently include an angular correlated uncertainty model in the retrieval algorithm (Sect. 2), evaluate the impact of correlated measurement uncertainty in retrieval uncertainties (Sect. 3), and estimate correlation strength by fitting residual analysis from both synthetic AirHARP and HARP2 data and real AirHARP measurements (Sect. 4). Discussion and conclusions are provided in Sect. 5. Although this study focuses on angular noise correlation, the conclusions on the impacts of correlations are also applicable to other instruments such as the hyperspectral measurements from both the SPEXone and the Ocean Color Instrument (OCI) that will be carried on NASA’s upcoming Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission (Werdell et al., 2019).
2.1 FastMAPOL retrieval algorithm
In this study, the FastMAPOL algorithm is used to retrieve aerosol and ocean optical properties from HARP measurements. The algorithm includes three main components: (1) a set of neuralnetworkbased radiative transfer forward models of the coupled atmosphere and ocean system (Gao et al., 2021a) and the corresponding analytical Jacobian matrix based on automatic differentiations derived from these neural networks (Gao et al., 2021b), (2) a multiangle cloud masking and data screening module (Gao et al., 2021b), and (3) an efficient uncertainty quantification component (Gao et al., 2022). Waterleaving signals in terms of remote sensing reflectance are derived with an additional neural network trained for the atmospheric correction process (Gao et al., 2021a).
The neuralnetwork forward models are trained for both reflectance and degree of linear polarization (DoLP) based on simulations from the successive orders of scattering radiative transfer model (RTSOS) developed by Zhai et al. (2009, 2010) and Zhai and Hu (2022). The atmosphere and ocean system is assumed to be a fourlayer system, and the radiative transfer interactions among them are fully considered in the RTSOS model. The bottom layer is the ocean water body in which an openocean biooptical model is used to parameterize scattering and absorption of ocean constituents based on the chlorophyll a concentration (Chl a; mg m^{−3}) (Gao et al., 2019, 2021a). The second layer is the ocean surface with its roughness parameterized by wind speed through a scalar Cox–Munk model (Cox and Munk, 1954). The third layer is an aerosol layer mixed with Rayleigh scattering. This layer extends from ocean surface to a height of 2 km with a uniform aerosol vertical distribution. The last layer contains atmospheric molecules from 2 km to the top of the atmosphere. The US standard atmospheric constituent profile is used to describe the molecular distributions (Anderson et al., 1986).
A total of nine parameters are used to describe the aerosol microphysical properties. There are four parameters for the complex refractive index of fine and coarse mode. Aerosol size distributions are parameterized by five volume densities for five size submodes with fixed effective radius and variance (Dubovik et al., 2006; Xu et al., 2016; Gao et al., 2018). Absorption by atmospheric gases is considered in the RTSOS simulation, with ozone density as the only variable. The radiant path geometries are represented by the solar and viewing zenith angles and the viewing azimuth angle relative to the solar direction. Therefore, a total of 15 parameters are used as forwardmodel input, with 11 of them defined as retrievable parameters. Details of the parameter ranges are listed in Appendix A and discussed in Gao et al. (2021a, 2022). To represent the forward model accurately and efficiently, the neuralnetwork (NN) architecture is optimized with an input layer of 15 parameters, followed by three hidden layers with 1024, 256, and 128 nodes and a final output layer with 4 nodes for each HARP band. The deeplearning Python library PyTorch is used for the training the NN (Paszke et al., 2019). The accuracy of the NN forward model is examined with an independent synthetic measurement dataset not used in training. An accuracy of better than 1 % for reflectance and better than 0.003 for DoLP has been achieved (Gao et al., 2021a). The uncertainties of the NN forward model are less than the instrument uncertainties of AirHARP and HARP2 (3 % in reflectance, 0.01 in DoLP for AirHARP, and 0.005 in DoLP for HARP2).
In this study, we will discuss the retrieval uncertainty and performance in aerosol properties, ocean surface wind speed, and Chl a in the ocean, as well as waterleaving signals based on the retrieval parameters. The waterleaving signal refers to the remote sensing reflectance (R_{rs}), which is the ratio of the upwelling waterleaving radiance and the downwelling solar irradiance just above the ocean surface (Mobley, 2022). R_{rs} can be estimated through the atmospheric correction process, which removes the contribution from the atmosphere and ocean surface from the total measurements at the sensor, and the additional BRDF (bidirectional reflectance distribution function) correction to reduce the dependency on the solar and viewing directions. Both atmospheric and BRDF corrections with their associated uncertainties are implemented using neural networks as discussed in Gao et al. (2021a, b) and followed by this study.
2.2 Retrieval cost function and uncertainty quantification
The maximum likelihood approach is used to retrieve the state parameters in FastMAPOL by minimizing a cost function that represents the difference between the measurements and the forwardmodel fitting (Rodgers, 2000)
where $\mathit{y}=\mathit{m}\mathit{f}\mathbf{\left(}\mathit{x}\mathbf{\right)}$ is the residual vector between measurement m and forward model f under retrieval parameters of x. Measurement vector m includes both reflectance (ρ_{t}) and DoLP (P_{t}), where the subscript t indicates the total signal measured by the instrument. The total number of measurements, N, at each pixel includes contributions from both reflectance and DoLP, which has been used in previous studies (Gao et al., 2021a, b).
The error covariance matrix S_{ϵ} in Eq. (1) specifies the uncertainties of each measurement and the correlation between different measurements at the same pixel, which is a symmetric matrix defined as
where i and j indicate the measurement at different angles and bands, and 𝔼 indicate the expectation values. To capture the angular uncertainty correlation, the autoregressive model of order of 1 (denoted as AR(1)) is used in the study of RSP data (Knobelspiesse et al., 2012) and adopted in this study for HARP data. AR(1) represents a linear Markov process with the error covariance matrix specified as
where σ_{t} is the total uncertainty, which includes both random noise and calibration uncertainty (σ_{c}). Only σ_{c} is assumed to be correlated between measurements at different viewing angles.
The ratios between random and calibration uncertainties may be different for reflectance and polarized signals (Knobelspiesse et al., 2019). The synthetic data are generated directly using the forward model; therefore, the contribution of forwardmodeling uncertainty is not considered for the synthetic data study. Δ_{θ} is the average angular grid size, which depends on the spectral bands. We model the correlation properties using the Δ_{θ} estimated from the viewing angles in the alongtrack direction to better represent the stripe filter characteristics used to conduct HARP angular measurement. The averaged Δ_{θ} is approximately 6.0^{∘} for AirHARP and 12^{∘} for HARP2 at 440, 550, and 870 nm bands and 2.0^{∘} for the 670 nm bands for both HARP instruments. r in Eq. (3) is the correlation parameter with a value between 0 and 1. For uncertainties with more complex structures, a general autoregression and moving average (ARMA) model can be used (Priestley, 1983; Brockwell and Davis, 1991). However, from our analysis based on the retrieval results from real AirHARP measurements, AR(1) works well for most cases. Detailed analysis can be found in Sect. 2.5 for theoretical basis, Sect. 4.3 for real data applications, and Sect. 5 for general discussions.
To better represent stronger correlations when it is close to 1, we define the correlation angles θ_{c} based on the correlation parameter r as
Therefore, θ_{c} indicates the angular range where magnitude in the correlation between angles is reduced by a factor of e. Similarly, correlation angles can be derived from r as
2.3 Uncertainty quantification
The pixelwise retrieval uncertainty can be quantified by mapping the measurement and forwardmodel uncertainties into retrieval parameter space (Rodgers, 2000):
where the Jacobian matrix K represents the partial derivatives of the measurements with respect to all the retrieval parameters. In this study, each retrieval parameter can only vary in a limited range, which imposes an implicit a priori constraint on the retrieval parameters. To capture its influence on retrieval uncertainties, we assume the a priori error matrix S_{a} in Eq. (6) to be diagonal with the a prior uncertainty for each state parameter approximated by its permitted range in retrievals (Gao et al., 2022). Both the parameter ranges and a priori values are listed in Table A1. The uncertainties are defined as the standard deviation (1σ) around the retrieval solution, which is estimated by the square roots of the diagonal elements of S. The uncertainties of variables which are a function of the retrieval parameters can also be derived from S and their derivatives. Due to the large number of retrieval parameters used in the retrieval, the evaluation of the retrieval uncertainties can be time consuming. The speed to compute uncertainties is improved using automatic differentiations based on neuralnetwork forward models (Gao et al., 2022). For example, the uncertainty of remote sensing reflectance (R_{rs}) can be derived using the automatic differentiation applied on the neural networks for BRDF correction and atmospheric correction components as discussed in Appendix A from Gao et al. (2021b).
The retrieval uncertainties estimated by error propagation (hereafter called theoretical retrieval uncertainty) as shown in Eq. (6) represent the optimal scenarios, with limitations such as the assumption that the retrieval parameters successfully converged to the global minima (more discussions in Sayer et al., 2020; Gao et al., 2022). When the retrievals converge to a local minimum, both the retrieval results and associated Jacobians can be less representative of the truth values and therefore lead to inaccurate error propagation and uncertainty estimation. To quantify the retrieval uncertainties based on actual retrieval results, the retrieval errors are defined as the difference between the retrieval results and the truth from synthetic data, which are then used to compute the retrieval uncertainty (hereafter called real retrieval uncertainty). The comparison between the theoretical and real uncertainties is useful to assess the optimal and actual performance of a retrieval algorithm. The Monte Carlo error propagation (MCEP) method is used in this study to conduct such comparison (Gao et al., 2022). MCEP samples the retrieval errors from theoretical retrieval uncertainties and then directly compares the error distributions between theoretical estimation and real retrievals. This method provides additional flexibility in analyzing their statistics. Multiple sets of random samples are generated from the theoretical uncertainties with their variations analyzed, which therefore provides a way to evaluate the impact of sample size in estimating uncertainties (Gao et al., 2022). This method is used to quantify the retrieval uncertainties with various correlation strengths in the next section. A summary of the terminology used for the error and uncertainties for measurement and retrieval results is provided in Table 2.
2.4 Eigenvector decomposition on error covariance matrix
The error covariance matrix with nondiagonal terms is challenging to implement efficiently in optimization algorithms, which typically operate in diagonal space with no correlation between measurements. The error covariance matrix also creates barriers to understand the retrieval uncertainties, as the input uncertainties are not for a single measurement but rather related to multiple measurements. To overcome these issues, we convert the measurements into a new space where the error covariance matrix is diagonalized. Therefore, conventional optimization and error analysis techniques can be readily used.
To achieve this goal, eigenvector decomposition is applied on the error covariance matrix (Rodgers, 2000) as
where D_{ϵ} is a positive diagonal matrix defined by the eigenvalues of S_{ϵ} and U is a unitary matrix. Examples of how the eigenvalues vary with correlation strength and the associated impact of those variations on Shannon information content and retrieval uncertainties are discussed in Appendix B. Based on Eq. (7), the cost function in Eq. (1) and the error propagation in Eq. (6) can be written as
where the original set of measurements, y, are converted to a new set of measurements, y^{′}, without any correlation, along with the respective Jacobian matrices:
Equations (8) and (9) are mathematical transformations that conveniently allow for working with diagonal matrices, with advantages and applications summarized below.

A clear expression of measurement uncertainty. The diagonal terms in matrix D_{ϵ} represent the uncertainties for the new measurement y^{′}, therefore providing insights on the accuracy of the measurements impacted by correlation. An example of this is shown in Appendix B.

Conducting minimization on the retrieval cost function. Equation (8) represents the cost function for noncorrelated measurement y^{′} and therefore can be used for conventional optimization algorithms such as the subspace trustregion interior reflective (STIR) algorithm (Branch et al., 1999) as used in current FastMAPOL algorithm (Gao et al., 2021a, b, 2022).

Generating correlated errors. To study and visualize angular uncertainty correlations, correlated errors need to be generated and then added to the synthetic data. To achieve this goal, we generated the errors in the space for y^{′} with random parameters sampled from a normal distribution assuming the eigenvalues in D_{ϵ} as its variance. These errors in y^{′} are then transformed back to the original space y through $\mathit{y}={\mathbf{U}}^{T}{\mathit{y}}^{\mathbf{\prime}}$. To demonstrate how angular correlations impact the errors, the correlated error samples with a correlation angle of ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{10}{}^{\circ}$ (r=0.9) and correlation angle of ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$ (r=0.98) are shown in Fig. 2a and c. A value of r=0.9 has been assumed in the study of RSP angular correlation by Knobelspiesse et al. (2012). With a larger θ_{c} the errors start to form a longer range of correlation with smoother variations. Note that the overall magnitude of the errors can vary within the full range as described by the calibration uncertainties. These errors are then added to the synthetic data and used to study retrieval. The fitting residuals from retrievals on the synthetic data with the added correlated errors are shown in Fig. 2c and d and discussed in the next section.
2.5 Correlation strength estimation using autocorrelation
Autocorrelation is a useful function to quantify correlation in a discrete data sequence and is defined as (Priestley, 1983)
where i and j are two indices of the datasets. Comparing with Eq. (2), the autocorrelation is equivalent to the autocovariance when 𝔼[y_{i}]=0. This method can be applied to the simulated noise generated in Sect. 2.4 and used to analyze the fitting residuals. However, the mean values and variance in the fitting residuals often vary with respect to the angular grids. This type of signal is classified as nonstationary and difficult to study by the AR models (Priestley, 1983). To overcome this issue, the original residual data y are processed by removing their mean and normalizing by their standard deviation. These normalized data are denoted as $\stackrel{\mathrm{\u0303}}{\mathit{y}}$. For the data within the same band and polarization state, the autocorrelation function on the normalized data is equal to their covariance as defined in Eq. (3),
We can estimate the correlation by analyzing the residuals between the measurement and forward model. The autocorrelation function is averaged over multiple pixels to reduce uncertainties for the analysis in both synthetic data and real retrieval residuals. The correlation parameter can then be derived as
Correlation angles, θ_{c}, are then computed based on r following the formula in Eq. (5). Furthermore, a partial autocorrelation function from a sequence of data can be computed, which removes the correlation due to lags higher than 1 (Priestley, 1983). If the AR(1) model is sufficient to describe the noise structure, only one additional term would be left besides the zeroorder term in the partial autocorrelation results. Therefore, partial autocorrelation can be used to validate our assumption in the uncertainty model. The Python packages, statsmodels (Seabold and Perktold, 2010) and SciPy (Virtanen et al., 2018), are used to conduct the autocorrelation analysis.
An example is shown in Fig. 3; the autocorrelation function and partial autocorrelation function are applied on the simulated errors and the retrieval residuals with examples from Fig. 2b and d. The autocorrelations are shown in Fig. 3a and c for the simulated errors and retrieval residuals in reflectance data. The partial autocorrelations for Fig. 3a and c are shown in Fig. 3b and d, respectively. For both cases in Fig. 3b and d, only the first order of data are prominent, which confirms that the data can be represented by the AR(1) process. If higher orders in the AR process are presented, more prominent data points will appear in Fig. 3b and d. The estimated correlation angles for the errors in Fig. 3a and b and residuals in Fig. 3c and d are approximately 30 and 15^{∘}, respectively, after converting correlation parameters to correlation angles but less than the actual correlation angle of 60^{∘}. The results show that autocorrelation can be a useful way to estimate correlation strength but with a tendency to underestimate due to the finite length of the data and overfitting of the retrievals (more discussion in the Sect. 4). Note that the Bartlett confidence interval that corresponds the autocorrelation of uncorrelated white noise can be calculated but is found to be much smaller than the results in Fig. 3 and is therefore not shown in the figure (Brockwell and Davis, 1991).
3.1 Synthetic data generated using the NN forward model
The neuralnetwork forward model discussed in Sect. 2.1 is used to generate 1000 sets of synthetic AirHARP data, and then the number of viewing angles at 440, 550 and 870 nm are downsampled to 10 to represent HARP2 data. A fixed solar zenith angle of 50^{∘} is used to represent the solar geometries of the AirHARP scenes over ocean from the ACEPOL field campaign (more information in Sect. 4.3). The aerosol properties, wind speed, and Chl a values are randomly sampled based on their allowed range, as discussed in Sect. 2 and Appendix A. The same sampling approach discussed in Gao et al. (2022) is conducted assuming that the aerosol optical depth (AOD) and finemode volume fraction are uniformly distributed within [0.01,0.5] and [0,1], respectively. A larger range of AOD values will be needed for applying this algorithm to cases of smoke and plume events. Realistic HARPlike viewing geometries are constructed by sampling the alongtrack and crosstrack viewing angles randomly and then converting to the actual viewing zenith and azimuth angles following the formulas provided in Gao et al. (2021b). Example viewing geometries for AirHARP and HARP2 for bands at 550 nm are provided in Fig. 4, with geometries at other bands constructed similarly.
To generate realistic measurements, correlated uncertainties with a correlation angle (θ_{c}) from 0, 1, 2, 5, 10, 20, 30, 60, and 120^{∘} are considered in this study, with corresponding correlation parameters of 0, 0.368, 0.607, 0.819, 0.905, 0.951, 0.967, 0.983, and 0.992, respectively. The corresponding correlated error samples are generated based on the error covariance matrix using the method discussed in Sect. 2.4. Examples of the correlated errors are already shown in Fig. 2 for AirHARP at 660 nm, with correlation angles of 10 and 60^{∘}, respectively.
Correlated errors for both the AirHARP and HARP2 instruments are generated according to the same 3 % uncertainty for reflectance but 0.01 in DoLP for AirHARP and 0.005 in DoLP for HARP2. These errors are added to the corresponding simulated reflectance and DoLP ($\sqrt{{Q}^{\mathrm{2}}+{U}^{\mathrm{2}}}/I$) for further studies. The reflectance is more likely to be dominated by systematic uncertainty like calibration (possibly correlated), while DoLP defined as the ratio between two measurements is more likely dominated by randomly generated uncertainty like shot noise (probably less correlated) (Knobelspiesse et al., 2012). Therefore, we assume two scenarios: (1) angular correlation only existed in reflectance measurement, not in DoLP measurement, and (2) both reflectance and DoLP have angular correlations with the same strength. Since the actual amount of correlation is not known, we designed our studies with the assumed correlation in the synthetic measurement but with either no information or full information on the correlation angle in the retrieval cost functions. Therefore, four scenarios are discussed in this study as summarized in Table 3 denoted by C1 to C4. We will discuss whether better retrieval results can be obtained if accurate correlation angles are considered in the retrieval cost function and whether we can estimate correlation from retrieval residuals in Sect. 4.
3.2 Retrieval uncertainties impacted by uncertainty correlation
Using the Monte Carlo error propagation (MCEP) method discussed in Gao et al. (2022), we compared both real and theoretical retrieval uncertainties with different correlation angle and testing scenarios as summarized in Table 3 for both synthetic AirHARP and HARP2 measurements. Figure 5 demonstrates the basic approach in MCEP by comparing the AirHARP retrieval results with ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$ for scenarios C3 and C4, where the real retrieval errors are sampled based on their retrieval uncertainties. The real uncertainties in both the root mean square error (RMSE) and the mean average error (MAE) are mostly larger when uncertainty is correlated (comparing b and a), with exceptions possibly due to statistical fluctuations in the Monte Carlo sampling. The theoretical uncertainties are similar because in both cases the correlation angles are assumed to be zero. After considering the same correlation angle in the retrieval cost function model as shown in Fig. 5c, both theoretical and real uncertainties are reduced. The real and theoretical uncertainties are similar to each other as shown in Fig. 5a, but the agreement degrades when correlation is considered (Fig. 5b, c). There are more real retrieval errors in the negative side as shown in Fig. 5c. This may be related to the convergence for cases with small AOD values as discussed in Gao et al. (2022). Although real uncertainties are generally larger than theoretical uncertainties, the differences are mostly associated with cases where AODs are small and underestimated relative to the truth values.
Both real errors and theoretical uncertainties have occasional outliers with large values possibly due to convergence to local minima instead of global minima, and this has large impacts on the RMSE values. For example, in Fig. 5b there are a few large estimated uncertainty values outside of the plot range that result in large RMSE for the theoretical error, almost comparable to the real uncertainty. Comparing MAE, however, shows that the theoretical value is much smaller than the real error, consistent with the histogram shape where the theoretical curve is narrower than the real curve. Given its robustness to outliers, MAE is used as the primary metric in this study.
Based on the MCEP method, we analyzed the retrieval uncertainties for synthetic AirHARP measurements from real errors and theoretical errors for various properties in Fig. 6, including the AOD, single scattering albedo (SSA), the real part of the refractive index (m_{r}), effective radius (r_{eff}), and effective variance (v_{eff}) for both fine and coarse modes and their combinations. fvf is the finemode volume fraction. Ocean surface wind speed and Chl a concentration are also retrieved. As recommended by Seegers et al. (2018), the uncertainty of Chl a is represented by a logtransformed metric:
where M is the total number of samples, which equals to the total number of synthetic measurement cases. The values of ${\mathrm{log}}_{\mathrm{10}}\left(\text{Chl}\phantom{\rule{0.25em}{0ex}}{a}_{i,\text{retrieval}}\right){\mathrm{log}}_{\mathrm{10}}\left(\text{Chl}\phantom{\rule{0.25em}{0ex}}{a}_{i,\text{truth}}\right)$ are sampled for both the real and theoretical uncertainties similar to Fig. 5, with detailed discussions provided in Gao et al. (2022). Furthermore, remote sensing reflectance at the four HARP wavelengths is derived after conducting the atmospheric correction using the retrieved aerosol and ocean properties (Gao et al., 2021a). Multiple sets of random errors are sampled and averaged to estimate the impact of sample size in estimating retrieval uncertainties and are shown in Fig. 6.
The real retrieval uncertainties for Scenario C3, in which correlation is considered in the simulated errors but not in the retrieval cost function, are found to be always increasing with the correlation angle. When the correct correlation angle is considered in the retrieval cost function (C4), the real retrieval uncertainty increases until θ_{c}, reaching 10 to 20^{∘}, and then slightly decreases at higher correlation angles, for most retrieval parameters. Similar behavior of the information content has been reported by Knobelspiesse et al. (2012) on the study of error correlation in RSP measurements. To understand how the correlation influences retrieval accuracy, we further analyze its impacts on the eigenvalues of the error covariance matrix with details discussed in Appendix B. The theoretical uncertainties for Scenario C4 are similar to the real uncertainties for ${\mathit{\theta}}_{\mathrm{c}}<\mathrm{10}{}^{\circ}$ except for refractive index, for which the theoretical values are almost half that of the real uncertainties. When ${\mathit{\theta}}_{\mathrm{c}}>\mathrm{10}{}^{\circ}$, the theoretical uncertainties then decrease much faster than the real uncertainties as θ_{c} increases. The difference mostly comes from the retrieval cases which underestimate the truth values as shown in Fig. 5. Note that the real retrieval uncertainties are mostly larger than the theoretical uncertainties without any correlation within the range of ${\mathit{\theta}}_{\mathrm{c}}<\mathrm{120}{}^{\circ}$. However, theoretical uncertainties predict that the retrieval uncertainties increase slightly with θ_{c} around 10^{∘} (or r=0.9) and then decrease. Its values can be smaller than those with zero correlation when θ_{c} is larger than around 20^{∘}.
The results for HARP2 are similar to that of AirHARP as shown in Fig. 7. The overall HARP2 retrieval uncertainties are slightly smaller than AirHARP retrieval uncertainties, but the difference is mostly within 20 % for the theoretical uncertainties and mostly within 10 % for the real uncertainties as shown in Fig. 7 for both ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$ and ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$. Although HARP2 measures fewer viewing angles at 440, 550, and 870 nm bands, its better DoLP accuracy still results in slightly smaller uncertainties in most cases. Note that the retrieval accuracies also depend on the total number of viewing angles used and the range of scattering angles as discussed in Gao et al. (2021b).
To understand more quantitatively how the correlation angle impacts both AirHARP and HARP2 retrievals, the ratios of the real uncertainties between scenarios C3 and C4 are represented as $R\left(\text{C4}\right)/R\left(\text{C3}\right)$ for each retrieval quantity as shown in Fig. 8. Note that both C3 and C4 considered uncertainty correlations in both reflectance and DoLP, but only C4 considered the same correlation in its retrieval cost function, and C3 assumes no correlation in its retrievals. Therefore, the ratio $R\left(\text{C4}\right)/R\left(\text{C3}\right)$ represents how much retrieval uncertainties can be reduced if the correct amount of correlations are known in the retrieval process. As shown in Fig. 8, the ratio is close to 1 for most parameters for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{10}{}^{\circ}$, with slightly larger impacts for effective variance for both fine and coarse mode as well as for wind speed uncertainty. Such ratios almost double for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$ for both AirHARP and HARP2 measurements. Note that sunglint has been removed in the retrieval as shown in Fig. 4, which may contribute to a larger wind speed retrieval uncertainty. When comparing the theoretical uncertainty T(C4) with R(C3), their ratio reduces to a value of 0.5 to 0.7 for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{10}{}^{\circ}$ for aerosol properties and further decreases to a value of 0.3 to 0.5 with ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$. The impacts on the remote sensing reflectance are generally smaller for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{10}{}^{\circ}$ but more significant for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$. This is because the theoretical uncertainty with a correlation angle decreases faster than the real uncertainties (see Fig. 6). These results demonstrate that there is potentially more room to reduce the retrieval uncertainty as predicted by the error propagation theory represented by Eq. (6), which warrants future investigations. The largest gap between the real and theoretical uncertainties is still for refractive index in both fine mode and coarse mode, which is the same as observed in Gao et al. (2022). The gaps between the real and theoretical uncertainties increase with increasing θ_{c}, which indicates degrading retrieval performance in the presence of correlated uncertainties.
To understand retrieval performance when the uncertainty correlation is only in reflectance for scenarios C1 and C2, similar ratios to Fig. 8 are plotted in Fig. 9. The ratios between the real uncertainties are almost always equal to a value of 1, with slightly larger impacts (less than 10 %) for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$. However, the theoretical errors are much smaller, with values ranging from approximately 0.75 for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{10}{}^{\circ}$ to 0.5 for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$ and with even smaller values observed for refractive index. This suggest that correlation in reflectance alone has the potential to be harvested to improve retrieval performance but is even harder to realize in real retrievals from current algorithms.
4.1 Cost function and fitting performance
As discussed in Gao et al. (2021b), the retrieval cost function can be characterized well by the χ^{2} distribution for synthetic data and for real data after removing anomalies such as cirrus cloud. Note that the cost function in Eq. (1) is called a χ^{2} function, but its histogram may not always follow a χ^{2} statistical distribution, which depends on how well the fitting residuals can represent the real uncertainty and the degree of freedom of the χ^{2} distribution (Gao et al., 2021a). Therefore, the histogram of successful retrievals is a useful indicator on how well the retrieval residuals compare with the assumed input uncertainty model (Rodgers, 2000). In this study, we found that the retrieval residuals can be represented well by the χ^{2} distribution with a degree of freedom equal to the total number of measurements used (N), which is twice that of the total number of viewing angles (N_{v}) when there is no correlation, as shown in Fig. 10. Note that the degree of freedom is defined for the χ^{2} distribution (Gao et al., 2021b).
However, the cost function histogram shifts to a smaller value for Scenario C3 (Fig. 10a), where correlated uncertainty is included in the synthetic data but no correlation is considered in the retrieval cost function. Smaller cost function values indicate smaller retrieval residual, which may be caused by overfitting of the data and also possibly lead to the larger real retrieval uncertainties as shown in Fig. 6. When the correct value of correlation is considered in the retrieval cost function as shown in Fig. 10b for Scenario C4, the histogram better approaches a χ^{2} distribution. For example, in the right panel in Fig. 10, the χ^{2} distributions with degrees of freedom of 40 and 20 is found to better fit the cost function histogram with ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{10}{}^{\circ}$ and ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$, which compares to results using all the measurement degrees of freedom (150). This suggests that the correlation in the uncertainty reduces its degree of freedom.
To understand how much overfitting impacts retrievals under a different strength of correlations, we compare the standard deviation of the retrieval residuals with the original simulated uncertainty for all 1000 cases. For reflectance, we consider the ratio between the simulated uncertainty and the reflectance as a convenient way to compare with the 3 % (or 0.03) uncertainty model for reflectance (Fig. 11). For the simulated uncertainty, the value of 0.03 is confirmed; however, the retrieval residuals become smaller with increasing θ_{c} for C3. The 670 nm band shows the largest reduction, where a ratio of 0.025 is observed for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{10}{}^{\circ}$, reaching 0.015 for ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{120}{}^{\circ}$, which may be due to the large number of angles (60) in the 670 nm bands versus the other bands (10 or 20). This behavior indicates overfitting, where the uncertainties are partially removed as real signals and lead to reduced residuals. When the correct correlation is considered in the model, the result is much like the assumption, with a very slight indication of overfitting (within 0.005 in most cases). Similarly for DoLP, the standard deviations of the fitting residuals at the 670 nm band are smaller than other bands for C3, which indicate larger overfitting at the 670 nm band (Fig. 11c). The overfitting is much reduced for C4 with a standard deviation close to 0.01 (Fig. 11d). The results for HARP2 are very similar, but with the assumed uncertainty for DoLP of 0.005 (not shown).
4.2 Correlation estimation results with synthetic data
Different amount of over fitting also partially removed the correlation in the fitting residuals. As shown in Fig. 12a and b, the estimated correlation is smaller than the truth. The estimated θ_{c} for 670 nm bands are smaller than other bands, which is consistent with previous studies where these bands overfit the data the most. Green bands seem to show the best results, which estimate ${\mathit{\theta}}_{\mathrm{c}}<\mathrm{10}{}^{\circ}$ (r=0.90) well but underestimate true ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{60}{}^{\circ}$ (r=0.983) as 20^{∘} (r=0.951) and ${\mathit{\theta}}_{\mathrm{c}}=\mathrm{120}{}^{\circ}$ (r=0.992) as 30^{∘} (r=0.967). Note that the value of correlation parameter r asymptotically approaches 1 and becomes harder to distinguish with a finite length of data. Meanwhile, we also computed the correlation angles from the added simulated errors, which also resulted in a smaller correlation angle compared to the truth but much closer than using the residual data. This may due to the finite length of the measurement, which is 90^{∘} after the glint is removed (total 120^{∘}). When the true correlation angle is considered in the model (C4), the estimated correlation angle is improved. The results of DoLP are slightly better where the retrieval residuals seem to estimate the results from simulated errors in a similar way for both C3 and C4 as shown in Fig. 12c and d. In real data, the correlation strength may be different for reflectance and DoLP; it would be useful to estimate the correlation angles for all the four bands and both reflectance and DoLP and analyze the difference by comparing it with the synthetic data study.
4.3 Correlation estimation results with real AirHARP data
The Aerosol Characterization from Polarimeter and Lidar (ACEPOL) field campaign was conducted from October to November of 2017 with the NASA's ER2 aircraft at a high altitude of approximately 20 km (Knobelspiesse et al., 2020). Measurements over a variety of scenes are conducted from four MAPs – AirHARP, AirMSPI, SPEX airborne, and RSP – and two lidar sensors – HSRL2 (Burton et al., 2015) and CPL (the Cloud Physics Lidar) (McGill et al., 2002). There is a total of five AirHARP ocean scenes available in the ACEPOL measurements with retrieval uncertainties studied by Gao et al. (2022) without considering angular uncertainty correlation. Gao et al. (2021a) reported that the retrieval results of both aerosol and ocean color signals have been found to be in good agreement with the AERONET Ocean Color site (Zibordi et al., 2009). However, the cost function histogram was much wider than expected due to the impacts of cirrus clouds. After removing the cirrus cloud impacts from the multipleangle measurement using an adaptive data screening method, the cost function histogram improved significantly, with much higher similarity with a χ^{2} distribution (Gao et al., 2021b). In this study, the same adaptive data screening methods are applied on all the five AirHARP ocean scenes, which removes cirrus clouds and other anomalies that could not be represented adequately by the current forward model. The resulting total number of measurements, including both reflectance and DoLP, are shown in Fig. 13, where the spatial distribution of pixels with many valid viewing angles is not uniform.
Correlation properties from real measurements are difficult to quantity. However, we show in Sect. 2.5 that, as also demonstrated by Sect. 4.2, the autocorrelation analysis on the retrieval residual can be used as a good estimator when correlation is not too strong with consistent behaviors to derive the correlation parameters. Therefore, using the ACEPOL data, we can estimate angular uncertainty correlation from retrieval residuals. From each scene we selected 200 pixels with the most available number of angles, which are clustered together around the region with the maximum number of measurements as shown in Fig. 14. The retrieval residual data were then normalized as discussed in Sect. 2.5 to remove impacts by the nonuniform mean and variance, which are often observed from the real data residuals. The autocorrelation and partial autocorrelation are calculated to assess whether the AR(1) model is sufficient, with examples shown in Fig. 14, for the fitting residuals of reflectance and DoLP from Scene 3 at the 670 nm band. Partial autocorrelation for reflectance showed similar results for the synthetic data in Fig. 3b, with only the firstorder term prominent, which suggests that the AR(1) model is sufficient to describe the fitting residual for reflectance, with higherorder contributions negligible. Figure 3d for DoLP shows that higherorder terms may also contribute to the uncertainty model; however, the overall correlation strength is small. From these plots, we can estimate the correlation parameters r following Eq. (14), where ${\stackrel{\mathrm{\u0303}}{R}}_{\mathrm{1}}$ corresponds to the firstorder point in the autocorrelation plots. The values of r are approximately 0.9 and 0.7 for reflectance and DoLP, respectively, corresponding to correlation angles of approximately 10 and 15^{∘} following Eq. (5).
We analyzed the fitting residuals for all the five AirHARP scenes at the four bands with results summarized in Fig. 15. The minimal number of available angles used in the correlation estimation at each band are shown in Fig. 15a. Note that different scenes have various numbers of angles removed due to the impacts of cirrus cloud or other anomalies. Those values are filled with zeros, which may reduce the strength of angular correlation. Therefore, the real correlation angle is likely larger than these values because only a subset of the total measurement data are used to estimate these angles. The estimated correlation angles for reflectance and DoLP vary mostly between 5 and 20^{∘}. The correlation angles are smaller for the 670 nm band, probably due to overfitting that partially removes the correlated errors as real retrieval signals, consistent with our observations based on synthetic data (Fig. 12). The correlation angles for the 440 and 550 nm bands are largest, which may be mostly close to the truth. The correlation angles are generally larger for reflectance, with a value between 10 and 20^{∘} compared to DoLP with a value between 5 and 15^{∘}, which suggests different amounts of correlation in the reflectance and DoLP data. Comparing the retrieval uncertainties for θ_{c} around 10 to 20^{∘} as shown in Fig. 6, the impact of the correlation to the real retrieval is small, but there is potential for large reduction of theoretical uncertainties by 25 % to 50 % if the correct correlation is considered in the retrieval cost function. This requires, however, that the retrieval algorithm is capable of achieving its optimal performance as described by its theoretical uncertainties. Furthermore, the information on the correlation properties is useful to parameterize realistic measurement uncertainties into synthetic data. The correlation angles for reflectance and DoLP are likely larger than 10 and 5^{∘}, respectively, which correspond to an estimated correlation parameter of 0.9 and 0.8.
In this study, we evaluate the impacts of angular correlation on the retrieval uncertainties for various aerosol microphysical and optical properties, ocean surface properties, and waterleaving signals. Theoretical uncertainties are derived based on error propagation, and the real uncertainties are obtained through the comparison of retrieved and true values. The theoretical and real uncertainties are compared and discussed. Only small angular correlation impacts are found on the real retrieval uncertainties unless the correlation strength is large (such as with a correlation angle larger than ${\mathit{\theta}}_{\mathrm{c}}>\mathrm{10}{}^{\circ}$). The impacts vary with different retrieval parameters. Theoretical uncertainties are more impacted by angular correlation, which suggests the retrieval algorithm may not always converge to the global minima and that there is potential room for algorithm improvement.
Studies on the fitting residuals from both synthetic data and real AirHARP measurements were conducted. Autocorrelation is useful to estimate the angular correlation, though it tends to underestimate when the correlation strength is strong, and thus overfitting of the measurements is likely. Analysis on the real data showed that the angular correlation is stronger in the reflectance data than DoLP, which makes sense because we expect that DoLP is less sensitive to systematic uncertainties that are more likely to be correlated. Partial autocorrelation analysis suggests that the uncertainty model considering a linear Markov process (AR(1)) is sufficient for reflectance but may need to be further studied for DoLP. From AirHARP retrieval residual analysis, the correlation angles for reflectance and DoLP are estimated to be larger than 10 and 5^{∘}, corresponding to correlation parameters larger than 0.9 and 0.8, respectively.
This work intends to provide basic methodology to analyze the measurement uncertainties with angular correlations, but the methods can also be applied in the spatial and spectral domains that may be more appropriate for other instruments. There are several remaining issues that need discussion in future works.

Application to real data. It is complex to analyze real data as we discussed in previous sections. The major challenges and possible issues that may impact uncertainty correlation estimation are summarized below.
 1.
The retrieval is based on a forward model which also has uncertainties, a portion of which may be correlated. This uncertainty will contribute to the fitting residuals and may impact correlation analysis, but it is difficult to quantify.
 2.
The fitting residuals are often not stationary with uniform mean and variance. To reduce this issue, the residuals are normalized, but it would be valuable to analyze how the mean value and variance depend on the angle, as this may provide insight into the modeling uncertainties.
 3.
Some residuals are not continuous with angle due to removed cirrus clouds, which may reduce the correlation.
 4.
Synthetic data analysis has demonstrated that the retrieval is likely to overfit the data when the correlation is strong.
 5.
The angular grids for HARP measurements are slightly nonuniform, which is likely to further reduce the correlation strength from autocorrelation analysis. To evaluate impacts of this feature, an uncertainty model considering the impact of the real angular grids needs to be built. But the variation of the angular grids is less than 1^{∘} (670 nm band) or 2^{∘} (other bands), which may impact more the cases with small correlation angles.
 1.

Lab calibration. Although the correlation strength is estimated from fitting residuals from real AirHARP measurements, it may be only used for qualitative discussions due to various issues discussed above. To obtain the actual correlation properties, lab characterizations are desired to separate measurement characteristics. Lab measurement signals may also be evaluated through autocorrelation and partial correlation functions, or more general ARMA models. It is also interesting to discuss the possible impacts in the uncertainty model due to the binning and collocation that happen in later processing steps.

Correlation strength as a fitting metric. Due to the limitation discussed above, our analysis on the fitting residuals may only provide the lower boundary for the correlation strength in terms of correlation parameter (r) or correlation angle (θ_{c}). Since the correlation strength represents properties of the fitting residuals, it can be also used as a metric to represent retrieval fitting performance, together with the cost function χ^{2} and variance of the fitting residuals as discussed in Sect. 4.1.

Signal correlation vs. uncertainty correlation. Nature as measured by the instrument and expressed in the forward model has inherent correlation, which becomes part of the retrieval process. The phenomena we observe tend to be only slowly changing with respect to view angle, and thus measurements at different view angles do not necessarily express retrieved parameters independently. This correlation is related to the actual signal in the measurement rather than its uncertainties. This type of correlation is captured by the Jacobian matrix. The overall information content of the measurements with respect to the set of retrieval parameters is determined by both the Jacobian matrix and the correlated uncertainty model as further discussed in Appendix B.

Future retrieval algorithm development. The pixelwise theoretical uncertainties achieve a reasonably good performance to represent real retrievals when no correlation is presented. Their performances on various retrieved geophysical properties are quantified by comparing with the real retrieval errors. The difference grows bigger when the angular correlation is stronger, which suggests convergence to local minima and indicates that more development is needed to improve the retrieval optimization. For the development of future algorithms with more retrieval parameters, such as aerosols with more complex shape and absorption properties and coastal waters with more complex biooptical properties, a better characterized error model, such as the one considering angular or spectral correlations, will be helpful to identify information useful for the retrievals and therefore improve retrieval performance and uncertainty assessment.
A total of 15 parameters are used as input of the forward model as discussed in Sect. 2.1 as listed in Table A1. The solar zenith (θ_{0}), viewing zenith angle (θ_{v}), viewing azimuth angle relative to solar direction (ϕ_{v}), and ozone density (${n}_{{\mathrm{O}}_{\mathrm{3}}}$) are assumed to be known input. All the other 11 parameters are retrieval parameters in the FastMAPOL algorithm, including the aerosol volume density for each submode (V_{i}), the real (m_{r}) and imaginary parts (m_{i}) of the refractive index for fine and coarsemode aerosols, ocean surface wind speed (w), and chlorophyll a concentration (Chl a).
Section 2.4 discussed that the error covariance matrix with nondiagonal terms can be diagonalized through eigendecomposition as shown in Eq. (7). The original measurements can be transformed into a new space without correlations, with uncertainty variance described by the eigenvalues in the diagonal matrix D_{ϵ} denoted as ${d}_{i}^{\mathrm{2}}$. To understand how the uncertainties vary with different correlation strength, the square root of the diagonal term in D_{ϵ} is used to represent the new measurement uncertainties. Results for different correlation angles are shown in Fig. B1.
The uncertainties in the new measurement space (y^{′}) show components both larger and smaller than the original uncertainties from different combinations of measurements. The corresponding Shannon information content (SIC) (Rodgers, 2000) is defined as
where the error covariance matrix S and the a priori matrix S_{a} in Fig. B2 are from Eq. (6). The Jacobian matrix and the error covariance matrix are converted into the diagonal space as shown in Sect. 2.4 and are used to represent the SIC in Eq. (B1). Different correlation strength will lead to a different unitary matrix, which is used to transform the Jacobian and error covariance matrix (Eq. 7). We analyzed the SIC and retrieval uncertainties for a sequence of correlation angles from 0 to 120^{∘}. As shown in Fig. B2, when the correlation strength is strong, the SIC is increasing with respect to correlation, and the theoretical retrieval uncertainties estimated from error propagation are decreasing, such as for finemode refractive index, wind speed, and Chl a. However, when the correlation is relatively weak, SIC decreases and uncertainties increase with the correlations. The correlation angles with maximum uncertainties are different with different retrieval parameters but generally fall within 30^{∘}. This behavior may relate how the measurement uncertainties are mapped to the retrieval parameters space through the Jacobian matrix. Similar behavior of the SIC has been reported by Knobelspiesse et al. (2012) for the RSP measurement. Understanding of the uncertainty properties such as correlation strength is useful to further exploit the information in the measurements.
The AirHARP data used in this study are available from the ACEPOL data portal (https://doi.org/10.5067/SUBORBITAL/ACEPOL2017/DATA001; ACEPOL Science Team, 2017).
MG, KK, BAF, PWZ, and BC formulated the study concept. MG generated the scientific data and wrote the original manuscript. PWZ developed the radiative transfer code used to train the NN models. XX and JVM provided and advised on the HARP data. All authors provided critical feedback and edited the manuscript.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors would like to thank the ACEPOL teams for conducting the field campaign. The numerical studies are conducted on the Poseidon supercomputer cluster at the NASA Ocean Biology Processing Group (OBPG). We thank the OBPG system team for supporting the highperformance computing. We thank Yunwei Cui, Zhonghuan Chen, Amir Ibrahim, Can Li, Xu Liu, Andy Sayer, and Jason Xuan for constructive discussions.
Meng Gao, Kirk Knobelspiesse, Bryan A. Franz, and Brian Cairns are supported by the NASA PACE project. PengWang Zhai is supported by NASA (grant no. 80NSSC20M0227). The ACEPOL campaign has been supported by the NASA Radiation Sciences Program, with funding from NASA (ACE and CALIPSO missions) and SRON (NWO/NSO project no. ALWGO/1609).
This paper was edited by Jian Xu and reviewed by two anonymous referees.
ACEPOL Science Team: Aerosol Characterization from Polarimeter and Lidar Campaign, NASA Langley Atmospheric Science Data Center (DAAC) [data set], https://doi.org/10.5067/SUBORBITAL/ACEPOL2017/DATA001, 2017. a
Anderson, G., Clough, S., Kneizys, F., Chetwynd, J., and Shettle, E.: AFGL Atmospheric Constituent Profiles (0–120 km), Air Force Geophysics Lab., Hanscom AFB, MA, USA, AFGLTR860110, 1986. a
Branch, M. A., Coleman, T. F., and Li, Y.: A Subspace, Interior, and Conjugate Gradient Method for LargeScale BoundConstrained Minimization Problems, SIAM J. Sci. Comput., 21, 1–23, https://doi.org/10.1137/S1064827595289108, 1999. a
Brockwell, P. J. and Davis, R. A.: Time Series: Theory and Methods, Springer, 2nd Edn., ISBN: 9781441903198, 1991. a, b
Burton, S. P., Hair, J. W., Kahnert, M., Ferrare, R. A., Hostetler, C. A., Cook, A. L., Harper, D. B., Berkoff, T. A., Seaman, S. T., Collins, J. E., Fenn, M. A., and Rogers, R. R.: Observations of the spectral dependence of linear particle depolarization ratio of aerosols using NASA Langley airborne High Spectral Resolution Lidar, Atmos. Chem. Phys., 15, 13453–13473, https://doi.org/10.5194/acp15134532015, 2015. a
Cairns, B., Russell, E. E., and Travis, L. D.: Research Scanning Polarimeter: calibration and groundbased measurements, Proc. SPIE, 3754, 186–196, https://doi.org/10.1117/12.366329, 1999. a, b
Chatterjee, S. and Simonoff, J. S.: Handbook of Regression Analysis, Wiley, ISBN: 9780470887165, 2012. a
Chowdhary, J., Cairns, B., Mishchenko, M., and Travis, L.: Retrieval of aerosol properties over the ocean using multispectral and multiangle Photopolarimetric measurements from the Research Scanning Polarimeter, Geophys. Res. Lett., 28, 243–246, https://doi.org/10.1029/2000GL011783, 2001. a
Cox, C. and Munk, W.: Measurement of the Roughness of the Sea Surface from Photographs of the Sun's Glitter, J. Opt. Soc. Am., 44, 838–850, 1954. a
Deschamps, P.Y., Breon, F.M., Leroy, M., Podaire, A., Bricaud, A., Buriez, J.C., and Seze, G.: The POLDER mission: instrument characteristics and scientific objectives, IEEE T. Geosci. Remote, 32, 598–615, https://doi.org/10.1109/36.297978, 1994. a
Diner, D. J., Beckert, J. C., Reilly, T. H., Bruegge, C. J., Conel, J. E., Kahn, R. A., Martonchik, J. V., Ackerman, T. P., Davies, R., Gerstl, S. A. W., Gordon, H. R., Muller, J. ., Myneni, R. B., Sellers, P. J., Pinty, B., and Verstraete, M. M.: Multiangle Imaging SpectroRadiometer (MISR) instrument description and experiment overview, IEEE T. Geosci. Remote, 36, 1072–1087, https://doi.org/10.1109/36.700992, 1998. a
Diner, D. J., Xu, F., Garay, M. J., Martonchik, J. V., Rheingans, B. E., Geier, S., Davis, A., Hancock, B. R., Jovanovic, V. M., Bull, M. A., Capraro, K., Chipman, R. A., and McClain, S. C.: The Airborne Multiangle SpectroPolarimetric Imager (AirMSPI): a new tool for aerosol and cloud remote sensing, Atmos. Meas. Tech., 6, 2007–2025, https://doi.org/10.5194/amt620072013, 2013. a
Diner, D. J., Boland, S. W., Brauer, M., Bruegge, C., Burke, K. A., Chipman, R., Girolamo, L. D., Garay, M. J., Hasheminassab, S., Hyer, E., Jerrett, M., Jovanovic, V., Kalashnikova, O. V., Liu, Y., Lyapustin, A. I., Martin, R. V., Nastan, A., Ostro, B. D., Ritz, B., Schwartz, J., Wang, J., and Xu, F.: Advances in multiangle satellite remote sensing of speciated airborne particulate matter and association with adverse health effects: from MISR to MAIA, J. Appl. Remote Sens., 12, 1–22, https://doi.org/10.1117/1.JRS.12.042603, 2018. a
Dubovik, O., Sinyuk, A., Lapyonok, T., Holben, B. N., Mishchenko, M., Yang, P., Eck, T. F., Volten, H., Muñoz, O., Veihelmann, B., van der Zande, W. J., Leon, J.F., Sorokin, M., and Slutsker, I.: Application of spheroid models to account for aerosol particle nonsphericity in remote sensing of desert dust, J. Geophys. Res.Atmos., 111, D11208, https://doi.org/10.1029/2005JD006619, 2006. a
Dubovik, O., Lapyonok, T., Litvinov, P., Herman, M., Fuertes, D., Ducos, F., Lopatin, A., Chaikovsky, A., Torres, B., Derimian, Y., Huang, X., Aspetsberger, M., and Federspiel, C.: GRASP: a versatile algorithm for characterizing the atmosphere, SPIE Newsroom, https://doi.org/10.1117/2.1201408.005558, 2014. a
Dubovik, O., Li, Z., Mishchenko, M. I., Tanré, D., Karol, Y., Bojkov, B., Cairns, B., Diner, D. J., Espinosa, W. R., Goloub, P., Gu, X., Hasekamp, O., Hong, J., Hou, W., Knobelspiesse, K. D., Landgraf, J., Li, L., Litvinov, P., Liu, Y., Lopatin, A., Marbach, T., Maring, H., Martins, V., Meijer, Y., Milinevsky, G., Mukai, S., Parol, F., Qiao, Y., Remer, L., Rietjens, J., Sano, I., Stammes, P., Stamnes, S., Sun, X., Tabary, P., Travis, L. D., Waquet, F., Xu, F., Yan, C., and Yin, D.: Polarimetric remote sensing of atmospheric aerosols: Instruments, methodologies, results, and perspectives, J. Quant. Spectrosc. Ra., 224, 474–511, https://doi.org/10.1016/j.jqsrt.2018.11.024, 2019. a, b
Dubovik, O., Fuertes, D., Litvinov, P., Lopatin, A., Lapyonok, T., Doubovik, I., Xu, F., Ducos, F., Chen, C., Torres, B., Derimian, Y., Li, L., HerrerasGiralda, M., Herrera, M., Karol, Y., Matar, C., Schuster, G. L., Espinosa, R., Puthukkudy, A., Li, Z., Fischer, J., Preusker, R., Cuesta, J., Kreuter, A., Cede, A., Aspetsberger, M., Marth, D., Bindreiter, L., Hangler, A., Lanzinger, V., Holter, C., and Federspiel, C.: A Comprehensive Description of MultiTerm LSM for Applying Multiple a Priori Constraints in Problems of Atmospheric Remote Sensing: GRASP Algorithm, Concept, and Applications, Frontiers in Remote Sensing, 2, 23, https://doi.org/10.3389/frsen.2021.706851, 2021. a
Fougnie, B., Marbach, T., Lacan, A., Lang, R., Schlüssel, P., Poli, G., Munro, R., and Couto, A. B.: The multiviewing multichannel multipolarisation imager – Overview of the 3MI polarimetric mission for aerosol and cloud characterization, J. Quant. Spectrosc. Ra., 219, 23–32, https://doi.org/10.1016/j.jqsrt.2018.07.008, 2018. a
Frouin, R. J., Franz, B. A., Ibrahim, A., Knobelspiesse, K., Ahmad, Z., Cairns, B., Chowdhary, J., Dierssen, H. M., Tan, J., Dubovik, O., Huang, X., Davis, A. B., Kalashnikova, O., Thompson, D. R., Remer, L. A., Boss, E., Coddington, O., Deschamps, P.Y., Gao, B.C., Gross, L., Hasekamp, O., Omar, A., Pelletier, B., Ramon, D., Steinmetz, F., and Zhai, P.W.: Atmospheric Correction of Satellite OceanColor Imagery During the PACE Era, Frontiers in Earth Science, 7, 145, https://doi.org/10.3389/feart.2019.00145, 2019. a
Gao, M., Zhai, P.W., Franz, B., Hu, Y., Knobelspiesse, K., Werdell, P. J., Ibrahim, A., Xu, F., and Cairns, B.: Retrieval of aerosol properties and waterleaving reflectance from multiangular polarimetric measurements over coastal waters, Opt. Express, 26, 8968–8989, https://doi.org/10.1364/OE.26.008968, 2018. a
Gao, M., Zhai, P.W., Franz, B. A., Hu, Y., Knobelspiesse, K., Werdell, P. J., Ibrahim, A., Cairns, B., and Chase, A.: Inversion of multiangular polarimetric measurements over open and coastal ocean waters: a joint retrieval algorithm for aerosol and waterleaving radiance properties, Atmos. Meas. Tech., 12, 3921–3941, https://doi.org/10.5194/amt1239212019, 2019. a
Gao, M., Zhai, P.W., Franz, B. A., Knobelspiesse, K., Ibrahim, A., Cairns, B., Craig, S. E., Fu, G., Hasekamp, O., Hu, Y., and Werdell, P. J.: Inversion of multiangular polarimetric measurements from the ACEPOL campaign: an application of improving aerosol property and hyperspectral ocean color retrievals, Atmos. Meas. Tech., 13, 3939–3956, https://doi.org/10.5194/amt1339392020, 2020. a
Gao, M., Franz, B. A., Knobelspiesse, K., Zhai, P.W., Martins, V., Burton, S., Cairns, B., Ferrare, R., Gales, J., Hasekamp, O., Hu, Y., Ibrahim, A., McBride, B., Puthukkudy, A., Werdell, P. J., and Xu, X.: Efficient multiangle polarimetric inversion of aerosols and ocean color powered by a deep neural network forward model, Atmos. Meas. Tech., 14, 4083–4110, https://doi.org/10.5194/amt1440832021, 2021a. a, b, c, d, e, f, g, h, i, j, k, l
Gao, M., Knobelspiesse, K., Franz, B. A., Zhai, P.W., Martins, V., Burton, S. P., Cairns, B., Ferrare, R., Fenn, M. A., Hasekamp, O., Hu, Y., Ibrahim, A., Sayer, A. M., Werdell, P. J., and Xu, X.: Adaptive Data Screening for MultiAngle Polarimetric Aerosol and Ocean Color Remote Sensing Accelerated by Deep Learning, Frontiers in Remote Sensing, 2, 46, https://doi.org/10.3389/frsen.2021.757832, 2021b. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q
Gao, M., Knobelspiesse, K., Franz, B. A., Zhai, P.W., Sayer, A. M., Ibrahim, A., Cairns, B., Hasekamp, O., Hu, Y., Martins, V., Werdell, P. J., and Xu, X.: Effective uncertainty quantification for multiangle polarimetric aerosol remote sensing over ocean, Atmos. Meas. Tech., 15, 4859–4879, https://doi.org/10.5194/amt1548592022, 2022. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p
Hannadige, N. K., Zhai, P.W., Gao, M., Franz, B. A., Hu, Y., Knobelspiesse, K., Werdell, P. J., Ibrahim, A., Cairns, B., and Hasekamp, O. P.: Atmospheric correction over the ocean for hyperspectral radiometers using multiangle polarimetric retrievals, Opt. Express, 29, 4504–4522, https://doi.org/10.1364/OE.408467, 2021. a
Hasekamp, O. P. and Landgraf, J.: Retrieval of aerosol properties over land surfaces: capabilities of multipleviewingangle intensity and polarization measurements, Appl. Optics, 46, 3332–3344, https://doi.org/10.1364/AO.46.003332, 2007. a
Hasekamp, O. P., Fu, G., Rusli, S. P., Wu, L., Noia, A. D., aan de Brugh, J., Landgraf, J., Smit, J. M., Rietjens, J., and van Amerongen, A.: Aerosol measurements by SPEXone on the NASA PACE mission: expected retrieval capabilities, J. Quant. Spectrosc. Ra., 227, 170–184, https://doi.org/10.1016/j.jqsrt.2019.02.006, 2019. a, b
Knobelspiesse, K., Cairns, B., Mishchenko, M., Chowdhary, J., Tsigaridis, K., van Diedenhoven, B., Martin, W., Ottaviani, M., and Alexandrov, M.: Analysis of finemode aerosol retrieval capabilities by different passive remote sensing instrument designs, Opt. Express, 20, 21457–21484, https://doi.org/10.1364/OE.20.021457, 2012. a, b, c, d, e, f, g
Knobelspiesse, K., Tan, Q., Bruegge, C., Cairns, B., Chowdhary, J., van Diedenhoven, B., Diner, D., Ferrare, R., van Harten, G., Jovanovic, V., Ottaviani, M., Redemann, J., Seidel, F., and Sinclair, K.: Intercomparison of airborne multiangle polarimeter observations from the Polarimeter Definition Experiment, Appl. Optics, 58, 650–669, https://doi.org/10.1364/AO.58.000650, 2019. a
Knobelspiesse, K., Barbosa, H. M. J., Bradley, C., Bruegge, C., Cairns, B., Chen, G., Chowdhary, J., Cook, A., Di Noia, A., van Diedenhoven, B., Diner, D. J., Ferrare, R., Fu, G., Gao, M., Garay, M., Hair, J., Harper, D., van Harten, G., Hasekamp, O., Helmlinger, M., Hostetler, C., Kalashnikova, O., Kupchock, A., Longo De Freitas, K., Maring, H., Martins, J. V., McBride, B., McGill, M., Norlin, K., Puthukkudy, A., Rheingans, B., Rietjens, J., Seidel, F. C., da Silva, A., Smit, M., Stamnes, S., Tan, Q., Val, S., Wasilewski, A., Xu, F., Xu, X., and Yorks, J.: The Aerosol Characterization from Polarimeter and Lidar (ACEPOL) airborne field campaign, Earth Syst. Sci. Data, 12, 2183–2208, https://doi.org/10.5194/essd1221832020, 2020. a
Kokhanovsky, A., Davis, A., Cairns, B., Dubovik, O., Hasekamp, O., Sano, I., Mukai, S., Rozanov, V., Litvinov, P., Lapyonok, T., Kolomiets, I., Oberemok, Y., Savenkov, S., Martin, W., Wasilewski, A., Di Noia, A., Stap, F., Rietjens, J., Xu, F., Natraj, V., Duan, M., Cheng, T., and Munro, R.: Spacebased remote sensing of atmospheric aerosols: The multiangle spectropolarimetric frontier, EarthSci. Rev., 145, 85–116, https://doi.org/10.1016/j.earscirev.2015.01.012, 2015. a
Lamquin, N., Mangin, A., Mazeran, C., Bourg, B., Bruniquel, V., and D’Ando, O. F.: OLCI L2 PixelbyPixelUncertainty Propagation in OLCI Clean Water Branch, ESA, Ref.: S3L2SD01C01ACRTN, 2013. a
Martins, J. V., FernandezBorda, R., McBride, B., Remer, L., and Barbosa, H. M. J.: The HARP hyperangular imaging polarimeter and the need for small satellite payloads with high science payoff for earth science remote sensing, in: IGARSS 2018 – 2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018, IEEE, 6304–6307, https://doi.org/10.1109/IGARSS.2018.8518823, 2018. a, b
McBride, B. A., Martins, J. V., Barbosa, H. M. J., Birmingham, W., and Remer, L. A.: Spatial distribution of cloud droplet size properties from Airborne HyperAngular Rainbow Polarimeter (AirHARP) measurements, Atmos. Meas. Tech., 13, 1777–1796, https://doi.org/10.5194/amt1317772020, 2020. a
McGill, M., Hlavka, D., Hart, W., Scott, V. S., Spinhirne, J., and Schmid, B.: Cloud Physics Lidar: instrument description and initial measurement results, Appl. Optics, 41, 3725–3734, https://doi.org/10.1364/AO.41.003725, 2002. a
Mishchenko, M. I. and Travis, L. D.: Satellite retrieval of aerosol properties over the ocean using polarization as well as intensity of reflected sunlight, J. Geophys. Res.Atmos., 102, 16989–17013, https://doi.org/10.1029/96JD02425, 1997. a
Mishchenko, M. I., Cairns, B., Kopp, G., Schueler, C. F., Fafaul, B. A., Hansen, J. E., Hooker, R. J., Itchkawich, T., Maring, H. B., and Travis, L. D.: Accurate Monitoring of Terrestrial Aerosols and Total Solar Irradiance: Introducing the Glory Mission, B. Am. Meteorol. Soc., 88, 677–692, https://doi.org/10.1175/BAMS885677, 2007. a
Mobley, C. D.: The Oceanic Optics Book, International Ocean Colour Coordinating Group (IOCCG), Dartmouth, NS, Canada, https://doi.org/10.25607/OBP1710, 2022. a
National Academies of Sciences, Engineering, and Medicine: Thriving on Our Changing Planet: A Decadal Strategy for Earth Observation from Space, The National Academies Press, Washington, DC, https://doi.org/10.17226/24938, 2018. a
Paszke, A., Gross, S., Massa, F., Lerer, A., Bradbury, J., Chanan, G., Killeen, T., Lin, Z., Gimelshein, N., Antiga, L., Desmaison, A., Kopf, A., Yang, E., DeVito, Z., Raison, M., Tejani, A., Chilamkurthy, S., Steiner, B., Fang, L., Bai, J., and Chintala, S.: PyTorch: An Imperative Style, HighPerformance Deep Learning Library, Advances in Neural Information Processing Systems 32, 8024–8035, http://papers.neurips.cc/paper/ 9015pytorchanimperativestylehighperformancedeeplearninglibrary.pdf (last access: 3 April 2023), 2019. a
Priestley, M. B.: Spectral Analysis and Time Series, Academic Press, ISBN: 9780125649223, 1983. a, b, c, d
Rodgers, C.: Inverse Methods for Atmospheric Sounding: Theory and Practice, World Scientific World Scientific Publishing, Singapore, https://doi.org/10.1142/3171, 2000. a, b, c, d, e
Sayer, A. M., Govaerts, Y., Kolmonen, P., Lipponen, A., Luffarelli, M., Mielonen, T., Patadia, F., Popp, T., Povey, A. C., Stebel, K., and Witek, M. L.: A review and framework for the evaluation of pixellevel uncertainty estimates in satellite aerosol remote sensing, Atmos. Meas. Tech., 13, 373–404, https://doi.org/10.5194/amt133732020, 2020. a
Seabold, S. and Perktold, J.: Statsmodels: Econometric and statistical modeling with python, in: 9th Python in Science Conference, Austin, Texas, https://doi.org/10.25080/Majora92bf1922011, 2010. a
Seegers, B. N., Stumpf, R. P., Schaeffer, B. A., Loftin, K. A., and Werdell, P. J.: Performance metrics for the assessment of satellite data products: an ocean color case study, Opt. Express, 26, 7404–7422, https://doi.org/10.1364/OE.26.007404, 2018. a
Smit, J. M., Rietjens, J. H. H., van Harten, G., Noia, A. D., Laauwen, W., Rheingans, B. E., Diner, D. J., Cairns, B., Wasilewski, A., Knobelspiesse, K. D., Ferrare, R., and Hasekamp, O. P.: SPEX airborne spectropolarimeter calibration and performance, Appl. Optics, 58, 5695–5719, https://doi.org/10.1364/AO.58.005695, 2019. a, b
Tanré, D., Bréon, F. M., Deuzé, J. L., Dubovik, O., Ducos, F., François, P., Goloub, P., Herman, M., Lifermann, A., and Waquet, F.: Remote sensing of aerosols by using polarized, directional and spectral measurements within the ATrain: the PARASOL mission, Atmos. Meas. Tech., 4, 1383–1395, https://doi.org/10.5194/amt413832011, 2011. a
Theys, N., Fioletov, V., Li, C., De Smedt, I., Lerot, C., McLinden, C., Krotkov, N., Griffin, D., Clarisse, L., Hedelt, P., Loyola, D., Wagner, T., Kumar, V., Innes, A., Ribas, R., Hendrick, F., Vlietinck, J., Brenot, H., and Van Roozendael, M.: A sulfur dioxide CovarianceBased Retrieval Algorithm (COBRA): application to TROPOMI reveals new emission sources, Atmos. Chem. Phys., 21, 16727–16744, https://doi.org/10.5194/acp21167272021, 2021. a
Virtanen, T. H., Kolmonen, P., Sogacheva, L., Rodríguez, E., Saponaro, G., and de Leeuw, G.: Collocation mismatch uncertainties in satellite aerosol retrieval validation, Atmos. Meas. Tech., 11, 925–938, https://doi.org/10.5194/amt119252018, 2018. a
Waquet, F., Riedi, J., Labonnote, L. C., Goloub, P., Cairns, B., Deuzé, J.L., and Tanré, D.: Aerosol Remote Sensing over Clouds Using ATrain Observations, J. Atmos. Sci., 66, 2468–2480, https://doi.org/10.1175/2009jas3026.1, 2009. a
Werdell, P. J., Behrenfeld, M. J., Bontempi, P. S., Boss, E., Cairns, B., Davis, G. T., Franz, B. A., Gliese, U. B., Gorman, E. T., Hasekamp, O., Knobelspiesse, K. D., Mannino, A., Martins, J. V., McClain, C. R., Meister, G., and Remer, L. A.: The Plankton, Aerosol, Cloud, Ocean Ecosystem Mission: Status, Science, Advances, B. Am. Meteorol. Soc., 100, 1775–1794, https://doi.org/10.1175/BAMSD180056.1, 2019. a
Xu, F., Dubovik, O., Zhai, P.W., Diner, D. J., Kalashnikova, O. V., Seidel, F. C., Litvinov, P., Bovchaliuk, A., Garay, M. J., van Harten, G., and Davis, A. B.: Joint retrieval of aerosol and waterleaving radiance from multispectral, multiangular and polarimetric measurements over ocean, Atmos. Meas. Tech., 9, 2877–2907, https://doi.org/10.5194/amt928772016, 2016. a
Xu, F., Diner, D. J., Dubovik, O., and Schechner, Y.: A Correlated MultiPixel Inversion Approach for Aerosol Remote Sensing, Remote Sens.Basel, 11, 746, https://doi.org/10.3390/rs11070746, 2019. a
Zhai, P.W. and Hu, Y.: An improved pseudo spherical shell algorithm for vector radiative transfer, J. Quant. Spectrosc. Ra., 282, 108132, https://doi.org/10.1016/j.jqsrt.2022.108132, 2022. a
Zhai, P.W., Hu, Y., Trepte, C. R., and Lucker, P. L.: A vector radiative transfer model for coupled atmosphere and ocean systems based on successive order of scattering method, Opt. Express, 17, 2057–2079, 2009. a
Zhai, P.W., Hu, Y., Chowdhary, J., Trepte, C. R., Lucker, P. L., and Josset, D. B.: A vector radiative transfer model for coupled atmosphere and ocean systems with a rough interface, J. Quant. Spectrosc. Ra., 111, 1025–1040, 2010. a
Zhang, M., Ibrahim, A., Franz, B. A., Ahmad, Z., and Sayer, A. M.: Estimating pixellevel uncertainty in ocean color retrievals from MODIS, Opt. Express, 30, 31415–31438, https://doi.org/10.1364/OE.460735, 2022. a
Zibordi, G., Mélin, F., Berthon, J.F., Holben, B., Slutsker, I., Giles, D., D’Alimonte, D., Vandemark, D., Feng, H., Schuster, G., Fabbri, B. E., Kaitala, S., and Seppälä, J.: AERONETOC: A Network for the Validation of Ocean Color Primary Products, J. Atmos. Ocean. Tech., 26, 1634–1651, https://doi.org/10.1175/2009JTECHO654.1, 2009. a
 Abstract
 Introduction
 Algorithm and methodology
 Retrieval uncertainties for AirHARP and HARP2
 Estimating correlation from residual analysis
 Discussions and conclusions
 Appendix A: Input parameters of the neuralnetwork forward model
 Appendix B: Eigenvector decomposition in the error covariance matrix
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Algorithm and methodology
 Retrieval uncertainties for AirHARP and HARP2
 Estimating correlation from residual analysis
 Discussions and conclusions
 Appendix A: Input parameters of the neuralnetwork forward model
 Appendix B: Eigenvector decomposition in the error covariance matrix
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References