Systematic Validation of Conductivity and Temperature from Ocean moored buoy data in the northern Indian Ocean with in situ ship based measurements

Similar documents
Inter-comparison of wave measurement by accelerometer and GPS wave buoy in shallow water off Cuddalore, east coast of India

MODELING INDIAN OCEAN CIRCULATION: BAY OF BENGAL FRESH PLUME AND ARABIAN SEA MINI WARM POOL

Country report - India NATIONAL INSTITUTE OF OCEAN TECHNOLOGY 07/10/2010

Atmospheric Forcing and the Structure and Evolution of the Upper Ocean in the Bay of Bengal

Effect of sea surface temperature on monsoon rainfall in a coastal region of India

Variability in the tropical oceans - Monitoring and prediction of El Niño and La Niña -

On some aspects of Indian Ocean Warm Pool

Long period waves in the coastal regions of north Indian Ocean

COMPARISON OF CONTEMPORANEOUS WAVE MEASUREMENTS WITH A SAAB WAVERADAR REX AND A DATAWELL DIRECTIONAL WAVERIDER BUOY

Variance-Preserving Power Spectral Analysis of Current, Wind and 20º Isothermal Depth of RAMA Project from the Equatorial Indian Ocean

PART I: DRAFT [PRACTICAL] GUIDELINES OF IOC, WITHIN THE CONTEXT OF UNCLOS, FOR THE COLLECTION OF OCEANOGRAPHIC DATA BY SPECIFIC MEANS

Climate Change Impacts on Sea Surface Temperature in the Eastern Mediterranean, Levantine Basin

ENSO and monsoon induced sea level changes and their impacts along the Indian coastline

Kathleen Dohan. Wind-Driven Surface Currents. Earth and Space Research, Seattle, WA

RECTIFICATION OF THE MADDEN-JULIAN OSCILLATION INTO THE ENSO CYCLE

Spatial Distribution of Sound Channel and Its Parameters in North Indian Ocean

Deep SOLO. Nathalie Zilberman, Dean Roemmich, and SIO float lab. 1. Deep SOLO float characteristics. 2. Deep SOLO float Deployment

Applications of Collected Data from Argos Drifter, NOAA Satellite Tracked Buoy in the East Sea

Validation of SST and Windspeed from TRMM using North Indian Ocean Moored Buoy Observations

LONG- TERM CHANGE IN PRE- MONSOON THERMAL INDEX OVER CENTRAL INDIAN REGION AND SOUTH WEST MONSOON VARIABILITY

Impact of Typhoons on the Western Pacific: Temporal and horizontal variability of SST cooling Annual Report, 2011 James F. Price

Temperature, Salinity, Dissolved Oxygen and Water Masses of Vietnamese Waters

EL NIÑO AND ITS IMPACT ON CORAL REEF ECOSYSTEM IN THE EASTERN INDIAN OCEAN

Goal: Develop quantitative understanding of ENSO genesis, evolution, and impacts

Data Analysis of the Seasonal Variation of the Java Upwelling System and Its Representation in CMIP5 Models

Subsurface Ocean Indices for Central-Pacific and Eastern-Pacific Types of ENSO

Interannual variation of northeast monsoon rainfall over southern peninsular India

TNA PROJECT PRELIMINARY REPORT 1 st Call of Proposals 12 January 3 April, 2012

THE northward-propagating intraseasonal oscillations

Seasonal variability of mixed layer depth (MLD) in the Bay of Bengal

Characteristics of Bay of Bengal Water mass in South Eastern Arabian Sea during

JCOMM Technical Workshop on Wave Measurements from Buoys

