North Atlantic dense water formation

[adele-morrison]

Bit of a left field idea: is it worth checking how North Atlantic dense water formation responds to changes in vertical resolution?

i.e. for an initial look, plot zonal average overturning in depth or density coordinates in the Atlantic.

[PaulSpence]

Definitely a good idea. We have ignored the NAtl overflows for too long! Assigned :)

[willaguiar]

Ok, I did a quick check of the Overturning in the last year of each simulation, to see how the N.A .respond to thickening the top cell.
Figures are...

Overturning contours:
Top: Overturning in the last year of control (1 meter top cell)
Middle: Anomaly (5mtop - control) on colors, lines are the same as top plot for the control
Bottom: Anomaly (halfmtop - control) on colors, lines are the same as top plot for the control

Line plot:
[a] maximum overturning at 26N, for $\sigma_2>1036$, i.e., NADW flow
[b] Absolute minimum overturning at 60S, for $\sigma_2>1037.05$, i.e., NADW flow

Interestingly it seems that the effect of dz in the North Atlantic overturning is reversed, i.e., as we increase the thickness of the top cell, we end up having more NADW, and less AABW
Any ideas on what could be the difference here?

Next steps is to analyze MLD in the North Atlantic and SWMT to see how/where NADW is formed in the model. Might give us a hint on what to look for.

possible hypothesis: DSW formation in the South is driven by salt fluxes, while in North Atlantic deep convection is likely driven by heat fluxes. So perhaps the sensitivity of the top cell depends on what is driving the mass transformation at surface?

[wghuneke]

Thanks for looking into this. Very intriguing that the result is opposite in the NA - I agree that checking whether salinity vs temperature is dominating the convection might help understanding what is going on.

Just noting here that the observed AMOC at 26N is about 17 Sv (McCarthy et al., 2015, 2004-2012 average) and the values here are much higher - possibly not the exact same calculation. Anyway, just saying that it might be the case that in the NH, the model produces more deep convection with a thicker uppermost grid cell and thereby overshoots. So the result would be consistent with our Antarctic analysis, i.e. a thinner uppermost grid cell is desirable.

[willaguiar]

Good Point Wilma..... That overturning line-plot is measured for the maximum below $\sigma_2>1036$, but that is not really the centre of the NADW overturning cell (which is actually at $\sigma_2>1036.5$), on the b contour. When I do the same line plot for $\sigma_2>1036.5$ indeed we get ~17SV for halfmtop, ~19SV for control by the end of the last year (Which supports your idea)

[dkhutch]

Yes very interesting that it's like a bipolar seesaw type response. Is it just good luck that we get closer to observations this way? In my mind 19 Sv in the model is a respectable answer for Amoc and a bit of ad hoc tuning could change it either way. It's interesting to see this though. Yes would be great to get some idea of whether temp or salt fluxes are more in control in each location. I think we saw on Friday that the expansion of MLD regions in the Antarctic margins tended to be on the periphery of the existing Maxima locations. Whereas when it contracted, it tended to contract in the middle of the Maxima, right? Wonder if it's similar pattern in the north.

…

On Thu, 7 Dec 2023, 4:32 pm Wilton Aguiar, ***@***.***> wrote: Good Point Wilma..... That overturning line-plot is measured for the maximum below $\sigma_2&gt;1036$, but that is not really the centre of the NADW overturning cell (which is actually at $\sigma_2&gt;1036.5$), on the b contour. When I do the same line plot for $\sigma_2&gt;1036.5$ indeed we get ~17SV for halfmtop, ~19SV for control by the end of the last year (Which supports your idea) — Reply to this email directly, view it on GitHub <#32 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEVNSK2QLAQCGMHTYK5VE6TYIFIIRAVCNFSM6AAAAAA6HWM2AKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQNBUGUYTCOBVGM> . You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

[willaguiar]

here it goes some considerations of the NADW formation on the North Atlantic (SWMT and Mass transports).
Apologies in advance for the loads of plots

Based on winter Mixed Layer Depth (DJFMA, to pinpoint deep convection), I defined the following regions and boundaries for SWMT and to calculate the mass transports.

SWMT and Mass flow in the Nordic Seas:
For reference, NADW seem to have $\sigma_0$ approx $27.8$, and $\sigma_2&gt;36.7$, in the model.
Nordic Seas dense-SWMT seem to increase in $\sigma_0&gt;28$, around the densest layers ([a-c]). In turn the Southward transport across the Nordic seas array increases in the densest classes ($\sigma_2&gt;36.8$), but is a very small increase (I guess it makes sense, as these denser waters would struggle to cross the Iceland ridge)

SWMT and Mass flow in the Northeast Atlantic:
Here the SWMT creates waters that are too light to be NADW ( dense transformation is in $\sigma_0&lt;27.8$). We also don’t see significant changes in SWMT between simulations. So likely no significant changes in NADW formation happen here

SWMT and Mass flow in the Labrador:
Here there is a huge dense-SWMT occurring ( more than 30 SV), so likely the bulk of the NADW formation in the model happens here [a]. This is very common in ocean models as far as I know, but it doesn’t match much observations e.g., OSNAP shows stronger MOC flows on East than west.
When thickening the top cell, it seems that dense-SWMT increases for $\sigma_0&gt;27.86$, and decreases for $\sigma_0&lt;27.86$. So the results is not as simple as we saw for the Antarctic Shelf SWMT where thickening the top cell weakens the SWMT mostly in the densest water classes.

