Tracking Air–Sea Exchange and Upper-Ocean Variability in the Indonesian–Australian Basin during the Onset of the 2018/19 Australian Summer Monsoon

: Sea surface temperatures (SSTs) north of Australia in the Indonesian–Australian Basin are significantly influenced by Madden–Julian oscillation (MJO), an eastward-moving atmospheric disturbance that traverses the globe in the tropics. The region also has large-amplitude diurnal SST variations, which may influence the air–sea heat and moisture fluxes, that provide feedback to the MJO evolution. During the 2018/19 austral summer, a field campaign aiming to better understand the influences of air–sea coupling on the MJO was conducted north of Australia in the Indonesian–Australian Basin. Surface meteorology from buoy observations and upper-ocean data from autonomous fast-profiling float observations were collected. Two MJO convective phases propagated eastward across the region in mid-December 2018 and late January 2019 and the second MJO was in conjunction with a tropical cyclone development. Observations showed that SST in the region was rather sensitive to the MJO forcing. Air–sea heat fluxes warmed the SST throughout the 2018/19 austral summer, punctuated by the MJO activities, with a 2°–3°C drop in SST during the two MJO events. Substantial diurnal SST variations during the suppressed phases of the MJOs were observed, and the near-surface thermal stratifications provided positive feedback for the peak diurnal SST amplitude, which may be a mechanism to influence the MJO evolution. Compared to traditionally vessel-based observation programs, we have relied on fast-profiling floats as the main vehicle in measuring the upper-ocean variability from diurnal to the MJO time scales, which may pave the way for using cost-effective technology in similar process studies.


