Review and suggestions for estimating particulate organic carbon and dissolved organic carbon inventories in the ocean using remote sensing data

Dissolved organic carbon（DOC） and particulate organic carbon（POC） are basic variables for the ocean carbon cycle.Knowledge of the distribution and inventory of these variables is important for a better estimation and understanding of the global carbon cycle.Owing to its considerable advantages in spatial and temporal coverage,remote sensing data provide estimates of DOC and POC inventories,which are able to give a synthetic view for the distribution and transportation of carbon pools.To estimate organic carbon inventories using remote sensing involves integration of the surface concentration and vertical profile models,and the development of these models is critical to the accuracy of estimates.Hence,the distribution and control factors of DOC and POC in the ocean first are briefly summarized,and then studies of DOC and POC inventories and flux estimations are reviewed,most of which are based on field data and few of which consider the vertical distributions of POC or DOC.There is some research on the estimation of POC inventory by remote sensing,mainly in the open ocean,in which three kinds of vertical profile models have been proposed：the uniform,exponential decay,and Gauss models.However,research on remote-sensing estimation of the DOC inventory remains lacking.A synthetic review of approaches used to estimate the organic carbon inventories is offered and the future development of methods is discussed for such estimates using remote sensing data in coastal waters.     

Heat oscillation in the upper ocean of the southern South China Sea

Data used in this study are temperature/depth profiles taken over the upper 400 m of the ocean in the southern South China Sea (4°-14° N, 106°-120° E) for the period 1961-1973. The data are analyzed on the grid 2 (latitude) by 2 (longitude) in space and bimonthly in time. The vertically averaged temperature (TAV) over the upper 100 m of the ocean is calculated as the estimate of the heat content in the upper ocean.The TAV is cooler in the northwest region of the study area and warmer in the southeast in the annual and seasonal mean figures. The first EOF (Empirical Orthogonal Function) of anomalous TAV accounts for 41 % of the total variance for the period 1961-1973. The time function associated with it displays a significant interannual changes in the heat content, with 2-4 a oscillation period and associated with the ENSO events. During ENSO event TAV increases with the tendency of increasing towards equator along the basin. This anomalous states also exist in the water layers below 100 m depth. T     

Global parameter estimation of the Cochlodinium polykrikoides model using bioassay data

Cochlodinium polykrikoides is a notoriously harmful algal species that inflicts severe damage on the aquacultures of the coastal seas of Korea and Japan. Information on their expected movement tracks and boundaries of influence is very useful and important for the effective establishment of a reduction plan. In general, the information is supported by a red-tide(a.k.a algal bloom) model. The performance of the model is highly dependent on the accuracy of parameters, which are the coefficients of functions approximating the biological growth and loss patterns of the C. polykrikoides. These parameters have been estimated using the bioassay data composed of growth-limiting factor and net growth rate value pairs. In the case of the C. polykrikoides, the parameters are different from each other in accordance with the used data because the bioassay data are sufficient compared to the other algal species. The parameters estimated by one specific dataset can be viewed as locally-optimized because they are adjusted only by that dataset. In cases where the other one data set is used, the estimation error might be considerable. In this study, the parameters are estimated by all available data sets without the use of only one specific data set and thus can be considered globally optimized. The cost function for the optimization is defined as the integrated mean squared estimation error, i.e., the difference between the values of the experimental and estimated rates. Based on quantitative error analysis, the root-mean squared errors of the global parameters show smaller values, approximately 25%–50%, than the values of the local parameters. In addition, bias is removed completely in the case of the globally estimated parameters. The parameter sets can be used as the reference default values of a red-tide model because they are optimal and representative. However, additional tuning of the parameters using the in-situ monitoring data is highly required.As opposed to the bioassay data, it is necessary because the bioassay data have limitations in terms of the in-situ coastal conditions.     

The validation of HY-2 altimetermeasurements of a significant wave height based on buoy data

