Estimating photosynthetically available radiation at the ocean surface from ADEOS-II global imager data

A simple, yet efficient and fairly accurate algorithm is presented to estimate photosynthetically available radiation (PAR) at the ocean surface from Global Imager (GLI) data. The algorithm utilizes plane-parallel radiation-transfer theory and separates the effects of the clear atmosphere and clouds, i.e., the planetary atmosphere is modeled as a clear atmosphere positioned above a cloud layer. PAR is computed as the difference between the incident 400–700 nm solar flux at the top of the atmosphere (known) and the solar flux reflected back to space by the atmosphere and surface (derived from GLI radiance), taking atmospheric absorption into account. Knowledge of pixel composition is not required, eliminating the need for cloud screening and arbitrary assumptions about sub-pixel cloudiness. For each GLI pixel, clear or cloudy, a daily PAR estimate is obtained. Diurnal changes in cloudiness are taken into account statistically, using a regional diurnal albedo climatology based on 5 years of Earth Radiation Budget Satellite (ERBS) data. The algorithm results are verified against other satellite estimates of PAR, the National Centers for Environmental Prediction (NCEP) reanalysis product, and in-situ measurements from fixed buoys. Agreement is generally good between GLI and Sea-viewing Wide Field-of-view Sensor (SeaWiFS) estimates, with root-mean-squared (rms) differences of 7.9 (22%), 4.6 (13%), and 2.7 (8%) Einstein/m²/day on daily, weekly, and monthly time scales, and a bias of only 0.8–0.9 (about 2%) Einstein/m²/day. The rms differences between GLI and Visible and Infrared Spin Scan Radiometer (VISSR) estimates and between GLI and NCEP estimates are smaller and larger, respectively, on monthly time scales, i.e., 3.0 (7%) and 5.0 (14%) Einstein/m²/day, and biases are 1.1 (2%) and −0.2 (−1%) Einstein/m²/day. The comparison with buoy data also shows good agreement, with rms inaccuracies of 10.2 (23%), 6.3 (14%), and 4.5 (10%) Einstein/m²/day on daily, weekly, and monthly time scales, and slightly higher GLI values by about 1.0 (2%) Einstein/m²/day. The good statistical performance makes the algorithm suitable for large-scale studies of aquatic photosynthesis.     

Reconstruction of radiation fluxes at the ocean surface from NOAA satellite data

Algorithms for obtaining operative information on the fluxes of the total solar radiation and radiation balance at the ocean surface are developed and realized on the basis of satellite observations within the spectral ranges 0·725–1·1 and 10·3–11·3 m at arbitrary scanning angles. For the 10-day and monthly values, the reconstruction errors remain at the level of the accuracy of the shipboard measurements. The atlas of charts of the above parameters (grid 0·5°×0·5°) is compiled using data obtained by the RVAkademik Vernadsky in the Atlantic Ocean during 1986–1989.Translated by Mikhail M. Trufanov.     

Overview of geostationary ocean color imager (GOCI) and GOCI data processing system (GDPS)

GOCI, the world??s first geostationary ocean color satellite, provides images with a spatial resolution of 500 m at hourly intervals up to 8 times a day, allowing observations of short-term changes in the Northeast Asian region. The GOCI Data Processing System (GDPS), a specialized data processing software for GOCI, was developed for real-time generation of various products. This paper describes GOCI characteristics and GDPS workflow/products, so as to enable the efficient utilization of GOCI. To provide quality images and data, atmospheric correction and data analysis algorithms must be improved through continuous Cal/Val. GOCI-II will be developed by 2018 to facilitate in-depth studies on geostationary ocean color satellites.     

Estimating the sea surface dynamic topography from Geosat altimetry data

Optimal interpolation method is applied to Geosat altimetry data both to remove orbit error and to separate temporal mean sea surface dynamic topography (SSDT) from temporal fluctuations around the mean. Loss of long-wavelength oceanic signals at orbit error reduction procedure is smaller in this method than that in conventional collinear methods, but the areal average height over the study domain is still removed as the orbit error. The fluctuation SSDT is quantitatively evaluated by sea level data from tide gauge stations at Japanese islands. The correlation coefficient of the two sea-level variations is 0.83 when the loss of the areal average is compensated by the seasonal variation of the areal average height determined from the climatological monthly-mean SSDT. In addition, the improvement of the geoid model by combined use of Seasat altimetry data and hydrographic data is validated through the estimated temporal mean SSDT. In a local area where hydrographic data contemporary with the Seasat mission exist, the geoid model has been significantly improved so that the absolute SSDT can be determined from combination of the altimetry data and geoid model; the absolute SSDT describes the onset event of a quasi-stationary large meander of the Kuroshio south of Japan very well. Outside this local area, however, errors of several tens of centimeters still remain in the improved geoid model.     