T
he Madden-Julian oscillation (MJO) is an eastward-moving disturbance of clouds, rainfall, winds, and pressure that traverses the globe in the tropics with a period in the range of 30-80 days.The MJO is the dominant mode of intraseasonal variability in the tropical Indo-Pacific warm pool (Zhang 2005).It modulates the timing and strength of monsoons and is an important influence for the development of some high-impact weather events such as tropical cyclones, flood, and drought in the Indo-Pacific region.The MJO offers promising perspectives to forecast tropical rainfall with a lead time of a couple of weeks; however, current state-of-art weather forecasting models have systematic biases and cannot reproduce the MJO well (e.g., Kim et al. 2014), possibly because they do not represent air-sea interaction processes associated with the MJO (DeMott et al. 2015).
During boreal summer, intraseasonal oscillations display prominent northwardnortheastward propagation and variability, extending much farther from the equator compared to the MJOs in austral summer and interacting with summer Asian monsoons (Wang and Xie 1997;Lee et al. 2013).After boreal summer, the first eastward-propagating MJO is often associated with the onset of the Australian monsoon (Hashiguchi et al. 1995;Duan et al. 2019), which normally takes place in early December (Kajikawa et al. 2010).The MJOs can also modulate the Australian monsoon rainfall and winds on intraseasonal time scales (Wheeler and Hendon 2004).
The strongest large-scale intraseasonal sea surface temperature (SST) variations in austral summer in the tropics are found north of Australia in the Indonesian-Australian Basin, forced by air-sea heat flux and wind-driven current anomalies associated with the MJO (Vialard et al. 2013;Marshall and Hendon 2014).In addition, large diurnal variation of SST (DV SST) has been modeled and observed in the region (Seo et al. 2014;Wang and Zhang 2017).The air-sea coupling at these two time scales can have significant influences on the MJO and its forecasts (DeMott et al. 2015), which has motivated a number of observation programs in the Indian Ocean: the Mirai Indian Ocean Cruise for the Study of the MJO (MISMO; Yoneyama et al. 2008), the air-sea interactions at several time scales in the Indian Ocean (Cirene; Vialard et al. 2009), Cooperative Indian Ocean Experiment on Intraseasonal Variability/Dynamics of the Madden-Julian Oscillation (CINDY-DYNAMO; Yoneyama et al. 2013;Moum et al. 2014), and Ocean Mixing and Monsoon/Air-Sea Interactions in the Northern Indian Ocean Regional Initiative (OMM-ASIRI; Wijesekera et al. 2016), mostly focused on the equatorial central Indian Ocean.Both observation and model studies show that the increase of daily mean SST can enhance the latent and sensible heat fluxes into the atmosphere, providing important feedbacks to the evolution of the MJO (Seo et al. 2014;Moum et al. 2014).Regional oceanography and changes of diurnal thermocline can affect the magnitude of diurnal SST variations nonlinearly, thus, exploring the diurnal variation of vertical mixing in the upper ocean is also critical to improve the MJO forecasting.
The Maritime Continent, a unique mixture of land and ocean straddling the equator between the Indian and Pacific Oceans, has been identified as a barrier/weakening mechanism for the eastward propagation of the MJOs by inhibiting convection development over the sea (Zhang and Ling 2017).There has been suggestion that warm SST anomalies in the Indonesian-Australian Basin, to the south of the Maritime Continent, may draw the MJO southward before reaching the Maritime Continent through air-sea interactions (Zhou and Murtugudde 2020).However, air-sea exchanges in the region during the MJO evolution has not been examined previously from field observations.
In November 2018, a field campaign was carried out in the region, under a bilateral collaboration between the Centre for Southern Hemisphere Oceans Research (CSHOR; CSHOR 2018), CSIRO, Australia, and the First Institute of Oceanography (FIO), Ministry of Natural Resources, China.This field campaign is also under the umbrella of the Years of the Maritime Continent (YMC) program during 2017-20 (www.jamstec.go.jp/ymc/).The goal of the YMC is to improve the understanding and prediction of local multiscale variability of the Maritime Continent weather-climate system and its global impact through observation and modeling.
In our campaign, one surface meteorology and ocean profiling buoy from the FIO, and eight profiling floats from CSHOR were deployed to the north of Australia, in a region that experiences relatively large MJO and diurnal SST variations (e.g., Duvel and Vialard 2007).The region also has large uncertainties in air-sea heat fluxes among different atmosphere reanalysis products (Yu et al. 2007;Hermes et al. 2019).For almost a decade, the site has been proposed as a high-priority extension of the tropical mooring array of the Indian Ocean [Research Moored Array for African-Asian-Australian Monsoon Analysis and Prediction (RAMA); McPhaden et al. 2009].The region will become a future focus of the Indian Ocean Observing System to improve the MJO and monsoon forecasting systems and a RAMA flux reference site has been reintroduced [Indian Ocean Observing System (IndOOS); Beal et al. 2019].
The aim of the field campaign is to obtain new insights into air-sea coupling in the Indonesian-Australian Basin using fast ocean profiling platforms to resolve the very nearsurface ocean temperature structures, as well as the key atmospheric variables such as humidity, air temperature, vector winds, and radiation parameters; such measurements in the region are unprecedented, to assess the upper-ocean response to the MJO forcing in the Indonesian-Australian Basin as compared to other regions in the Indian Ocean.The data will be used to understand the coupled model performance in capturing the scale, strength, and atmospheric signature of the MJO and diurnal warming events and to improve model parameterizations for the simulation of upper-ocean processes on these time scales.The data can also be used to better understand oceanography processes during the transition between the Australian monsoon and the trade wind regimes.This paper reviews the field campaign and highlights some key research results to be delivered under the CSHOR-FIO collaboration.

