Responses of the Southern Ocean mixed layer depth to the eastern and central Pacific El Niño events during austral winter

Based on the Ocean Reanalysis System version 5(ORAS5)and the fifth-generation reanalysis datasets derived from European Centre for Medium-Range Weather Forecasts(ERA5),we investigate the different impacts of the central Pacific(CP)El Niño and the eastern Pacific(EP)El Niño on the Southern Ocean(SO)mixed layer depth(MLD)during austral winter.The MLD response to the EP El Niño shows a dipole pattern in the South Pacific,namely the MLD dipole,which is the leading El Niño-induced MLD variability in the SO.The tropical Pacific warm sea surface temperature anomaly(SSTA)signal associated with the EP El Niño excites a Rossby wave train propagating southeastward and then enhances the Amundsen Sea low(ASL).This results in an anomalous cyclone over the Amundsen Sea.As a result,the anomalous southerly wind to the west of this anomalous cyclone advects colder and drier air into the southeast of New Zealand,leading to surface cooling through less total surface heat flux,especially surface sensible heat(SH)flux and latent heat(LH)flux,and thus contributing to the mix layer(ML)deepening.The east of the anomalous cyclone brings warmer and wetter air to the southwest of Chile,but the total heat flux anomaly shows no significant change.The warm air promotes the sea ice melting and maintains fresh water,which strengthens stratification.This results in a shallower MLD.During the CP El Niño,the response of MLD shows a separate negative MLD anomaly center in the central South Pacific.The Rossby wave train triggered by the warm SSTA in the central Pacific Ocean spreads to the Amundsen Sea,which weakens the ASL.Therefore,the anomalous anticyclone dominates the Amundsen Sea.Consequently,the anomalous northerly wind to the west of anomalous anticyclone advects warmer and wetter air into the central and southern Pacific,causing surface warming through increased SH,LH,and longwave radiation flux,and thus contributing to the ML shoaling.However,to the east of the anomalous anticyclone,there is no statistically significant impact on the MLD.

The Effect of Boreal Summer Intraseasonal Oscillation on Mixed Layer and Upper Ocean Temperature over the South China Sea

Intraseasonal oscillation of the mixed layer and upper ocean temperature has been found to occur over the South China Sea(SCS)in the summer monsoon season based on the multiple reanalysis and observational data in this study.The method of composite analysis and an upper ocean temperature equation assisted the analysis of physical mechanisms.The results show that the mixed layer depth(MLD)in the SCS has a significant oscillation with a 30-60 d period over the SCS region,which is closely related to boreal summer intraseasonal oscillation(BSISO)activities.The MLD can increase(decrease)during the positive(negative)phase of the BSISO and usually lags behind by approximately one-eighth of the lifecycle(5 days)of the BSISO-related convection.The BSISO may cause periodic anomalies at the air-sea boundary,such as wind stress and heat flux,so it can play a dominant role in modulating the variation in MLD.There also are significant intraseasonal seawater temperature anomalies in both the surface and subsurface layers of the SCS.In addition,during the initial phase of the BSISO,the temperature anomaly signals of the thermocline are obviously opposite to the sea surface temperature(SST),especially in the southern SCS.According to the results from the analysis of the temperature equation,the vertical entrainment term caused by BSISO-related wind stress is stronger than the thermal forcing during the initial stage of convection,and it is more significant in the southern SCS.

Influence of the upper mixed layer depth on Langmuir turbulence characteristics

