Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Storyboard #37

Open
wghuneke opened this issue Sep 9, 2021 · 78 comments
Open

Storyboard #37

wghuneke opened this issue Sep 9, 2021 · 78 comments

Comments

@wghuneke
Copy link
Collaborator

wghuneke commented Sep 9, 2021

Suggested figures that go into paper. Add example figure, short description, suggestions for panels/what needs to be changed...

@wghuneke
Copy link
Collaborator Author

wghuneke commented Sep 9, 2021

Figure 1: Experimental design (enhancing/relaxation of climatological wind field)
a) climatological wind field in JRA
b) historical trends in JRA and CMIP
c) projected trends in CMIP
d) applied perturbation

@adele-morrison
Copy link
Owner

I've also added a figure summary page here. I suggest we keep the discussion here in this issue, but can put the result of this discussion on that figure summary page, to keep a cleaner version and somewhere to have the figures in order.

@wghuneke
Copy link
Collaborator Author

wghuneke commented Sep 9, 2021

Sounds good. Thought it's easier to add figures in an issue.

@julia-neme
Copy link
Collaborator

Figure 1: Experimental design (enhancing/relaxation of climatological wind field)
a) climatological wind field in JRA
b) historical trends in JRA and CMIP
c) projected trends in CMIP
d) applied perturbation

For c), trends for the period 2015-2100, or perhaps 2100-2080 mean minus 2015-2005 mean?

@adele-morrison
Copy link
Owner

You can add figures to the summary page just the same as issues, I've added an example one in now:
https://github.com/adele157/easterlies-collaborative-project/blob/master/Figure_outline.md

This issue is good for discussion though, because we'll get email notifications and I'm guessing there will be lots of discussion to be had!

@adele-morrison
Copy link
Owner

Should we decide what journal we're aiming for as well? i.e. Are we aiming for 4 or 10 figures?

@adele-morrison

This comment has been minimized.

@matthew-england-unsw
Copy link
Collaborator

Yes agreed Adele - good call.

@julia-neme

This comment has been minimized.

@matthew-england-unsw

This comment has been minimized.

@matthew-england-unsw
Copy link
Collaborator

Should we decide what journal we're aiming for as well? i.e. Are we aiming for 4 or 10 figures?

I'd vote for 10 I think. Unless we make 10 and decide that 4 tell the essential story, and then we can move the other 6 to Suppl. Info?

@StephenGriffies
Copy link
Collaborator

I find it easier to write the long manuscript first to be sure the story is solid and complete. Then, could cut back to shorter GRL size if warranted/desired.

@wghuneke
Copy link
Collaborator Author

I suggest we need figures that show (somewhat in this order):

  • two-time scale response (e.g. using SSH): quick Ekman response in first year, reversal of SSH signal within a couple of years (we think due to sea ice export increasing the sea ice/DSW formation)
  • increase in DSW formation (e.g. bottom age, bottom temp/salt, map of surface water mass transformation rates)
  • sea ice changes #10, increase in formation, increase in export

@AndyHoggANU
Copy link
Collaborator

I've been thinking a bit about this figure order. I think we have two possibilities:

  1. Something linear, where we build up a picture gradually -- starting with temp and salt evolution, SSH response, sea ice evolution and export, and culminating in the DSW and off-shelf age changes. (Note that we may need a second set of figures showing different experiments to distinguish between hypotheses.)
  2. Show the headline result first (well, as figure 3, assuming a CMIP 6 winds + model configuration figure). This, I think, would be the DSW change/age change figure. Then, go back to investigating other quantities to disentangle hypotheses.

I quite like the second option for this paper, accepting that it is a little unusual. Any thoughts?

Also, I am not quite 100% convinced that we have fully disentangled the Ekman from the sea ice export hypothesis. Do we have a good experiment that could actually separate these two possibilities? That might be important for the storyline I am proposing...

@adele-morrison
Copy link
Owner

I like this suggestion Andy, but I would be tempted to show the SSH time series before the DSW plot. That way we could start off by saying that we might naively to expect the easterlies to pile up water on the shelf, and indeed we do see that initially, but it is quickly overwhelmed by a larger opposing response by the DSW. Then go into DSW mechanisms. I agree I don’t think we’ve distangled the mechanisms completely yet, let’s discuss today what analysis / simulations we need to do that properly.

@matthew-england-unsw
Copy link
Collaborator