Design of the campaign
Buoy and profiling float measurements.A miniature version of the Bailong buoy system from the FIO (Cole et al. 2011) was used in the program (SB1), which has a surface meteorological package to measure the surface air temperature, humidity, pressure, winds, and shortwave and longwave radiation at a 10-min interval, and a subsurface (within the upper 500 m) oceanographic package to measure upper-ocean temperature, salinity, pressure, oxygen, and current speed and direction at a 30-min interval.The details of the instrumentation are listed in Tables ES1 and ES2.We chartered an offshore support vessel, Limitless, from the Offshore Unlimited Pty. Ltd. for the period of 19-24 November 2018 to deploy the buoy and the profiling platforms.After a regional bathymetry survey on 20 November, the buoy system was deployed on 21 November at 16°51.432΄S, 115°13.314΄E,on a plateau about 300 nautical miles (~555 km) north of Exmouth, Western Australia, just outside of the border of the Australian Exclusive Economic Zone (Fig. 1, Fig. SB1).The real-time data from the buoy were transmitted to the FIO laboratory every 10 min via satellites.
On 22 November, six Air Launched Autonomous Micro Observer (ALAMO) profiling floats (Jayne and Bogue 2017) from MRV Systems [one with SBE-41 conductivitytemperature-depth (CTD) sensor and five with RBR CTD sensors], and two Electromagnetic Autonomous Profiling Explorer (EM-APEX) floats were deployed near the buoy, at a radius of approximately 5 km.Two ALAMO floats failed within 2 days of deployment.The ALAMO floats profiled the temperature and salinity only during ascending phase, and the EM-APEX floats profiled temperature, salinity and current velocity during both ascending and descending phases (Sanford et al. 2005;Hsu et al. 2017).The vertical resolution of ALAMO float measurements using the RBR sensors was 0.1 m in the upper 5 m and 1 m down to 20 m, and that using the Seabird sensors was 1 m in the upper 20 m.The RBR CTD sensors allow temperature and salinity measurements within a few centimeters below the sea surface.The vertical resolution of EM-APEX float measurements was about 3.5 m (Sanford et al. 2011).The nominal temporal resolution between profiles on ALAMO floats and EM-APEX floats was about 3 h and 30 min, respectively.
Other datasets used in the study.Several other datasets are used to help interpret atmosphere and ocean processes during the field experiment.The NOAA extended reconstructed SST v5 is used to identify SST anomalies in the Indo-Pacific Ocean, as infrared The CSHOR-FIO field experiment off north Australia The CSHOR-FIO experiment was designed to use collocated surface meteorology buoy and fast-profiling floats, as a cost-effective observation strategy, to improve our understanding the air-sea coupling at the diurnal cycle and MJO time scales (Fig. SB1).In this study, the surface meteorology records and the sea surface temperature measurements at 1-m depth from the buoy were used to derive airsea fluxes.COARE3.0 bulk formula (Fairall et al. 2003) was used to calculate the latent and sensible heat fluxes using the buoy data.Net radiation fluxes were calculated by assuming the Stefan-Boltzmann constant of 5.6697 × 10 -8 W m -2 K -4 , longwave emissivity of 0.985, and surface albedo for shortwave radiation of 6%.SST measurements may significantly underestimate the SST cooling under convective systems (e.g., Duvel and Vialard 2007).The propagation of MJO convection ("General conditions during the campaign" section) is identified using daily outgoing longwave radiation (OLR) data from NOAA and precipitation data from the near-real-time legacy product of Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis (3B42RT), version 7 (Huffman and Bolvin 2018).The measurements of SST from Himawari-8 satellite and sea surface height anomalies (SSHA) from Integrated Marine Observing System (IMOS) data portal (http://imos.aodn.org.au/webportal) are used to compare with the float measurements ("Featured research studies in the experiment" section).More details of these datasets are described in the online supplementary text (https://doi .org/10.1175/BAMS-D-19-0278.2).

General conditions during the campaign
The global ocean was under the influence of a weak El Niño event in 2018/19 austral summer.Cool SST anomalies of up to −1.5°C were observed in the southeast Indian Ocean and the Indonesian-Australian Basin, whereas the central equatorial Pacific and the western Indian Ocean experienced warmer than usual conditions (Fig. 1, Fig. SB2).The cool SST anomalies in the southeast Indian Ocean started to fade around January 2019.The evolution of SST anomalies is consistent with a composite analysis of El Niño influences in the region (Zhang et al. 2017), that is, SST in the region is typically cooler (warmer) than normal by ~0.2°C (~0.5°C) before (after) an El Niño peak.
The intraseasonal SST and OL R a noma l ie s nor t h of Australia have standard deviations of about 0.4°C and Sea surface temperature variability in the Indonesian-Australian Basin and El Niño Intense SST variability exists in the region to the north of Australia.Warm SSTs are usually observed during and after the peak of an El Niño event in austral summer, such as during the 1997/98 and 2015/16 events (Fig. SB2), due to increased solar radiation and reduced latent heat release (Zhang et al. 2017).During the evolution of some of the El Niño events, SSTs are cooler than normal before the El Niño peak.and Hendon 2004) (www.bom.gov.au/climate/mjo/), or when the deep convection of the MJOs was located at the Maritime Continent.During the active phase of MJO1, air humidity has large fluctuations between the dry and wet states, and the wet state was fully established toward the end of December 2018 (Fig. 3, Fig. ES2).A low-humidity period during 16-20 December might be influenced by the dry southerly winds from the Australian continent, as shown in the Analysis Chart Archive from the Bureau of Meteorology (www.bom.gov.au/australia/charts/archive/index.shtml).There was a lowsurface-air-pressure signal near the end of December, but it was not associated with any local synoptic weather patterns, according to the Bureau of Meteorology chart.The surface wind measured at the buoy was generally from the west to southwest, except around 17 January 2019 at the onset of MJO2 convection, when easterly wind was observed, and during 25-28 January under the influence of TC Riley (SB3).The active phases of the MJOs were generally associated with low sea level pressure and high wind speed, which were preceded with lowwind periods of less than 5 m s -1 .After the passage of MJO2, the wind speed weakened until the end of February 2019.
Significant diurnal variations of SST and surface air temperature were observed during the low-wind periods, in the absence of wind stirring, the feature of MJO's suppressed phase (Fig. 3).Both the surface air temperature and SST steadily rose prior to the onset of the MJO1 convection, and then experienced a sharp drop of ~2°C during the MJO1 convective phase.The SST rose again at the buoy site after the passage of the MJO1 active phase, with peak SST reaching 30°C just before the arrival of the MJO2 active phase.In other words, the accumulation of heat (and associated moisture) preceded the onset of the MJO active phase.MJO2 and TC Riley then drove a ~3°C SST drop at the buoy site.