The upper mixed layer depth(h)has a significant seasonal variation in the real ocean and the low-order statistics of Langmuir turbulence are dramatically influenced by the upper mixed layer depth.To explore the influence of the upper mixed layer depth on Langmuir turbulence under the condition of the wind and wave equilibrium,the changes of Langmuir turbulence characteristics with the idealized variation of the upper mixed layer depth from very shallow(h=5 m)to deep enough(h=40 m)are studied using a non-hydrostatic large eddy simulation model.The simulation results show that there is a direct entrainment depth induced by Langmuir turbulence(h_(LT))within the thermocline.The normalized depthaveraged vertical velocity variance is smaller and larger than the downwind velocity variance for the ratio of the upper mixed layer to a direct entrainment depth induced by Langmuir turbulence h/h_(LT)<1 and h/h_(LT)>1,respectively,indicating that turbulence characteristics have the essential change(i.e.,depth-averaged vertical velocity variance(DAVV)DADV for Langmuir turbulence)between h/h_(LT)<1 and h/h_(LT)>1.The rate of change of the normalized depth-averaged low-order statistics for h/h_(LT)<1 is much larger than that for h/h_(LT)>1.The reason is that the downward pressure perturbation induced by Langmuir cells is strongly inhibited by the upward reactive force of the strong stratified thermocline for h/h_(LT)<1 and the eff ect of upward reactive force on the downward pressure perturbation becomes weak for h/h_(LT)>1.Hence,the upper mixed layer depth has significant influences on Langmuir turbulence characteristics.

Mixed layer warming by the barrier layer in the southeastern Indian Ocean

The southeastern Indian Ocean is characterized by the warm barrier layer(BL)underlying the cool mixed layer water in austral winter.This phenomenon lasts almost half a year and thus provides a unique positive effect on the upper mixed layer heat content through the entrainment processes at the base of the mixed layer,which has not been well evaluated due to the lack of proper method and dataset.Among various traditional threshold methods,here it is shown that the 5 m fixed depth difference can produce a reliable and accurate estimate of the entrainment heat flux(EHF)in this BL region.The comparison between the daily and monthly EHF warming indicates that the account for high-frequency EHF variability almost doubles the warming effect in the BL period,which can compensate for or even surpass the surface heat loss.This increased warming is a result of stronger relative rate of the mixed layer deepening and larger temperature differences between the mixed layer and its immediate below in the daily-resolving data.The interannual EHF shows a moderately increasing trend and similar variabilities to the Southern Annular Mode(SAM),likely because the mixed layer deepening under the positive SAM trend is accompanied by enhanced turbulent entrainment and thus increases the BL warming.

Seasonal variability of the mixed layer depth determined using an improved maximum angle method in the Arctic basins

To investigate the spatiotemporal variations in the mixed layer depth(MLD)in the Arctic basins,a new criterion to determine the MLD,called the improved maximum angle method(IMAM),was developed.A total of 45123 potential density profiles collected using Ice-Tethered Profilers(ITPs)in the Arctic basins during 2005-2021 were used to demonstrate the method’s effectiveness.By comparing the results obtained by the fixed threshold method(FTM),percentage threshold method(PTM),and maximum gradient method(MGM)for profiles in the Canada Basin,Makarov Basin,and Eurasian Basin,we determined that the quality index(1.0 for perfect identification of the MLD)of the IMAM regarding the assessment of the MLD determination method reached 0.94,which is much greater than those of other criteria.Moreover,two types of the density profiles were identified based on the mixed layer development stage.The MLDs of the typical profiles determined using the IMAM were found to have better consistency with the original definition.By utilizing the new mixed layer criterion,the seasonal variations and regional differences in the MLD in the Arctic basins were analyzed.Spatially,the summer and winter MLDs in the Canada Basin were the shallowest(13.55 m in summer,26.76 m in winter)than those in the Makarov(29.51 m in summer,49.08 m in winter)and Eurasian(20.36 m in summer,46.81 m in winter)basins due to the stable stratification in the upper ocean and the subsequent small effects of dynamic and thermodynamic processes(wind-driven stirring and brine rejection)in the Canada Basin.Seasonally,in the three Arctic basins,the average MLD was shallowest(22.77 m)in summer;it deepened through autumn and reached a winter maximum(41.12 m).

On the subtropical Northeast Pacific mixed layer depth and its influence on the subduction