I think that’s a great idea Andy if the aim is to make the paper headline about the DSW response. But over the wider ocean there’s a largely coherent circumpolar (and counter-intuitive) response that is pretty impressive too, that remains when we pulled away the DSW formation region wind anomalies. To me that’s the headline story, and the DSW changes are kind of secondary or at least consequent on that. But that could just be me, DSW enthusiasts would be taken in by the storyline you propose. :-)

@wghuneke
Copy link
Collaborator Author

wghuneke commented Sep 23, 2021

Thanks for updating Fig 1 @julia-neme, looks good!

Just a few minor things: Can you change the colorbar ticks (control) to a spacing of 1? Or maybe even 2 and for the anomaly 0.5 and then you can increase the font size for all the ticks and the other text. And maybe show a few less arrows, I find it looks slightly busy, but that might just be me.

@wghuneke
Copy link
Collaborator Author

Also adding a short summary on what we discussed today (mostly for @StephenGriffies ):

We continued our discussion about the paper outline, Adele and Andy updated the Figure_outline.md document accordingly.

We have three hypotheses for different mechanisms that could explain the observed response: A local wind (katabatics) mechanism, B Upwelling mechanism, C Sea ice mechanism. We can rule out A with the additional perturbation experiment we did. The next steps are to do a bit more analysis/run additional simulations to disentangle B and C. Atm, we think C is more relevant.

Adele is currently running two more simulations for the UP case: a "meridional winds only" and a "zonal winds only" case, where only one component is perturbed and the other is kept as in the control run. The idea is that it will change the curl, hence Ekman upwelling, but not so much the sea ice advection.

@matthew-england-unsw
Copy link
Collaborator

Thanks @wghuneke great summary from today. 👍

@matthew-england-unsw
Copy link
Collaborator

PS: @julia-neme given that the wind experiments are simply +10% and -10% south of the purple line, the 2 x anomaly plots on the right are maybe not needed, and so instead you could use the full 'real estate' of the page width to devote to a single circumpolar wind vector plot....? I'm also curious how this looks without uniform vectors (which can be deceiving in the 'blank' low wind regions). e.g. if taking up the full page width and using a proper vector scale, we will maybe see clearly the full wind field and then knowing the expts are +/- 10% anomalies on this mean field could be sufficient?

@adele-morrison
Copy link
Owner

Yes agreed Matt, I think we could lose the 2 right panels. Julia, what do you think about using the model land/sea mask too, which cuts out some more of the Ross and Weddell?

@julia-neme
Copy link
Collaborator

Yeah! Will update based on all these comments :) Thanks!

@julia-neme
Copy link
Collaborator

julia-neme commented Sep 27, 2021

Fig 1 updated

With input from @adele157 and @matthew-england-unsw, how about this Fig. 1:

fig1_experimentalsetup_v2_streamplot (1)
a) Wind streamlines colorcoded by wind velocity, 1000m isobath in black contour, zero divide line poleward of which we apply the perturbation in dark blue and 1000m, 3000m depth contours in light grey. The boxes enclose the masked out regions for the additional experiment. b) Zonal and c) meridional components of the wind, with 1000m isobath in black.

I've opted for streamlines since the quiver plots are a bit messy and difficult to interpret. The advantage of having the 3 panels is that a) makes it easier to understand the wind field, b) and c) allow to focus on the zonal/meridional components separately.

@julia-neme
Copy link
Collaborator

I had a go at Figure 2:
figure2
Top panel: monthly SSH anomalies. Bottom panel: monthly magnitude of the zonally integrated overturning streamfunction at 55S.