Featured research studies in the experiment
Strong DV SST events captured by fast profiling floats.A diurnal cycle of SST is observed when the absorption of insolation exceeding the outgoing heat near the ocean surface, in the absence of wind stirring (Fig. SB4).There were two strong DV SST events that occurred in December 2018, starting on 2 and 30 December 2018, respectively, with dSST > 1.5°C (Figs.3-5).Both events were

Tropical Cyclone Riley
The tropical disturbance Riley was identified by BoM in the Indonesian-Australian Basin around 15 January 2019 (www.bom.gov.au/cyclone/history/index.shtml).It became a category 2 tropical cyclone on 26 January 2019 (Fig. SB3).The storm initially moved southward, and then moved westward after 24 January, without making landfall in Australia.On 26 January, the storm passed the Bailong buoy, with a translation speed of ~6 m s -1 .The closest distance of the storm eye to the buoy was ~60 km, with a maximum wind speed of ~33.4 m s -1 , and a radius of maximum wind of ~40 km.Satellite imagary clearly shows the rainbands around TC Riley.The precipitation rate within ±3° of the FIO buoy position is measured to be more than 4 mm h −1 at its peak.Tropical cyclones can induce significant wind-driven current, causing strong shear instability and cooling of the upper ocean.associated with low wind speed and low cloud cover (Gentemann et al. 2003;Matthews et al. 2014), during the suppressed phases of the two MJO events.The DV SST was also high during the second half of February 2019, after the passage of MJO2; however, most of the profiling floats had stopped sampling.
On 3 December 2018, the significant DV SST event in the Indonesian-Australian Basin was captured by the Himawari-8 satellite measurements (Fig. 4).DV SST in regions covered by persistent clouds cannot be detected.Nevertheless, substantial spatial variations of DV SST amplitude were observed.The DV SST was high along the Australian coastal regions and extended to the northwest, covering the field work region; low DV SST regions existed to the east.Most of the dSST near the profiling floats were >1.5°C, consistent with the buoy (Fig. 3) and profiling floats (Fig. 5) measurements.Thus, the profiling floats sampled a region with relatively high DV SST in the Indonesian-Australian Basin on 3 December 2018.
The profiling float data not only accurately captured the DV SST variability, but also revealed the evolution of ocean mixed layer, which is critical for the DV SST development.(The base of the surface mixed layer is defined as the depth where difference between potential density ρ and the mean ρ in the upper 3 m is equal to or greater than 0.1 kg m -3 .)The mixed layer depth was about 50 m before the onset of the strong DV SST event in early December 2018 (Fig. 5).Shallowing of the mixed layer and strong DV SST developed within 2 days from 2 December 2018.The DV SST magnitude was more than 2°C, with peak SST reaching 29°C, from a background mixed layer temperature of just above 26°C (Fig. 4).After the first strong DV SST event, three more ALAMO floats failed before 21 December 2018, so that the second strong DV SST event of more than 2°C was only captured by a single ALAMO float, 9207 (Fig. 5).
Buoyancy frequency N 2 , which is proportional to the vertical gradient of potential density, was