Air temperature at ocean surface derived from surface-level humidity

A new method deriving surface air temperature from specific humidity is proposed. Surface atmospheric pressure and relative humidity in addition to specific humidity are necessary in order to derive surface air temperature. Assuming effects of variation of atmospheric pressure and relative humidity are small, climatological values are used for those values. Derived surface air temperature is compared with in situ surface air temperature. A cross-correlation coefficient is high and the rms error is small. However, the agreement between them varies spatially. The errors are largest in the eastern equatorial region and high-latitudes. The former may be caused by a large sampling error and remarkable internannual variation related to ENSO phenomena. On the other hand, the latter may be related to sensitivity of saturated vapor curve to air temperature. Sensible heat fluxes are estimated by using derived surface air temperature and compared with that by in situ data. For the whole North Pacific, a cross-correlation coefficient, a mean error and an rms difference are 0.89 W m^–2, 0.58 W m^–2 and 8.03 W m^–2, respectively.     

Surveillance of waste disposal activity at sea using satellite ocean color imagers: GOCI and MODIS

Korean Geostationary Ocean Color Imager (GOCI) and Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua observations of the variation in ocean color at the sea surface were utilized to monitor the impact of nutrient-rich sewage sludge disposal in the oligotrophic area of the Yellow Sea. MODIS revealed that algal blooms persisted in the spring annually at the dump site in the Yellow Sea since year 2000 to the present. A number of implications of using products of the satellite ocean color imagers were exploited here based on the measurements in the Yellow Sea. GOCI observes almost every hour during the daylight period, every day since June 2011. Therefore, GOCI provides a powerful tool to monitor waste disposal at sea in real time. Tracking of disposal activity from a large tanker was possible hour by hour from the GOCI timeseries images compared to MODIS. Smaller changes in the color of the ocean surface can be easily observed, as GOCI resolves images at smaller scales in space and time in comparison to polar orbiting satellites, e.g., MODIS. GOCI may be widely used to monitor various marine activities in the sea, including waste disposal activity from ships.     

Estimating the underwater light field from remote sensing of ocean color

We present a new approach that incorporates two models to estimate the underwater light field from remote sensing of ocean color. The first employs a series of analytical, semi-analytical, and empirical algorithms to retrieve the spectrum of inherent optical properties (IOPs), including the absorption and the backscatter coefficients, from the spectrum of remote sensing reflectance. The second model computes the profile of photosynthetically available radiation E _0,PAR (z) for a vertically homogeneous water column using the information of the retrieved IOPs and the ambient optical environment. This computation is based on an improved look-up table technology that possesses high accuracy, comparable with the full solution of the radiative transfer equation, and meets the computational requirement of remote sensing application. This new approach was validated by in situ measurements and an extensive model-to-model comparison with a wide range of IOPs. We successfully mapped the compensation depth by applying this new approach to process the SeaWiFS imagery. This research suggests that E _0,PAR (z) can be obtained routinely from ocean-color data and may have significant implications for the estimation of global heat and carbon budget.     

Estimating the potential of ocean wave power resources

The realistic assessment of an ocean wave energy resource that can be converted to an electrical power at various offshore sites depends upon many factors, and these include estimating the resource recognizing the random nature of the site-specific wave field, and optimizing the power conversion from particular wave energy conversion devices. In order to better account for the uncertainty in wave power resource estimates, conditional probability distribution functions of wave power in a given sea-state are derived. Theoretical expressions for the deep and shallow water limits are derived and the role of spectral width is studied. The theoretical model estimates were compared with the statistics obtained from the wave-by-wave analysis of JONSWAP based ocean wave time-series. It was shown that the narrow-band approximation is appropriate when the variability due to wave period is negligible. The application of the short-term models in evaluating the long-term wave power resource at a site was illustrated using wave data measured off the California coast. The final example illustrates the procedure for incorporating the local wave data and the sea-state model together with a wave energy device to obtain an estimate of the potential wave energy that could be converted into a usable energy resource.     

An analysis of the second-order Doppler return from the ocean surface

A theoretical analysis of the scattering of high-frequency (HF) electromagnetic waves from rough surfaces is proposed. The analysis, when applied to a model for the ocean surface, leads to new predictions of its second-order backscattered radar cross section in addition to those provided by existing theories. Some of the predictions of the theory have already been verified experimentally.     