I'm not sure who wins the bet based on panel 2! Perhaps this is not the best way to measure bottom water transport? I'm currently working on trying to separate thermosteric from halosteric contributions (#44) to add to the top panel.

@adele-morrison
Copy link
Owner

Yep, I think we may need to look at DSW transport across the 1000m isobath instead! I guess the transport change just isn't making it that far north in this time. Is this 55S?

@AndyHoggANU
Copy link
Collaborator

OK, or maybe 2000m isobath?

@matthew-england-unsw
Copy link
Collaborator

(Thanks Julia!!). Ummm.... the values look a bit strange? e.g. there's a strong semi-annual cycle from what I can tell? I count 8 cycles per 4 yrs. And that range from max to min within a year is min O(10 Sv) shooting up to O(30-40 Sv) max values? But not from summer to winter: it's 2 peaks per annum..... ?

@wghuneke
Copy link
Collaborator Author

I started working on the bottom age plot (paper Fig 3) and created a list with items we need to agree upon:

  • Which projection do we want to use? Stereo is nice, but I find it hard to see details on the continental slope. (compare my two version of the age plot)
  • 1 km isobath: We should use the same definition
  • Naming of experiments
  • Which years do we want to show (if it's not a time series)? Average over the last year or last 5 years? (I found it weird to average the age over 5 years, but maybe that's ok.)

@AndyHoggANU
Copy link
Collaborator

Yeah, I was thinking that they would just be too similar, but perhaps we can do it this way (I think it is self-explanatory, but I can explain if you want).
Figure5-2

@matthew-england-unsw
Copy link
Collaborator

Thanks Andy. I think if it’s that similar it can also just be noted in the caption / ms text.

@adele-morrison
Copy link
Owner

adele-morrison commented Nov 7, 2021 via email

@AndyHoggANU
Copy link
Collaborator

Yep, that works. I will calculate the percentage.

@julia-neme
Copy link
Collaborator

First take at figure 8:

figure-8-v1

It becomes really clear how the increase in DSW comes from the meridional component of the wind!!

@adele-morrison
Copy link
Owner

adele-morrison commented Nov 7, 2021 via email

@julia-neme
Copy link
Collaborator

That one has got 6 month rolling mean, I could do either 12 or annual averages??

@adele-morrison
Copy link
Owner

Yeah, maybe give them a try, see what it looks like? Alternatively, we could show anomalies for the DSW export (I know we had a discussion on this on Thursday, can't remember if we decided that was a good or bad idea? @AndyHoggANU?), then the seasonality could be removed by subtracting the control climatology.

@julia-neme
Copy link
Collaborator

I think we had talk about showing the control as well? Anyways, here are the two options, 12 month running mean and anomalies from control:

figure-8-v2
figure-8-v4

@adele-morrison
Copy link
Owner

I like the 12 month running mean.

@matthew-england-unsw
Copy link
Collaborator

Same! :-) PS: interesting to see that the seasonal cycle in DSW anomalies in the WIND+zonal and WIND+meridional cases are ~6 months out of phase, and then approximately cancel out in the total WIND+ case. Could indicate places where zonal winds are cross-shelf I guess. But it's kinda curious that this would be seasonally opposite to the meridional case... ( a side issue of course)

@AndyHoggANU
Copy link
Collaborator

Yeah, I see the point, but maybe it is interesting to include the seasonal variation of DSW export -- partly because the meridional and zonal are out of sync!! ANyone else find that interesting?

@matthew-england-unsw
Copy link
Collaborator

Haha yes I also saw that and thought it was interesting, see my comment above. :-)

@AndyHoggANU
Copy link
Collaborator

OK, if you look at the first draft of figure 8 (not plotted as an anomaly) you can see that what is happening here is that the control has a seasonal peak in summer; this peak is delayed in the zonal wind case and onset earlier in the meridional wind case. So, it might be interesting, but maybe it will be hard to explain succinctly ... so I am now leaning towards saying we should just focus on the 12-month running mean.

@PaulSpence
Copy link
Collaborator

PaulSpence commented Nov 9, 2021

Hi Folks, attempt at Fig 6 - calculation of Ekman pumping binned into surface density bins.

Screen Shot 2021-11-10 at 9 48 32 am

I separated the area integration into on shelf, and within perturbation zone off shelf to reveal upwelling in CDW off shelf. The WIND+x and WIND+y are the zonal and meridional wind components of WIND+. Clearly WIND+x dominates the upwelling/downwelling. It is difficult to isolate the upwelling of CDW using surface rho classes since the surface rho changes a lot seasonally and in spatially, as shown below.

@adele-morrison
Copy link
Owner

Pretty plots Paul! Is something wrong with the WIND+y lines? I'm thinking it should at least be similar to the control (i.e. there should still be upwelling over the open ocean, just no more upwelling perhaps than the control).

Maybe a few north-south transects of temperature with these rho0 isopycnals overlaid would help in identifying on what isopycnal the warm CDW is located?

@PaulSpence
Copy link
Collaborator

Re: Is something wrong with the WIND+y lines? I'm thinking it should at least be similar to the control (i.e. there should still be upwelling over the open ocean, just no more upwelling perhaps than the control).

WIND+y is just the dx_tauyf component of WIND+, not the meridional wind perturbation expt. So WIND+x + WIND+y =WIND+

@adele-morrison
Copy link
Owner

Got it! Sorry!

@wghuneke
Copy link
Collaborator Author