Diurnal variation of SST and upper-ocean temperature
Before sunrise, net outgoing longwave radiation, latent heat flux, and sensible heat flux cool the near sea surface waters, and thereby the nighttime convective mixing exists above the seasonal thermocline (Fig. SB4).After sunrise, the absorption of insolation increases the SST to a peak (SST max ) in the afternoon.The penetrating solar radiation also increases temperature gradient near the surface and forms a diurnal warm layer (DWL).The DWL with a high gradient of potential density can stabilize the upper ocean and inhibit the turbulent mixing (Fairall et al. 2003;Brainerd and Gregg 1993).The nighttime convective mixing will repeat and erode the DWL after the sunset.The magnitude of DV SST is evaluated by computing the difference between daily maximum SST, SST max , and the foundation SST, SST fnd , which is defined using the mean SST before the sunrise (e.g., Bernie et al. 2005;Zhang et al. 2016).The window of averaging between 0000 and 0530 LT has negligible effects to the estimates of SST fnd (Karagali and Høyer 2014).calculated to explore the density stratification during the strong DV SST events (Fig. 5).High N 2 can inhibit the turbulent mixing.After sunrise, the heating from solar insolation started to form diurnal stratification (diurnal warm layer) near the ocean surface.The thermal stratification  prevented the mixing of ocean surface warm water with colder water in the subsurface layer, so that the shallow warm layer was insulated from the thermocline water below to provide a positive feedback to eventuate the strong DV SST event.The high vertical sampling resolution of the ALAMO floats was able to depict the fine warm-layer structure in the top meter of the upper ocean.The surface mixed layer shoaled by more than 40 m within 2 days during the two strong DV SST events, starting from 1 and 29 December 2018, respectively.Compared to the SST before the strong DV SST events, the daily mean SST at most floats increased by at least 1°C, and the higher SST was sustained for more than a week, presumably contributed by the increasing stratification in the upper ocean, prohibiting vertical turbulent mixing.Understanding the dynamics of the warm layer is critical for simulating the DV SST in numerical forecasting models.
Upper-ocean responses to the MJOs and the tropical cyclone.During December 2018-February 2019, there was net heat input into the ocean via the air-sea interface on most days except during the two short periods of the active MJOs, in mid-December 2018 and again in late January 2019 (Fig. 3).The average net heat flux into the ocean during December 2018-February 2019 was 121 W m -2 , which would warm a 50-m surface mixed layer by more than 4°C during a 3-month period.Based on the buoy observations, the surface temperature rose by less than 2°C.Thus, the heat input might be partly balanced by a northward wind-driven transport of cooler waters from the south, due to the prevailing westerly winds.The net heat flux into the ocean during the 2018/19 austral summer was somewhat higher than the climatology derived from OAFlux, a daily air-sea heat flux product (Fig. ES1; Yu et al. 2008), which might be associated with the cooler SST that enhanced ocean heat uptake.
The two active phases of the MJOs were associated with large latent heat releases from the ocean (Fig. 3).MJO1 induced a peak net heat loss of over 100 W m -2 on 16 December 2018.MJO2, in conjunction with the tropical cyclone, caused a heat loss over a 4-day period, with peak heat loss over 270 W m -2 for 2 days on 24-25 January 2019.There was ~2°C drop in SST in MJO1 and a 3°C drop in MJO2 in response to the surface heat losses (Fig. 3).A 1.5°C drop of SST was observed under the influence of a strong MJO event in conjunction with a tropical cyclone during a field campaign in the western Indian Ocean in early 2008 (Vialard et al. 2008(Vialard et al. , 2009)).Thus, the SSTs in the eastern Indian Ocean warm pool appear to be more susceptible to MJO disturbances, which is supported by the (~50%) higher ratio between intraseasonal SST and OLR variability in the Indonesian-Australian Basin, especially the southern half of the basin, compared with the western Indian Ocean (Fig. ES3; SST data are based on a high-resolution multisensor SST analysis; Beggs et al. 2011).Upper-ocean temperatures rose again after MJO2 toward the end of February 2019, which was captured by both the buoy and one ALAMO float.Prior to or during the active phase of the MJOs, precipitation and horizontal advection anomalies may freshen the ocean surface layer (Fig. 6).The fresher water may then be subducted underneath the relatively high-salinity water at ~60-m depth after the passage of the MJO convection.
The passage of the MJO convections brought strong westerly w inds to force eastward surface currents, which interacted with the rich mesoscale eddy field in the Indonesian-Australian Basin (Figs. 7 and 8).A cyclonic eddy persisted to the northeast of the deployment site in late November-early December 2018, inducing a predominant northward drift of the profiling floats (Fig. 7).The eddy weakened and shifted slightly to the west from early December.In the meantime, a high sea level anomaly center approached the area from the east, north of 15°S, inducing southeastward surface geostrophic current anomalies.In late December 2018, the cyclonic eddy seemed to reemerge, centered around 115°E.During the active phase of MJO1 from 13 to 20 December 2018, the EM-APEX float em8487 was located at the south of one anticyclonic eddy (Fig. 7).Colder water below the surface mixed layer was entrained into the surface layer during the MJO convective phase, likely due to the enhanced mixing induced by shear instability (Fig. 6), as the surface mixed layer was deepened from 10 to 50 m.This basic description of the oceanographic conditions during the field campaign will provide background information for future ocean mixing and air-sea coupling studies utilizing the dataset.
Inertial motions and internal tides.Vertical turbulence mixing plays a crucial role in the SST evolution.Ocean processes affecting the turbulence mixing include ocean stratification, internal waves (oscillations of thermocline depth at tidal or other frequencies), inertial waves (interior fluid motions using Earth rotation as restoring force, typically driven by surface winds), and background current shears.Both internal waves and inertial waves can propagate from their generation sites and cause enhanced turbulence mixing when they break.The roles of different processes in the SST budget on the MJO time scale have been assessed in other Indian Ocean regions (e.g., Cuypers et al. 2013); however, mixed layer processes have not been well observed in the Indonesian-Australian Basin, where the strong internal tidal energy may be an important contributing factor to the vertical mixing.Strong internal (tide) wave influence and water mass interleaving features were also observed in the float data (Fig. 6).
During the CSHOR-FIO field program, the two EM-APEX floats, em8487 and em8488, were deployed on 22 November 2018, and continuously obtained temperature, salinity, and current velocity profiles in the upper 300 m for 1.5 months (Fig. 7), providing a rich dataset to study the interplay among MJO forcing, mesoscale eddy, inertial waves, and internal waves in affecting the surface mixed layer processes in the Indonesian Australian Basin.The strong cold-core cyclonic eddy located to the east of the floats' initial positions drove the floats northward in the first 10 days.Being entrained into the cyclonic eddy, the drifting speed of float em8488 was much faster than em8487.After the cyclonic eddy dissipated in about 10 days, the floats were then affected by the positive sea level anomalies associated with the warm-core anticyclonic eddy in the north and started to drift eastward.The float motion in general agreed with the geostrophic current due to the low wind speed (i.e., weak wind-driven Ekman current) during most parts of the field experiment (Fig. 3).Most measured current velocity from em8487 was less than 0.3 m s -1 in the upper 50 m (red arrows in Figs.7 and 8), except before and during the active convection phase of MJO1, when the measured surface velocity steered to the north compared with the satellite-estimated geostrophic current (Fig. 7b).
The current velocity measurement from em8487 in successive ±2-day windows are used in a harmonic analysis (Fig. 8), which assumes the ocean current constituted by one inertial wave (period ≈ 43.5 h), one diurnal K 1 tide (period = 24 h) and one semidiurnal M 2 tide (period = 12.42 h).The diurnal (K 1 ) and semidiurnal (M2) tidal motion results from the gravity of the moon.The magnitude of inertial motions was generally less than 0.15 m s -1 .Besides the inertial waves in the upper 50 m generated by the local wind, internal waves below 100-m depth (i.e., below the seasonal thermocline) may be due to vertical propagations of internal waves cascaded from inertial waves generated from other locations.
The semidiurnal tides are significant, up to 0.15 m s -1 .The magnitude of semidiurnal tides was generally higher than the diurnal tides.The observation of stronger semidiurnal tides in the region was consistent with the previous observation and model results (Holloway et al. 2001;Rayson et al. 2011).Similar with the internal waves at the inertial frequency, the EM-APEX float measurements of semidiurnal internal tides may be used for exploring the tidal energy flux (e.g., Lien et al. 2013) around the Maritime Continent.It will also be interesting to study the breaking of internal waves, which would enhance the vertical mixing in the upper ocean (e.g., Cuypers et al. 2017).Vertical mixing due to wind-generated internal gravity waves could drive an upper-ocean heating rate up to 0.42°C month -1 , as observed in the western Indian Ocean (Cuypers et al. 2013).