HY-2 has been launched by China on August 16, 2011 which assembles multi-microwave remote sensing payloads in a body and has the ability of monitoring ocean dynamic environments. The HY-2 satellite data need to be calibrated and validated before being put into use. Based on the in-situ buoys from the Nation- al Data Buoy Center （NDBC）, Ku-band significant wave heights （SWH, hs） of HY-2 altimeter are validated. Eleven months of HY-2 altimeter Level 2 products data are chose from October 1, 2011 to August 29, 2012. Using NDBC 60 buoys yield 902 collocations for HY-2 by adopting collocation criteria of 30 min for tempo- ral window and 50 km for a spatial window. An overall RMS difference of the SWH between HY-2 and buoy data is 0.297 m. A correlation coefficient between these is 0.964. An ordinary least squares （OLS） regression is performed with the buoy data as an independent variable and the altimeter data as a dependent vari- able. The regression equation of hs is hs （HY-2）=0.891 × hs （NDBC）＋0.022. In addition, 2016 collocations are matched with temporal window of 30 rain at the crossing points of HY-2 and Jason-2 orbits. RMS difference of Ku-band SWH between the two data sets is 0.452 m.     

Global air-sea surface carbon-dioxide transfer velocity and flux estimated using ERS-2 data and a new parametric formula

Using data from the European remote sensing scatterometer（ERS-2） from July 1997 to August 1998,global distributions of the air-sea CO₂ transfer velocity and flux are retrieved.A new model of the air-sea CO₂ transfer velocity with surface wind speed and wave steepness is proposed.The wave steepness（5） is retrieved using a neural network（NN） model from ERS-2 scatterometer data,while the wind speed is directly derived by the ERS-2 scatterometer.The new model agrees well with the formulations based on the wind speed and the variation in the wind speed dependent relationships presented in many previous studies can be explained by this proposed relation with variation in wave steepness effect.Seasonally global maps of gas transfer velocity and llux are shown on the basis of the new model and the seasonal variations of the transfer velocity and llux during the 1 a period.The global mean gas transfer velocity is 30 cm/h after area-weighting and Schmidt number correction and its accuracy remains calculation with in situ data.The highest transfer velocity occurs around 60°N and 60°S,while the lowest on the equator.The total air to sea CO₂ llux（calculated by carbon） in that year is 1.77 Pg.The strongest source of CO₂ is in the equatorial east Pacific Ocean, while the strongest sink is in the 68°N.Full exploration of the uncertainty of this estimate awaits further data.An effectual method is provided to calculate the effect of waves on the determination of air-sea CO₂ transfer velocity and fluxes with ERS-2 scatterometer data.     

A new merged dataset of global ocean chlorophyll a concentration with higher spatial and temporal coverage

Understanding the ocean's role in the global carbon cycle and its response to environmental change requires a high spatio-temporal resolution of observation.Merging ocean color data from multiple sources is an effective way to alleviate the limitation of individual ocean color sensors(e.g.,swath width and gaps,cloudy or rainy weather,and sun glint) and to improve the temporal and spatial coverage.Since the missions of Sea-Viewing Wide Field-of-View Sensor(Sea Wi FS) and Medium-spectral Resolution Imaging Spectrometer(MERIS) ended on December 11,2010 and May 9,2012,respectively,the number of available ocean color sensors has declined,reducing the benefits of the merged ocean color data with respect to the spatial and temporal coverage.In present work,Medium Resolution Spectral Imager(MERSI)/FY-3 of China is added in merged processing and a new dataset of global ocean chlorophyll a(Chl a) concentration(2000–2015) is generated from the remote sensing reflectance(Rrs(λ)) observations of MERIS,Moderate-resolution imaging spectra-radiometer(MODIS)-AQUA,Visible infrared Imaging Radiometer(VIIRS) and MERSI.These data resources are first merged into unified remote sensing reflectance data,and then Chl a concentration data are inversed using the combined Chl a algorithm of color index-based algorithm(CIA) and OC3.The merged data products show major improvements in spatial and temporal coverage from the addition of MERSI.The average daily coverage of merged products is approximately 24% of the global ocean and increases by approximately 9% when MERSI data are added in the merging process.Sampling frequency(temporal coverage) is greatly improved by combining MERSI data,with the median sampling frequency increasing from 15.6%(57 d/a) to 29.9%(109 d/a).The merged Chl a products herein were validated by in situ measurements and comparing them with the merged products using the same approach except for omitting MERSI and Glob Colour and MEa SUREs merged data.Correlation and relative error between the new merged Chl a products and in situ observation are stable relative to the results of the merged products without the addition of MERSI.Time series of the Chl a concentration anomalies are similar to the merged products without adding MERSI and single sensors.The new merged products agree within approximately 10% of the merged Chl a product from Glob Colour and MEa SUREs.     