The present climate simulations of the mixed layer depth（MLD） and the subduction rate in the subtropical Northeast Pacific are investigated based on nine of the CMIP5 models. Compared with the observation data,spatial patterns of the MLD and the subduction rate are well simulated in these models. The spatial pattern of the MLD is nonuniform, with a local maximum MLD（〉140 m） region centered at（28°N, 135°W） in late winter. The nonuniform MLD pattern causes a strong MLD front on the south of the MLD maximum region, controls the lateral induction rate pattern, and then decides the nonuniform distribution of the subduction rate. Due to the inter-regional difference of the MLD, we divide this area into two regions. The relatively uniform Ekman pumping has little effect on the nonuniform subduction spatial pattern, though it is nearly equal to the lateral induction in values. In the south region, the northward warm Ekman advection（–1.75×10–7 K/s） controls the ocean horizontal temperature advection（–0.85×10–7 K/s）, and prevents the deepening of the MLD. In the ensemble mean, the contribution of the ocean advection to the MLD is about –29.0 m/month, offsetting the sea surface net heat flux contribution（33.9 m/month）. While in the north region, the southward cold advection deepens the MLD（21.4 m/month） as similar as the heat flux（30.4 m/month）. In conclusion, the nonuniform MLD pattern is dominated by the nonuniform ocean horizontal temperature advection. This new finding indicates that the upper ocean current play an important role in the variability of the winter MLD and the subduction rate.

Response of the mixed layer depth and subduction rate in the subtropical Northeast Pacific to global warming

The response of the mixed layer depth(MLD)and subduction rate in the subtropical Northeast Pacific to global warming is investigated based on 9 CMIP5 models.Compared with the present climate in the 9 models,the response of the MLD in the subtropical Northeast Pacific to the increased radiation forcing is spatially nonuniform,with the maximum shoaling about 50 m in the ensemble mean result.The inter-model differences of MLD change are non-negligible,which depend on the various dominated mechanisms.On the north of the MLD front,MLD shallows largely and is influenced by Ekman pumping,heat flux,and upper-ocean cold advection changes.On the south of the MLD front,MLD changes a little in the warmer climate,which is mainly due to the upper-ocean warm advection change.As a result,the MLD front intensity weakens obviously from 0.24 m/km to0.15 m/km(about 33.9%)in the ensemble mean,not only due to the maximum of MLD shoaling but also dependent on the MLD non-uniform spatial variability.The spatially non-uniform decrease of the subduction rate is primarily dominated by the lateral induction reduction(about 85%in ensemble mean)due to the significant weakening of the MLD front.This research indicates that the ocean advection change impacts the MLD spatially non-uniform change greatly,and then plays an important role in the response of the MLD front and the subduction process to global warming.

Simulation and future projection of the mixed layer depth and subduction process in the subtropical Southeast Pacific

The present climate simulation and future projection of the mixed layer depth(MLD)and subduction process in the subtropical Southeast Pacific are investigated based on the geophysical fluid dynamics laboratory earth system model(GFDL-ESM2 M).The MLD deepens from May and reaches its maximum(>160 m)near(24°S,104°W)in September in the historical simulation.The MLD spatial pattern in September is non-uniform in the present climate,which shows three characteristics:(1)the deep MLD extends from the Southeast Pacific to the West Pacific and leads to a"deep tongue"until 135°W;(2)the northern boundary of the MLD maximum is smoothly near 18°S,and MLD shallows sharply to the northeast;(3)there is a relatively shallow MLD zone inserted into the MLD maximum eastern boundary near(26°S,80°W)as a weak"shallow tongue".The MLD nonuniform spatial pattern generates three strong MLD fronts respectively in the three key regions,promoting the subduction rate.After global warming,the variability of MLD spatial patterns is remarkably diverse,rather than deepening consistently.In all the key regions,the MLD deepens in the south but shoals in the north,strengthing the MLD front.As a result,the subduction rate enhances in these areas.This MLD antisymmetric variability is mainly influenced by various factors,especially the potential-density horizontal advection non-uniform changes.Notice that the freshwater flux change helps to deepen the MLD uniformly in the whole basin,so it hardly works on the regional MLD variability.The study highlights that there are regional differences in the mechanisms of the MLD change,and the MLD front change caused by MLD non-uniform variability is the crucial factor in the subduction response to global warming.

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor- Yamada level 2.5 turbulence closure scheme, was implemented to investi- gate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic ener- gy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corre- sponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the tempera- ture gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.

Simulation and Exploration of the Mechanisms Underlying the Spatiotemporal Distribution of Surface Mixed Layer Depth in a Large Shallow Lake

