2.1. Data sets used
The India Meteorological Department (IMD) best-tracked cyclone data (obtained from http://www.imd.gov.in) for the period 1993–2019 is utilized to identify the tracks and dates of the evolution of cyclones over the Bay of Bengal region. This data includes the maximum three-minute sustained wind speed (in knots) and minimum sea-level pressure (in mb). Classification of tropical disturbances based on wind speed is as given by IMD. The high-resolution (25 km or 0.25°) Optimum Interpolated Sea Surface Temperature (OISST, Reynolds et al. 2007) daily data is obtained from the website (http://www.esrl.noaa.gov/psd). A high resolution (0.083 x 0.083 degree, horizontal) Ocean reanalysis product from Marine Copernicus (available at https://resources.marine.copernicus.eu/, product ID: GLOBAL_MULTIYEAR_PHY_001_030) is used to compute different terms involved in mixed layer heat budget equation (Section 2.2). This product is available from January 1993 to December 2019. This Ocean reanalysis product assimilates real/near real-time ocean in-situ and satellite observations and is forced by high-resolution ERA-Interim and ERA5 atmospheric reanalysis, hence capable to represent TC associate/induced Ocean variability (eg. Avila-Alonso et al., 2021). More information and documentation of the data are available fromhttps://catalogue.marine.copernicus.eu/documents/PUM/CMEMS-GLO-PUM-001-030.pdf. These ocean data are also used to compute the ocean mixed layer depth (MLD) and depth of 26°C isotherm (D26). Ideally, the pre-monsoon season doesn’t show the formation of a barrier layer over the Bay of Bengal (Vissa et al. 2012, 2013); hence we have adopted temperature-based criteria for defining mixed layer. The MLD is defined as the extent of the ocean surface up to which the temperature decreases by 0.8°C from the surface temperature (Vissa et al. 2013; Thangaprakash et al. 2016). To study the influence of the Viyaru cyclone on basin-scale physical properties of BoB, we have examined the daily OISST; MLD and D26 obtained from Marine Copernicus, for the period 1993–2019. The TropFlux (Kumar et al. 2012, 2013; Girishkumar et al. 2013) data are used for computing different terms involved in Eq. 2 (section 2.2) of mixed layer heat budget, and the data are available from www.incois.gov.in/tropflux. All the data that are utilized in the present study is obtained for the period 1993 to 2019, and the same period is considered in computing climatological values.
2.2. Method - the Mixed Layer Heat Budget analysis
The Ocean mixed layer heat budget analysis is carried out to identify the dominant processes responsible for the surface ocean changes (up to the MLD) associated with the Viyaru cyclone. To understand the mechanisms and dominant processes associated with changes in the ocean mixed layer (i.e. during the active cyclone period, pre and post of the cyclone period) the mixed layer heat budget analysis is utilized. The MLD heat budget is computed using the method described in prior research (Thangaprakash et al. 2016; KBRR et al. 2020; Vijith et al. 2020). The temperature of the mixed layer is susceptible to surface heat fluxes, horizontal and vertical heat transport, and heat diffusion. However, heat diffusion via molecular processes is relatively minor in comparison to the other components, so it is ignored for ease of analysis. The mixed layer budget equation (Thangaprakash et al. 2016; KBRR et al. 2020; Vijith et al. 2020) is provided hereunder.
$$\frac{\partial T}{\partial t}=\frac{{Q}_{net}}{{\rho c}_{p}h}-\left[\left(u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}\right)+ H\left(Wh + \frac{dh}{dt}\right)\frac{(T-Th)}{h}\right] +\text{R}\text{e}\text{s}\text{i}\text{d}\text{u}\text{a}\text{l}$$
1
The temperature tendency term is given in the LHS of Eq. (1). The net surface heat flux, horizontal advection, heat flux due to vertical processes, and the residual are the separate components in the R.H.S of the equation. The net surface heat flux is indicated by Qnet.
Qnet = Qshortwave − (Qpen+ Qlongwave+ Qlatent+ Qsensible) (2)
Qshortwave= Net shortwave radiation,
Qpen= Shortwave radiation that penetrates below the ML,
Qlongwave= Net longwave radiation,
Qlatent= Latent heat flux,
Qsensible= Sensible heat flux.
where ρ (1,024 kg m–3) is the density of seawater, Cp (3,993 J kg–1K–1) is the specific heat capacity of seawater, h (m) is the mixed layer depth (MLD), u is the zonal current velocity (ms–1) and v is the meridional current velocity (m s–1), T (°C) is then mixed layer temperature, ∂T/∂x and ∂T/∂y is the horizontal gradient of temperature, t is time (in days), Wh is the vertical advection (m day–1). Th is the temperature just below the MLD (Girishkumar et al. 2013).
Section 2.1 describes the various data sets utilised in the mixed-layer heat budget study. Except for the climatological values, all of the variables in the mixed layer budget equation are given that are presented in this study are de-trended and removed for the seasonal cycle, to quantify the changes occurring due to the Viyaru cyclone.