How does the Indian Ocean subtropical dipole trigger the tropical Indian Ocean dipole via the Mascarene high？

The variation in the Indian Ocean is investigated using Hadley center sea surface temperature(SST)data during the period 1958–2010.All the first empirical orthogonal function(EOF)modes of the SST anomalies(SSTA)in different domains represent the basin-wide warming and are closely related to the Pacific El Ni o–Southern Oscillation(ENSO)phenomenon.Further examination suggests that the impact of ENSO on the tropical Indian Ocean is stronger than that on the southern Indian Ocean.The second EOF modes in different domains show different features.It shows a clear east-west SSTA dipole pattern in the tropical Indian Ocean(Indian Ocean dipole,IOD),and a southwest-northeast SSTA dipole in the southern Indian Ocean(Indian Ocean subtropical dipole,IOSD).It is further revealed that the IOSD is also the main structure of the second EOF mode on the whole basin-scale,in which the IOD pattern does not appear.A correlation analysis indicates that an IOSD event observed during the austral summer is highly correlated to the IOD event peaking about 9 months later.One of the possible physical mechanisms underlying this highly significant statistical relationship is proposed.The IOSD and the IOD can occur in sequence with the help of the Mascarene high.The SSTA in the southwestern Indian Ocean persists for several seasons after the mature phase of the IOSD event,likely due to the positive wind–evaporation–SST feedback mechanism.The Mascarene high will be weakened or intensified by this SSTA,which can affect the atmosphere in the tropical region by teleconnection.The pressure gradient between the Mascarene high and the monsoon trough in the tropical Indian Ocean increases(decreases).Hence,an anticyclone(cyclone)circulation appears over the Arabian Sea-India continent.The easterly or westerly anomalies appear in the equatorial Indian Ocean,inducing the onset stage of the IOD.This study shows that the SSTA associated with the IOSD can lead to the onset of IOD with the aid of atmosphere circulation and also explains why some IOD events in the tropical tend to be followed by IOSD in the southern Indian Ocean.     

世界海洋蒸发波导时空统计规律研究

The statistical features of the evaporation duct over the global ocean were comprehensively investigated with reanalysis data sets from the National Centers for Environmental Prediction.These data sets have time and spatial resolutions of 1 h and 0.313°×0.312°,respectively.The efficiency of the analysis was evaluated by processing weather buoy data from the Pacific Ocean and measuring propagation loss in the Yellow Sea of China.The distribution features of evaporation duct height(EDH) and the related meteorological factors for different seas were analyzed.The global EDH is generally high and demonstrates a latitudinal distribution for oceans at low latitudes.The average EDH is approximately 11 m over oceans beside the equator with a latitude of less than 20°.The reasons for the formation of the global EDH features were also analyzed for different sea areas.     

The influence of ENSO on sea surface temperature variations in the China seas