The aquatic eco-environment is significantly affected by temporal and spatial variation of the mixed layer depth （MLD） in large shallow lakes. In the present study, we simulated the three-dimensional water temperature of Taihu Lake with an unstructured grid with a finite-volume coastal ocean model （FVCOM） using wind speed, wind direction, short-wave radiation and other meteorological data measured during 13-18 August 2008. The simulated results were consistent with the measurements. The temporal and spatial distribution of the MLD and the possible relevant mechanisms were analyzed on the basis of the water temperature profile data of Taihu Lake. The results indicated that diurnal stratification might be established through the combined effect of the hydrodynamic conditions induced by wind and the heat exchange between air and water. Compared with the net heat flux, the changes of the MLD were delayed approximately two hours. Furthermore, there were significant spatial differences of the MLD in Taihu Lake due to the combined impact of thermal and hydrodynamic forces. Briefly, diurnal stratification formed relatively easily in Gonghu Bay, Zhushan Bay, Xukou Bay and East Taihu Bay, and the surface mixed layer was thin. The center of the lake region had the deepest surface mixed layer due to the strong mixing process. In addition, Meiliang Bay showed a medium depth of the surface mixed layer. Our analysis indicated that the spatial difference in the hydrodynamic action was probably the major cause for the spatial variation of the MLD in Taihu Lake.

A study of the mixed layer of the South China Sea based on the multiple linear regression

Multiple linear regression （MLR） method was applied to quantify the effects of the net heat flux （NHF）, the net freshwater flux （NFF） and the wind stress on the mixed layer depth （MLD） of the South China Sea （SCS） based on the simple ocean data assimilation （SODA） dataset. The spatio-temporal distributions of the MLD, the buoyancy flux （combining the NHF and the NFF） and the wind stress of the SCS were presented. Then using an oceanic vertical mixing model, the MLD after a certain time under the same initial conditions but various pairs of boundary conditions （the three factors） was simulated. Applying the MLR method to the results, regression equations which modeling the relationship between the simulated MLD and the three factors were calculated. The equations indicate that when the NHF was negative, it was the primary driver of the mixed layer deepening; and when the NHF was positive, the wind stress played a more important role than that of the NHF while the NFF had the least effect. When the NHF was positive, the relative quantitative effects of the wind stress, the NHF, and the NFF were about i0, 6 and 2. The above conclusions were applied to explaining the spatio-temporal distributions of the MLD in the SCS and thus proved to be valid.

Interannual variability of mixed layer depth and heat storage of upper layer in the tropical Pacific Ocean

By using the upper layer data(downloaded from the web of the Scripps Institution of Oceanography),the interannual variability of the heat storage of upper layer(from surface to 400 m depth) and the mixed layer depth in the tropical Pacific Ocean are investigated. The abnormal signal of the warm event comes from the central and west Pacific Ocean, whereas it is regarded that the abnormal signal of the warm event comes from the east Pacific Ocean in the popular viewpoint. From the viewpoint on the evolution of the interannual variability of the mixed layer depth and the heat storage of the whole upper layer, the difference between the two types of El Nino is so small that it can be neglected. During these two El Nino/La Nina events(1972/1973 and 1997/1998), other than the case of the heat storage or for the mixed layer depth, the abnormal signal propagates from the central and west Pacific Ocean to the east usually by the path along the equator whereas the abnormal signal propagates from the east to the west by the path northern to the equator. For the interannual variability, the evolution of the mixed layer depth corresponds to that of the heat storage in the upper layer very well. This is quite different from the evolution of seasonality.

Wave breaking on turbulent energy budget in the ocean surface mixed layer

As an important physical process at the air-sea interface, wave movement and breaking have a significant effect on the ocean surface mixed layer (OSML). When breaking waves occur at the ocean surface, turbulent kinetic energy (TKE) is input downwards, and a sublayer is formed near the surface and turbulence vertical mixing is intensively enhanced. A one-dimensional ocean model including the Mellor-Yamada level 2.5 turbulence closure equations was employed in our research on variations in turbulent energy budget within OSML. The influence of wave breaking could be introduced into the model by modifying an existing surface boundary condition of the TKE equation and specifying its input. The vertical diffusion and dissipation of TKE were effectively enhanced in the sublayer when wave breaking was considered. Turbulent energy dissipated in the sublayer was about 92.0% of the total depth-integrated dissipated TKE, which is twice higher than that of non-wave breaking. The shear production of TKE decreased by 3.5% because the mean flow fields tended to be uniform due to wave-enhanced turbulent mixing. As a result, a new local equilibrium between diffusion and dissipation of TKE was reached in the wave-enhanced layer. Below the sublayer, the local equilibrium between shear production and dissipation of TKE agreed with the conclusion drawn from the classical law-of-the-wall (Craig and Banner, 1994).

