Effect of heat on stratification hindering convection/ DSW formation

[willaguiar]

I was looking at profiles of $\sigma_0$ and $\delta \sigma_0 / \delta z$ for the 1mtop and 5mtop simulation, and noticed that at the upper 50 m, 1mtop is less stratified than 5mtop (below).

Perhaps that stronger stratification in 5mtop requires more energy to break, hindering convection and DSW formation on 5mtop?

That would add up to the fact that 5mtop has lower surface density, so it would have even lower potential energy to break stratification.

I cant explain why 5mtop would be more stratified tho.

[adele-morrison]

It would be interesting to see how these differences develop, eg what do these plots look like in the first month? Does the difference start from the surface and propagate downward?

[willaguiar]

Good point. I was just looking at that. Below are delta values (5m - 1m). each color is a month from IC, 1=Jan:

Looking at the $\Delta \sigma_0$.... perhaps a thinner cell is more efficiently stirred/mixed laterally by wind, making the surface denser. In turn, the thicker cell has move volume, is less efficiently mixed laterally by winds and some of its high density signal travels down to the bottom of the mixed layer ….. making the 5mtop more stratified?

[AndyHoggANU]

OK, so within the first month there is significantly lighter water in the upper 15m of the 5mtop case. And maybe more dense water in 15-50m range? That seems more like a difference in the vertical transport. I hate to say it, but a @PaulSpence -inspired salt budget over this first month (maybe at daily resolution?) might be informative??

[adele-morrison]

That's really interesting to see! It definitely starts at the surface in summer in the first month. It propagates down all of a sudden in May - I wonder if that's when some winter convection starts to happen normally in the control and in the 5mtop case the convection is weakened due to the already enhanced surface stratification.

Maybe going even closer to the start would be interesting, i.e. repeating the plot on the right of density profile anomalies using daily data for the first 10 days or so? I wonder if we'd see it starting constrained even more to the top layers?

[willaguiar]

Just updating here the same plots for panan 5mtop - 1mtop. The behavior is fairly similar, but perhaps more concentrated on the surface, and with a smaller intensity

[willaguiar]

I was thinking.......
As sea ice melts, it would add freshwater and volume on the surface cell of the model. If I remember correctly from Andy 's $z*$ presentation, increases in SSH are distributed along $z*$ levels as changes in $dzt$, proportionally to the mean cell thickness. Perhaps that means that 5m cell can hold bigger $dzt$ increase and bigger increase in FW content than 1m cell.

Complementing that..... If I remember correctly, if the $z^*$ cell grows too much, then the extra grow relative to the maximum $z^{target}$ is used to calculate the vertical flux/velocity. I imagine it is easier to reach that maximum volume increase in the 1m cell, and thus its vertical velocities are likely bigger too ( exporting more of the extra FW to deeper levels and not a allowing to much FW to build up at surface).

This is a guess made knowing very little of the $z*$ vertical coords, so one question is....Does that make any sense?

[willaguiar]

Based on the density anomaly/stratification plots above, We thought it would be useful to see how the salinity over the shelf evolves in the few first months of the simulation + wether the freshwater transport to the shelf was altered between the simulations.

For clarification, I calculated the freshwater transport as the residual of the mass transport (Qx, and Qy) minus the salt content , as in

where S is the absolute salinity. And after that I extracted the freshwater transport across the 1km isobath.

First, we can see that there is a continuous decrease in salinity over the top 10 m meters of the shelf in the 5m model [c] in 90 days.
When we compute the freshwater transport across the 1km isobath, we can see below that right at the surface the freshwater transport of the 5m case increases a lot ([b] compared to [a]). Even when integrated over an average Ekman Layer depth of 35 m, the 5m model has stronger southward freshwater transport [d].

This only clarifies that DSW shutdown is likely due to stronger Southward freshwater transport in the 5m. But that could still be either by differences in Ekman pumping/transport or differences in ice melting I guess.

[willaguiar]

Update on Freshwater/Salt effect- part 1

Context: We know the MOM5_5m simulation has reduced DSW formation relative to MOM5_1m, because MOM5_5m shelf is fresher at surface, then resisting DSW formation by buoyancy fluxes. This is because MOM5_5m shelf gets very fresh/light at the surface in the first 30-60 days prior to form DSW (plot below).

The goal to evaluate the salt budget in the next comment, is to figure out what causes this freshening at the surface.

---