Positive SST anomalies usually appear in remote ocean such as the China seas during an ENSO event.By analyzing the monthly data of HadISST from 1950 to 2007,it shows that the interannual component of SST anomalies peak approximately 10 months after SST anomalies peak in the eastern equatorial Pacific.As the ENSO event progresses,the positive SST anomalies spread throughout the China seas and eastward along the Kuroshio extension.Atmospheric reanalysis data demonstrate that changes in the net surface heat flux entering into the China seas are responsible for the SST variability.During El Ni o,the western north Pacific anticyclone is generated,with anomalous southwester lies prevailing along the East Asian coast.This anticyclone reduces the mean surface wind speed which decreases the surface heat flux and then increases the SST.The delays between the developing of this anticyclone and the south Indian Ocean anticyclone with approximately 3–6 months cause the 2–3 months lag of the surface heat flux between the China seas and the Indian Ocean.The northwestern Pacific anticyclone is the key process bridging the warming in the eastern equatorial Pacific and that in the China seas.     

The Relationships among the Various Particulars of Japanese Fishing Boat

The leading aim of this paper is to explain with plain equations the relationships among the various particulars of Japanese fishing boats and then using these equations, determine a method for calculating the various particulars necessary in the basic design stage.Almost all the equations shown in tables 3 and 4 take tha form ( 1 ) Y=a + bx.In special cases.however, forms ( 2 ) Y=a + b In X and ( 3 ) Y=aebx are also assumed.Regression coefficients a and b arc obtained by regression analysis and r2xy for all the equations are larger than 0.9.The procedure for calculating the various particularsof fishing boats is shown by the arrow diagram in Fig.1.First of all, calculate the full load displacement △f and fish hold capacity FH from gross tonnage, then using Af and FH as initial conditions, determine two sets of various particulars. The average values got form the two sets of data are the various particulars of fishing boat. The comparison between the average calculated value and that which is currently     

Japanese Ocean Flux Data Sets with Use of Remote Sensing Observations (J-OFURO)

We have constructed ocean surface data sets using mainly satellite data and called them Japanese Ocean Flux data sets with Use of Remote sensing Observations (J-OFURO). The data sets include shortwave radiation, longwave radiation, latent heat flux, sensible heat flux, and momentum flux etc. This article introduces J-OFURO and compares it with other global flux data sets such as European Centre for Medium Range Weather Forecasting (ECMWF) and National Center for Environmental Prediction (NCEP) reanalysis data and da Silva et al. (1994). The usual ECMWF data are used for comparison of zonal wind. The comparison is carried out for a meridional profile along the dateline for January and July 1993. Although the overall spatial variation is common for all the products, there is a large difference between them in places. J-OFURO shortwave radiation in July shows larger meridional contrast than other data sets. On the other hand, J-OFURO underestimates longwave radiation flux at low- and mid-latitudes in the Southern Hemisphere. J-OFURO latent heat flux in January overestimates at 10°N–20°N and underestimates at 25°N–40°N. Finally, J-OFURO shows a larger oceanic net heat loss at 10°N–20°N and a smaller loss north of 20°N in January. The data of da Silva et al. in July show small net heat loss around 20°S and large gain around 20°N, while the NCEP reanalysis (NRA) data show the opposite. The da Silva et al. zonal wind speed overestimates at low-latitudes in January, while ECMWF wind data seem to underestimate the easterlies. This revised version was published online in August 2006 with corrections to the Cover Date.     

利用HOAPS资料研究南海海气界面热通量时空分布

基于第二版本HOAPS(Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite data)潜热、感热和海表温度(SST)3个参量的15 a(1988~2002年)逐月平均资料,利用经验正交方法分解分析了这3个参量在南海的时空分布.结果表明,在夏季模态,潜热表现为南高北低,感热表现为中间低两边高,两者主要都是海洋向大气输送热量,但大气有时也向南海中部输送感热;在冬季模态,潜热和感热的高值区都在南海北部,东北部有一强中心,该中心主要是由风场引起的;夏季SST的变化导致全年SST呈准半年周期变化.冬季SST的变化滞后于潜热变化1个月;除夏季和冬季模态外,冬夏转换季节模态也十分明显;HOAPS与NCEP(National Center of Environment Prediction)资料相比,两者3个参量的时空分布大体一致,区别在于HOAPS资料能更好地反映参量的一些细微特征.     

Validation of gridded data set of global surface wind/wind-stress vector field