Role of Horizontal Density Advection in Seasonal Deepening of the Mixed Layer in the Subtropical Southeast Pacific

The mechanisms behind the seasonal deepening of the mixed layer （ML） in the subtropical Southeast Pacific were investigated using the monthly Argo data from 2004 to 2012. The region with a deep ML （more than 175 m） was found in the region of （22°-30°S, 105°-90°W）, reaching its maximum depth （-200 m） near （27°-28°S, 100°W） in September. The relative importance of horizontal density advection in determining the maximum ML location is discussed qualitatively. Downward Ekman pumping is key to determining the eastern boundary of the deep ML region. In addition, zonal density advection by the subtropical countercurrent （STCC） in the subtropical Southwest Pacific determines its western boundary, by carrying lighter water to strengthen the stratification and form a ＂shallow tongue＂ of ML depth to block the westward extension of the deep ML in the STCC region. The temperature advection by the STCC is the main source for large heat loss from the subtropical Southwest Pacific. Finally, the combined effect of net surface heat flux and meridional density advection by the subtropical gyre determines the northern and southern boundaries of the deep ML region： the ocean heat loss at the surface gradually increases from 22~S to 35~S, while the meridional density advection by the subtropical gyre strengthens the strat- ification south of the maximum ML depth and weakens the stratification to the north. The freshwater flux contribution to deepening the ML during austral winter is limited. The results are useful for understanding the role of ocean dynamics in the ML formation in the subtropical Southeast Pacific.

Parameterization for the Depth of the Entrainment Zone above the Convectively Mixed Layer

It has been noted that when the convective Richardson number Ri* is used to characterize the depth of the entrainment zone, various parameterization schemes can be obtained. This situation is often attributed to the invalidity of parcel theory. However, evidence shows that the convective Richardson number Ri* might be an improper characteristic scaling parameter for the entrainment process. An attempt to use an innovative parameter to parameterize the entrainment-zone thickness has been made in this paper. Based on the examination of the data of water-tank experiments and atmospheric measurements, it is found that the total lapse rate of potential temperature across the entrainment zone is proportional to that of the capping inversion layer. Inserting this relationship into the so-called parcel theory, it thus gives a new parameterization scheme for the depth of the entrainment zone. This scheme includes the lapse rate of the capping inversion layer that plays an important role in the entrainment process. Its physical representation is reasonable. The new scheme gives a better ordering of the data measured in both water-tank and atmosphere as compared with the traditional method using Ri*. These indicate that the parcel theory can describe the entrainment process suitably and that the new parameter is better than Ri*.

Contrasting dynamic characteristics of shear turbulence and Langmuir circulation in the surface mixed layer

Large eddy simulation （LES） is used to investigate contrasting dynamic characteristics of shear turbulence （ST） and Langmuir circulation （LC） in the surface mixed layer （SML）. ST is usually induced by wind forcing in SML. LC can be driven by wave-current interaction that includes the roles of wind, wave and vortex forcing. The LES results show that LC suppresses the horizontal velocity and greatly modifies the downwind velocity profile, but increases the vertical velocity. The strong downweUing jets of LC accelerate and increase the downward transport of energy as compared to ST. The vertical eddy viscosity Km of LC is much larger than that of ST. Strong mixing induced by LC has two locations. They are located in the 26s-36s （Stokes depth scale） and the lower layer of the SML, respectively. Its value and position change periodically with time. In contrast, maximum Km induced by ST is located in the middle depth of the SML. The turbulent kinetic energy （TKE） generated by LC is larger than that by ST. The differences in vertical distributions of TKE and Krn are evident. Therefore, the parameterization of LC cannot be solely based on TKE. For deep SML, the convection of large-scale eddies in LC plays a main role in downward transport of energy and LC can induce stronger velocity shear （S2） near the SML base. In addition, the large-scale eddies and Sz induced by LC is changing all the time, which needs to be fully considered in the parameterization of LC.