Recent developments in ocean data telemetry at Woods HoleOceanographic Institution

The Woods Hole Oceanographic Institution is developing techniques for telemetering oceanographic data from the deep ocean to the laboratory in near real time. Three general approaches that provide a link between subsurface instruments and surface buoys equipped with satellite transmitters are being pursued. These approaches are: cabled systems that use electromechanical cables to connect subsurface instruments to a central controller; high data rate acoustic modems to transfer information between multiple remote units and a central controller; and inductive modems that use standard mechanical mooring lines as the transmission medium between instruments deployed on the mooring and a central controller. These telemetry systems are targeted for general use by the oceanographic community and are designed to be power efficient, low in cost, and capable of integration with most oceanographic data collection systems     

Coupling winds to ocean surface currents over the global ocean

A Wind stress–Current Coupled System (WCCS) consisting of the HYbrid Coordinate Ocean Model (HYCOM) and an improved wind stress algorithm based on Donelan et al. [Donelan, W.M., Drennan, Katsaros, K.B., 1997. The air–sea momentum flux in mixed wind sea and swell conditions. J. Phys. Oceanogr. 27, 2087–2099] is developed by using the Earth System Modeling Framework (ESMF). The WCCS is applied to the global ocean to study the interactions between the wind stress and the ocean surface currents. In this study, the ocean surface current velocity is taken into consideration in the wind stress calculation and air–sea heat flux calculation. The wind stress that contains the effect of ocean surface current velocity will be used to force the HYCOM. The results indicate that the ocean surface velocity exerts an important influence on the wind stress, which, in turn, significantly affects the global ocean surface currents, air–sea heat fluxes, and the thickness of ocean surface boundary layer. Comparison with the TOGA TAO buoy data, the sea surface temperature from the wind–current coupled simulation showed noticeable improvement over the stand-alone HYCOM simulation.     

Behavior of fresh water injected at the surface of a uniformly rotating ocean

The behavior of low density fresh water injected at the surface of a uniformly rotating saline water was investigated on the basis of a tank experiment. The injected water mass shows a clockwise circulation and grows gradually with an axisymmetric convex shape, until it breaks into two vortices at a critical size. An experimental formula for the change of radius of the water mass with time for the axisymmetric stage is obtained. It is shown that within our experimental range of values the radius of the water mass increases almost in proportion tot ^1/2, wheret is the elapse time, while the inviscid theory indicates that the radius should increase in proportion tot ^1/4. The dependence of the radius on elapse time is essential for forecasting the extent of discharged waters. The position of the maximum azimuthal velocity is fixed at \(V = - ge^{ - a^2 q^2 } \) within our experimental range of values wherer is the radial coordinate,f the Coriolis parameter,v the viscosity coefficient andQ the flow rate of injection, respectively. This radius corresponds to the radial scale derived by Gillet al. (1979). The steadiness of the position of the maximum azimuthal velocity may be essential in partition of the water mass into inner and outer regions and in the understanding the derived experimental formula. The critical radius for breaking is also investigated. The radius is shown to be independent ofQ and to be almost proportional to (Δ _ρ / _ρ )^1/2 f ^-1 whereρ is the density of the saline water andΔρ the density difference between the saline and injected waters. Even after the water supply is cut off in the axisymmetric stage, the radius of the water mass increases at almost the same rate as before, while its thickness decreases. The behavior after supply cut-off is discussed in the Appendix.     

Estimating effective shale area distributions from subsurface data

Estimation of shale bed continuity is needed in reservoir modeling and simulation as these constitute permeability barriers and baffles. However, efforts to characterize continuity of shale beds from 2D outcrop or well log sections lack critical 3D geometry and areal dimensions. A different approach is to estimate shale bed continuity by measuring the “effective shale area”. While it is a more indirect measure of shale bed continuity as it relies upon the sealing capacity of shales in selected, well-documented field examples, it eliminates the problem of the unknown but critical third dimension. Analyses of effective shale area distributions are made within structurally simple fields where hydrocarbon contacts were directly controlled by the sealing shales being measured, and faulting is limited. Datasets are exclusive to siliciclastic reservoirs formed in paleoenvironments ranging from non-marine to shallow water marine.Results show that effective shale areas are separated into two distinct distributions of fully marine versus non-marine/marginal marine origin. Shales deposited in non-marine to marginal marine paleoenvironments display smaller areas where effective sealing capacity is reduced by channel body incision. In fields where hydrocarbon columns are sealed by fully marine shale, larger areas are supported, often but not always, approaching structural spill points. In the case of reservoirs capped by these marine flooding shales, one can logically infer an abrupt decrease in channelization and resulting sizeable increase in continuity of shale bed seals.Similar shale types in traps formed in different structural regimes (e.g. extensional versus contractional domains) show comparable areal distributions, supporting use of data from different locations where the controls on hydrocarbon seal areas are well known. All shales measured here are areally large and would probably be accommodated in reservoir simulation deterministically, but use of the dimensions provided here would clearly be an aid in situations where a limited number of field appraisal wells are available. Discriminating between deterministic shales and stochastic shales, which are several orders of magnitude smaller in size, is also important for reservoir modeling and simulation. Results shown here should also be of great utility as input into discovery appraisal and field development scenario construction.     