Summary
The CSHOR-FIO atmosphere-ocean field campaign in the Indonesian-Australian Basin during the austral summer of 2018/19 collected rich surface meteorology and ocean profiling data from diurnal to intraseasonal time scales, to study the interplay among the MJO forcing, mesoscale eddy, inertial waves, and internal waves in affecting the surface mixed layer processes in the Indonesian-Australian Basin.The dataset will provide validations of high-resolution numerical models to study the coupled processes associated with the evolution of the MJO in the region.The success of the field campaign was built on a bilateral collaboration between CSHOR and FIO.Whereas the surface meteorology buoy provided realtime measurements of key air-sea exchange processes associated with the MJOs, tropical cyclones, and the transition between the Australian monsoon and the trade wind regimes, the influences from two MJO events and one tropical cyclone on the upper ocean were captured during the core period of the campaign (December 2018-February 2019) by most of the autonomous profiling instruments.The campaign was a unique combination of traditional buoy measurement and the new cost-effective profiling technology to a process study.The buoy system lasted well beyond the full transition from premonsoon onset through the end of the monsoon of 2018/19 austral summer.The data will be crucial to evaluate reanalysis products for the Indonesian-Australian Basin where large uncertainties of air-sea net heat flux exist (e.g., Yu et al. 2007).
The field campaign provides unique observations of simultaneous air-sea exchanges and upper-ocean variability in the Indonesian-Australian Basin for the first time, building upon the process studies in the Indo-Pacific warm pool started from the Coupled Ocean-Atmosphere Research Program (COARE) back in 1992/93 (Webster and Lukas 1992).It appears that SSTs in the region may be more susceptible to the MJO associated atmospheric forcing.The data can be used to explore multiscale processes from diurnal to MJO variability, from warm layer-mixed layer processes to their interactions with strong internal tide activities in the region.The collected data will be used to validate weather and climate forecast models, because previous studies often overlook the air-sea coupling processes in the Indonesian-Australian Basin, which has stronger MJO response compared to other tropical regions.The measurements can also improve our understanding on the oceanographic conditions for the relatively weak MJO variability, and the delayed onset of the Australian summer monsoon rainfall in the tropical north region during the 2018/19 austral summer.The cool SST anomalies may stabilize the atmosphere and inhibit convection in the Indonesian-Australian Basin.The profiling data from ALAMO and EM-APEX floats, which are the main vehicle for the upper-ocean measurements, provide new tools for exploring mixed layer dynamics.This may pave the way of using cost-effective autonomous technology in process studies in other parts of the World Ocean.R/V Investigator revisited the air-sea coupling processes in the eastern part of the Indonesian-Australian Basin during the 2019/20 austral summer, cocontributed from the CSHOR project, with additional atmospheric sounding measurements and more ocean profiling including direct turbulence measurements.The combination of these measurements would allow us to improve the understanding of the coupled ocean-atmosphere processes in the Indonesian Australian Basin.