The spatiotemporal variation of the wind-induced near-inertial energy flux in the mixed layer of the South China Sea

On the basis of the QSCAT/NCEP blended wind data and simple ocean data assimilation （SODA）, the wind-induced near-inertial energy flux （NIEF） in the mixed layer of the South China Sea （SCS） is estimated by a slab model, and the model results are verified by observational data near the Xisha Islands in the SCS. Then, the spatial and temporal variations of the NIEF in the SCS are analyzed. It is found that, the monthly mean NIEF exhibits obvious spatial and temporal variabilities, i.e., it is large west of Luzon Island all the year, east of the Indo-China Peninsula all the year except in spring, and in the northern SCS from May to Septem- ber. The large monthly mean NIEF in the first two zones may be affected by the large local wind stress curl whilst that in the last zone is probably due to the shallow mixed layer depth. Moreover, the monthly mean NIEF is relatively large in summer and autumn due to the passage of typhoons. The spatial mean NIEF in the mixed layer of the SCS is estimated to be about 1.25 mW/m2 and the total wind energy input from wind is approximately 4.4 GW. Furthermore, the interannual variability of the spatial monthly mean NIEF and the Nifio3.4 index are negatively correlated.

Factors influencing the climatological mixed layer depth in the South China Sea:numerical simulations

The mixed layer depth （MLD） in the upper ocean is an important physical parameter for describing the upper ocean mixed layer. We analyzed several major factors influencing the climatological mixed layer depth （CMLD）, and established a numerical simulation in the South China Sea （SCS） using the Regional Ocean Model System （ROMS） with a high-resolution （1/12~x 1/12~） grid nesting method and 50 vertical layers. Several ideal numerical experiments were tested by modifying the existing sea surface boundary conditions. Especially, we analyzed the sensitivity of the results simulated for the CMLD with factors of sea surface wind stress （SSWS）, sea surface net heat flux （SSNHF）, and the difference between evaporation and precipitation （DEP）. The result shows that of the three factors that change the depth of the CMLD, SSWS is in the first place, when ignoring the impact of SSWS, CMLD will change by 26% on average, and its effect is always to deepen the CMLD; the next comes SSNHF （13%） for deepening the CMLD in October to January and shallowing the CMLD in February to September; and the DEP comes in the third （only 2%）. Moreover, we analyzed the temporal and spatial characteristics of CMLD and compared the simulation result with the ARGO observational data. The results indicate that ROMS is applicable for studying CMLD in the SCS area.

Mixed layer in the Sea of Japan:numerical simulation andlong-term data analysis

Seasonal variation and topography of the mixed layer in the Sea of Japan are studied by comparison of results from long-term observation data analysis and from numerical simulation with the MHI oceanic model (Shapiro. 1998. Marine Hydrophysical Journal, 6: 26-40). The data are retrieved from Oceanographic A tlas of the Bering Sea, Okhotsk Sea, and Japan/East Sea (Rostov, Rostov, Dmitrieva, et al. 2003. Pacific Oceanography, 1(1):70-72). The simulated and long-term patterns are compared. An impact of surface buoyancy flux, wind, and convergence/divergence of surface currents upon the mixed layer in the Sea of Japan is analyzed.

Changes in Mixed Layer Depth and Spring Bloom in the Kuroshio Extension under Global Warming

The mixed layer is deep in January-April in the Kuroshio Extension region. This paper investigates the response in this region of mixed layer depth （MLD） and the spring bloom initiation to global warming using the output of 15 models from CMIP5. The models indicate that in the late 21st century the mixed layer will shoal and the MLD reduction will be most pronounced in spring at about 33~N on the southern edge of the present deep-MLD region. The advection of temperature change in the upper 100 m by the mean eastward flow explains the spatial pattern of MLD shoaling in the models. Associated with the shoaling mixed layer, the onset of spring bloom inception is projected to advance due to the strengthened stratification in the warming climate.