Next Article in Journal
Implementing Machine Learning Algorithms to Predict Particulate Matter (PM2.5): A Case Study in the Paso del Norte Region
Previous Article in Journal
PM2.5-Related Health Risk during Chinese Spring Festival in Taizhou, Zhejiang: The Health Impacts of COVID-19 Lockdown
Previous Article in Special Issue
What Made the Sustained Intensification of Tropical Cyclone Fani in the Bay of Bengal? An Investigation Using Coupled Atmosphere–Ocean Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Role of Phytoplankton Biomacromolecules in Controlling Ocean Surface Roughness

by
Amadini Jayasinghe
1,*,
Scott Elliott
1,
Georgina A. Gibson
2 and
Douglas Vandemark
3
1
Los Alamos National Laboratory, Los Alamos, NM 87545, USA
2
International Arctic Research Center, University of Alaska, Fairbanks, AK 99775, USA
3
Ocean Process Analysis Laboratory, University of New Hampshire, Durham, NH 03824, USA
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(12), 2101; https://doi.org/10.3390/atmos13122101
Submission received: 1 November 2022 / Revised: 2 December 2022 / Accepted: 7 December 2022 / Published: 14 December 2022

Abstract

:
Satellite altimetric data routinely map sea surface topography by measuring the ocean return signal. One source of altimeter measurement contamination occurs when the radar ocean backscatter becomes unusually large, a situation termed a Sigma-0 bloom. Past research suggests Sigma-0 blooms are associated with weak wind and natural surface slick conditions where capillary waves at the air–sea interface are suppressed. To date, no explicit connection between these conditions and Sigma-0 bloom presence has been provided. Using a series of simplified equations, our reduced model determines capillary wave heights from estimates of planktonic carbon concentrations and regional wind speed. Our results suggest that the radar signal reflection increases as capillary wave height decreases. This relationship depends on surfactant concentration, surfactant composition, and wind speed. Model sensitivity analysis indicates that the interface reflectivity depends on biological activity and wind speed. Our proposed simplified model provides a method to identify potential Sigma-0 bloom regions. We conclude that because of the demonstrated impact of biological surfactants on ocean roughness, it is necessary to consider the biological activity, i.e., phytoplankton bloom events, when interpreting signals from radar altimetry and when developing ocean hydrology models.

1. Introduction

1.1. Measuring Sea Surface Height with Satellites

Changes in the surface topography of the ocean are now routinely used to monitor ocean circulation dynamics and sea level rise. Satellite altimeters provide these data by measuring the return signal from a radar pulse [1]. Radar signals in the microwave range are transmitted toward the ocean, and the time taken for them to travel from the satellite antenna to the water surface and back is calculated. The return signal (waveform) also provides a time evolution of reflected power, from which the range parameter (i.e., satellite to surface distance) is extracted. The waveform shape and amplitude depend strongly on the local roughness of the ocean surface. Radar altimetry thus also provides an estimate for the backscattering cross-section (Sigma-0, σ0), or the proportion of the radiation that is returned in the direction of the incident signal. The strength of the return depends on the roughness of the surface and the wavelength of incidence [2].
At times, across the global ocean, regions have been found to exhibit unusually high values for Sigma-0. These localized phenomena are referred to as Sigma-0 blooms. Altimeter Sigma-0 blooms, and their high backscattering tendency, indicate that the sea surface has become smooth and is thus acting more as a specular reflector than a scattering boundary with respect to the incident electromagnetic radar signal. These Sigma-0 bloom cases can be thought of as a contamination of the radar return signal because they can cause a breakdown in the basic assumptions required to estimate surface height from the return power signal [3,4]. The phenomenon impacts a significant portion of satellite-derived wave height data, and thus the overall success of the radar altimetry method. During the historic TOPEX/Poseidon joint altimeter mission approximately 5–10% of the altimeter data were affected by the presence of Sigma-0 blooms. Yet, this potential error source is often not considered during analysis [3,5]. Efforts to explore mechanisms and reasons behind the reflected Sigma-0 blooms might, at the very least, help to mitigate such contamination in sea surface height mapping [6].

1.2. Sigma-0 Blooms and Capillary Waves

A likely general explanation for Sigma-0 bloom events is that they are associated with seawater surfaces with damped capillary waves, i.e., smooth interfaces [7]. Even the largest planetary ocean wind waves start as tiny capillary waves in the ocean microlayer, the upper 1000 micrometers of the water column. These waves, of very short wavelength, arise when surface wind stress induces small vortices in a completely flat sea. As these small capillary ripples propagate and develop turbulent eddies, the waves grow and become gravity waves. Capillary waves act as the “fastening system” of the sea-air interface, providing the opportunity for wind to exert friction. As such, they govern momentum transfer from atmospheric forcing. However, what is most critical to this study is the fact that capillary wave dynamics are also controlled by the surface tension of the ocean’s microlayer.

1.3. Biological Impact on Surface Tension

The existence of organic surfactants at the sea-air interface modulates physical properties such as surface tension and surface pressure. Small wave structures are supported by turbulence and microscopic energy flow in an interplay of the rising-spreading motions, which generate area and the push-pull on water molecules and organics. Van der Waals interactions among hydrophobic tails of the surfactant molecules and the Coulombic attraction, associated with the embedded head groups, balance to determine the surface tension [8,9,10,11].
The height of capillary waves and surface ripples is suppressed through turbulent energy losses into a planar Van der Waals sink, derived from intermolecular forcing that reduces the electromagnetic attraction of the water molecules and thus the surface tension [8,11]. Surface tension, and hence ocean roughness, are controlled to a great extent by amphiphilic molecules that partially coat the oceanic surface [8,12].
The biomacromolecules that comprise the ocean surface films are the polymers of life, i.e., carbohydrates, lipids, proteins, and nucleic acids [11]. Chained carbon compounds from detritus are a major source of such surfactants, as are proteins, generated both as the catalysts of phytoplanktonic metabolic activity and as structural components [9,13] within the phytoplankton cells. When effects are extensive and visible to the naked eye, the thin films of bio-macromolecules are referred to as slicks. In addition to slicks from relatively fresh organic material, the presence of oil spills and related organic carbon compounds has also been found to partially alter ripple size at low wind speed [7,14]. Elliott et al. (2018) hypothesized that ripple suppression occurs strongly and preferentially in regions where slicks develop [8]. Thus, a possible explanation for the highly reflective Sigma-0 regions is the presence of biogenic surface slicks [3].

1.4. The Role of Phytoplankton in Sigma-0 Blooms