Indian Ocean Seasonal Cycle Jérôme Vialard (IRD) LOCEAN Paris France From Schott & McCreary (Prog. Oc.

Sea Surface Temperature of the Bay of Bengal Derived from the TRMM Microwave Imager*,+

N.V. Nanda Kumar, A. Nagarjuna and D.C. Reddy

Buoy data assimilation to improve wave height assessment in Bay of Bengal during monsoon seasons

Mechanistic links between the tropical Atlantic and the Indian monsoon in the absence of El Nino Southern Oscillation events

CHANGE OF THE BRIGHTNESS TEMPERATURE IN THE MICROWAVE REGION DUE TO THE RELATIVE WIND DIRECTION

The 2008 North Atlantic Bloom Experiment. Calibration Report #3

QC for Hydrographic Data GEOSECS to GO-SHIP. Susan M. Becker Scripps Institution of Oceanography Oceanographic Data Facility

Modification of the Stratification and Velocity Profile within the Straits and Seas of the Indonesian Archipelago

Temperature, salinity, density, and the oceanic pressure field

Oxygen measurements: sensors accuracy and scientific needs

An evaluation of XBT depth equations for the Indian Ocean

SEASONAL SEA LEVEL VARIABILITY AND ANOMALIES IN THE SINGAPORE STRAIT

Wave Energy Atlas in Vietnam

3 The monsoon currents in an OGCM

FIU Digital Commons. Florida International University. Henry O. Briceño Florida International University,

COMPARISON OF CONTEMPORANEOUS WAVE MEASUREMENTS WITH A SAAB WAVERADAR REX AND A DATAWELL DIRECTIONAL WAVERIDER BUOY

CORIOLIS, A FRENCH PROJECT FOR IN SITU OPERATIONAL OCEANOGRAPHY. S. Pouliquen, A. Billant, Y. Desaubies, G. Loaec, F. Gaillard, G.

INDIA METEOROLOGICAL DEPARTMENT (MINISTRY OF EARTH SCIENCES) SOUTHWEST MONSOON-2010 END OF SEASON REPORT

National Report of China

Argo Float Pressure Offset Adjustment Recommendations. 1: IORGC/JAMSTEC, 2-15, Natsushima-cho, Yokosuka, , Japan

South Bay Coastal Ocean Observing System California Clean Beaches Initiative

TRIAXYS Acoustic Doppler Current Profiler Comparison Study

Lecture 33. Indian Ocean Dipole: part 2

WEBB RESEARCH CORPORATION

7.4 Temperature, Salinity and Currents in Jamaica Bay

Lecture 20. Active-weak spells and breaks in the monsoon: Part 1

Ebenezer Nyadjro NRL/UNO. Collaborators: Dr. George Wiafe University of Ghana Dr. Subrahmanyam Bulusu University of South Carolina 1

Response of tropical cyclone potential intensity over the north Indian Ocean to global warming

Fine-Scale Survey of Right and Humpback Whale Prey Abundance and Distribution

Energy from seas and oceans

Acoustic communication for Maya Autonomous Underwater Vehicle - performance evaluation of acoustic modem.

.y..o ~ - \ o ~ ~~~I bl:..ill & ~j.a,_,.,ui J-1 ~4 b~

Seasonal Variation of Stratification in the Gulf of Thailand

Wintertime intraseasonal SST variability in the tropical South Indian Ocean and Role of Ocean Dynamics in the MJO Initiation

An IOOS Operational Wave Observation Plan Supported by NOAA IOOS Program & USACE

WHOTS Mooring Subsurface Instrumentation

Sea-Bird Scientific Carol Janzen, Ph.D. Oceanography David Murphy, M.S.E.E. Dan Quittman, M.S.E.E.

SEA SURFACE TEMPERATURE RETRIEVAL USING TRMM MICROWAVE IMAGER SATELLITE DATA IN THE SOUTH CHINA SEA

Lightning distribution with respect to the monsoon trough position during the Indian summer monsoon season

Stability of Temperature and Conductivity Sensors of Argo Profiling Floats

Lecture 13 El Niño/La Niña Ocean-Atmosphere Interaction. Idealized 3-Cell Model of Wind Patterns on a Rotating Earth. Previous Lecture!

Tidal regime along Vietnam coast under impacts of sea level rise

Investigation of Common Mode of Variability in Boreal Summer Intraseasonal Oscillation and Tropospheric Biennial Oscillation

Measured broadband reverberation characteristics in Deep Ocean. [E.Mail: ]

PERSPECTIVES FOR ACCURACY AND QUALITY ASSURANCE OF CTD & BIOGEOCHEMICAL DATA STREAMS FROM AUTONOMOUS PLATFORMS

CHAPTER 7 Ocean Circulation

Wave-Current Interaction in Coastal Inlets and River Mouths

ValidatingWindProfileEquationsduringTropicalStormDebbyin2012

MOTUS Wave Buoys. Powered By the Aanderaa MOTUS Directional Wave Sensor

COMPARISON OF DEEP-WATER ADCP AND NDBC BUOY MEASUREMENTS TO HINDCAST PARAMETERS. William R. Dally and Daniel A. Osiecki

Air-Sea Interaction Spar Buoy Systems

SCIENCE OF TSUNAMI HAZARDS

Vertical in situ profiles of nitrate and oxygen in the northern Japan Sea

Figure 4, Photo mosaic taken on February 14 about an hour before sunset near low tide.

Validation of Tsunameter in laboratory environment

MEASUREMENT COMMUNICATION INSIGHT

Effects of very harsh 2012 winter conditions on the deep layer of the Southern Adriatic Sea

Available online at ScienceDirect. Procedia Engineering 116 (2015 )

Subsurface Ocean Temperature Indices for Central-Pacific and Eastern-Pacific Types of El Niño and La Niña Events

The events associated with the Great Tsunami of 26 December 2004 Sea Level Variation and Impact on Coastal Region of India

Available online at ScienceDirect. Procedia Engineering 116 (2015 )

CHAPTER 6 DISCUSSION ON WAVE PREDICTION METHODS

Analysis of 2012 Indian Ocean Dipole Behavior

Thermohaline Front at the Mouth of Ise Bay

BUYER S GUIDE AQUAlogger 530WTD

Increasing trend of break-monsoon conditions over India - Role of ocean-atmosphere processes in the Indian Ocean

TUNA RESEARCH IN INDIA

Transcription:

Indian Journal of Geo-Marine Sciences Vol. 45(2), February 2016, pp. 224-229 Systematic Validation of Conductivity and Temperature from Ocean moored buoy data in the northern Indian Ocean with in situ ship based measurements Simi Mathew *, V. R. Shamji 1, G. Vengatesan, M. ArulMuthiah & R. Venkatesan Ocean Observation Systems, National Institute of Ocean Technology, Pallikaranai P.O., Chennai 600100, India 1 Faculty of Maritime Studies, King Abdulaziz University, P. O. Box 80401 - Jeddah 21589, Kingdom of Saudi Arabia. *[E-mail: simi@niot.res.in] Received 29 Octoberber 2014; revised 20 May 2015 A continuous data on the surface meteorological, surface and subsurface data on temperature, salinity and current is available online with the deployment of twelve moored buoys; Ocean Moored Network for the Northern Indian Ocean (OMNI) with seven buoys in the Bay of Bengal and five in the Arabian Sea. This network of OMNI buoys has been providing data which are of great relevance to the climate research community to constantly monitor the seasonal, intra-seasonal, annual and inter-annual variations in the northern Indian Ocean. In the long run the data accuracy is another important part of the program, to ensure the quality of the data delivered from these buoys, especially for the sub-surface data where satellite cannot reach. Conductivity-Temperature-Depth profile is done adjacent to the buoy locations during various phases of service of the buoys. This operation carried out in three phases in the Bay of Bengal has been compared with the buoy data, in order to ensure the quality of the conductivity-temperature measurements taken by the OMNI buoys. The high correlation between both the data sets ensures the quality of the data delivered by the Conductivity- Temperature sensors attached along the mooring line of the OMNI buoys. [Keywords: OMNI, CTD, BoB, SST] Introduction The Ministry of Earth Sciences, Government of India has initiated the National Data Buoy Program (NDBP) in 1997 under National Institute of Ocean Technology (NIOT), Chennai to make real time surface meteorological and near-surface oceanographic measurements. This has been later evolved as comprehensive ocean observation network program involving moored buoys which can provide continuous real time data on subsurface temperature, salinity and currents, named as the OMNI buoys. At present twelve OMNI buoys are being operational and they are deployed at seven specific locations in the Bay of Bengal (BoB) and five specific locations in the Arabian Sea (AS) as shown in figure 1. In addition to that, NIOT still maintains seven met-ocean buoys mostly for the coastal waters. In the northern BoB, two OMNI buoys are deployed adjacent to each other for the uninterrupted flow of data from this location, due its significance for various research interests. The lifetime of the buoy depends on the sensor performance, maintenance on timely basis and vandalism issues. The continuous inflow of data from these buoys has Fig. 1 The OMNI buoy locations in the Arabian Sea and Bay of Bengal made a leap in the studies of ocean currents, evolution of warm pool in the AS and BoB 1,2, intra-seasonal oscillations 3, barrier layer formation, propagation of planetary waves, fresh water influx 4 and study on the genesis and responses of the cyclones in the both the seas 5,6.

SIMI et al.: SYSTEMATIC VALIDATION OF CONDUCTIVITY AND TEMPERATURE FROM OMNI BUOY DATA 225 Materials and Methods OMNI buoy measures conductivity and temperature upto 500m depth with the help of Conductivity Temperature (CT) sensors attached to the mooring line at discrete depths at 5, 10, 15, 20, 30, 50, 75,100, 200 and 500m as shown in figure 2. Data is collected every one hour and transmitted through satellite to the land station once in every three hours. However as the dynamic ranges of these ocean parameters are very large and thereby sea truth validation assumes high priority. It is in this context the Conductivity Temperature Depth (CTD) surveys were carried out in a proximate location to the buoy to provide the sea truth data to validate the sub-surface CT sensor data. The CTD surveys were carried out under the deployment and maintenance of data buoy program by the NIOT. This was done in several phases, in this paper we are trying to elucidate the inter-comparison exercise conducted during three phases in the BoB; one during the deployment of the new four OMNI buoys during May-June, 2011in the BoB. The second phase was carried out during the maintenance of RAMA buoys in the BoB, during Fig. 2 Schematic diagram of the OMNI buoy. August-September, 2011. Third one was during April-May, 2012 during the service of the OMNI buoys after one year of performance in the sea. Similar kind of exercise was done in the AS too; one during the deployment of new buoys during October, 2012 and another one during the maintenance of the buoys during September- October, 2013. Ship based CTD samplings were carried out in close proximity to the buoy location, mostly within 5km of the selected OMNI buoy location. Ship-borne CTD, onboard SagarKanya and SagarNidhi was released mostly upto 3500m depth, depending on the prevailing weather conditions. Data were logged in with high resolution and processed on board ship lab using Sea Bird software. Processed data consists of depth, conductivity, temperature, salinity, dissolved oxygen, chlorophyll etc. The conductivity and temperature data from the OMNI buoy is validated through intercomparison with in situ CTD profiler data. Sensor types used in the both are of the same make (Sea Bird); SBE37-IM is used in the OMNI buoys and SBE-9 is the one used in the ship. Results and Discussions The Inter-comparison During May-June, 2011 a major cruise operation was conducted in the BoB for the deployment of the four new OMNI buoys. The CT sensors attached to the mooring of the OMNI buoys are tested and calibrated before deployment. CTD records data with very high vertical resolution of even less than 1m, depending on the speed with which it is lowered from surface to the desired depth. During the first phase which was during the onset of the Southwest Monsoon (SWM) season; the temperature and conductivity profiles at all the four buoy locations showed well mixed upper layer. The CTD sampling was carried out within one hour after deployment of BD13 at 11.02 N, 86.49 E on 5th June 2011. Ship CTD as well as the buoy recorded a well mixed upper layer with an isothermal layer depth (ILD) of around 50m. Both the data were in good agreement even up to 500m depth as shown in figure 3. CTD measurements taken at the other OMNI buoy locations too showed very good correlation of more than 0.99. OMNI buoys BD10 and BD11 located in the north-central and western BoB, respectively; showed double thermocline structure, which were reflected in both the datasets. BD08 buoy which is located 370 km

226 INDIAN J. MAR. SCI., VOL. 45, NO. 2 FEBRUARY 2016 south of the Ganges-Brahmaputra river-mouth recorded shallow ILD of around 30 m depth, with a sub-surface maximum in conductivity, which was in agreement with the CTD measurements too. Non-linearity in the temperature and salinity measurements can be expected in regions of strong river discharges especially for normal discharge years 7. Good correlation in the data ensures the quality of the OMNI buoy data which provides valid information on the impact of the strong river discharges especially in the northern BoB. temperature and salinity as the SWM season progress and the same was reflected in the buoy data too. But there was a mismatch in the conductivity values along the halocline. Since in the halocline region where the salinity varies rapidly with depth, even small vertical displacement of the mooring, due to the mooring scope of S type mooring will result in large variation in the conductivity values. Best option in such cases will be to go for daily average values wherein the pressure recorded at 500m depth by the Conductivity Temperature Depth (CTD) sensor is within a desired limit. Strong freshwater influx during the SWM season has reduced the ILD with respect to phase 1 measurements. CTD survey at south location, near BD 13 buoy showed an increase in the ILD compared to the other buoys. A sub-surface maximum in conductivity were also recorded by both the measurements. Inter-comparison showed good correlation even after three months of the buoy s performance at this specific location. A sharp gradient in the thermocline and halocline at the BD11 location during phase 2 of the cruise replaced the double thermocline and halocline present during the phase 1 as shown in figure 4. Both the CTD and CT sensors were in good agreement with each other with an over-all correlation of 0.99. Fig. 3 Inter-comparison between the ship CTD and the OMNI buoy data during the first phase of the cruise. The second phase of CTD survey was carried out during August-September, 2011; the peak of SWM season and on completion of 3 months of successful functioning of buoys at sea. Measurements obtained near to all the four buoy locations were compared, except for the BD10 location whose CT sensors malfunctioned. A ship based CTD sampling conducted near the BD08 location recorded the decrease in surface Fig. 4 Inter-comparison between the ship CTD and the OMNI buoy data during the second phase of the cruise.

SIMI et al.: SYSTEMATIC VALIDATION OF CONDUCTIVITY AND TEMPERATURE FROM OMNI BUOY DATA 227 Table 1 The drag in conductivity and temperature as obtained from calibration sheets of CT sensors retrieved from BD12 OMNI buoy. Depth of CT(m) Serial Number Drift since last calibration ( C/Year) - as received Third phase of inter-comparison exercise was carried out in the BoB during April-May, 2012; almost after one year of successful operation of the buoys in the BoB. During this cruise, CTD operation was carried out near BD08, BD10 and BD12 buoy locations. The good correlation between CTD and buoy data proved the good quality of data delivered from the buoys even after one year of operation in the sea. High correlation at BD08 location, confirms the mismatch in the conductivity data found at 20m and 30m depth during the second phase of intercomparison exercise was due to the vertical displacement of the mooring along with the sharp halocline structure present at that time. The CT sensors of the BD10 OMNI buoy were working fine after maintenance operation during August-September, 2011. So an intercomparison was carried out Drift since last calibration ( C/Year) -final Drift since last Drift since last calibration calibration (psu/year) (psu/year) - final - as received 5 SBE37SIP 60648-8051 -0.00048-0.00027 0.004 NA 10 SBE37IM 59729-7980 -0.00011-0.00012 0.0057 0.0062 15 SBE37IM 59729-7981 -0.00013-0.00005 0.0028 0.0035 20 SBE37IM 59729-7982 -0.00028-0.00028 0.0041 0.0053 30 SBE37IM 59729-7983 -0.00044-0.00048 0.0063 0.0081 50 SBE37IM 59729-7984 -0.00035-0.00031 0.0057 0.007 75 SBE37IM 59729-7985 -0.00027-0.00017 0.0026 0.0037 100 SBE37IM 59729-7986 -0.0002-0.00026 0.0022 0.0028 200 SBE37IM 59729-7987 -0.00024-0.0002-500 SBE37IM 59729-7991 -0.00016 - -0.0002 - Fig. 5 Inter-comparison between the ship CTD and the OMNI buoy data in the Bay of Bengal during the service of the buoys during April-May, 2012.

228 INDIAN J. MAR. SCI., VOL. 45, NO. 2 FEBRUARY 2016 during this cruise and the correlation was pretty high as shown in figure 5, even the fine features in the upper 75m column was very well captured. At BD12 location the sensors were replaced with new sensors during this cruise, since the sensors have completed two years of life after the first calibration. Retrieved sensors from BD12 were sent back to the manufacturers for re-calibration. Drift values noted by the manufacturer for the set of sensors used in BD12 are given in Table 1. The drift values are only in the third decimal place for conductivity measurements while in the fourth and fifth decimal place for temperature measurements. The choice whether to in-corporate those drift values in the measurements solely depends on the accuracy level with which the user aimed at, depending on the project. Inter-comparison near BD12 location were carried out with two sets of data, one which has completed its two years life and also with the new set of CT sensors attached to the BD12 buoy during this cruise as represented by the last two panels in the figure 5. CT sensor at 5m was not working during preservice and it was replaced with new one during the service. The only CT sensor which was showing deviation of more than 2 ms/cm in conductivity value was the one at 75m depth. It has nothing to do with the drift values since the drift value noted for the particular sensor at 75m depth with serial numbersbe-37 IM 7985 is very low as listed in Table 1. Difference in conductivity value must have arised due to the time difference between the CTD observations and buoy data prior to the service. The buoy data before service corresponds to 22 GMT on 29thApril and the CTD data is for 13GMT on 30th April. The time series data from 75m from the BD12 buoy shows an increase in the conductivity value for a short period and it is reflected in the intercomparison plot too. On the other side the time difference between CTD data and post-service sub-surface CT sensor data from the BD12 location are only 1 hour and both the data correlated very well with correlation coefficient of 0.99. Even the sub-surface increase in the conductivity values were very well captured inboth the measurements prior to and after service of the BD12 buoy as shown in the last two panels of figure 5. Discussion Seven OMNI buoys are placed in the BoB in order to closely monitor the temporal and spatial variations in the temperature and salinity structure of the BoB. These parameters are very crucial for studying the driving forces behind the seasonal, intra-seasonal and inter-annual variability in the BoB. The only source for subsurface data is either through Argo, moored buoys or the data collected during major experiments. Validation operation of the subsurface CT sensors of the OMNI buoys has been carried out regularly with an interval of 6 to 12 months time in order to ensure the good quality of the data. Drag values in the CT sensors are also obtained from the manufacturers while the sensors are re-calibrated after one year of operation in the sea. The CT sensor seems to be perfectly match with the ship based CTD observations even after one full year of operation in the sea which gives good confidence over the sub-surface CT sensor data delivered by the OMNI buoys. Acknowledgement We thank Ministry of Earth Sciences (MoES) for extending funding support to this program.we express our sincere thanks to Dr. M. A. Atmanand, Director, N.I.O.T. for his valuable guidance. We sincerely acknowledge the team efforts of NIOT technical and administrative staff and FugroOceanor, Norway who provided the OMNI buoys. Special thanks are due to Dr. DamodarShenoy, N.I.O., Goa for his valuable suggestions. Thanks are due to MoES, INCOIS and members of expert committee for the support of this project and for extending the ship time. References 1 Sengupta, D., P. K. Ray and G. S. Bhat, Spring warming of the eastern Arabian Sea and Bay of Bengal from buoy data. Geophys. Res. Lett., 29 (15) (2002), 1734, doi: 10.1029/2002GL015340. 2 J. Vimala, R. Venkatesan, G. Latha and R. R. Rao, Observed buildup and collapse of warm pool in eastern Arabian Sea and Bay of Bengal from moored buoy SST records from 1998-2008, The International Journal of Ocean and Climate Systems, 13-22, Vol. 5(2014), doi:10.1260/1759-3131.5.1.13. 3 Sengupta, D. and Ravichandran, M., Oscillations of Bay of Bengal Sea surface temperature during the 1998 summer monsoon. Geophys. Res. Lett., 28 (10) (2001), 2033-2036.

SIMI et al.: SYSTEMATIC VALIDATION OF CONDUCTIVITY AND TEMPERATURE FROM OMNI BUOY DATA 229 4 Chaitanya A.V.S., F. Durand, Simi Mathew, V. V. Gopalakrishna, F. Papa, M. Lengaigne, J. Vialard, Ch. Krantikumar and R. Venkatesan. Observed year-toyear sea surface salinity variability in the Bay of Bengal during 2009-2014 period, Ocean Dynamics, (2014), DOI 10.1007/s10236-014-0802-x. 5 Jossia Joseph, K., Balchand, A. N., Hareeshkumar, P. V., and Rajish, G., Inertial oscillation forced by the September 1997 cyclone in the Bay of Bengal. Current Sci., 92(2007), 790 794. 6 R. Venkatesan, Simi Mathew, J. Vimala, G. Latha, M. Arul Muthiah, S. Ramasundaram, R. Sundar, Lavanya R. And M. A. Atmanand, Signatures of very severe cyclonic storm Phailin in met-ocean parameters observed by moored buoy network in the Bay of Bengal. Current Sci., Vol. 107(2014), 589-595. 7 Fabien Durand, Fabrice Papa, Atiqur Rahman and Sujit Kumar Bala, 2011: Impact of Ganges- Brahmaputra interannual discharge variations on Bay of Bengal salinity and temperature during 1992-1999 period, J. Earth Syst. Sci., 120, No. 5 (2011), pp. 859-872.