Gridded data of global surface wind/wind-stress vectors, called J-OFURO v2, were obtained from satellite scatterometer (QuikSCAT/SeaWinds) data for the decade from August 1999 to July 2009. The data were validated by comparing with (1) in situ measurements from moored buoy observations, (2) a data product from the same scatterometer using a different gridding procedure (IFREMER), and (3) data products based on numerical models (NCEP-1 and NCEP-2). The results on averaging all buoy data revealed lower mean differences, lower root-mean-square (RMS) differences, and higher correlations for the wind/wind-stress components for the J-OFURO v2 product than the other data products. The statistical values calculated for each buoy station showed tendencies of decreasing reliability with increasing latitude in the mid-latitude region, while the reliabilities in the equatorial areas were low. We performed intercomparisons between the J-OFURO v2 data and the other data sets to clarify discrepancies among different wind products in open ocean regions with few moored buoys. We determined that the meridional wind components from the NCEP products exhibited poleward deviations compared to data from the J-OFURO v2 product. Relatively high mean differences, high RMS differences, and low correlations were found in the equatorial ocean for the NCEP products. Striped features were spatially correlated with buoy locations in the equatorial Pacific, which suggested that the reliability of the NCEP products was governed by buoy locations in this region.     

西北太平洋海气界面热通量时空分布特征研究

基于第三版本HOAPs (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data)海表面温度、潜热通量、感热通量、海表面空气比湿以及海表面风场5个参量的18 a(1988～2005年)逐月平均资料,利用经验正交函数和奇异值分解方法分析了异常潜热和感热通量场在西北太平洋的时空分布特征及造成这种分布的主要影响因素.EOF的分析结果表明,异常潜热通量场主要体现为第一第二两个模态的变化,第一模态显示整个海域呈同相变化且在时间上呈准年周期变化,第二模态则描述了分别位于10°N,25°N和40°N的3个极值中心并伴随多年振荡,由因子载荷分布可知热带太平洋是第二模态的行为中心,因此该模态可能与ENSO事件相关.异常感热通量场则主要表现为第一模态的变化,在时间上呈准年周期变化并伴随有多年时间尺度的振荡.奇异值分解方法的分析结果表明异常海表面风场是异常潜热和感热通量场时空变化的重要影响因素.     

1985-2005年东海海表温度时空变化特征分析

利用1985-2005年间AVHRR辐射计获取的东海海表温度资料,研究其时空变化特征.1985-2005年间东海海表平均温度总体呈上升趋势.对东海海温距平(SSTA)进行经验正交函数(EOF)分解,得到其2个主要变化模态,EOF1占总方差变化的34.8％,主要变化周期为3～5 a,表示1985-2005年中东海绝大部分...     

渤海冬季风生环流的年际变化特征及机制分析

为探究渤海冬季风生环流的年际变化及机制,在近35年NCEP CFSR大气强迫下,利用ROMS海洋模式对渤海冬季流场进行了高分辨率的数值模拟。模式结果基本重现了已知的渤海冬季风生环流的主要特征。渤海冬季深度平均流场的EOF结果显示:首先,第一模态中风向的改变引起了一个环绕整个渤海的环流结构变化;其次,第二模态反映了渤海冬季环流存在的线性增强趋势与风场增强有关,尤其在1995年之后。从模式结果中发现,风场对渤海冬季流场存在的年际变化起到重要作用,风向的偏转对环流年际特征的影响强于风速变化。     

全球百年海表面温度年际和年代际变化特征分析