When conceived, the radar altimeter term “Sigma-0 bloom” was only meant as an analogy to phytoplankton blooms; however, here we demonstrate the direct connections between these biological and physical features. The bio-macromolecules, including fatty acids, amino acids, and polysaccharide-protein complexes, which enrich the microlayer, are often amphiphilic structures and are mainly produced during food web cycling—primarily by phytoplankton producers and zooplankton grazers. They are generated as metabolic by-products released as individual cells are broken during consumption [9,10,11,15]. A significant fraction of marine dissolved organic matter, which primarily originates from phytoplankton, contributes to global surface activity [10]. There are additional contributions to the surface bio-macromolecules from atmospheric deposition and river runoff, although these are mainly restricted to coastal regions.

2. Analysis and Results

We start by presenting our reduced model approach (Section 2.1) and then discuss the impact of surface slicks and wind speed on wave geography (Section 2.2), demonstrate the model utility by computing surface pressure for the North Atlantic from lower trophic level dynamics (Section 2.3), show how roughness length and capillary wave height varies spatially (Section 2.4), and how this relates to Sigma-0 frequency (Section 2.5). We then explore the sensitivity of surfactant concentration and surface pressure threshold to capillary wave height (Section 2.6) and examine the impact of surfactant concentration on capillary wave heights (Section 2.7) and the impact of surface pressure threshold on capillary wave height (Section 2.8). Finally, we assess model performance using regional ship track measurements (Section 2.9).

2.1. Reduced Model Approach

Previously established relationships can be used to explore potential linkages between biological surfactants at the ocean’s surface and global Sigma-0 blooms. By down-scaling extremely complex, dynamic, and pervasive air–sea interactions linking aqueous and atmospheric boundary layers, we can derive compact geochemical steady-state equations for predicting ripple size from surface biology. We assume that ripple size or elemental height (E0, Equation (1)) is based on the proportional relationship established in a micrometeorological study [16,17], which relates it to the roughness length (z0). This key parameter connects capillary wave height (E0) to the chemical adsorption of the surfactants at the air–sea interface:
E 0 = 10 z 0
The roughness length can be determined empirically (Equation (2), [12]). This formula can be used to determine roughness length values for different surfactant concentrations (z0,surf). When there is absolutely no adsorption at the sea-air interface, a condition we assume occurs at a surfactant concentration of zero, we represent the roughness length value as z0,pure. The roughness length value for the maximum adsorption surfactant concentration is denoted by (z0,surf.100%). The term Log10(z0) is equal to [Log10(z0,pure) − Log10(z0,surf.100%)], the logarithmic difference between the maximum and minimum roughness length values [12].
L o g 10 z 0 , s u r f = L o g 10 z 0 , p u r e π π t h r e s h Δ L o g 10 z 0
This equation also uses the interfacial surface pressure (π), i.e., the difference in interfacial tension of the microlayer with and without surfactants [11,12]. Due to the adsorption of the surfactants at the interface, a change to the surface tension force occurs, leading to the tangential generation of surface pressure. The proportionality between the surfactant concentration and surface pressure enables us to use surface pressure (π) as the parameter that represents the effect of the surfactant concentration on the roughness length. A threshold surface pressure (πthresh) is an arbitrary maximum value that estimates the surfactant influence on capillary wave geometries (height). Threshold surface pressure is assumed to be the arbitrary 2D surface pressure at which we will obtain a saturated roughness length value where no more reduction will occur [12,18].
π = i = 1 2 θ i . π i , m a x
Where πi,max represents the maximum 2D surface pressure value of ith surface-active macromolecule (i.e., i = 1 or 2, 1 = protein, 2 = lipid, etc.) and θi represents the fraction of interfacial coverage by each surface-active macromolecule. While numerous heteropolycondensates could be considered in these calculations, for simplicity, here we restrict our analysis to proteins and lipids, arguably the two most important in terms of surface activity. The fraction of interfacial coverage by each macromolecule at the sea-air interface can be calculated using the empirical adsorption isotherm Equation (4) [12], which accounts for the surface adsorption concentration of each macromolecule type and its surface adsorption equilibrium.
θ i = ( 1 / K i , ) n i . ( a i C i ) n i 1 + i = 1 2 ( 1 / K i   ) n i . ( a i C i ) n i   w h e r e   i = 1 p r o t , 2 l i p
Where Ci is the micromolar carbon atom concentration of the ith macromolecule, Ki is the half-saturation micromolar carbon atom concentration for the ith macromolecule (μM carbon), and ni (nprot = 1, nlip = 8) is the effective coverage power, a parameter unique to the macromolecular type which defines the shape of the adsorption curve. Finally, ai is the uncertainty parameter that allows exploration of the uncertainties created by the absolute concentration difference of measured and calculated concentrations. The micromolar carbon atom concentration of each of the macromolecules (Ci) can be calculated from the observed planktonic carbon concentration (Equation (5)).
C i = g Z P c K + P c 1 γ P i % τ i
Where Pc is the planktonic carbon atom concentration, Z is the zooplankton concentration, g represents zooplankton grazing on the phytoplankton, K is the half-saturation concentration for ingestion, and γ is the zooplankton assimilation efficiency. τi represents the lifetime of the ith macromolecule, and Pi% is the fraction (0–1) of the ith macromolecule in a typical planktonic cell (inside the cytosol). Given Equation (5), the zonal protein and lipid concentrations at a steady state are readily calculated. With this set of five equations, we can now relate capillary wave height to the ecosystem dynamics of a region. In Section 2.3, we illustrate the utility of this approach through a case study of the North Atlantic. A summary of symbols, descriptions, unit values, and sources is provided in Table 1.

2.2. Impact of Surface Slicks and Wind Speed on Wave Geography