New go on Figure 3 with bottom age and salinity anomalies for WIND+ and WIND-. Figure shows the symmetry between the up and down case and introduces the changes in DSW formation.

I plot the annual average of year 10 (2159), could change the salinity to an average over the last 5 years. (?)
Fig3_Bottom_Age_Salt_Anomaly_cartesian

@wghuneke
Copy link
Collaborator Author

Question on Figure 9 on temperature changes: What depth levels and which experiments do we want to show? The temperature response is not really part of the mechanisms we discuss in the paper. The motivation to include a temperature plot was more to inform the community about what all of that means for the temperature and melting...

Here is a go of bottom temperature for the WIND+ and WIND- cases only. Any suggestions?
Fig9_Temp_anomaly_cartesian

@adele-morrison
Copy link
Owner

Looks good. Yeah, I reckon salinity averaged over the last 5 years is good. Perhaps could try enlarging both figures too (e.g. stretch vertical axis of Fig 3, and for Fig 9 put one underneath the other, so they're bigger?).

@PaulSpence
Copy link
Collaborator

PaulSpence commented Nov 11, 2021

Some updated Ekman/rho plots for the Amundsen, Denman and Wedell region are here:
#20 (comment)
Perhaps rho>27.7 for CDW classification?

@matthew-england-unsw
Copy link
Collaborator

Question on Figure 9 on temperature changes: What depth levels and which experiments do we want to show? The temperature response is not really part of the mechanisms we discuss in the paper. The motivation to include a temperature plot was more to inform the community about what all of that means for the temperature and melting...

Here is a go of bottom temperature for the WIND+ and WIND- cases only. Any suggestions? Fig9_Temp_anomaly_cartesian

Very nice Wilma! I think surface and bottom temperature responses would be good to plot. Like you said, not because it's core to the mechanism, more for community interest. Not necessarily on same plot either; maybe we want T and S changes for both WIND+ and WIND- on same plot. So one Fig. for surface, one for bottom.

@adele-morrison
Copy link
Owner

Haha Matt, I was literally writing this at the same time. :)

For the bottom temperature plot, how well can we explain these changes?

The warming regions I think can be easily linked to the changing vertical heat fluxes associated with the change in the lower overturning cell. i.e. Increased winds -> enhancement of lower cell -> warming in upper ocean near and downstream of DSW production regions.

For the cooling (i.e West Antarctica and much of eastern East Antarctica), I'm less sure of the mechanism. I think it may just be that the isopycnals are being pushed down against the coast. It would be interesting to see if this cooling is coming from the zonal, not meridional component of the change. If it is a heaving signal, do we want to show that by plotting a transect somewhere, or calculating the heave component of the temperature change?

@wghuneke
Copy link
Collaborator Author

wghuneke commented Nov 12, 2021

I prepared the temperature anomaly (at the surface and bottom) figure for WIND+ and WIND-:

  • I'm not sure if we need both experiments as there are only very few places where the response is not symmetric (and we show the symmetry throughout the paper, e.g. Fig 3)
  • Maybe we can replace the WIND- panels with some figure that explains the cooling response (see Adele's comment above). A transect and/or a map of the depth change of some representative isotherm?
  • The bottom salinity is at the moment already shown in Fig 3 which is why I didn't prepared the same figure for salinity. And the surface salinity would be a mirror of the sea ice plot, so maybe not that informative?

Fig9_Temp_surface_bottom_anomaly

@AndyHoggANU
Copy link
Collaborator

Here is a revised Fig 8 that I have sent to Julia to review. Others might have comments. If yes, please post them here. Or comment in the pull request. Whatever ...
Unknown-4

@julia-neme
Copy link
Collaborator

Looks great Andy, thank you! I think we had decided to show 12 month rolling mean though?

@AndyHoggANU
Copy link
Collaborator

Oh yeah, you’re right. Well, the code is there now — I just didn’t have time to extend to the other 3 potential versions. Did you want to do that, or should I?

@PaulSpence
Copy link
Collaborator

PaulSpence commented Nov 22, 2021

Figure of sea ice velocity and thickness anomalies:
Screen Shot 2021-11-22 at 6 02 22 pm

Figure of sea ice transport (volume*velocity) and thickness anomalies:
Screen Shot 2021-11-22 at 6 12 15 pm

I think the velocity offers more clarity than the transport.

@adele-morrison
Copy link
Owner

adele-morrison commented Nov 22, 2021 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
Development

No branches or pull requests

7 participants