利用100 a(1903—2002年)HadiSST的逐月资料,将全球海表面温度异常(SSTA)作为整体进行经验正交函数分解(EOF分解),提取了控制各大洋SSTA的主导模态和各大洋之间的联合模态,分析了相应的空间分布和时间序列。研究结果表明:SSTA变化最剧烈的海区是赤道中东太平洋、西北大西洋湾流海区和北太平洋黑潮延伸体海区。热带太平洋厄尔尼诺和南方涛动是主导模态并具有2~7 a周期的年际变化;SSTA变化第二主模态和第三主模态都是以大约70 a为周期的年代际变化为主的跨大洋联合模态。第二主模态的空间分布主要表现为中纬度北太平洋和北大西洋反位相、热带太平洋与大西洋反位相的太平洋-大西洋双偶极子型分布。SSTA变化的第三模态主要呈现南北半球海洋反位相的特征,为北太平洋-北大西洋-南大洋联合模态。第四模态基本上是反映各个不同海域特有的局地海洋-大气相互作用模态,该模态的时间序列具有1~4 a周期的年际变化和约9 a周期的年代际变化。     

热带太平洋-南海-印度洋海面风与海面温度年际变化整体耦合的主模态

利用COADS资料和相关分析、经验正交函数分解、奇异值分解和小波分析等方法对热带太平洋-南海-印度洋的海表面温度(SST)及海面风场(SSW)进行了研究,确定了热带太平洋、南海以及印度洋作为1个整体时海面风场及海面温度场相互作用的主模态。在把热带太平洋、南海和印度洋作为整体的情况下,研究了该模态的时空分布特征。各种方法研究均表明,该模态在热带太平洋上为ENSO模态;在热带印度洋上为对应于ENSO态的印度洋“单极型”模态;在南海则表现为整个海盆尺度上的一致型。     

赤道潜流变化及其与ENSO的关系

应用TOGA/TAO实测海流和SODA全球海洋再分析资料,分析研究了热带太平洋赤道潜流的变化特征及其与ENSO循环的关系,初步探讨了赤道潜流变化对两种不同类型(中部型和东部型)ENSO事件的影响。对赤道潜流距平场的EOF分解表明,第一模态为"东太平洋潜流模态",特征向量主要反映东太平洋赤道潜流的变化情况;第二模态为"中太平洋潜流模态",特征向量主要反映中太平洋赤道潜流的变化情况;这两个主要模态基本可以反映赤道潜流距平场的主要信息。这两个模态对应的时间系数有明显的年际变化特征,并且与NINO指数都有较好的负相关,最大相关均通过99.9%的信度检验。相关分析表明潜流第一模态的变化滞后于NINO指数变化,而第二模态的变化则提前于NINO指数变化,即ENSO事件对东太平洋赤道潜流的变化有一定影响,而中太平洋的潜流变化又会对ENSO事件产生影响。通过个例分析表明,中太平洋赤道潜流的变化对两种不同类型ENSO事件的发生有重要影响。在东部型ENSO事件发生前,中太平洋赤道潜流增强,这样西太平洋暖池的异常海温在潜流的引导下快速向东传播,直接到达东太平洋形成东部型ENSO事件;中部型ENSO事件发生前,中太平洋赤道潜流明显减弱,西太平洋的异常海温不能迅速东传速,在中太平洋堆积上升到达海面,使得中部型ENSO事件爆发。     

热带太平洋气候变率的三类模态及与ENSO强度变化的关系

基于Kaplan等重建的1856—2001年全球海表水温距平(SSTA)资料，用小波变换分析了热带太平洋SSTA的气候变率，对热带太平洋SSTA分别进行2—8、8—30和30—100a带通滤波，然后进行EOF分解。结果发现，ENSO模态具有5a左右的年际变化和15a左右的年代际变化2种显著周期，当二者位相相同时，ENSO事件加强，当二者位相相反时，ENSO事件减弱，当年际变化不明显时，显著的年代际变化也可单独导致ENSO事件；热带太平洋SSTA气候态变率以西太平洋暖池和赤道两侧的热带中东太平洋明显海温同号异常为主要特征，具60a左右的周期，其相位变化与气候跃变及El Nifio事件的类型有密切联系；长期增温倾向加大了El Nifio事件的振幅。文章最后指出，ENSO事件强度变化是由年际、年代际和气候态等3类模态变率共同作用的结果，在ENSO预报模式中考虑并引入年代际和气候态变化对ENSO循环的影响，是提高ENSO预测水平的有效途径之一。