In addition to the macromolecules in the surface slick, capillary wave height is also dependent on wind velocity. In order to understand the dual impact of the surface pressure and the wind velocity on the capillary wave height, we represent the relationship visually (Figure 1) using the dynamics outlined in Equations (1) and (2). Assuming a surface pressure ratio (π/πthresh) whose influence increases incrementally to reach a threshold value of 1.0 due to the presence of a surface slick and applying Equations (1) and (2), we show that while element height increases with increasing wind speed, the length is greatly reduced by increases in the surface pressure ratio (Figure 1). When slicking is introduced to varying degrees, the capillary wave height is demonstrably reduced for all wind speeds but has a maximal effect at intermediate velocities. At π/πthresh = 0 there is the assumption of no surfactants to dampen the surface roughness, and thus the capillary wave height is only dependent on the wind speed. Conversely, at a π/πthresh = 1, when the microlayer is assumed completely covered in surfactant, a minimum wind speed of 5 ms−1 is required for the growth of any capillary waves. Figure 1 reflects the relationship between capillary wave height (E0) vs. windspeed at π/πthresh = 0 and unity, representing absent and maximal slick conditions, respectively initially developed by Barger (1970) [16]. However, here we have additionally shown how the relationship could change at intermediate π/πthresh values. At π/πthresh ≥ 0.1 and below a wind speed of about 6 ms−1 the capillary wave height is < 1 cm. This portion of the parameter space readily generates capillary wave heights below the wavelengths relevant to microwave radar. We can therefore hypothesize a smooth reflective interfacial barrier under such conditions, with relatively high radar return. Capillary wave heights approach the centimeter scale when surface slicks are absent, thus falling within the altimetric diffuse reflection window (microwave frequency range). To zeroth order, we propose that scattering would be expected under these conditions.

2.3. Computing Surface Pressure for the North Atlantic from Lower Trophic Level Dynamics

To demonstrate the utility of our reduced model approach in identifying Sigma-0 blooms in real world data, we apply it to satellite chl-a data for the North Atlantic region. Though marine organo-chemical measurement sets remain relatively rare, the North Atlantic is relatively data-rich, making this region a suitable area for our case study. The Northern Atlantic has been previously divided into sub-zones with distinct environmental and biogeochemical conditions (Figure 2, [19]), i.e., the Atlantic Arctic Province (ARCT), the North Atlantic Drift Province (NADR), the North Atlantic Subtropical Gyral Province East (NASE), the North Atlantic Subtropical Gyral Province West (NASW), the North Atlantic Tropical Gyral Province (NATR), and the Western Tropical Atlantic Province (WTRA). Applying our model to the observed chl-a concentration in each of these unique Longhurst sub-regions allows us to demonstrate how the surface pressure predictions will vary with environmental conditions (Table 2).
For each sub-zone in the North Atlantic (Figure 2), we extracted average summer (June/July/August) chl-a data from NASA Aqua Modis satellite observations for a fifteen-year period, (2002–2016).
Chlorophyll can be converted to phytoplanktonic carbon concentration (Pc), assuming a C: chl-a ratio of 50 mgC m−3:1 mgchl-a m−3 [20]. Standard macromolecular cytosol percentages (Pprot% = 0.6, Plip% = 0.2, [12]) can then be applied to determine the fraction of the living carbon that is a cellular amphiphilic protein or lipid mass. The ecodynamic parameter values for g, K, Z, and γ can be approximate for each North Atlantic province, and they closely follow those used in the early but highly effective Fasham differential equations as described in Sarmiento et al. (1993) [15,21]. Surfactant lifetimes (τi) are assumed to depend on microbial and photochemical degradation and, along with the uncertainty factors (ai), were taken from a global assessment of macromolecule surface activity [11].
Using model Equations (3)–(5), we computed surface pressure from our estimates of regional chl-a concentration, and by assuming a π threshold value of 3 mNm−1, we computed the surface pressure ratios. We show that the higher the regional chl-a concentration, the greater the surface coverage by protein. Correspondingly, the surface pressure value and surface pressure ratio are similarly high. At 0.34 mNm−1 the ARTC has the highest surface pressure value, whereas the NASE region has the lowest computed surface pressure value of 0.07 mNm−1 (Table 2).

2.4. Spatial Variability of Roughness Length and Capillary Wave Height

As shown in Section 2.3, the combined effect of surface tension and wind speed determines the capillary wave height. Using the surface pressure values in each region of the North Atlantic (Figure 3, Table 3) along with the average summer wind speed values, we can compute the regional variation in the roughness length and the capillary wave height. By applying Equation (2) and our plot of elemental height vs. wind speed (Figure 1), the roughness length values are calculated for the relevant wind speed values (Table 3). The capillary wave heights can then be calculated using the direct proportionality shown in Equation (1).
Due to their differences in average wind and chl-a concentration, the North Atlantic regions have different capacities for capillary wave height. Increased capillary wave height values are achieved with low surface pressure and high wind speed values, and the low capillary wave height values are more likely at low wind speed and high surface pressure values. Although the highest surface pressure ratio value (0.11) was calculated in the ARCT region, the regionally averaged wind speed value is low (2 ms−1), resulting in the lowest computed capillary wave height (0.45 cm). Conversely, even though a low surface pressure ratio (0.04) was in the WTRA region, the higher regionally averaged wind speed (4 ms−1) led to a capillary wave height value of 1.96 cm. The NATR region has the highest reported wind speed values (6 ms−1) with a low surface pressure ratio of 0.03, resulting the largest simulated capillary wave height of 2.91 cm.

2.5. Capillary Wave Height and Sigma-0 Frequency

The relationship between capillary wave height and the reflectivity of the surface is examined by plotting the fraction of a Longhurst region that has a Sigma-0 bloom occurrence at any point in the time series from 2001 to 2019 (Sigma-0 frequency, σf) against the calculated capillary wave height for each of the North Atlantic ecological regions considered. The annual average σf values (Table 4) were calculated by averaging estimates for σf from Jason 1, 2 and 3 satellite mission data for 2001–2019. We impose seasonality on the annual frequency with a scaling based on older (1992–1999) seasonal TOPEX/Poseidon data to synthesize the previously observed peaks and valleys in the data in summer and winter [3].

2.6. Sensitivity of Surfactant Concentration and Surface Pressure Threshold to Capillary Wave Height

The sensitivity of our reduced model predicting the connection of sea-air interface reflectivity to Sigma-0 blooms was tested by varying the surfactant concentrations and the surface pressure threshold. Negative relationships were observed between the capillary wave height and Sigma-0 frequency (Figure 4 and Figure 5) indicating that the lower the wave height the higher Sigma-0 frequency.
To provide a reference framework and assist the reader in visualizing the potential for real-world correlations at the basin scale, we have plotted annual and seasonal Sigma-0 frequency against the capillary wave heights calculated in the absence of a surface slick (Figure 4A). This reference panel shows the distribution of friction-inducing wave heights for a pristine, abiotic ocean. No surfactant monolayer or film has been allowed and the Sigma-0 frequency values were permitted to range freely. Sigma-0 frequency values for the lower latitude WTRA and NATR regions were considered constant and a-seasonal because they lie along, or near, the equator where the solar angle and the biological activity in the upper water column can be assumed to be relatively constant [19]. Our other regions in the North Atlantic boast strong seasonal variability in both solar input and biological activity; as such, we computed both annual and summer Sigma-0 values for these regions.