Salinity changes due to salt or freshwater content
Salinity (S) can change either due to changing freshwater content, or due to changing salt content. With the DS/Dt output (salt tendency), I can calculate the daily salt content in each cell ($SC_{shelf}$). The daily mass (M) of each cell will be $M = SC_{shelf} \times$ $(S*1e-3)^{-1}$ , and the daily freshwater content of each cell will be the difference $FW_{shelf} = M - SC_{shelf}$ . With these values, I can calculate the salinity changes due to changes in salt content ($S^{salt}$) or changes in freshwater content ($S^{FW}$) as:

A comparison of these three could give us a hint on whether salt or freshwater content changes determines the surface freshening in MOM5_5m. The biggest freshening on the shelf happens in the upper 10 m of the water column (below), so I focused on this layer.

Below is $S^{salt}$ and $S^{FW}$ for the upper 10m, and $S^{salt}$ explains almost all the salinity changes. Therefore, the salt content budget might help us figure out the reason for the surface shelf freshening in MOM5_5m.

[willaguiar]

Sal budget on the shelf

For the salt budget I output the following terms, and summed them fluxes on the top 10m of the shelf:

With $F_{surf}$ being the sum of the outputs: 'sfc_salt_flux_ice' + 'sfc_salt_flux_restore'

Below is the salt budget terms for MOM5_1m (a,b) and MOM5_5m (b,c), and the 5m-1m difference (d,e) The dominant term in salt removal in MOM5_1m is the advection, which is approx. balanced by salt gain by the implicit diffusion (a,b). Both advective salt loss and the salt gain by implicit diffusion decreases in MOM5_5m (b,e), but the compensation is not perfect (perhaps due to increased stratification), resulting in stronger salt removal in the first 30 days.

---

• Bias warning: I can’t seem to close the salt budget, i.e., DS/Dt do not match the right side of the equation. For context, below is Fsurf calculated from these flux diagnostics, and from the residual of the equation in the 1m case, and the residual is ~20 times smaller than the surface fluxes, so be mindful of that.But, the surface forcing difference between simulations still very small and and can’t explain salinity changes, independent if I compare the diagnostics or residual term. Meanwhile, any suggestion on how to close the salt budget is welcome.

[adele-morrison]

For dS/dt are you using the online tendency term?

[willaguiar]

For dS/dt are you using the online tendency term?

Yes - for dS/dt I used the salt_tendency output

[adele-morrison]

FYI, we solved the budget closure (salt_vdiffuse_impl already includes the surface fluxes).

[willaguiar]

Perhaps this plot below helps explain a little more our 5m signal…

Below are the change in SSH(a), Ekman pumping (b), vertical salt fluxes in the upper 20m (positive upwards, c,d), and salinity anomalies [3], all for the shelf only ( south of 1km isobath). All that for the absolute wind experiments. Blue/orange contours in d and e are 1m and 5m isopycnals respectively.

Here a interesting story seem to appear. In the 5m case we have a fast development of positive anomalies in SSH near the Antarctic Coast (a, isobath shallower than 400m), suggesting increase in Ekman Pumping. The calculated Ekman pumping in the same region seem to increase from 1m to 5m cases. The downwards salt advection also increases south of the 400m isobath (d,e) where the Ekman pumping increased, decreasing the salinity in the upper shelf, and making the shelf densities decrease (contours in d,e).

This suggest that an increase in Ekman pumping, pumps out salt out of the surface, making these waters too fresh and DSW resistant.

[willaguiar]

When we calculate the Ekman transport across the 1km isobath ( offline, assuming boussinesq, + southward), we see that the Ekman transport actually increases in the first month, which would increase the SSH at the coast, and increase the Ekman pumping… I think this ties up with our hypothesis than changes in Ekman pumping make the surface waters fresh and DSW resistant. The Ekman transport changes are small, but since we only change the top 5m of the grid, the pumping change likely happens only in this thin 5m slice ( ~15% of Ekman layer depth). So they might be important, even tho small.

---

Some things to consider in the explanations above are:
1- We still don’t know why Ekman transport would increase in the 5m case with absolute winds. Especially since the upepr cell is slower in the 5m case( I have some guesses for it)
2- This is a simplistic story, as we see that there is a balance between Implicit diffusion and advection in the top 20m of the ocean ( check comment with salt balance), so advection alone might not be the best approach to explain things. I am remaking these calculations for a thicker layer ( perhaps about 50m or 100m) where implicit diffusion might not be as important. If the relationship holds, then perhaps we can reserve our explanation to only advection