Fig. SB1 .
Fig. SB1.A map of the CSHOR-FIO experiment.The offshore support vessel Limitless was chartered from Exmouth, Western Australia, during 19-26 Nov 2018, to carry out the field experiment centered at 16°51.432'S, 115°13.314'E(purple dot).As shown in the inset, one surface meteorology and ocean profiling buoy from the FIO was deployed at the center (red dot), and six ALAMO floats (black dots) and two EM-APEX floats (green dots) were deployed surrounding the buoy.

Fig. 2 .
Fig. 2. (a) Daily OLR anomalies and (b) TRMM 3B42RT, version 7, precipitation rates averaged over 20°S-10°N in the Indo-Pacific Ocean during November 2018-February 2019.The solid red lines denote the longitude of the Bailong buoy system.The ellipses denote the active phases of the two MJOs.(c) The time series of satellite measurements averaged within 3° of the buoy position.

Fig. 3 .
Fig. 3. (a) Sea level pressure, (b) wind speed and direction at 4-m height above the sea surface, (c) air temperature at 3-m height above the sea surface and ocean temperatures at 1-and 10-m depth below the sea surface, (d) specific humidity at 3-m height above the sea surface (blue line) and at sea surface derived from SST (red line), and (e) shortwave radiation, latent heat flux, and net air-sea heat flux, measured on the Bailong buoy.The net heat flux is derived from the Bailong buoy surface data (based on COARE algorithm; Fairall et al. 2003).The sampling frequency is 30 min and in (e), the net heat flux is averaged daily.The gray bars denote the convective periods of the two MJO events.