In maps with mean of the SWMT for $\sigma_0&gt;27.86$ in the simulations (a-f), and of anomalies in MLD (g), we can see that the region where we have increased dense transformation in the Labrador sea approximately matches where we have increased deep convection (a-b, mld in g). And that this is mostly driven by heat fluxes (c-d, in contrast with e-f, scale in e-f is 15% of c-d). So increasing the top cell thickness, increases and concentrates the heat-flux induced convection in that red blob in the Labrador sea, increasing the formation of the densest NADW.

Finally, the plot below is an Atlantic-only sum of the northward transports,. Southward NADW transport is negative in green in b ($\sigma_2&gt;36.7$), and in c anomalies are as 5mtop-1mtop over control transport contour lines. Yellow line 36.98. We can see a decrease in UNADW southward transport ($\sigma_2&lt;36.98$, red), and an increase in LNADW ($\sigma_2&gt;36.98$,blue). The sum of the transport in 45N and 26N (d), show that the total NADW transport actually doesn’t change that much (~ 1SV), but the proportion of UNADW vs LNADW does.

Summary:

Most of the NADW formation in the model happens through deep convection in the Labrador Sea. Thickening the model top cell causes stronger heat-flux-induced convection in a concentrated region in the Labrador sea, which creates denser NADW. Dense NADW ( + mixing likely?) increases NADW inflow slightly ( ~ 1 SV).

One question that remains… why thickening the top cell causes an increase in denser water transformation in the North, and decrease in the South? Perhaps this is related to transformation in the North being heat-driven, and in the south being salt-driven.

Ps:Van Roekel experiment with the + kpp bias in OSBL during free convection with thick top cells uses heat fluxes to drive the convection. So perhaps the issue here is the sensitivity of the kpp to heat fluxes that is not present for salt fluxes, because salt diffusion is slower.

Code

[PaulSpence]

Interesting. I think some people would be unimpressed with the lack of WMT/convection in the GIN sea.

[willaguiar]

Interesting. I think some people would be unimpressed with the lack of WMT/convection in the GIN sea.

I agree... probably the forcing might not be the best in the North Atlantic (perhaps not much of a model issue but more of a forcing issue??)

[wghuneke]

Great analysis. @willaguiar , did we run an experiment with different KPP parameters? Do you think it would be a helpful addition? (I forgot if we had a discussion around this already and what the arguments were for/against this experiment.)

[willaguiar]

did we run an experiment with different KPP parameters? Do you think it would be a helpful addition?

Thanks for the suggestion! I think we gave up on the idea after having similar results for MOM5-kpp and MOM6-epbl for DSW/AABW. But it might be worth to reconsider it, especially if we further want to explain what happens in the N.A.

[adele-morrison]

Nice to see Wilton! And curious that it’s so different to the Antarctic. Did you check that the Atlantic SWMT is heat not salt dominated?

if it is heat dominated: for a thick top cell, perhaps if it’s the top cell equilibrating with the atmosphere (i.e losing heat), can you do that more efficiently with a larger volume top cell? e.g. the model doesn’t have to wait another time step to advect more heat upwards to the top cell. Just brainstorming, not sure if that’s the correct way to think about it. @AndyHoggANU thoughts on what’s going on here?

[dkhutch]

I think what @adele-morrison suggests is certainly plausible, i.e. a thickening of the top cell could somehow enable more heat exchange, while the freshwater/salt exchange is more limited in the south. I think Wilton was saying on Friday that the northern DSW formation is heat dominated.

[willaguiar]

Good tip to look at heat and Salt transformation…:
Below is the comparison of SWMT from heat and salt between models. Left are the salt, heat and total components for the control experiments (1mtop), middle are the salt components for the 1mtop and 5mtop experiments, and right are the same as middle but for heat fluxes.
Top is Antarctic shelf panan (MOM6-EPBL), middle is Antarctic Shelf OM2 (KPP), and bottom is Labrador Sea OM2.

Salt and heat transformation response to top cell thickening in Antarctic Shelf is the same (decreases transformation, a to f). Antarctic heat and salt transformation response to top cell thickening is not sensitive to the mixing scheme (a,b,c compared to d,e,f). Still, increasing the top cell thickness in Labrador Sea favor (salt and heat driven) denser deep water transformation (and hinder lighter deep water transformation). So probably it isn’t the shift from salt to heat forced SWMT that causes the differing hemispheric responses.

Perhaps the different continent geographies is causing contrasting responses instead?

I.e., In the Labrador Sea convection site, thickening the top cell might lead to locations with lower heat loss (lower latitudes green blob below) to cease the deep convection. But then, because there isn’t a continent north of Labrador, the convection can still happen more at North, where heat loss would be more intense and capable of surpassing the top cell thickening effect. The only physical constrain here would be the presence of sea ice, and indeed the sea ice edge seem border the location where MLD increases (white line & orange blob, figs below). In these higher latitudes heat transformation now happens over thicker cells and with stronger heat loss, which could (possibly) convect more and denser waters locally.
In turn, In the Southern Hemisphere, the Antarctic continent prohibits convection to move southwards where heat loss would be stronger. Add to that: SWMT on OM2 and MOM6 Antarctic Shelf is driven by salt fluxes, so heat loss wouldn’t be a constraining factor anyways?