2.7. Impact of Surfactant Concentration on Capillary Wave Heights

Our model results show that, in the absence of surfactants, there was no relationship between Sigma-0 frequency and capillary wave height (Figure 4A). When surfactant concentrations were increased, the capillary wave heights tended to be reduced, most notably in the ARCT and NADR regions due to the high surfactant concentration resultant by the high biological activity in those regions. Assuming the surfactants consisted solely of proteins (protx3), the summer Sigma-0 frequency negatively correlated (r = −0.31, Figure 4B) with the capillary wave height. In contrast, the annual Sigma-0 frequencies were weakly positively correlated (r = 0.16, Figure 4B), with the capillary wave height. Adding amphiphilic lipids to the surfactants with the primary multiplier of ten (to bring the concentration closely into an agreement with observations, Ogunro et al., 2015 [11]), did not notably change the relationship between the capillary wave height and the Sigma-0 frequency (r = −0.31, Figure 4C), indicating that lipid (×10) has only a small impact on the capillary wave height. At extreme lipid concentrations (×100 there was a weakly negative correlation between wave height and the Sigma-0 frequency (r = −0.13, Figure 4D), while on an annual basis the relationship was positive (r = 0.36, Figure 4D). Under such high lipid concentrations, the reduction in wave height in the ARCT and NADR regions is large relative to the other areas. Returning to a more moderate lipid concentration but increased protein concentration (×10) increased the relative strength of the negative correlation (r = −0. 27, Figure 4E), indicating surfactant behavior of proteins is high in diminishing the wave height, leading to a smooth surface which is more visible in the NASE region. Finally, with moderate protein but extreme lipid, the relative strength of the negative correlation was reduced only slightly (r = −0.22, Figure 4F), with respect to Figure 4E. The regions ARCT and NADR shift to the maximum wave height reduction, as the surface pressure reaches to its saturation point with the increased lipid concentration (Figure 4F) as these regions have the highest biological activity (ARCT = 0.85 mgm−3, NADR = 0.58 mgm−3). Additionally, the relative wave height reduction of the WTRA regions is considerable as it’s reported the third highest biological activity (0.20 mgm−3) in combination with the wind speed (4 ms−1). The regions (ARCT, NADR and WTRA) with higher biological activity are more sensitive to wave height reduction.

2.8. Impact of the Surface Pressure Threshold on Capillary Wave Height

The surface pressure threshold varies in a range based on the surfactant adsorption at the sea-air interface. For the surfactant concentration combination protein(×3) + lipid(×10), an additional series of sensitivity tests was carried out adjusting the surface pressure threshold value utilized in Equation (2) from the preliminary estimate of 3 mNm−1, which is the maximum surface pressure value that exists in polar and coastal capillary waves seasonally [8,12], downward to 1 mNm−1 to represents film pressure of natural sea water [22], and upward to 10 mNm−1 to represent extreme surface pressure values [23]. The surface pressure threshold values could be changed based on a regions tendency to support slicks. The area and chemical composition of natural slicks are dependent on phytoplankton speciation, nutrient availability, and physical processes of the region. For the seasonal Sigma-0 frequency, negative correlations continue to exist with capillary wave height with both the reduced surface pressure threshold (r = −0.26, Figure 5A) and the increased surface pressure threshold (r = −0.43, Figure 5B). Conversely, the relationship between annual Sigma-0 frequencies and capillary wave height was weakly positively correlated with both low (r = 0.18, Figure 5A) and high (r = 0.10, Figure 5B) surface pressure thresholds. The damping of the capillary waves was more prominent at the 1 mNm−1 surface pressure threshold as capillary wave heights were reduced to 1.0 × 10−5 (ARCT), 1.6 × 10−4 (NADR), 1.3 × 10−1 (NASE), 1.5 (NATR), 4.5 × 10−2 (WTRA) cm (summer, Figure 5A). In contrast, the capillary wave heights are relatively high 5.4 × 10−1 (ARCT), 7.0 × 10−1 (NADR), 1.4 (NASE), 2.9 (NATR), 2.0 (WTRA) cm when surface pressure threshold was increased to 10 mNm−1 (summer, Figure 5B). A reduced pressure threshold is thus equivalent to the accumulation of several surfactants. Different region-wise changes to the capillary wave heights arose when implementing the different surface pressure thresholds of 1 mNm−1 and 10 mNm−1. With a reduced surface pressure, we find that the relative reduction in capillary wave height was ARCT > NADR > WTRA > NASE > NATR (Figure 5A) while with the increased surface pressure the change of the wave height was ARCT > NADR > NASE > WTRA > NATR (Figure 5B). In both instances, the regions where the highest chl-a concentrations were reported (ARCT, NADR, Table 2) were the regions where the capillary wave heights decreased more prominently. In the other regions (NASE, NATR and WTRA) both wind speed and surfactant concentration played a role in determining the resultant capillary wave height.

2.9. Regional Ship Track Comparisons

Observations of surface macromolecules are sparse, but one recent basin-scale field measurement dataset does focus specifically on surfactants within the open water microlayer. Using these observations, Sabbaghzadeh et al. (2017) [24] provided calibrated true ambient surface adsorption concentration (surface activity) for the entire Atlantic based on a well-known benchmark. We compare our surfactant concentrations determined from monthly averaged chl-a satellite observations extracted along the location of the ship track to their observed microlayer surfactant activity expressed as an equivalent adsorption effect of Triton-X-100(T-X-100)—a reference surfactant comprising a polyethylene oxide chain attached to an aromatic lipid. The surface adsorption concentration patterns along the two cruise tracks were markedly different (Figure 6), likely reflecting the different paths they took as they traversed the North Atlantic Ocean [24]. Between 10–30 °N the T-X-100 concentrations observed during the AMT25 cruise (~0.1–0.3 mgL−1) were generally at least half the concentration observed at the same latitude on the later AMT24 cruise (~0.4–1.0 mgL−1). The AMT24 cruise found increased concentrations between 5–10 °N that were not observed during AMT25. The concentration of T-X-100 on both cruises increased to a peak of ~1.7–1.8 mgL−1 from 35–40 °N. Protein and lipid concentrations for the North Atlantic computed from satellite chl-a observations (Equation (5), Figure 3) broadly reproduce these observed patterns, with lower values below 35 °N. The modeled protein and lipid concentrations began to increase around 40 °N and were highest at 50 °N, the northernmost limit of the cruise track. The model also shows a dip-pattern from the maximum in equatorial waters (0–5 °N) falling to a minimum in the gyral regime (10–30 °N) and then rising again in the subtropics (30–50 °N) that is broadly consistent with observations.

3. Summary and Discussion

Ocean altimetric data are known to carry contamination due to the Sigma-0 bloom phenomenon, which leads to erroneous results in calculations of the sea surface height. Sigma-0 blooms (high backscattering values) are known to occur when the ocean surface is smooth and absent capillary waves [14,25]; usually in the presence of surface slicks and at low wind conditions [25,26]. However, to date, a valid qualitative or quantitative explanation for the blooms that connect the extent of surface slicking and wind speed to the capillary wave height, and thereby to the reflectivity of the sea-air interface, has not been provided. Here, a quantitative and qualitative biogeochemical approach was used to explain the origin of the surface slicks and how surfactant concentration could work with wind speed in determining capillary wave heights. Using a reduced bio-geo-chemical model we demonstrate that in the North Atlantic, regions of high biological activity and low wind speed will have capillary waves damped to heights less than typical radar altimeter wavelengths, resulting in a smooth sea-air interface with high radar backscattering values.
Our results suggest that the fraction of highly reflective ocean surface tends to increase as capillary wave amplitude decreases. We show that this is because as the capillary wave heights decrease the fraction of wave heights below the incident wavelength increases, which means that the sea-air interface presents as a smoother surface, and thus causes an increase in backscattering of a radar return signal. The amount of directional return depends on the roughness of the surface [2]). Radar altimeters transmit in the microwave region, i.e., at the centimeter scale. Capillarity of the global ocean interface takes the form of millimeter-centimeter scale objects, so the reflection versus scattering of an incoming radar signal will be highly dependent on damping [3,7]. This inverse relationship between reflectivity and wave heights is more obvious in low wind summertime conditions. During the summer, the biological activity and the injection of organic molecules are relatively high, which leads to a larger reduction of capillary wave height and, thus a smoother surface which is further facilitated by the low wind conditions. We can argue that an increase of the Sigma-0 bloom frequency caused by the biological activity would mean an increase in the backscattering value, i.e., the reflected power.
The increase in Sigma-0 frequency with the reduction in capillary wave heights depends on both the concentration of surfactant macromolecules, and on wind speed. The lysis of phytoplankton cells due to zooplankton grazing brings biomacromolecules to the surface ocean [9,13]. A fraction of these biomacromolecules contains chemical structures that tend to stay at the sea-air interface since the macromolecules contain both hydrophilic and hydrophobic parts. Adsorption of these macromolecules (mostly the lipids and proteins) coats at the sea-air interface with organic carbon chains and thereby changes the intermolecular interactions at the sea-air interface —readily depleting the micro-fastening system of the sea. We have demonstrated, with our modeled results, that in regions with lower regionally averaged wind speed (i.e., 2 ms−1), (i.e., ARCT, NADR, and NASE), the surfactant concentration resulting from high biological activity will play a dominant role in damping the capillary wave height. Under low wind conditions the region (ARCT) that has the highest biological activity tends to have the lowest capillary wave height (0.45 cm). In contrast, the region with the lowest surfactant concentration (NASE) tends to have the highest capillary wave height (1.33 cm). Even though both NASE and NATR had relatively low biological activity, wave heights in the two regions were quite different. The capillary wave height increase in the WATR region was facilitated by the stronger regional wind speed (4 ms−1). The combined effect of these biological and physical variables is also visible in the NATR region which has the highest regional wind speed (6 ms−1). The relatively large capillary waves that one would expect were not observed. We argue that this is due to the surfactants from the moderate biological activity (indicated by the regional chl-a concentrations), damping the development of capillary waves.
Our study suggests that the impact of surfactants on capillary wave heights is sensitive to both the concentration and composition of surfactants. This sensitivity was demonstrated through the change to capillary wave height in response to changes in the surfactants. The relative heights of the capillary waves in each region shifted from ARCT < NADR < NASE < WTRA < NATR to ARCT < NADR < WTRA < NASE < NATR with the increase in the protein surfactant concentration by a factor of 10.
The adsorption isotherms for proteins and lipids are different due to their structural divergence. Under the same surface pressure condition and concentrations, lipids tend to adsorb more strongly onto the interface than the proteins. The adsorption process also depends on the concentrations of surface-active materials [27]. We assume that during the cell disruption the protein: lipid injection ratio would be 60%:20%. and protein had a longer lifetime (10 days) than the lipid (2 days). As such, modeled lipids concentrations were lower at the sea-air interface than proteins and the protein adsorption contributed more to the damping of capillary waves. If the lipid concentrations are same as proteins, the damping effect by lipids would be more prominent. The adsorption isotherm in the presence of both protein and lipid biomacromolecules is a combination of the individual protein and lipid isotherms, although the weight of the protein isotherm on the combined isotherms is relatively high.
Knowledge of biological distributions will allow the identification of regions where radar altimetry ocean reflection models should not assume a clean surface. Our research has provided evidence of the extent to which biological surface activity can impact the structure of the surface ocean. Consideration of this biological activity is critical when interpreting radar return signals from satellite altimeters used to map sea surface topography. By considering satellite derived chlorophyll concentration (our proxy for biological activity) and regional wind speed data our simple model provides estimates of expected capillary wave height and thereby predicts the potential for a reflective surface. Our model shows interconnected functionality between biological, chemical, and physical parameters for the capillary height damping calculations. Our simplified modeling approach could be used to identify the potential distribution of Sigma-0 bloom events. More extensive Sigma-0, chlorophyll, and wind speed data, ideally on a finer temporal and spatial scale, would be needed to further validate our predictions, and resolve the occurrence of Sigma-0 bloom and correct for them.
The relative impact of the surfactant concentration and surface pressure could be quantified by the relative reduction percentage [(E0,no surfactantE0,with surfactant)100/E0,no surfactant]. For regions with high biological activity (ARCT, NADR), the relative reduction percentage for all surfactant concentrations at 3 mN/m lies between 95.0–99.9%, indicating the damping is almost 100%. For low biological productivity regions NAST and NATR damping lie in the 59–95% and 21–54%, respectively. The larger ranges in these regions are due to variability in surfactant composition.
The drag coefficient is an important parameter in the transfer of wind momentum to the ocean [28], thus impacting the ocean circulation. As the magnitude of the drag partially depends on the height of the capillary waves, our research suggests that momentum flux in a region will be impacted by the presence of surface biology. The drag coefficient of the surface ocean has been previously speculated to have a dependency on the macromolecular carbon concentration [12]. Here, using our reduced model, we have reinforced this finding and demonstrated that the presence of surfactants at the sea-air interface could impact surface ocean roughness by as much as 21–99% in highly biologically active regions. To our knowledge, this link between biology, chemistry, and the physical domain is not currently represented in any Earth System Model and the drag coefficient is usually set to a constant value on the global scale. We argue that ignoring the impact of biological macromolecules on capillary wave height could result in understanding the drag coefficient and, thus, bias model estimates of ocean dynamics.
Several global marine biogeochemical models currently represent dissolved organic carbon, but only in a simple generic form [29,30]. We advocate an extension of this bio-chemical variable to explicitly represent the dominant macromolecules (proteins and lipids) such that a spatially explicitly biologically responsive drag coefficient can be computed. In our validation process, we observed that the surfactant adsorption data from two different cruise tracks did not closely resemble the calculated protein and lipid concentrations from our simple model. We speculate that this could be due to a miss-match in the timing of the satellite data which was a summer average (June–August) with the time of ship-board sampling (Sep–Oct). Furthermore, satellites are limited by optical depth so usually ‘see’ deeper (~5–10 m) than the very ocean surface. Nevertheless, our simple model was able to replicate the broad patterns and thus shows promise as a useful tool for altimetric data correction.