Fig
Fig. SB3.(top) Himawari-8 satellite image on 25 Jan 2019, showing Tropical Cyclone Riley off the northwest coast of Australia.(bottom) The trajectory of TC Riley (time interval between dots is 1 day) off northern Australia.The storm passed the buoy (red asterisk) on 26 Jan 2019, with a maximum wind speed V max (colors) of more than 30 m s -1 .The satellite image was obtained from the Japan Meteorology Agency website: https://www.data.jma.go.jp/mscweb/data/himawari/index.html.

Fig. SB4 .
Fig. SB4.Schematic plot of DV SST (red line) and potential density profiles in the upper ocean (blue lines).During a diurnal cycle, peak SST (SST max = mean SST from 1200 to 1700 LT) typically occurs in the local afternoon.The magnitude of DV SST is quantified by computing the dSST, the difference between the daytime peak and nighttime foundation SST (SST fnd ) during a diurnal cycle.In this study, we have used the mean SST from 1200 to 1700 LT as SST max , and from 0000 to 0500 LT as SST fnd .

Fig. 4 .
Fig. 4. Trajectories of six ALAMO floats (blue lines in the inset) and two EM-APEX floats (green lines in the inset) from 22 Nov to 3 Dec 2018.The background map shows dSST on 3 Dec, estimated using the Himawari-8 satellite observations.SST in regions covered by clouds cannot be detected by the infrared measurement.

Fig. 6 .
Fig. 6.Measured (a) temperature and (b) salinity from ALAMO float 9207.The black dashed lines denote the surface mixed layer depth.

Fig. 7 .
Fig. 7. Trajectories of the two EM-APEX floats (interval between dots = 1 day), em8487 (blue line), and em8488 (cyan line) on maps of 10-day-averaged sea surface height anomalies (shading) and geostrophic current anomalies (black arrows) derived from satellite altimeter data (obtained from IMOS) during November-December 2018.The red arrows are measured current velocity from the EM-APEX floats averaged in the upper 50-m depth.

Fig. 8 .
Fig. 8. (a),(e) Measurements of current velocity taken by float em8487.The measured current velocity is decomposed with harmonic analysis within a successive window of ±2 days, assuming the current is constituted by (b),(f) one inertial wave, (c),(g) one diurnal tide, and (d),(h) one semidiurnal tide.Columns show the (left) zonal and (right) meridional components.The black dashed lines and dots denote the mixed layer depth and missing measurements, respectively.