Daily variability of dissolved inorganic radiocarbon at three sites in the surface ocean

We report radiocarbon measurements of dissolved inorganic carbon (DIC) in surface water samples collected daily during cruises to the central North Pacific, the Sargasso Sea and the Southern Ocean. The ranges of Δ¹⁴C measurements for each cruise (11–30‰) were larger than the total uncertainty (7.8‰, 2-sigma) of the measurements. The variability is attributed to changes in the upper water mass that took place at each site over a two to four week period. These results indicate that variability of surface Δ¹⁴C values is larger than the analytical precision, because of patchiness that exists in the DIC Δ¹⁴C signature of the surface ocean. This additional variability can affect estimates of geochemical parameters such as the air–sea CO₂ exchange rate using radiocarbon.     

Heat balance of the surface layer of the sea at ocean weather station T

Heat balance of the upper 200 m of the sea south of Japan is studied, by the use of marine meteorological and oceanographic data at Ocean Weather Station T (29°N, 135°E), intensively obtained from June 1950 to November 1953. Local time change of the heat content in the surface layer and the net heat flux through the air-sea interface are calculated directly from these data, and the heat convergence in the sea is estimated from their residuals. Regarding the relative importance of one- and three-dimensional processes, it is found that, on a time scale of a few days to one month, the variation of heat content depends on heat convergence in the sea, while on a seasonal time scale, the heat content is determined primarily by the heat flux through the sea surface in December through February, by heat convergence within the sea from March to May, and by both processes from June to November. It is inferred that the heat convergence in the sea is caused by advection of water masses which are bounded by sharp fronts. Spectral analysis of sea surface temperature indicates that they typically take 2 to 3 days to pass the station, and their typical size is estimated as around 20 km by assuming the typical advection velocity of water masses to be 10 cm s^?1.     

Microwave radiometric determination of wind speed at the surface of the ocean during BESEX

Microwave radiometric measurements were made at wavelengths ranging from 0.8 to 2.8 cm at altitudes from 0.16 to 11 km under well documented meteorological conditions over the Bering Sea. It is shown that determinations of wind speed at the ocean surface and liquid water content of the clouds may be made from such data. Determinations were made from two simultaneous but independent sets of radiometric measurements. The wind speeds and liquid water contents made from these two sets showed remarkable agreement. Independent estimates of these parameters made from in situ measurements showed reasonable agreement as well.     

Estimating climatic temperature change in the ocean with syntheticacoustic apertures

An acoustic tomography simulation is carried out in the eastern North Pacific ocean to assess whether climate trends are better detected and mapped with mobile or fixed receivers. In both cases, acoustic signals from two stationary sources are transmitted to ten receivers. Natural variability of the sound-speed field is simulated with the Naval Research Laboratory (NRL) layered-ocean model. A sequential Kalman-Bucy filter is used to estimate the sound speed field, where the a priori error covariance matrix of the parameters is estimated from the NRL model. A spatially homogeneous climate trend is added to the NRL fluctuations of sound speed, but the trend is not parameterized in the Kalman filter. Acoustic travel times are computed between the sources and receivers by combining sound speeds from the NRL model with those from the unparameterized climate trend. The effects of the unparameterized climate trend are projected onto parameters which eventually drift beyond acceptable limits. At that time, the unparameterized trend is detected. Mobile and fixed receivers detect the trend at about the same time. At detection time, however, maps from fixed receivers are less accurate because some of the unparameterized climate trend is projected onto tile spatially varying harmonics of the sound-speed field. With mobile receivers, the synthetic apertures suppress the projection onto these harmonics. Instead, the unparametrized trend is correctly projected onto the spatially homogeneous portion of the parameterized sound-speed field