4. Concluding Remarks

In conclusion, radar reflectivity of the ocean surface depends on the damping of the capillary wave dimension, which we have shown is determined by both surfactant and wind speed. Because of the impact of surfactants on ocean roughness, and thus the drag coefficient, we believe biological surface activity should be considered when interpreting radar altimeter data.

Author Contributions

Conceptualization, S.E., D.V. and A.J.; methodology, S.E. and G.A.G.; formal analysis, S.E. and A.J.; writing—original draft preparation, G.A.G. and A.J.; writing—review and editing, S.E. and D.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Regional and Global Model Analysis (RGMA) program of the US Department of Energy’s Office of Science as a contribution to the HiLAT-RASM project. This work was also supported by the New Mexico Institute of Mining and Technology (NMT) and NOAA grants NA16NOS0120020 and NA19NES4320002, the latter was a grant to the Univ. of Maryland Cooperative Institute for Satellite Earth System Studies (CISESS).

Institutional Review Board Statement

This is not applicable.

Informed Consent Statement

This is not applicable.

Data Availability Statement

Data would be available by contacting the authors.

Acknowledgments

This research used resources provided by the Los Alamos National Laboratory Institutional Computing Program. Los Alamos National Laboratory is operated by Triad National Security, LLC, for the National Nuclear Security Administration of the U.S. Department of Energy (Contract No. 89233218CNA000001). The authors also would like to thank Gayan Rubasinghege and Jeff Altig at the Department of Chemistry of NMT.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brown, G.S. The average impulse response of a rough surface and its application. IEEE Trans. Antenna Propag. 1977, 25, 67–74. [Google Scholar] [CrossRef]
  2. Wright, J. Backscattering from capillary waves with application to sea clutter. IEEE Trans. Antennas Propag. 1966, 14, 749–754. [Google Scholar] [CrossRef]
  3. Mitchum, G.T.; Hancock, D.W., III; Hayne, G.S.; Vandemark, D.C. Blooms of σ0 in the TOPEX Radar Altimeter Data. J. Atmos. Ocean. Technol. 2004, 21, 1232–1245. [Google Scholar] [CrossRef]
  4. Quartly, G.D. Optimizing sigma σ0 Information From the Jason-2 Altimeter. IEEE Geosci. Remote Sens. Lett. 2009, 6, 398–402. [Google Scholar] [CrossRef]
  5. Thibaut, P.; Amarouche, L.; Zanife, O.Z.; Steunon, N.; Vincent, P.; Raizonville, P. Jason-1 altimeter ground processing look up correction tables. Mar. Geod. 2004, 27, 409–431. [Google Scholar] [CrossRef]
  6. Gómez-Enri, J.; Vignudelli, S.; Quartly, G.D.; Gommenginger, C.P.; Cipollini, P.; Challenor, P.G.; Benveniste, J. Modeling Envisat RA-2 waveforms in the coastal zone: Case study of calm water contamination. IEEE Geosci. Remote Sens. Lett. 2010, 7, 474–478. [Google Scholar] [CrossRef] [Green Version]
  7. Cox, C.S.; Zhang, X.; Duda, T.F. Suppressing breakers with polar oil films: Using an epic sea rescue to model wave energy budgets. Geophys. Res. Lett. 2017, 44, 1414–1421. [Google Scholar] [CrossRef]
  8. Elliott, S.; Burrows, S.; Cameron-Smith, P.; Hoffman, F.; Hunke, E.; Jeffery, N.; Liu, Y.; Maltrud, M.; Menzo, Z.; Ogunro, O. Does marine surface tension have global biogeography? Addition for the OCEANFILMS Package. Atmosphere 2018, 9, 216. [Google Scholar] [CrossRef] [Green Version]
  9. Kujawinski, E.B.; Farrington, J.W.; Moffett, J.W. Evidence for grazing-mediated production of dissolved surface-active material by marine protists. Mar. Chem. 2002, 77, 133–142. [Google Scholar] [CrossRef]
  10. Liss, P.; Watson, A.; Bock, E.; Jahne, B.; Asher, W.; Frew, N.; Hasse, L.; Korenowski, G.; Merlivat, L.; Phillips, L. Report Group 1-Physical processes in the microlayer and the air-sea exchange of trace gases. In The Sea Surface and Global Change; Cambridge University Press: Cambridge, UK, 1997; pp. 1–33. [Google Scholar]
  11. Ogunro, O.O.; Burrows, S.M.; Elliott, S.; Frossard, A.A.; Hoffman, F.; Letscher, R.T.; Moore, J.K.; Russell, L.M.; Wang, S.; Wingenter, O.W. Global distribution and surface activity of macromolecules in offline simulations of marine organic chemistry. Biogeochemistry 2015, 126, 25–56. [Google Scholar] [CrossRef]
  12. Elliott, S.; Menzo, Z.; Jayasinghe, A.; Allen, H.C.; Ogunro, O.; Gibson, G.; Hoffman, F.; Wingenter, O. Biogeochemical equation of state for the sea-air interface. Atmosphere 2019, 10, 230. [Google Scholar] [CrossRef] [Green Version]
  13. Ẑutić, V.; Ćosović, B.; Marčenko, E.; Bihari, N.; Kršinić, F. Surfactant production by marine phytoplankton. Mar. Chem. 1981, 10, 505–520. [Google Scholar] [CrossRef]
  14. Cheng, Y.; Tournadre, J.; Li, X.; Xu, Q.; Chapron, B. Impacts of oil spills on altimeter waveforms and radar backscatter cross section. J. Geophys. Res. Ocean. 2017, 122, 3621–3637. [Google Scholar] [CrossRef] [Green Version]
  15. Sarmiento, J.L.; Slater, R.D.; Fasham, M.J.R.; Ducklow, H.W.; Toggweiler, J.R.; Evans, G.T. A seasonal three-dimensional ecosystem model of nitrogen cycling in the North Atlantic euphotic zone. Glob. Biogeochem. Cycles 1993, 7, 417–450. [Google Scholar] [CrossRef]
  16. Barger, W.R.; Garrett, W.D.; Mollo-Christensen, E.L.; Ruggles, K.W. Effects of an artificial sea slick upon the atmosphere and the ocean. J. Appl. Meteorol. 1970, 9, 396–400. [Google Scholar] [CrossRef]
  17. Peixoto, J.P.; Oort, A.H. Observed mean state of the atmosphere. In Physics of Climate; Springer: New York, NY, USA, 1992; pp. 131–165. [Google Scholar]
  18. Garrett, W.D.; Bultman, J.D. Capillary-wave damping by insoluble organic monolayers. J. Colloid Sci. 1963, 18, 798–801. [Google Scholar] [CrossRef]
  19. Longhurst, A.R. Ecological Geography of the Sea; Academic Press: San Diego, CA, USA, 1998. [Google Scholar]
  20. Fasham, M.; Sarmiento, J.L.; Slater, R.; Ducklow, H.; Williams, R. Ecosystem behavior at Bermuda Station “S” and Ocean Weather Station “India”: A general circulation model and observational analysis. Glob. Biogeochem. Cycles 1993, 7, 379–415. [Google Scholar] [CrossRef]
  21. Fasham, M.J.; Ducklow, H.W.; McKelvie, S.M. A nitrogen-based model of plankton dynamics in the oceanic mixed layer. J. Mar. Res. 1990, 48, 591–639. [Google Scholar] [CrossRef]
  22. Modini, R.; Russell, L.; Deane, G.; Stokes, M. Effect of soluble surfactant on bubble persistence and bubble-produced aerosol particles. J. Geophys. Res. Atmos. 2013, 118, 1388–1400. [Google Scholar] [CrossRef]
  23. Barger, W.R.; Daniel, W.H.; Garrett, W.D. Surface chemical properties of banded sea slicks. Deep. Sea Res. Oceanogr. Abstr. 1974, 21, 83–89. [Google Scholar] [CrossRef]
  24. Sabbaghzadeh, B.; Upstill-Goddard, R.; Beale, R.; Pereira, R.; Nightingale, P.D. The Atlantic Ocean surface microlayer from 50 °N to 50 °S is ubiquitously enriched in surfactants at wind speeds up to 13 m s−1. Geophys. Res. Lett. 2017, 44, 2852–2858. [Google Scholar]
  25. Tournadre, J.; Chapron, B.; Reul, N.; Vandemark, D. A satellite altimeter model for ocean slick detection. J. Geophys. Res. Ocean. 2006, 111, 4. [Google Scholar] [CrossRef]
  26. Carswell, J.R.; Donnelly, W.J.; McIntosh, R.E.; Donelan, M.A.; Vandemark, D.C. Analysis of C and Ku band ocean backscatter measurements under low-wind conditions. J. Geophys. Res. Ocean. 1999, 104, 20687–20701. [Google Scholar] [CrossRef] [Green Version]
  27. Jarvis, N.; Garrett, W.; Scheiman, M.; Timmons, C. Surface chemical characterization of surface-active material in seawater. Limnol. Oceanogr. 1967, 12, 88–96. [Google Scholar] [CrossRef]
  28. Peixoto, J.P.; Oort, A.H. Exchange processed between the earth’s surface and the atmosphere. In Physics of Climate; Springer: New York, NY, USA, 1992; pp. 131–165. [Google Scholar]
  29. Manizza, M.; Menemenlis, D.; Zhang, H.; Miller, C.E. Modeling the recent changes in the Arctic Ocean CO2 sink (2006–2013). Glob. Biogeochem. Cycles 2019, 33, 420–438. [Google Scholar] [CrossRef]
  30. Seitzinger, S.; Harrison, J.; Dumont, E.; Beusen, A.H.; Bouwman, A. Sources and delivery of carbon, nitrogen, and phosphorus to the coastal zone: An overview of Global Nutrient Export from Watersheds (NEWS) models and their application. Glob. Biogeochem. Cycles 2005, 19, 4. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Elemental height (cm) as a function of 10 m wind speed and varying surface pressure ratios. Shaded area indicates parameter space where E0 < 1 cm.
