World Ocean Thermocline Weakening and Isothermal Layer Warming

This paper identifies world thermocline weakening and provides an improved estimate of upper ocean warming through replacement of the upper layer with the fixed depth range by the isothermal layer, because the upper ocean isothermal layer (as a whole) exchanges heat with the atmosphere and the deep layer. Thermocline gradient, heat flux across the air–ocean interface, and horizontal heat advection determine the heat stored in the isothermal layer. Among the three processes, the effect of the thermocline gradient clearly shows up when we use the isothermal layer heat content, but it is otherwise when we use the heat content with the fixed depth ranges such as 0–300 m, 0–400 m, 0–700 m, 0–750 m, and 0–2000 m. A strong thermocline gradient exhibits the downward heat transfer from the isothermal layer (non-polar regions), makes the isothermal layer thin, and causes less heat to be stored in it. On the other hand, a weak thermocline gradient makes the isothermal layer thick, and causes more heat to be stored in it. In addition, the uncertainty in estimating upper ocean heat content and warming trends using uncertain fixed depth ranges (0–300 m, 0–400 m, 0–700 m, 0–750 m, or 0–2000 m) will be eliminated by using the isothermal layer. The isothermal layer heat content with the monthly climatology removed (i.e., relative isothermal layer heat content) is calculated for an individual observed temperature profile from three open datasets. The calculated 1,111,647 pairs of (thermocline gradient, relative isothermal layer heat content) worldwide show long-term decreasing of the thermocline gradient and increasing of isothermal layer heat content in the global as well as regional oceans. The global ocean thermocline weakening rate is (−2.11 ± 0.31) × 10−3 (°C m−1 yr−1) and isothermal layer warming rate is (0.142 ± 0.014) (W m−2).