TY - JOUR AU - Nayak,, Shailesh AB - Abstract Radiometric observations are carried out to estimate the net shortwave solar radiation in the upper ocean during the Arabian Sea Monsoon Experiment Phase II (ARMEX-II). Radiant heating rates in the Arabian Sea Warm Pool (ASWP) region from in situ radiometric measurements are presented for the first time. The estimated shortwave fluxes (300–700 nm) from the observed data are 101, 62, 40 and 29 W m-2 at 10, 20, 30 and 40 m depths respectively. About 5% of the surface light reaches below 50 m. Heating of the water column by penetrating solar radiation is substantial and estimated to be 0.18, 0.11, 0.08 and 0.07 °C day-1 in the upper 10, 20, 30 and 40 m respectively. Simultaneous observations from space-based sensors (SeaWiFS and IRS-P4 Ocean Colour Monitor) show an increase in solar absorption in regions with enhanced concentration of biologically active constituents. ARMEX, Arabian Sea Warm Pool, radiant heating rate, chlorophyll 1 Introduction Study of evolution of the Arabian Sea Warm Pool (ASWP) prior to the onset of the southwest (summer) monsoon is one of the key objectives of the Arabian Sea Monsoon Experiment (ARMEX). It is expected to eventually lead to a better understanding of monsoon variability. The role of shortwave radiant heating in conjunction with biologically active constituents in the thermodynamics of the ASWP has been studied for the first time during ARMEX Phase II (ARMEX-II). The chlorophyll-a (Chl-a) pigments in phytoplankton lead to strong optical attenuation (Smith and Baker 1978, Morel 1988). Therefore they can play a significant role in the upper ocean heat budget (Godfrey and Lindstrom 1989). An effort has been made to understand the relationship between the sub-surface temperature of the ASWP, radiant heating rate due to shortwave flux and Chl-a concentration over this domain. This paper demonstrates significant changes in the upper ocean heating rates and possible relationship with optically active bio-constituents using observations from cruise 190 of the research ship ORV Sagar Kanya (SK190) and Chl-a concentration from remote sensing satellites (SeaWiFS and OCM). 2 Data and methods Data from the SK190 cruise of ARMEX-II, conducted in the southeastern Arabian Sea (figure 1) from 14 March to 10 April 2003, are used to study the evolution of the pre-monsoon warm pool. In situ measurements of upper-ocean optical, physical, biological and meteorological parameters were made from ORV Sagar Kanya during this cruise. The optical profiles were used to estimate the radiant heating rates in terms of net penetrating radiation over the 300–700 nm wavelength domain. During SK190, optical measurements could be made only from 14 to 23 March 2003 because of a telemetry problem between the profiler and deck unit. However, the spatial coverage of sampling during this short period is good (figure 1, table 1) and the data can be treated as representative of the ASWP waters. Table 1 Station locations of concurrent bio-optical and physical measurement locations during SK190 (cf figure 1). St. no Position (°E/°N) Date St. no Position (°E/°N) Date 01 74.00/12.66 14 March 09 72.66/08.29 18 March 02 72.70/12.07 15 March 10 72.66/08.28 18 March 03 72.20/11.87 15 March 11 74.45/09.19 19 March 04 71.75/09.98 16 March 12 75.33/09.63 19 March 05 72.65/08.28 17 March 13 76.00/07.89 20 March 06 72.66/08.28 17 March 14 74.50/09.21 21 March 07 72.65/08.28 17 March 15 74.50/09.21 23 March 08 72.65/08.28 18 March St. no Position (°E/°N) Date St. no Position (°E/°N) Date 01 74.00/12.66 14 March 09 72.66/08.29 18 March 02 72.70/12.07 15 March 10 72.66/08.28 18 March 03 72.20/11.87 15 March 11 74.45/09.19 19 March 04 71.75/09.98 16 March 12 75.33/09.63 19 March 05 72.65/08.28 17 March 13 76.00/07.89 20 March 06 72.66/08.28 17 March 14 74.50/09.21 21 March 07 72.65/08.28 17 March 15 74.50/09.21 23 March 08 72.65/08.28 18 March Table 1 Station locations of concurrent bio-optical and physical measurement locations during SK190 (cf figure 1). St. no Position (°E/°N) Date St. no Position (°E/°N) Date 01 74.00/12.66 14 March 09 72.66/08.29 18 March 02 72.70/12.07 15 March 10 72.66/08.28 18 March 03 72.20/11.87 15 March 11 74.45/09.19 19 March 04 71.75/09.98 16 March 12 75.33/09.63 19 March 05 72.65/08.28 17 March 13 76.00/07.89 20 March 06 72.66/08.28 17 March 14 74.50/09.21 21 March 07 72.65/08.28 17 March 15 74.50/09.21 23 March 08 72.65/08.28 18 March St. no Position (°E/°N) Date St. no Position (°E/°N) Date 01 74.00/12.66 14 March 09 72.66/08.29 18 March 02 72.70/12.07 15 March 10 72.66/08.28 18 March 03 72.20/11.87 15 March 11 74.45/09.19 19 March 04 71.75/09.98 16 March 12 75.33/09.63 19 March 05 72.65/08.28 17 March 13 76.00/07.89 20 March 06 72.66/08.28 17 March 14 74.50/09.21 21 March 07 72.65/08.28 17 March 15 74.50/09.21 23 March 08 72.65/08.28 18 March Figure 1 Open in new tabDownload slide Locations (squares) for bio-optical measurements taken during ARMEX-II (14–23 March 2003). Figure 1 Open in new tabDownload slide Locations (squares) for bio-optical measurements taken during ARMEX-II (14–23 March 2003). 2.1 Optical data in the 400–700 nm range Optical radiometer/profilers, the SPMR (SeaWiFS Profiler Multi-channel Radiometer, or the profiler) and the SMSR (SeaWiFS Multi-channel Surface Reference, or the surface unit) deployed in the present experiment were designed and developed by Satlantic Inc., Canada. The profiler measures the upwelling radiance, Lu(z, λ), and the downwelling irradiance, Ed(z, λ), at seven discrete wavelengths (centred approximately at 412, 443, 490, 510, 555, 670 and 780 nm). The surface unit measures the downwelling irradiance, Es(0+, λ), just above the sea surface at these wavelengths. The nominal half-power bandwidth for each channel is 10 nm except for the last two (20 nm). Typical saturation of the radiance and irradiance signals are 20 µW cm-2 nm-1 sr-1 and 300 µW cm-2 nm-1, respectively. The detectors, made of Si-photodiodes, have spectrally corrected cosine response (for irradiance sensors) with noise equivalent irradiance (NEI) ∼5 × 10-5 W cm-2 nm-1 and noise equivalent radiance (NER) ∼1 × 10-6 W cm-2 nm-1 sr-1. The system data rate in free-fall mode is nearly 10 Hz with an optimized profiler speed of 0.4–0.7 m s-1 where it collects approximately 10–12 samples per metre depending on the free-fall speed of the profiler. Both systems are equipped with internal tilt and roll sensors for data quality checking during the casts. Optical data from individual profiles were excluded when the profiler tilt was greater than 5°. This reduces the uncertainty that could arise from the wave action and multiple reflections in the near surface waters. The radiometer units were deployed away from the vessel to avoid ship-induced perturbations and self-shadowing in the light field (Mueller and Austin 1993). 2.2 Shortwave flux correction for 300–400 nm The climatological mean (from the Earth Radiation Budget Experiment, ERBE) incident shortwave radiation at the surface during March is 250–280 W m-2. The solar radiation incident at the surface lies in the wavelength range 300–3000 nm. For a clear sky tropical region, the net shortwave flux (over 300–3000 nm) reaching the Earth’s surface is about twice the radiative flux in the 300–700 nm range (Baker and Frouin 1987). The solar radiation in the wavelength range 300–700 nm, however, is important in the upper ocean heat balance because it penetrates into the sub-surface layers. The total net shortwave flux at different layers beneath the surface must include the irradiances from 300 to 700 nm. Although the contribution of radiation in the wavelength range greater than 700 nm is negligibly small below the first few centimetres from the sea surface, the contribution from the 300 to 400 nm wavelength range is significant (figure 2). Since the Satlantic radiometer has only seven discrete bands in the ∼400–800 nm range, we performed a statistical analysis using hyper-spectral data (LI-COR, LI-1800) over the Arabian Sea (centred at 16°23′N/71°22′E) from a comparison in April 1997 to obtain the relationship between the 300–400 nm flux in relation to the total 300–700 nm flux. LI-COR is a portable spectro-radiometer, which measures hemispheric spectral irradiance (with cosine response) from 300 nm to 900 nm at 1 nm wavelength intervals. After examining each profile (eight profiles at depths 02, 05, 10, 20, 30 and 40 m) it was clear that the irradiance values beyond 700 nm are negligible below the sea surface due to a strong attenuation of longer wavelengths. The per cent contributions of irradiances over the 300–400 nm range (per cent fractional flux) are estimated as a function of depth based on linear regressions within the photic zone. Equations (1.1), (1.2) and (1.3) are used effectively to compute the per cent contribution of the downwelling irradiance profile, ⁠, the upwelling irradiance profile, ⁠, and the surface downwelling irradiance, ⁠, respectively: 1.1 1.2 1.3 So, the directional shortwave flux over 300–700 nm can be calculated as 2EN may be substituted by Es (i.e., z = 0+), Ed or Eu profiles from the Satlantic observations. The Satlantic measured upwelling radiances, Lu(z, λ), were multiplied by a Q-factor to obtain the upwelling irradiance profiles (Morel and Gentili 1991, 1993). All integrations were performed using the MATLAB® trapezoidal scheme. Prior to integration, all SPMR/SMSR irradiance spectra were interpolated to 1 nm wavelength using cubic splines (de Boor 1978). Figure 2 Open in new tabDownload slide Sample LI-COR (a) downwelling, (b) upwelling and (c) surface irradiance spectra taken during April 1997. The fractional irradiance over 300–400 nm (especially underwater) contributes significantly to the net shortwave (400–700 nm) irradiance. Panel (d) renders a typical Satlantic profile showing changes in the downward irradiance spectrum with depth, on 18 March 2003 around 1300 hrs LST. Correction for the shadowed region is applied to all Satlantic profiles in this study. Note the rapid attenuation of energy in wavelengths greater than 600 nm below 10 metre depth (panels (a) and (d)). Figure 2 Open in new tabDownload slide Sample LI-COR (a) downwelling, (b) upwelling and (c) surface irradiance spectra taken during April 1997. The fractional irradiance over 300–400 nm (especially underwater) contributes significantly to the net shortwave (400–700 nm) irradiance. Panel (d) renders a typical Satlantic profile showing changes in the downward irradiance spectrum with depth, on 18 March 2003 around 1300 hrs LST. Correction for the shadowed region is applied to all Satlantic profiles in this study. Note the rapid attenuation of energy in wavelengths greater than 600 nm below 10 metre depth (panels (a) and (d)). 2.3 Ancillary measurements The profiler also carried ancillary sensors—a WetSTAR fluorometer that measures the Chl-a fluorescence, a thermistor to measure the water temperature, and a Viatran pressure transducer to estimate the pressure converted depth. A significant underestimation of the Chl-a concentration measured by the WetSTAR fluorometer has been found (Raman 2003). Nevertheless, the relative accuracies of concentrations are adequate for the purpose of this study. To overcome the absolute accuracies, all Chl-a data were normalized to the cruise maximum. 3 Results Radiant heating of the upper ocean due to the absorption of sunlight by seawater and biogenous matter can be determined from the net irradiance spectrum (Ivanoff 1977, Morel 1988, Seigel et al 1995), which is defined as 3 where En, Ed and Eu, respectively, are the net irradiance, downwelling (measured) and upwelling (estimated) irradiances. The apparent changes in the net irradiance below the water surface is observed as rapid attenuation occurs in wavelengths greater than 600 nm. Below a depth of 10 m, the contribution of longer (red) wavelengths to the net irradiance is negligible (less than 5%). Thus the total net irradiance at greater depths can be attributed to the blue–green regions of the solar spectrum (with a peak at ∼490 nm). The surface flux and net solar radiation flux at each depth from the Satlantic data are calculated by integrating the net irradiance over the wavelength range 400–700 nm as 4 The net penetrating solar flux from the net irradiance spectrum over 300–700 nm, ⁠, was then corrected using equation (2). The transmission function, Tr(z), for penetrating solar radiation within the upper ocean, which gives the amount of energy available at any depth for a given surface energy, is the total net solar flux normalized by the total incident flux at the surface 5 Figure 3 represents the SK190 cruise mean transmission function, 〈Tr(z)〉, which follows an ideal near-exponential decay of the net flux with depth (Jerlov 1976). More than 50% of the surface flux is absorbed within the top 5 m. Thus, if the incident flux and the transmission function are known, the net solar flux profile can be retrieved from equation (5). Figure 3 Open in new tabDownload slide Cruise mean transmission function, Tr(z); the thick grey line shows the exponential fit to the mean transmission curve. Figure 3 Open in new tabDownload slide Cruise mean transmission function, Tr(z); the thick grey line shows the exponential fit to the mean transmission curve. Radiant heating rate for an upper layer ocean (from surface to a depth z) in the upper ocean is a function of the net solar flux just below the surface, ⁠, and the net flux at that depth, ⁠. Relating the rate of energy absorption to the rate of change of temperature, the radiant heating rate may be calculated from the first law of thermodynamics 6 where, ρw is typical seawater density (∼1025 kg m-3) and Cp is the specific heat of seawater at constant pressure (4100 J kg-1 °C-1). The net flux just beneath the surface layer can be approximately obtained from the net incident surface flux as ⁠, where the sea surface albedo (α) is 0.06 (Payne 1972). Our discrete observations during the entire cruise span almost equally from morning till evening (figure 4). Hence the cruise mean net penetrating flux counts for the diurnal cycle of the intensity of the solar energy. The mean day time incident flux at the ocean surface during SK190, 〈Es(0+)〉, was 308.5 W m-2. Table 2 summarizes the cruise mean transmission function at different depths and radiant heating rate of different layers in the upper ocean in the ASWP region. The unprecedented radiant heating rates for the ASWP domain, computed from the in situ observations, are comparable to those found by Seigel et al (1995) in the western Pacific warm pool. Table 2 Cruise mean transmission function, and radiant heating rates at different depths during the operation period. Depth (m) 〈Tr(z)〉 (%) 〈En(z)〉 (W m-2) Layer (m) 〈[∂Tz/∂t]rad〉 (°C day-1) 10 35.84 101.435 0–10 0.179 18 20 21.97 062.170 10–20 0.112 80 30 14.07 039.835 20–30 0.084 12 40 08.96 025.365 30–40 0.070 03 50 05.32 015.059 40–50 0.058 44 60 03.15 008.922 50–60 0.053 36 Depth (m) 〈Tr(z)〉 (%) 〈En(z)〉 (W m-2) Layer (m) 〈[∂Tz/∂t]rad〉 (°C day-1) 10 35.84 101.435 0–10 0.179 18 20 21.97 062.170 10–20 0.112 80 30 14.07 039.835 20–30 0.084 12 40 08.96 025.365 30–40 0.070 03 50 05.32 015.059 40–50 0.058 44 60 03.15 008.922 50–60 0.053 36 Table 2 Cruise mean transmission function, and radiant heating rates at different depths during the operation period. Depth (m) 〈Tr(z)〉 (%) 〈En(z)〉 (W m-2) Layer (m) 〈[∂Tz/∂t]rad〉 (°C day-1) 10 35.84 101.435 0–10 0.179 18 20 21.97 062.170 10–20 0.112 80 30 14.07 039.835 20–30 0.084 12 40 08.96 025.365 30–40 0.070 03 50 05.32 015.059 40–50 0.058 44 60 03.15 008.922 50–60 0.053 36 Depth (m) 〈Tr(z)〉 (%) 〈En(z)〉 (W m-2) Layer (m) 〈[∂Tz/∂t]rad〉 (°C day-1) 10 35.84 101.435 0–10 0.179 18 20 21.97 062.170 10–20 0.112 80 30 14.07 039.835 20–30 0.084 12 40 08.96 025.365 30–40 0.070 03 50 05.32 015.059 40–50 0.058 44 60 03.15 008.922 50–60 0.053 36 Figure 4 Open in new tabDownload slide Above surface net flux over the entire operation at different times of the day during the entire cruise. Figure 4 Open in new tabDownload slide Above surface net flux over the entire operation at different times of the day during the entire cruise. Time variation of the surface Chl-a and temperature is shown in figure 5. During the study period, a deep chlorophyll maximum (DCM) was found at around 60–80 m depth which is shallower than the summer monsoon (June–September) mixed layer depth (MLD) in the Arabian Sea (Prasad 2004). The deepening of the MLD after the warm pool event is due to the strong wind-driven vertical mixing. There was no significant change in the Chl-a concentration during this short observation except for stations 3 and 12–14, where the DCM was found at shallower depths. There was a significant fall in SST by 0.9 °C (figure 5) between stations 1 and 2 accompanied by cooling down to ∼100 m depth, associated with an increase in the Chl-a concentration. There was also a corresponding decrease in the MLD (estimated approximately from the temperature profile) from station 1 to station 2. The cooling and the increase in the Chl-a is probably because the ship moved into a region with an increased Chl-a concentration and a different thermal structure (see below). The OCM snapshot on 15 March reveals an eddy and traces of high chlorophyll concentration at its periphery (figure 6). This picture coincides with the in situ surface chlorophyll high on 15 March at station 3 as the Chl-a concentration in OCM images are around 0.5–0.7 µg l-1, which is higher than normal in open ocean (case 1) waters. Hence cooling and enhanced Chl-a concentration may be attributed to the eddy-induced vertical upwelling and injection of subsurface nutrients near station 3. An increase in the chlorophyll concentration is responsible for a reduction in the transmission of net solar radiation to depth (figure 5). The increase in the upper layer of the Chl-a concentration resulted in a decrease in the penetrative radiation at stations 2–4, which allows ∼10% of the surface energy down below the MLD (∼30 m) (table 2). Figure 7 shows the time evolution of the spatially averaged (70°–77°E/6°–14°N) satellite derived Chl-a concentration and light attenuation coefficients at 490 nm (K490) around the study area from November 2002 to July 2003. A clear increasing trend in the surface chlorophyll and light attenuation are seen from March until May 2003. Although in situ measurements were obtained for a short period, the results are significant, showing the possible coupling between the upper layer biological and physical processes and most importantly a baseline for the ASWP heating rates. Figure 5 Open in new tabDownload slide Surface and depth profiles of temperature (a), (b), chlorophyll-a (c), (d), and transmission function (e) for the entire cruise period. The chlorophyll values are normalized with the maximum value obtained during the entire cruise. Figure 5 Open in new tabDownload slide Surface and depth profiles of temperature (a), (b), chlorophyll-a (c), (d), and transmission function (e) for the entire cruise period. The chlorophyll values are normalized with the maximum value obtained during the entire cruise. Figure 6 Open in new tabDownload slide OCM derived Chl-a concentration over the southeastern Arabian Sea. The right panel shows the blow-up region (box in the left panel) and the point in the open circle marks the location of station 3 (72.20E/11.87N). Figure 6 Open in new tabDownload slide OCM derived Chl-a concentration over the southeastern Arabian Sea. The right panel shows the blow-up region (box in the left panel) and the point in the open circle marks the location of station 3 (72.20E/11.87N). Figure 7 Open in new tabDownload slide Time evolution of spatially averaged Chlorophyll-a concentration and light attenuation at 490 nm (K490) over the ASWP domain. The lines in the graph represent monthly values taken from the SeaWiFS standard product. The solid diamonds and squares correspond to the OCM derived Chl-a concentration and K490 for the cruise period. Figure 7 Open in new tabDownload slide Time evolution of spatially averaged Chlorophyll-a concentration and light attenuation at 490 nm (K490) over the ASWP domain. The lines in the graph represent monthly values taken from the SeaWiFS standard product. The solid diamonds and squares correspond to the OCM derived Chl-a concentration and K490 for the cruise period. 4 Discussion and conclusion This experiment was targeted to the understanding of the upper layer thermodynamics of the Arabian Sea Warm Pool and to quantify with the aid of novel optical instruments the radiant heating rates during the pre-onset of the southwest monsoon. About 14% of the total incident shortwave flux penetrates into the mixed layer (30 m), which eventually results in a heating rate of 0.08 °C day-1 in the 20–30 m water slab. The warming rate in the surface 10 m is about 0.18 °C day-1. In addition to this, satellite ocean colour data are quite useful for understanding the bio-physical coupling over large dynamic areas. Since remotely-sensed ocean colour information is not just a surface manifestation (the pigment concentration retrieved by means of remote ocean colour sensors is a value representative of the upper few metres of the water column i.e., the first attenuation depth), they can be treated as important components of upper ocean dynamics. This work also suggests that biogenic matter can have a significant influence on the absorption of sunlight (in 300–700 nm wavelength domains) as a function of depth in the upper ocean. The trapped solar energy eventually enhances radiant heating rates within the upper layers of the ASWP developing a shallow MLD before the onset of summer monsoon period (transition), when the impact of the wind-induced mixing is minimum. So the enhanced chlorophyll concentration could increase the heating rates until May although inter-annual variability cannot be ruled out. The present study demonstrates the importance of solar penetration and warm pool thermodynamics based on a small dataset collected during ARMEX-II. Nevertheless, further field observations of subsurface solar radiation together with optical and biological measurements are needed to understand the interrelationship between radiant heating and phytoplankton. Acknowledgments The authors thank the Department of Science and Technology, New Delhi for support of ARMEX, and the Department of Ocean Development, New Delhi for ship time on ORV Sagar Kanya. SeaWiFS data were provided by the SeaWiFS Project, NASA/Goddard Space Flight Center and ORBIMAGE. OCM data were obtained from the NRSA Data Centre, Department of Space (India). YP thanks Sam Lavender for helpful discussions. References Baker K S , Frouin R . , 1987 Relation between photosynthetically available radiation and total insolation at the ocean surface under clear skies , Limnol. Oceanogr. , vol. 32 (pg. 1370 - 7 ) Google Scholar Crossref Search ADS WorldCat de Boor C . , 1978 , A Practical Guide to Splines New York Springer Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Godfrey J S , Lindstrom E J . , 1989 The heat budget of the equatorial western Pacific surface mixed layer , J. Geophys. Res. , vol. 94 (pg. 8007 - 17 ) Google Scholar Crossref Search ADS WorldCat Ivanoff A . , 1977 Oceanic absorption of solar energy , Modelling and Prediction of the Upper Layer of the Ocean Krauss E B . New York Pergamon (pg. 47 - 71 ) Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Jerlov N G . , 1976 , Marine Optics New York Elsevier pg. pp 229 Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Morel A . , 1988 Optical modeling of the upper ocean in relation to its biogenous matter content (case I waters) , J. Geophys. Res. , vol. 93 (pg. 10749 - 68 ) Google Scholar Crossref Search ADS WorldCat Morel A , Gentili B . , 1991 Diffuse reflectance of oceanic waters: its dependence on sun angle as influenced by the molecular scattering contribution , Appl. Opt. , vol. 30 (pg. 4427 - 38 ) Google Scholar Crossref Search ADS PubMed WorldCat Morel A , Gentili B . , 1993 Diffuse reflectance of oceanic waters II. Bidirectional aspects , Appl. Opt. , vol. 32 (pg. 6864 - 79 ) Google Scholar Crossref Search ADS PubMed WorldCat Mueller J L , Austin R W . , 1993 Ocean optics protocols , NASA Tech. Memo 104566 , vol. 5 pg. 45 OpenURL Placeholder Text WorldCat Payne R E . , 1972 Albedo of the sea surface , J. Atmos. Sci. , vol. 29 (pg. 959 - 70 ) Google Scholar Crossref Search ADS WorldCat Prasad T G . , 2004 A comparison of mixed layer dynamics between the Arabian Sea and Bay of Bengal: one-dimensional model results , J. Geophys. Res. , vol. 109 pg. C03035 Google Scholar Crossref Search ADS WorldCat Raman M . , 2003 Personal communication Seigel D A , Ohlmann J C , Washburn L . , 1995 Solar radiation, phytoplankton pigments and the radiant heating of the equatorial Pacific warm pool , J. Geophys. Res. , vol. 100 (pg. 4885 - 91 ) Google Scholar Crossref Search ADS WorldCat Smith R C , Baker K S . , 1978 Optical classification of natural waters , Limnol. Oceanogr. , vol. 23 (pg. 260 - 7 ) Google Scholar Crossref Search ADS WorldCat © 2005 Nanjing Institute of Geophysical Prospecting TI - Radiant heating rates and surface biology during the Arabian Sea Monsoon Experiment JF - Journal of Geophysics and Engineering DO - 10.1088/1742-2132/2/1/003 DA - 2005-03-01 UR - https://www.deepdyve.com/lp/oxford-university-press/radiant-heating-rates-and-surface-biology-during-the-arabian-sea-7nDMWTETMi SP - 16 VL - 2 IS - 1 DP - DeepDyve ER -