Figure 1. Elemental height (cm) as a function of 10 m wind speed and varying surface pressure ratios. Shaded area indicates parameter space where E0 < 1 cm.
Atmosphere 13 02101 g001
Figure 2. Location of the North Atlantic sub-ecological regions as defined by Longhurst (1998) [19]. Region acronyms are defined explicitly in the text.
Figure 2. Location of the North Atlantic sub-ecological regions as defined by Longhurst (1998) [19]. Region acronyms are defined explicitly in the text.
Atmosphere 13 02101 g002
Figure 3. Summer (June/July/August) average chl-a concentration in the North Atlantic. The average was computed from Aqua MODIS data for 2002–2016. The outline of the Longhurst ecological regions identified in Figure 1 is shown.
Figure 3. Summer (June/July/August) average chl-a concentration in the North Atlantic. The average was computed from Aqua MODIS data for 2002–2016. The outline of the Longhurst ecological regions identified in Figure 1 is shown.
Atmosphere 13 02101 g003
Figure 4. Annual (•) and summer (Atmosphere 13 02101 i001) log Sigma-0 frequency vs. the log of capillary wave height when the protein (Prot) and lipid (Lip) concentrations are scaled by varying amounts. (A) no surfactants (B) baseline protein ×3 but no lipids, (C) baseline protein ×3 and baseline lipid ×10, (D) baseline protein ×3 and baseline lipid ×100, (E) baseline protein ×10 and baseline lipid ×10, (F) baseline protein ×10 and baseline lipid ×100. All computations were performed at the surface pressure threshold value of 3 mNm−1.
Figure 4. Annual (•) and summer (Atmosphere 13 02101 i001) log Sigma-0 frequency vs. the log of capillary wave height when the protein (Prot) and lipid (Lip) concentrations are scaled by varying amounts. (A) no surfactants (B) baseline protein ×3 but no lipids, (C) baseline protein ×3 and baseline lipid ×10, (D) baseline protein ×3 and baseline lipid ×100, (E) baseline protein ×10 and baseline lipid ×10, (F) baseline protein ×10 and baseline lipid ×100. All computations were performed at the surface pressure threshold value of 3 mNm−1.
Atmosphere 13 02101 g004
Figure 5. Log Sigma-0 frequency (σf) arrayed against the log of capillary wave height for (A) πthresh = 1 mNm−1 and (B) πthresh = 10 mNm−1. All computations were performed at the surfactant concentration of protx3+lipx10.
Figure 5. Log Sigma-0 frequency (σf) arrayed against the log of capillary wave height for (A) πthresh = 1 mNm−1 and (B) πthresh = 10 mNm−1. All computations were performed at the surfactant concentration of protx3+lipx10.
Atmosphere 13 02101 g005
Figure 6. Comparison of observed North Atlantic surface activity measured as T-X-100 concentration and estimated surface activity (surface adsorption concentration) as determined by modeled protein and lipid concentrations.
Figure 6. Comparison of observed North Atlantic surface activity measured as T-X-100 concentration and estimated surface activity (surface adsorption concentration) as determined by modeled protein and lipid concentrations.
Atmosphere 13 02101 g006
Table 1. Model parameter descriptions and values. The use of alternate values for aprot and alip is discussed in Section 2.7 and the choice of alternative values for πthresh is discussed in Section 2.8.
Table 1. Model parameter descriptions and values. The use of alternate values for aprot and alip is discussed in Section 2.7 and the choice of alternative values for πthresh is discussed in Section 2.8.
SymbolDescriptionValueUnits
PcPlanktonic carbon atom
concentration
Calculated by zoneμM Carbon
gZooplanktonic grazing rate1.0d−1
ZZooplankton concentration0.5μM Carbon
KHalf saturation for ingestion7.0μM Carbon
γ Assimilation efficiency0.5-
τprotLifetime of proteins10d
τlipLifetime of lipids2d
Pprot%Percentage of the protein content in a typical planktonic cell (inside the cytosol)60 so 0.6 as ratio-
Plip%Percentage of the lipid content in a typical planktonic cell (inside the cytosol)20 so 0.2 as ratio-
CiMicromolar carbon atom concentration of the ith macromoleculeCalculated by zoneμM Carbon
KprotHalf saturation micromolar carbon atom concentration for the protein30μM Carbon
KlipHalf saturation micromolar carbon atom concentration for the lipid2μM Carbon
θiEffective fractional interfacial coverage of the ith macromoleculeCalculated by zone-
aprotScale factor for protein concentration uncertainties3, 10-
alipScale factor for lipid concentration uncertainties10, 100-
πprot, maxMaximum 2D surface pressure of protein20mNm−1 (mJm−2)
πlip, maxMaximum 2D surface pressure of lipid50mNm−1 (mJm−2)
πthreshThreshold surface pressure defined regionally1, 3, 10mNm−1
π/πthreshSurface pressure ratio for logarithmic scalingCalculated by zone-
z0z0,surf = Roughness length with the slick
z0,pure = Roughness length without the slick
Calculated by zonecm
E0Capillary wave heightCalculated by zonecm
uWind speed for the Longhurst zoneClimatologyms−1
σfSigma-0 frequencyRemotely sensed (radar altimetry)-
Table 2. Computing surface pressure from chl-a in the North Atlantic ecological regions assuming a threshold surface pressure (πthresh) of 3 mNm−1.
Table 2. Computing surface pressure from chl-a in the North Atlantic ecological regions assuming a threshold surface pressure (πthresh) of 3 mNm−1.
Zone
Longhurst
Summer Average chl-a
(mgm−3)
Planktonic Carbon (Pc) (μM C)[Protein] (Cprot)
(μM C)
Protein Coverage Fraction (Θprot)[Lipid] (Clip)
(μM C)
Lipid Coverage Fraction (Θlip)π
(mNm−1)
π/πthreshold
ARCT0.853.540.520.0170.038.34 × 10−150.340.11
NADR0.582.420.400.0130.031.01 × 10−150.260.09
NASE0.130.520.110.0040.013.11 × 10−200.070.02
NATR0.130.560.120.0040.015.37 × 10−200.080.03
WTRA0.200.850.170.0060.011.08 × 10−180.110.04
Table 3. Roughness length and elemental height values calculated for each Longhurst region.
Table 3. Roughness length and elemental height values calculated for each Longhurst region.
Zone LonghurstWind Speed (ms−1)π/πthresholdz0,surf
(cm)
E0,surf
(cm)
ARCT20.110.040.45
NADR20.090.060.61
NASE20.020.131.33
NATR60.030.292.91
WTRA40.040.201.96
Table 4. Annual and summer (June/July/August) Sigma-0 frequency values on a base ten logarithmic scale. Seasonal values are computed from these annual values based on seasonal TOPEX results shown in Mitchum et al. (2004) [3].
Table 4. Annual and summer (June/July/August) Sigma-0 frequency values on a base ten logarithmic scale. Seasonal values are computed from these annual values based on seasonal TOPEX results shown in Mitchum et al. (2004) [3].
Zone LonghurstLog10f)
AnnualSeasonal
ARCT−1.2−0.72
NADR−1.4−0.92
NASE−0.75−0.27
NATR−1.5−1.5
WTRA−0.8−0.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jayasinghe, A.; Elliott, S.; Gibson, G.A.; Vandemark, D. The Role of Phytoplankton Biomacromolecules in Controlling Ocean Surface Roughness. Atmosphere 2022, 13, 2101. https://doi.org/10.3390/atmos13122101

AMA Style

Jayasinghe A, Elliott S, Gibson GA, Vandemark D. The Role of Phytoplankton Biomacromolecules in Controlling Ocean Surface Roughness. Atmosphere. 2022; 13(12):2101. https://doi.org/10.3390/atmos13122101

Chicago/Turabian Style

Jayasinghe, Amadini, Scott Elliott, Georgina A. Gibson, and Douglas Vandemark. 2022. "The Role of Phytoplankton Biomacromolecules in Controlling Ocean Surface Roughness" Atmosphere 13, no. 12: 2101. https://doi.org/10.3390/atmos13122101

APA Style

Jayasinghe, A., Elliott, S., Gibson, G. A., & Vandemark, D. (2022). The Role of Phytoplankton Biomacromolecules in Controlling Ocean Surface Roughness. Atmosphere, 13(12), 2101. https://doi.org/10.3390/atmos13122101

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop