Effect of Wind Speed on Moderate Resolution Imaging Spectroradiometer ( MODIS ) Aerosol Optical Depth over the North Pacific

The surface-wind speed influences on aerosol optical depth (AOD), derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua daily observations over the central North Pacific during the period 2003–2016, have been investigated in this study. The cloud coverage is relatively low over the present investigation area compared to other marine areas, which favors AOD derived from passive remote sensing from space. In this study, we have combined MODIS AOD with 2 m wind speed (U2m) on a satellite-pixel basis, which has been interpolated from National Centers for Environmental Prediction (NCEP) reanalysis. In addition, daily averaged AOD derived from Aerosol Robotic Network (AERONET) measurements in the free-troposphere at the Mauna Loa Observatory (3397 m above sea level), Hawaii, was subtracted from the MODIS column AOD values. The latter was to reduce the contribution of aerosols above the planetary boundary layer. This study shows relatively strong power-law relationships between MODIS mean AOD and surface-wind speed for marine background conditions in summer, fall and winter of the current period. However, previous established relationships between AOD and surface-wind speed deviate substantially. Even so, for similar marine conditions the present relationship agrees reasonable well with a power-law relationship derived for north-east Atlantic conditions. The present MODIS retrievals of AOD in the marine atmosphere agree reasonably well with ground-based remote sensing of AOD.


Introduction
Oceans cover approximately two-thirds of the Earth's surface and are the major source of natural aerosol mass.In terms of mass, sea spray produced at the ocean surface dominate global aerosol flux with an annual sea salt production of 0.3 Pg to 30 Pg [1] and organic matter of 8-50 Tg C year −1 [2].Due to high hygroscopicity of sea salt and regional abundance, the aerosol provides a major contribution to the scattering of solar radiation, resulting in a cooling effect [3][4][5][6][7].Moreover, serving as an efficient cloud condensation nuclei the sea salt aerosol exerts a strong influence on cloud formation, properties and lifetime [8,9].The importance of sea-spray aerosol's impact on the radiative balance, both on global and regional scales, has led to extensive research.This has aimed for a better understanding of emissions, composition, transformation, transport and removal, as well as interactions with radiation and clouds, (e.g., [1,[10][11][12][13][14]).To represent oceanic aerosols explicitly in global climate models it is necessary to describe primary marine aerosol emission fluxes quantitatively, with respect to their size and chemical composition [13].
The main physical parameters that determine sea-spray aerosol production are: (i) sea ice cover ( [10]; (ii) wind speed [10,15]; and (iii) surface water temperature [11,16,17].Sea ice acts as a lid over the sea and over open leads between ice floes, and the ice limits the fetch and prevents wave breaking, which is needed for sea-spray formation.At an ice cover of about 80-90% the sea-spray emissions are reduced with one order of magnitude over the patches of open water, and on average by two orders of magnitude over larger areas of the same ice fraction, compared to open sea.The wind stresses the ocean surface and when waves break air is entraining into the water.Air bubbles rise to the surface, creating white colored areas (whitecaps), where they eventually burst and release aerosol droplets into the atmosphere [1,[18][19][20].The wind effect on sea spray production is often included in the source parameterizations through an expression of the white cap coverage, but it can also be done through the amount of entrained air [17].The effect of increasing wind is highly non-linear and results in about a factor 10 increase in sea-spray emission per 10 m/s increase.The water temperature changes the bubble spectra through the combined influence of surface tension, viscosity, coalescence, rise velocity and oxygen saturation [11,16].The net effect is a decreasing sea-spray production with an increasing water temperature.Although there have been many achievements in understanding sea-spray aerosol production, the uncertainty in existing parametrizations remains sufficiently large [1,13].Recently proposed parametrizations of the production flux result in two orders of magnitude higher aerosol concentrations simulated than observed, raising questions about the validity of such estimates and methods to attain them [13].
Unlike the continental atmospheric boundary layer, the marine atmospheric boundary layer has mostly a near-neutral stratification and, thus, we do not expect to find a well-mixed layer.In addition, decoupled layers frequently occur over the ocean [21][22][23][24][25] that create vertical gradients in sea-salt aerosol concentrations in the lower troposphere [12,14,[26][27][28][29][30].Relatively strong vertical gradients have also been observed and modeled for coarse mode aerosols (>0.5 µm radius), and modeled for submicron aerosols, in otherwise well-mixed marine boundary layers (MBLs) without thermodynamic evidence of decoupling [12].These findings are explained by the gravitational deposition velocity of the coarse-mode particles.This implies that in situ measurements of sea-salt concentration at one altitude do not represent the entire marine boundary layer.
Since winds near the surface are an important factor in the formation of breaking waves and bubbles that burst [10,13,31], which results in the release of sea-spray aerosol, this supports investigations of the relationship between aerosol optical depth (AOD) and surface-wind speed.The other two major physical drivers, sea ice and surface-water temperature, change much slower than the wind speed, for which the variation is found on a time scale smaller than about five days (the synoptic time scale).It is therefore possible to study the effect of surface-wind speed on AOD separated from temperature and sea ice even at high latitudes.The wind dependency allows for the establishment of an empirical relationship between AOD and wind speed that can be used in global and regional model simulations and/or, more importantly, for evaluating model simulations of AOD in the marine atmosphere.
The AOD-wind speed relationship has been investigated in several previous studies, although there is no consensus between the findings.Some studies suggest a linear relationship [32][33][34][35], while others a power-law relationship [14,36,37] or an exponential relationship [38].The discrepancy between different relationships is probably attributed to different instrumentations, platforms, locations and periods used in the investigations.
In the present study, we investigate how wind speed influences AOD that is derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) daily observations over the North Pacific during the period 2003-2016.The relationship estimated between AOD and surface-wind speed is compared to results obtained in several previous studies.In addition, we examine how surface-wind speed influences the contribution of coarse and fine-mode aerosols to the AOD.The satellite-derived AOD is validated against sun photometer measurements carried out at stationary ground-based sites and during ship cruises.

Moderate Resolution Imaging Spectroradiometer (MODIS) Aerosol Retrievals
MODIS on the Aqua satellite measures scattered sunlight at 36 spectral frequencies and provides daily global aerosol products.From 700 km above sea level and with nadir views, between the latitudes 0 • and 55 • , MODIS observes the planet with a swath of about 2330 km.MODIS Aqua makes 14-15 orbits per day with local equatorial crossing time approximately 13:30 (ascending node).All MODIS data are organized into 5-min intervals or granules.With a 10 × 10 km spatial resolution at nadir, the output grid results in 135 × 203 pixels for one granule.
MODIS retrievals of AOD over the land and ocean are based on the dark-target algorithm, thoroughly described in the literature [39][40][41][42][43][44].The algorithm calculates AOD over an expected dark surface, while over the ocean the effects of sun glint and the presence of foam (whitecaps) on the water need to be accounted for.The retrievals of AOD, both over land and ocean, are based on look-up tables (LUTs), where normalized satellite-detected scattered radiance are compared with pre-computed radiative transfer, with a set of likely aerosol situations as well as prescribed Rayleigh scattering and surface reflection.The ocean algorithm consists of nine tropospheric aerosol model types: four fine modes (with effective radius (r e ) < 0.25 µm) and five coarse modes (r e > 1.0 µm), where each type is a single mode of log-normal size distribution.To relate the measured reflectance to AOD, the algorithm finds a solution of weighted combinations of fine and coarse modes that approximates best the observed spectral reflectance [43].MODIS data is presented in different collections consisting of data products that were generated by similar but not necessarily the same versions of the retrieval algorithm.The AOD that is derived by the MODIS land and ocean retrievals at 555 nm has expected error envelopes of ∆AOD = ±0.05± 0.15•AOD [42] and In this study, we use MODIS Aqua Collection 6 (C6) level 2 standard 10-km products for best-quality retrievals (quality flag = 3) of AOD, mainly over the ocean surface.Data were taken from the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center's Atmosphere Archive and Distribution System (http://ladsweb.nascom.nasa.gov).In previous MODIS collections, the calculation of sun glint and foam on the water surface were based on a single value of wind speed (6 m/s) and foam fraction (0.16%) [43].For better representation of the surface roughness, additional wind speeds (2 m/s, 10 m/s and 14 m/s) and corresponding foam fractions (0.01%, 1% and 3%) were added to MODIS C6.C6 also includes 2 m wind speed (U 2m ) on a satellite-pixel basis, which has been interpolated from National Centers for Environmental Prediction (NCEP) reanalysis (1 • × 1 • ) GDAS 10 m wind data as a required ancillary field.In addition to total AOD and U 2m , aerosol size parameters such as Ångström exponent and contribution of fine-mode aerosols to the AOD, have also been investigated.The MODIS Ångström exponent (α MOD ) has been calculated based on AOD retrieved at the two wavelengths 555 nm and 858 nm:

Aerosol Optical Depth (AOD) from Sun Photometer Measurements
The Aerosol Robotic Network (AERONET) provides ground-based measurements of aerosol optical properties.AERONET is equipped with Cimel sun photometers that measure direct sun radiance at eight spectral bands.Measurements are taken every 15 min during daylight hours and are available in three categories: level 1.0 (unscreened), level 1.5 (cloud screened) and level 2.0 (quality assured).In this study, level 2.0 data are used.AERONET-derived AOD are expected to be accurate within <±0.01 for wavelengths longer than 440 nm, (e.g., [45]).The component of AERONET, Maritime Aerosol Network (MAN), provides ship-borne AOD measurements with Microtops II sun photometers calibrated at the NASA Goddard Space Flight Center.The estimated uncertainty of the MAN AOD does not exceed ±0.02 [46].
To validate MODIS AOD we compared the retrievals against AOD derived from AERONET observations at the three ground-based stations: Lanai (20 m, 20.7 • N, 156.9 • W), Coconut Island (0 m, 21.4 • N, 157.8 • W) and Midway Island (0 m, 28.2 • N, 177.4 • W).Observations during several cruises predominant in the North Pacific (Table 1) have also been used in the validation.AOD from AERONET and MAN is available at the 500 nm wavelength.We adjusted the AERONET and MAN AOD values to correspond to the MODIS AOD wavelength of 555 nm using the Ångström power law: where α AER is the AERONET Ångström exponent derived from the wavelengths 440 nm and 675 nm.

Investigation Area: Central North Pacific
Figure 1 shows the investigation area (12.0 • N to 35.0 • N, 145.0 • W to 168.0 • W) in the central North Pacific and an example of a MODIS Aqua scene with AODs for 18 September 2014.All MODIS Aqua AODs that are available within the black box and corresponding to the period 2003-2016 are included in the investigation.This results in 26,210 granules of C6 data.This area was chosen in order to minimize the influence from continental aerosols on the air masses investigated, This region is highly influenced by a persistent sub-tropical high-pressure system, which shifts north during the summer months and south during the winter season.The cloud coverage is relatively low over the central North Pacific compared to other marine areas, which favors satellite retrievals of AOD.The north-eastern trade-winds over this region are very consistent.They prevail around Hawaii between 85% and 95% of the time during summer and around 50% to 80% of the time during the winter season [47].

FLEXible TRAjectories (FLEXTRA) Trajectory Model
FLEXible TRAjectories (FLEXTRA) were used to calculate mean wind trajectories operating directly on data from the European Centre for Medium-Range Weather Forecasts (ECMWF).Meteorological-level fields from the ECMWF T213 L31 model were used.FLEXTRA has an ensemble option to assess the accuracy of the trajectories, and different options for the interpolation of the wind data.For more details of the FLEXTRA model, see [48,49].In this study, we analyzed trajectories that run 7 days back in time from Mauna Loa mountain station (19.54 • N, 155.58 • W, 3000 m above sea level, Figure 1), Hawaii, with starting times 00, 06, 12 and 18 UTC each day within the period 2003-2010.

Validation of MODIS AOD against Ground-Based Measurements
Figure 2a,b shows comparisons of MODIS AOD against AOD derived from AERONET and MAN measurements, respectively, carried out at coastal ground-based stations and during ship cruises, respectively.We used a box of 3 × 3 pixels (30 × 30 km) around each ground-based station or ship, although only cases with an availability of at least three AOD values were included in the comparison.In addition, satellite overpass windows of 30 min and 2.5 h were used here for AERONET and MAN, respectively.Previous studies have shown that MODIS retrievals of AOD over coastal sites are associated with significant uncertainties (Section 4), although when using the 30-min time window reasonable agreement between MODIS and AERONET was obtained (see below).The relative (normalized) root mean square deviation (relative RMSD) was determined for collocated satellite and ground-based AOD values, averaged according to the procedure described in [50].

Validation of MODIS AOD against Ground-Based Measurements
Figure 2a,b shows comparisons of MODIS AOD against AOD derived from AERONET and MAN measurements, respectively, carried out at coastal ground-based stations and during ship cruises, respectively.We used a box of 3 × 3 pixels (30 × 30 km) around each ground-based station or ship, although only cases with an availability of at least three AOD values were included in the comparison.In addition, satellite overpass windows of 30 min and 2.5 h were used here for AERONET and MAN, respectively.Previous studies have shown that MODIS retrievals of AOD over coastal sites are associated with significant uncertainties (Section 4), although when using the 30-min time window reasonable agreement between MODIS and AERONET was obtained (see below).The relative (normalized) root mean square deviation (relative RMSD) was determined for collocated satellite and ground-based AOD values, averaged according to the procedure described in [50].
Figure 2a shows that MODIS AOD agree reasonably well with sun photometer measurements carried out at the three marine ground-based sites.Most AOD values are below 0.2.In the validation against AERONET, the percentage of the MODIS values that is within the expected uncertainty is 71%.For ship-born measurements, 75% of the collocated points lie within the expected uncertainty (right figure).The vertical and horizontal error bars in Figure 2 show that some of the MODIS and AERONET/MAN AODs are associated with large variability.
In addition, satellite overpass windows of 30 min and 2.5 h were used here for AERONET and MAN, respectively.Previous studies have shown that MODIS retrievals of AOD over coastal sites are associated with significant uncertainties (Section 4), although when using the 30-min time window reasonable agreement between MODIS and AERONET was obtained (see below).The relative (normalized) root mean square deviation (relative RMSD) was determined for collocated satellite and ground-based AOD values, averaged according to the procedure described in [50].

Evaluation of MODIS-Derived AOD versus Surface Wind Speed
Fourteen years (2003-2016) of MODIS daily column AOD, at the wavelength 555 nm, have been matched with surface (2-m) wind speed over the central North Pacific.However, daily averaged AERONET AOD measured in the free-troposphere at the Mauna Loa Observatory (3397 m above sea level, 19.54 • N, 155.58 • W), Hawaii, was subtracted from the MODIS column AOD values.The latter was undertaken to reduce the contribution of aerosols above the planetary boundary layer, which are expected to be weak related to surface-wind speed.For the few days with missing AERONET data we used AOD averaged over the two neighboring days.The adjusted AOD values were binned with respect to surface wind speed, with bins of 0.5 m/s within the range 0-10 m/s. Figure 3 shows an example of all MODIS-adjusted AOD values on a pixel basis, versus surface-wind speed, derived from observations over the present investigation area during summer 2007.Cases with AOD values around and below 0.1 suggest that natural background aerosol was dominating, while the higher values may be due to hazier conditions and/or the effects of cloud contamination.However, the median absolute deviation suggests that the latter values are relatively few with respect to all AODs. Figure 3 shows that the median AOD starts to rise at a wind speed about 3.5 m s −1 , although a decrease appears for the highest wind speeds.
Furthermore, Figure 4 presents only the binned data on a seasonal basis, although for all years investigated here.The winter, spring, summer and fall seasons are represented by the months December-February, March-May, June-August and September-November, respectively.In summer, fall and winter Figure 4 shows that the median AOD increases roughly by a factor of 2 within the wind speed range about 3.5-9 m/s.In spring the increase in AOD is similar in absolute terms, however, the background AOD is substantially higher than in other seasons.This is likely due to the influence of continental aerosols (Section 4).Moreover, as in Figure 3, MODIS AOD do not exhibit an increase at wind speeds higher than ~9 m/s for any of the seasons (further discussed in Section 3.7).
Figure 4 shows that the MODIS-adjusted median AOD for the lower wind speeds in spring is about two times higher than in other seasons, and about 25% higher for the highest wind speeds.The figure also shows that the mean AODs (orange solid line) estimated for low wind speeds for winter are somewhat higher than AODs corresponding to summer and fall (more clearly shown in Figure 6b.The rate of change in AOD with increasing wind speed is smaller for winter than the other seasons.
speed, derived from observations over the present investigation area during summer 2007.Cases with AOD values around and below 0.1 suggest that natural background aerosol was dominating, while the higher values may be due to hazier conditions and/or the effects of cloud contamination.However, the median absolute deviation suggests that the latter values are relatively few with respect to all AODs. Figure 3 shows that the median AOD starts to rise at a wind speed about 3.5 m s −1 , although a decrease appears for the highest wind speeds.Furthermore, Figure 4 presents only the binned data on a seasonal basis, although for all years investigated here.The winter, spring, summer and fall seasons are represented by the months December-February, March-May, June-August and September-November, respectively.In summer, fall and winter Figure 4 shows that the median AOD increases roughly by a factor of 2 within the wind speed range about 3.5-9 m/s.In spring the increase in AOD is similar in absolute terms, however, the background AOD is substantially higher than in other seasons.This is likely due to the influence of continental aerosols (Section 4).Moreover, as in Figure 3, MODIS AOD do not exhibit an increase at wind speeds higher than ~9 m/s for any of the seasons (further discussed in Section 3.7).Figure 4 shows that the MODIS-adjusted median AOD for the lower wind speeds in spring is about two times higher than in other seasons, and about 25% higher for the highest wind speeds.The figure also shows that the mean AODs (orange solid line) estimated for low wind speeds for winter are somewhat higher than AODs corresponding to summer and fall (more clearly shown in Figure 6b.The rate of change in AOD with increasing wind speed is smaller for winter than the other seasons.
Figure 5a,b show seasonal wind speed probability density function for MODIS derived AOD, normalized according to the total number of derived AOD corresponding to a season and year, respectively.The majority of the AOD values were derived at wind speeds in the range between about 3.5 m/s and 8 m/s; thus, rather few AOD values correspond to lower and higher wind speeds.are somewhat higher than AODs corresponding to summer and fall (more clearly shown in Figure 6b.The rate of change in AOD with increasing wind speed is smaller for winter than the other seasons. Figure 5a,b show seasonal wind speed probability density function for MODIS derived AOD, normalized according to the total number of derived AOD corresponding to a season and year, respectively.The majority of the AOD values were derived at wind speeds in the range between about 3.5 m/s and 8 m/s; thus, rather few AOD values correspond to lower and higher wind speeds.For the most common wind speeds, Figure 5b shows that summer is associated with the highest number of available AODs.However, with the exception of 2011, 2014, 2015 and 2016 the winter season is associated with the highest number of available AODs over the entire wind speed range (not shown).In addition, more available AODs at the highest wind speeds in winter may be attributable to a higher For the most common wind speeds, Figure 5b shows that summer is associated with the highest number of available AODs.However, with the exception of 2011, 2014, 2015 and 2016 the winter season is associated with the highest number of available AODs over the entire wind speed range (not shown).In addition, more available AODs at the highest wind speeds in winter may be attributable to a higher frequency of low pressures in the westerly flow over the North Pacific (Section 3.4).The minimum number of observations is found for spring, with the exception of 2014 (not shown).
To identify whether the change in adjusted AOD with increased surface-wind speed is caused by sea-spray aerosol, within a particular size-range, we investigated the MODIS-derived Ångström exponent and the contribution of fine mode aerosols to the AOD. Figure 6a shows the seasonally averaged Ångström exponent according to the period 2003-2016 versus surface-wind speed.Since the extinction efficiency of larger-sized aerosol particles is nearly spectrally neutral, the decrease in the Ångström exponent indicates a shift towards larger aerosol particles at higher wind speeds.At lower wind speeds, the production of sea spray is low and it is reasonable to assume that the presence of secondary sub-micron aerosols such as ammonium sulfate, and maybe also organic compounds, dominates over the central North Pacific.The decrease in the Ångström exponent is evident for the current wind-speed range, although for winter the decrease is less salient compared to the other seasons.
Figure 6b shows the contribution of MODIS fine-mode adjusted AOD (for aerosols smaller than 0.5 µm in diameter) to the total adjusted AOD at 555 nm plotted against surface-wind speed and averaged over the period 2003-2016.Here, we focus on the results obtained for summer, fall and winter, for which we expect influences mainly from marine background aerosols.As shown in Figure 6b, the relationship between AOD and wind speed consists of two visible phases: AOD associated mainly with fine-mode aerosols at lower wind speeds and a distinct increase of larger aerosols that starts to occur around 3.5 m/s.Note, it is suggested that that white caps start to form at wind speeds of around 3.5 m s −1 [51][52][53][54].Even so, the parameterizations derived in the first three studies for low wind speeds, 3-5 m s −1 , include whitecap fractions that differ by more than one order in magnitude [53].Figure 6 shows that the lowest contribution of fine-mode aerosols to the total AOD at low wind speeds is occurring in winter, which is in line with the results of the Ångström exponent.This finding could be consistent with a contribution from secondary aerosol formation to fine aerosol, which should be smaller in winter due to a lower biological production.
extinction efficiency of larger-sized aerosol particles is nearly spectrally neutral, the decrease in the Ångström exponent indicates a shift towards larger aerosol particles at higher wind speeds.At lower wind speeds, the production of sea spray is low and it is reasonable to assume that the presence of secondary sub-micron aerosols such as ammonium sulfate, and maybe also organic compounds, dominates over the central North Pacific.The decrease in the Ångström exponent is evident for the current wind-speed range, although for winter the decrease is less salient compared to the other seasons.Figure 6b shows the contribution of MODIS fine-mode adjusted AOD (for aerosols smaller than 0.5 μm in diameter) to the total adjusted AOD at 555 nm plotted against surface-wind speed and averaged over the period 2003-2016.Here, we focus on the results obtained for summer, fall and winter, for which we expect influences mainly from marine background aerosols.As shown in Figure 6b, the relationship between AOD and wind speed consists of two visible phases: AOD associated mainly with fine-mode aerosols at lower wind speeds and a distinct increase of larger aerosols that starts to occur around 3.5 m/s.Note, it is suggested that that white caps start to form at wind speeds of around 3.5 m s −1 [51][52][53][54].Even so, the parameterizations derived in the first three studies for low wind speeds, 3-5 m s −1 , include whitecap fractions that differ by more than one order in magnitude [53].Figure 6 shows that the lowest contribution of fine-mode aerosols to the total AOD at low wind speeds is occurring in winter, which is in line with the results of the Ångström exponent.This finding could be consistent with a contribution from secondary aerosol formation to fine aerosol, which should be smaller in winter due to a lower biological production.

Aerosol Robotic Network (AERONET) Free Tropospheric AOD
Aerosols in the entire atmospheric column, between the sea surface and the top of the atmosphere, contribute to the AOD derived from MODIS observations.Quantifying the relative contributions of the boundary-layer and free-tropospheric aerosols to the total column AOD is not possible with passive remote sensing from the space and ground.Without deep convection that can effectively transport aerosols from the boundary layer to higher elevations in the free troposphere, long-range transport is expected to have a greater role in determining the number of aerosols in the free

Aerosol Robotic Network (AERONET) Free Tropospheric AOD
Aerosols in the entire atmospheric column, between the sea surface and the top of the atmosphere, contribute to the AOD derived from MODIS observations.Quantifying the relative contributions of the boundary-layer and free-tropospheric aerosols to the total column AOD is not possible with passive remote sensing from the space and ground.Without deep convection that can effectively transport aerosols from the boundary layer to higher elevations in the free troposphere, long-range transport is expected to have a greater role in determining the number of aerosols in the free troposphere.The AERONET observation site at mount Mauna Loa is mostly above the boundary layer and, therefore, occasionally affected by long-range transport of anthropogenic aerosols.
Figure 7 shows the annual cycle of AOD that was measured at the Mauna Loa site during the period 2003-2016.The most evident feature is the springtime AOD, which is about 2-3 times higher than in the rest of the year.Furthermore, the figure indicates that the averaged AOD derived from observations at Mauna Loa make up about 20% of the total column AOD (Figure 6b) for the central North Pacific in summer, fall and winter.troposphere.The AERONET observation site at mount Mauna Loa is mostly above the boundary layer and, therefore, occasionally affected by long-range transport of anthropogenic aerosols.
Figure 7 shows the annual cycle of AOD that was measured at the Mauna Loa site during the period 2003-2016.The most evident feature is the springtime AOD, which is about 2-3 times higher than in the rest of the year.Furthermore, the figure indicates that the averaged AOD derived from observations at Mauna Loa make up about 20% of the total column AOD (Figure 6b) for the central North Pacific in summer, fall and winter.

Air Mass Transport
Analyzing a large number of back-trajectories helps in evaluating influences of long-range transport of aerosols on the observed annual cycle of the MODIS and AERONET AOD.In Sections 3.2 and 3.3 it was demonstrated that AOD peaks in spring-in the free troposphere as well as in the boundary-layer during both low and high wind speeds.Figure 8

Air Mass Transport
Analyzing a large number of back-trajectories helps in evaluating influences of long-range transport of aerosols on the observed annual cycle of the MODIS and AERONET AOD.In Sections 3.2 and 3.3 it was demonstrated that AOD peaks in spring-in the free troposphere as well as in the boundary-layer during both low and high wind speeds.Figure 8 displays seasonal density of the FLEXTRA 7-day back-trajectories for the period 2003-2016.The trajectory density reflects the probability of a trajectory crossing any grid cell (regardless of the trajectory altitude) and highlights the source regions of the air masses in each season.The figure suggests that air masses approach Mauna Loa from the north-eastern sector mainly in summer.It also shows that the area covered by the 7-day trajectories is larger and extends further to the west in winter and spring.Thus, it is more likely that the air masses in these seasons to a larger extent arrived from southern and central Asia before reaching the Hawaii region.

Power-Law Relationship between MODIS AOD and Surface-Wind Speed
The relationship between AOD and surface-wind speed is often represented by a power-law function.Figure 9 shows annual mean MODIS-adjusted AOD (black circles) and corresponding one standard deviation for the period 2003-2016, although excluding the spring season due to substantially higher AOD values compared to the other seasons.The derived power law relationship, AOD = 0.001 × U 1.86 + 0.05, indicate a relatively strong wind dependency.The constant zero offset term (AOD ~ 0.05) reflects a marine background situation, where we assume the presence mainly of secondary aerosols such as ammonium sulfate that are formed independently of the local wind speed.However, the relatively long turn-over time for the submicron sea-salt particles (27) means that they may contribute to AOD also during situations with low wind speeds.Even so, we believe a minor contribution from these aerosols, since the local wind speed is probably a good substitute for the Lagrangian wind speed over the accumulation-mode lifetime (10).

Comparison with Previous Studies
Figure 10 and Table 2 presents relationships between AOD and surface wind speed from various studies.These relationships established for marine areas deviate substantially.Some of the studies suggest a linear relationship, (e.g., [35,55]), others a power-law relation ( [36,37,56] current study) or an exponential relationship [38].Despite the large differences, almost all relationships show an

Power-Law Relationship between MODIS AOD and Surface-Wind Speed
The relationship between AOD and surface-wind speed is often represented by a power-law function.Figure 9 shows annual mean MODIS-adjusted AOD (black circles) and corresponding one standard deviation for the period 2003-2016, although excluding the spring season due to substantially higher AOD values compared to the other seasons.The derived power law relationship, AOD = 0.001 × U 1.86 + 0.05, indicate a relatively strong wind dependency.The constant zero offset term (AOD ~0.05) reflects a marine background situation, where we assume the presence mainly of secondary aerosols such as ammonium sulfate that are formed independently of the local wind speed.However, the relatively long turn-over time for the submicron sea-salt particles (27) means that they may contribute to AOD also during situations with low wind speeds.Even so, we believe a minor contribution from these aerosols, since the local wind speed is probably a good substitute for the Lagrangian wind speed over the accumulation-mode lifetime (10).

Comparison with Previous Studies
Figure 10 and Table 2 presents relationships between AOD and surface wind speed from various studies.These relationships established for marine areas deviate substantially.Some of the studies suggest a linear relationship, (e.g., [35,55]), others a power-law relation ( [36,37,56] current study) or an exponential relationship [38].Despite the large differences, almost all relationships show an increase in AOD, roughly with a factor of 2 over the wind-speed range 0-10 m/s.For several of the non-linear relationships, AOD begins to increase at a wind speed of about 4 m s −1 .The main differences between the relationships presented in Figure 10 and Table 2 are likely caused by differences in one or several of the following factors used in the studies: methodology, range of wind speeds, location and time period.AOD 500 = 4.9 × 10 −5 × U 3 − 3.7 × 10 −5 × U 2 + 0.017 3rd degree polynomial n/a November-December 1995 SP g [57] a NP = North Pacific, b S = Simulation, c MH = Mace Head, d NA = North-east Atlantic, e MI = Midway, Island, f AS = Arabien Sea, g SP = South Pacific, h MBL = marine boundary layer.
Concerning power-law relationships, based AOD derived from Precision Filter Radiometer (PFR) measurements at Mace Head, Ireland, Mulcahy et al. [36] found a power exponent of ~2.2 (R 2 = 0.97).Glantz et al. [37] found a power exponent of 2.3 (R 2 = 0.98) for AOD derived from SeaWiFS observations, while O'Dowd et al. [56] reported a value ~2.5 based on MODIS AOD.These three studies are based on different investigation periods: Glantz et al. [37] used data only for September 2001, O'Dowd et al. [56] only for December 2006, and Mulcahy et al. [36] carried out measurements during winter in the period 2002-2004.In comparison to the previous studies the present study is actually more comprehensive as it covers three seasons and many years.Even so, the AOD dependence on wind speed estimated in the present study is in best agreement with the parameterization derived by Mulcahy et al. [36].
Kiliyanpilakkil and Meskhidze [35] obtained logistic regression based on five years of Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) measurements onboard the CALIPSO satellite.In the latter study, they used only single-layer AODs below an altitude of 2 km over 15 different ocean regions.AOD values substantially lower than 0.05 were found for low and intermediate wind speeds, which is likely explained by the method used: based on CALIOP data they were able to select purely the maritime wind-induced component of AOD.This is in line with the results obtained with the COSMO-ART model by Lundgren et al. (2013) where only the sea salt component was included in the simulation.Furthermore, Shinozuka et al. [57] found a cubic polynomial regression fit between AOD and surface-wind speed based on aircraft measurements over the south-western Pacific Ocean, south of Australia in November-December 1995.In their study, the calculation of total AOD over the altitudes up to 2000 m was based on ambient scattering coefficient derived from optical particle counter (OPC) measurements.The relationship estimated from this study is also established only based on sea-salt contribution to AOD.AOD500 = 0.00055 × U 2.195    Concerning power-law relationships, based AOD derived from Precision Filter Radiometer (PFR) measurements at Mace Head, Ireland, Mulcahy et al. [36] found a power exponent of ~2.2 (R 2 = 0.97).Glantz et al. [37] found a power exponent of 2.3 (R 2 = 0.98) for AOD derived from SeaWiFS observations, while O'Dowd et al. [56] reported a value ~2.5 based on MODIS AOD.These three studies are based on different investigation periods: Glantz et al. [37] used data only for September 2001, O'Dowd et al. [56] only for December 2006, and Mulcahy et al. [36] carried out measurements during winter in the period 2002-2004.In comparison to the previous studies the present study is actually more comprehensive as it covers three seasons and many years.Even so, the AOD dependence on wind speed estimated in the present study is in best agreement with the parameterization derived by Mulcahy et al. [36].2.
Huang et al. [34], O'Dowd et al. [56] and Satheesh et al. [38] found, instead, relatively high AOD values (around 0.1) at low wind speeds.When investigating AOD over the South Pacific, using data from the Advanced Along-Track Scanning Radiometer (AATSR) on board the European Space Agency (ESA) ENVISAT satellite, Huang et al. [34] found a linear relationship.However, a warning of a possible bias in the implemented algorithm was also reported.O'Dowd et al. [56] used MODIS-derived AOD over the North Pacific and explained elevated AOD levels at low wind speeds with increased aerosol scattering due to high relative humidity and/or secondary marine aerosol production and/or influence of long-range transport of anthropogenic aerosols.In addition, Satheesh et al. [38] analyzed MODIS-derived AOD over the Arabian Sea for the period 2001-2003 and noted a substantial contribution of mineral dust from land which, indeed, contributed to the high AODs.
Furthermore, for low wind speeds (0-4 m/s) the present AOD values are similar to ground-based measurements carried out by Smirnov et al. [55].The linear increase of AOD with higher wind speed in the latter study can, however, be questionable, since sea-spray aerosol flux (or rather white cap coverage) dependency on surface wind speed is suggested to be proportional to U 3.5  [31] or U 3.41  (e.g., [51]) or U 3.75  [58].

Comparison between MODIS Collections 5 and 6
MODIS aerosol products are organized according to different collections consisting of data that were generated by different versions of the retrieval algorithms.When comparing Collection 5 (C5) and C6, a significant difference in the AOD-wind speed relationship appears.Figure 11 shows C5 and C6 AOD versus surface wind speed, averaged with respect to autumn (September, October and November) of the period 2003-2012.AOD appears to be more sensitive to surface-wind speed in C5 and it does not show the decrease in AOD at about 9 m/s as with C6 and all the relationships obtained in the previous studies.One likely reason for this is that the modified algorithm overestimate the effect of foam on the water surface at the higher wind speeds.Overall, the present results are consistent with the relationship found between MODIS AOD and wind speed by Levi et al. [43], who highlighted that the C6 algorithm retrieves lower (higher) values of AOD than the C5 when wind speed is higher (lower) than 6 m/s.Thus, the actual relationship between AOD and surface wind speed is likely somewhat stronger that the one obtained here with the C6 algorithm.On the other hand, since MODIS AOD is in this study related to surface wind speed at 2 m, the relationship will be weaker when transforming the relationship to wind speeds at 10 m.The latter can be performed by applying the wind-profile power law: where U is the wind speed (m/s) at height z (m); U r is the known wind speed at a reference height z r ; and α is an empirically derived coefficient, which is ~0.11 for open water.

Comparison between MODIS Collections 5 and 6
MODIS aerosol products are organized according to different collections consisting of data that were generated by different versions of the retrieval algorithms.When comparing Collection 5 (C5) and C6, a significant difference in the AOD-wind speed relationship appears.Figure 11 shows C5 and C6 AOD versus surface wind speed, averaged with respect to autumn (September, October and November) of the period 2003-2012.AOD appears to be more sensitive to surface-wind speed in C5 and it does not show the decrease in AOD at about 9 m/s as with C6 and all the relationships obtained in the previous studies.One likely reason for this is that the modified algorithm overestimate the effect of foam on the water surface at the higher wind speeds.Overall, the present results are consistent with the relationship found between MODIS AOD and wind speed by Levi et al. [43], who highlighted that the C6 algorithm retrieves lower (higher) values of AOD than the C5 when wind speed is higher (lower) than 6 m/s.Thus, the actual relationship between AOD and surface wind speed is likely somewhat stronger that the one obtained here with the C6 algorithm.On the other hand, since MODIS AOD is in this study related to surface wind speed at 2 m, the relationship will be weaker when transforming the relationship to wind speeds at 10 m.The latter can be performed by applying the wind-profile power law: where U is the wind speed (m/s) at height z (m); Ur is the known wind speed at a reference height zr; and α is an empirically derived coefficient, which is ~0.11 for open water.

Discussion
AOD derived from MODIS Aqua observations over the central North Pacific has been validated against AOD obtained from sun-photometer measurements carried out at three marine ground-based stations and during several ship cruises.MODIS 555 nm AOD was found to vary within the expected uncertainties of the MODIS retrievals over the ground-based sites and ship tracks for 71% and 75%, respectively, of the compared cases.Figure 2 shows that many of the MODIS AOD values are substantially

Discussion
AOD derived from MODIS Aqua observations over the central North Pacific has been validated against AOD obtained from sun-photometer measurements carried out at three marine ground-based stations and during several ship cruises.MODIS 555 nm AOD was found to vary within the expected uncertainties of the MODIS retrievals over the ground-based sites and ship tracks for 71% and 75%, respectively, of the compared cases.Figure 2 shows that many of the MODIS AOD values are substantially higher than the expected error envelopes.The largest deviation is found in the comparison between MODIS and MAN AOD, where MODIS AOD in some cases exceeds AERONET AOD by almost a factor of 2. When investigating the largest deviation in Figure 2b we found that the AOD values outside the estimated uncertainties ranges belong to nine ship cruises of the 15 in Table 1.Thus, there is not a ship track that stands out concerning the deviation.Another possible cause for the deviation is cloud contamination.Several studies [28,44,[59][60][61] found that MODIS tends to overestimate AOD compared to AERONET over ocean areas under cloudy conditions.Zhang et al. [62] analyzed cloud artifacts over remote oceans and found that cloud contamination explains 60-90% of the correlation between MODIS AOD and cloud fraction.Similar deviation as the one found here in MODIS AOD derived from observations over the current coastal sites has also been reported in previous studies [43,[63][64][65][66]. Mélin et al. [64] suggest that further improvements of atmospheric correction in coastal waters may require an additional level of detail to describe the complex aerosol mixtures more accurately, and specific water bio-optical models for an appropriate representation of the boundary condition.
The radiance leaving the ocean surface may also cause overestimation in MODIS AOD.An extreme situation is sun glint phenomenon that occurs when sunlight reflects off the surface of water at the same angle than the satellite sensor views it.Although a screening of this effect is included in the MODIS algorithm, in any case it may induce uncertainties in the final product.White foam (whitecap) formations, especially during higher wind-speed conditions, may also induce uncertainties in MODIS-derived AOD.Levi et al. [43] state that despite the improvement in the MODIS Collection 6 algorithm, the errors caused by surface brightness may still have a large impact on the AOD estimated, particularly for cases with relatively low AOD values such as in the marine background atmosphere.Even so, the variability in AOD for the different wind speeds, shown in Figure 3, is in any case relatively constant around the mean values.
Substantially higher MODIS and AERONET (Mauna Loa) AODs were found in the present study for spring compared to the other seasons investigated.These results seem to be consistent with previous studies, e.g., [67][68][69][70].They all found that trans-Pacific transport of air pollutants from China is most efficient in spring due to active cyclonic activity.Frequent convection in spring lifts particulate matter from the Asian continent to the free troposphere, where prevailing westerlies rapidly transport aerosols across the ocean.When continental air reaches the high-pressure system in the North Pacific, it subsides in the middle troposphere where aerosol residence time is much longer than in the boundary layer.A number of studies [71][72][73][74][75] have also found mineral dust transport from Asia over the North Pacific in spring.In addition, the present back-trajectory analysis suggests a frequent transport of air masses from Asia to the investigation area also in winter, actually most obvious during this season (Figure 8).However, the vertical transport of aerosols from the Asian boundary layer to the free troposphere in winter is expected to be rather limited due to suppressed convection.This means that the transport of aerosols over the North Pacific in winter mainly takes place in the boundary layer, where efficient wet deposition of the aerosols is expected [71].
The AOD dependence on wind speed estimated in the current study is in reasonable agreement with the parameterization derived by Mulcahy et al. [36], although the latter relationship is slightly stronger.AOD derived from the PFR measurements is assumed to be more accurate than AOD derived from MODIS observations.Furthermore, the analysis by Mulcahy et al. [36] was done for selected days that fall into the category of "stable wind conditions", which are defined as days with the lowest variability in daily wind speeds (less than 2 m/s).At the end, only 14 days were left for analyzing the AOD-wind speed dependence.It is remarkable that the present averaged AOD-wind speed relationship derived from MODIS Aqua observations over the central North Pacific over 14 years are relatively close to the relationship found at Mace Head, since the two studies focused on different oceans with large differences in wind patterns (trade versus westerly winds) and sea-surface temperatures.Note however that the relationship found at Mace Head by Mulcahy et al. (2008) is obtained for substantially higher wind speeds (obtained with 2-min time resolution) than the present relationship.Even so, the majority of the MODIS AOD values over the central North Pacific were derived at a wind speed of around 6 m s −1 and there are probably not many places on Earth that have a mean wind speed substantially higher than 10 m s −1 .Although the surface wind speed is the most important factor for determining sea spray aerosol production over ice-free water surfaces, the water temperature plays an important role as well [11,16,17].Therefore, using the present relationship between AOD and surface-wind speed over other regions on Earth, where the surface-water temperature is substantially lower, induces uncertainties in the estimates of AOD

Conclusions
This study shows relatively strong relationships between MODIS mean AOD and surface-wind speed over the central North Pacific for marine background conditions in summer, fall and winter in the period 2003-2016.In addition, the enhanced AOD, due to increased wind speeds, is mainly associated with coarse-mode aerosols.Thus, AOD associated with particles larger than 0.25 µm in radius increases substantially for wind speeds higher than about 3.5 m/s, while only a minor increase is found for the fine-mode aerosols.This is consistent with a fine mode of secondary aerosols that form independent of the wind speed.
To summarize, the current study shows a distinct increase in MODIS AOD for higher wind speeds.The present results show that MODIS Aqua retrievals of AOD in the marine atmosphere agree reasonably well with ground-based remote-sensing observations.To our knowledge, the established parameterization is derived for the longest period of MODIS Aqua observations over the ocean.It can serve as a promising reference for future investigations of the AOD-wind speed relationship.

Atmosphere 2018, 9 , 19 Figure 1 .
Figure 1.The investigation area (12° N-35° N, 145° W-168° W) in the central North Pacific is shown in (b) and denoted with a black box in (a).Aerosol optical depths (AODs) were derived based on MODIS Aqua observations on 18 September 2014 (right).Aerosol Robotic Network (AERONET) ground-based data were obtained from the stations at Coconut Island (21.4°N, 157.8°W) and Lanai (20.7°N, 156.9°W) on Hawaii, as well as from Midway (28.2°N, 177.4°W) west of the investigation area (orange solid circle in (a)).

Figure 1 .
Figure 1.The investigation area (12 • N-35 • N, 145 • W-168 • W) in the central North Pacific is shown in (b) and denoted with a black box in (a).Aerosol optical depths (AODs) were derived based on MODIS Aqua observations on 18 September 2014 (right).Aerosol Robotic Network (AERONET) ground-based data were obtained from the stations at Coconut Island (21.4 • N, 157.8 • W) and Lanai (20.7 • N, 156.9 • W) on Hawaii, as well as from Midway (28.2 • N, 177.4 • W) west of the investigation area (orange solid circle in (a)).

Figure 2 .
Figure 2. (a) Comparison between Moderate Resolution Imaging Spectroradiometer (MODIS) and AERONET AOD; and (b) MODIS and MAN AOD at 555 nm.The solid and dashed lines represent the 1-to-1 line and estimated uncertainties of the MODIS retrievals over ocean (ΔAOD = +0.04+ 0.1•AOD, ΔAOD = −0.02− 0.1•AOD), respectively.Error bars correspond to the MODIS/AERONET and MODIS/MAN AODs' comparison and denote one standard deviation.

Figure 3 .
Figure 3. MODIS-adjusted AOD (grey dots) on pixel basis versus surface wind speed for summer 2007.The grey solid circles denote the median AOD and corresponding median absolute deviation for the binned data.The black solid and red dashed lines denote a power fit and corresponding 95% confidence interval, respectively.The R 2 and RMSE values were calculated based on the power fit of the median values.

Figure 3 .
Figure 3. MODIS-adjusted AOD (grey dots) on pixel basis versus surface wind speed for summer 2007.The grey solid circles denote the median AOD and corresponding median absolute deviation for the binned data.The black solid and red dashed lines denote a power fit and corresponding 95% confidence interval, respectively.The R 2 and RMSE values were calculated based on the power fit of the median values.Atmosphere 2018, 9, x FOR PEER REVIEW 7 of 19

Figure 4 .
Figure 4. MODIS-adjusted median AOD (555 nm) versus surface wind speed, subdivided with respect to season, for each year of the period 2003-2016.The orange solid line in each figure denote MODIS mean AOD, averaged over the investigation period.

Figure 4 .
Figure 4. MODIS-adjusted median AOD (555 nm) versus surface wind speed, subdivided with respect to season, for each year of the period 2003-2016.The orange solid line in each figure denote MODIS mean AOD, averaged over the investigation period.

Figure
Figure 5a,b show seasonal wind speed probability density function for MODIS derived AOD, normalized according to the total number of derived AOD corresponding to a season and year, respectively.The majority of the AOD values were derived at wind speeds in the range between about 3.5 m/s and 8 m/s; thus, rather few AOD values correspond to lower and higher wind speeds.

Figure 5 .
Figure 5. Frequency distributions of all AOD, derived from observations over the Northern Pacific during the period 2003-2016, which have been normalized with respect to the total number of derived AOD corresponding to the (a) season and (b) year.

Figure 5 .
Figure 5. Frequency distributions of all AOD, derived from observations over the Northern Pacific during the period 2003-2016, which have been normalized with respect to the total number of derived AOD corresponding to the (a) season and (b) year.

Figure 6 .
Figure 6.(a) Mean Ångström exponent, derived from MODIS Aqua observations over the period 2003-2016, versus surface wind speed and subdivided according to season; (b) MODIS Aqua fine mode adjusted median AOD (thin lines) and total adjusted median AOD (bold lines) obtained for different seasons of the period 2003-2016 versus surface wind speed.

Figure 6 .
Figure 6.(a) Mean Ångström exponent, derived from MODIS Aqua observations over the period 2003-2016, versus surface wind speed and subdivided according to season; (b) MODIS Aqua fine mode adjusted median AOD (thin lines) and total adjusted median AOD (bold lines) obtained for different seasons of the period 2003-2016 versus surface wind speed.

Figure 7 .
Figure 7. Bold and thin solid line denote daily mean AOD and corresponding one standard deviation, respectively, obtained from yearly time series observed at the AERONET Mauna Loa station during the period 2003-2016.
displays seasonal density of the FLEXTRA 7-day back-trajectories for the period 2003-2016.The trajectory density reflects the probability of a trajectory crossing any grid cell (regardless of the trajectory altitude) and highlights

Figure 7 .
Figure 7. Bold and thin solid line denote daily mean AOD and corresponding one standard deviation, respectively, obtained from yearly time series observed at the AERONET Mauna Loa station during the period 2003-2016.

Atmosphere 2018, 9 , 19 Figure 8 .
Figure 8. FLEXible TRAjectories (FLEXTRA) ECMWF 7-day back-trajectories starting at Mauna Loa (19.54°N, 155.58°W), Hawaii, at a height of 3000 m above sea level, subdivided with respect to (a) winter; (b) spring; (c) summer and (d) fall.Color shading denotes the probability of occurrence for air masses passing over a particular grid-cell before they arrived at Mauna Loa (see explanation in the text).

Figure 8 .
Figure 8. FLEXible TRAjectories (FLEXTRA) ECMWF 7-day back-trajectories starting at Mauna Loa (19.54 • N, 155.58 • W), Hawaii, at a height of 3000 m above sea level, subdivided with respect to (a) winter; (b) spring; (c) summer and (d) fall.Color shading denotes the probability of occurrence for air masses passing over a particular grid-cell before they arrived at Mauna Loa (see explanation in the text).

AOD500
North Pacific, b S = Simulation, c MH = Mace Head, d NA = North-east Atlantic, e MI = Midway, Island, f AS = Arabien Sea, g SP = South Pacific, h MBL = marine boundary layer.

Figure 9 .
Figure 9. MODIS-adjusted mean AOD (black solid circle) with corresponding one standard deviation, averaged with respect to the median annual curves obtained for summer, fall and winter in Figure 4.The blue solid line is a power-law fit of the mean AOD values.

Figure 9 .
Figure 9. MODIS-adjusted mean AOD (black solid circle) with corresponding one standard deviation, averaged with respect to the median annual curves obtained for summer, fall and winter in Figure 4.The blue solid line is a power-law fit of the mean AOD values.Atmosphere 2018, 9, x FOR PEER REVIEW 12 of 19

Figure 10 .
Figure 10.The present relationship between adjusted AOD and surface-wind speed in comparison to relationships obtained in previous studies.More information is given in Table 2. Kiliyanpilakkil and Meskhidze [35] obtained logistic regression based on five years of Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) measurements onboard the CALIPSO satellite.In the latter study, they used only single-layer AODs below an altitude of 2 km over 15 different ocean regions.AOD values substantially lower than 0.05 were found for low and intermediate wind speeds, which is likely explained by the method used: based on CALIOP data they were able to select purely the maritime wind-induced component of AOD.This is in line with the results obtained with the COSMO-ART model by Lundgren et al. (2013) where only the sea salt component was included in the simulation.Furthermore, Shinozuka et al. [57] found a cubic polynomial regression fit between AOD and surface-wind speed based on aircraft measurements over the south-western Pacific Ocean, south of Australia in November-December 1995.In their study, the calculation of total AOD over the altitudes up to 2000 m was based on ambient scattering coefficient derived from optical particle

Figure 10 .
Figure 10.The present relationship between adjusted AOD and surface-wind speed in comparison to relationships obtained in previous studies.More information is given in Table2.

Figure 11 .
Figure 11.MODIS C5 and C6 AOD (black dashed and blue solid lines, respectively) versus surface wind speed, averaged with respect to September, October and November for the years 2003-2012.

Figure 11 .
Figure 11.MODIS C5 and C6 AOD (black dashed and blue solid lines, respectively) versus surface wind speed, averaged with respect to September, October and November for the years 2003-2012.

Table 1 .
List of Maritime Aerosol Network (MAN) cruises, investigation areas, and number of measurement days.

Table 2 .
Current and previous parametrizations of AOD as a function of surface wind speed.