<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Martha<div class=""><br class=""></div><div class="">Never tried a density equivalent of a temperature threshold — nice idea!</div><div class=""><br class=""></div><div class="">Thanks!</div><div class="">Matt</div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On May 2, 2017, at 8:11 AM, Martha Buckley <<a href="mailto:marthabuckley@gmail.com" class="">marthabuckley@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><div class="gmail_extra">Thanks Matt! ECCO v4 does output DRHODR. I'd be curious about getting your code at some point, but for my present analysis, which involves comparing to other products (without DRHODR available) and observations, I can't really take that route. Unfortunately, I think my only choices are (1) use whatever arbitrary MLD definition is output by the various ocean data products I am using or (2) recompute MLD from some standard definition for all products using monthly mean T/S, neither of which is preferable. </div><div class="gmail_extra">With that said, the focus of my work is not on mixed layer depths, it is on predictability of upper ocean heat content, so the accuracy of the MLD definition may not be key. </div><div class="gmail_extra"><br class=""></div><div class="gmail_extra">Regarding the different criteria in different regions, I thought that using a density equivalent of a temperature threshold took care of this to some degree since the thermal expansion coefficient is smaller at low temperatures. </div><div class="gmail_extra"><br class=""></div><div class="gmail_extra">Thanks!</div><div class="gmail_extra">Martha</div><div class="gmail_extra"><br class=""></div><div class="gmail_extra"><br class=""></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Mon, May 1, 2017 at 11:41 PM, Matthew Mazloff<span class="Apple-converted-space"> </span><span dir="ltr" class=""><<a href="mailto:mmazloff@ucsd.edu" target="_blank" class="">mmazloff@ucsd.edu</a>></span><span class="Apple-converted-space"> </span>wrote:<br class=""><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-color: rgb(204, 204, 204); border-left-style: solid; padding-left: 1ex;"><div style="word-wrap: break-word;" class="">Hello<div class=""><br class=""></div><div class="">To get MLD from a model solution one should really look at first and second vertical derivatives of density. Does ECCO v4 output DRHODR diganostic? If so I can share code to get MLD.</div><div class=""><br class=""></div><div class="">Observations are too noisy to look at vertical derivatives of density so people typically use thresholds. But an appropriate threshold is regionally variable so there are somewhat complicated algorithms for the global problem (e.g.<span class="Apple-converted-space"> </span><a href="http://mixedlayer.ucsd.edu/" target="_blank" class="">http://mixedlayer.ucsd.edu</a>). So to compare to data it is appropriate to use a threshold criteria.</div><div class=""><br class=""></div><div class="">I have looked at MLD from a density threshold criteria, from the Holte algorithm, and from DRHODR in the MITgcm. Really they aren’t qualitatively different, but there are quantitative differences with the threshold criteria tending to give deeper values by missing some restratification events. I like the DRHODR method, but since it can’t be compared to obs the threshold is nice</div><div class=""><br class=""></div><div class="">I would suggest using a potential density cutoff is fine for most applications, just tune it to your location. E.g. 0.3 may be fine in subtropics but may want 0.2 in higher lats. MLD is very nonlinear, but isn’t as bad as mixing<span class="Apple-converted-space"> </span><i class="">layer</i><span class="Apple-converted-space"> </span>depth, which is the instantaneous depth of mixing and is very non-linear and should be calculated on the fly from vertical diffusivity or KPP_HBL or some other method. MLD is sort of the integrated mixing layer depth….so maybe even taking it from pot dense calculated from monthly mean T and S is OK for an approximate average value, though obviously not preferable!</div><div class=""><br class=""></div><div class="">Hope this quick brain dump helps!</div><div class=""><br class=""></div><div class="">Matt</div><div class=""><div class="h5"><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""><div class=""><blockquote type="cite" class=""><div class="">On May 1, 2017, at 7:44 PM, Fenty, Ian G (329C) <<a href="mailto:Ian.Fenty@jpl.nasa.gov" target="_blank" class="">Ian.Fenty@jpl.nasa.gov</a>> wrote:</div><br class="m_2197159868532103759Apple-interchange-newline"><div class=""><div class="">Martha,<br class=""><br class="">I wouldn't say that there is a single MLD definition for ECCO v4 state estimates. Yes, there is a MLD output field that uses the 0.8 deg C criteria but I believe that MLD has no prognostic effect on the solution. Model output can be used to calculate alternative MLDs. For example, Gael provides three alternative MLD criteria (Kara, Suga, Boyer) to show the wintertime and summertime MLD plots as part of the standard analysis package. <br class=""><br class="">See pp 140-142 of<span class="Apple-converted-space"> </span><br class=""><a href="ftp://snowwhite.jpl.nasa.gov/data4/Release3//doc/standardplots.pdf" target="_blank" class="">ftp://snowwhite.jpl.nasa.gov/<wbr class="">data4/Release3//doc/<wbr class="">standardplots.pdf</a><br class=""><br class="">In the Labrador Sea I used the following criteria to identify the base of mixed layer: the shallowest depth level with a potential density exceeding the profile minimum by 13 g m^−3.<br class=""><br class="">Note that ECCOv4 uses the GGL90 package not KPP so there is no KPPhbl field.<br class=""><br class="">If you want exact monthly-mean MLDs from the state estimate using an alternative criteria (e.g., potential density cutoff) then the forward model would have to be rerun since the MLD diagnosed from monthly-mean T and S are almost certainly not equal to the monthly mean of MLD because of nonlinearities.<br class=""><br class="">-Ian<br class=""><br class=""><br class=""><br class="">On 5/1/17, 7:18 PM, "<a href="mailto:ecco-support-bounces@mit.edu" target="_blank" class="">ecco-support-bounces@mit.edu</a><span class="Apple-converted-space"> </span>on behalf of Menemenlis, Dimitris (329C)" <<a href="mailto:ecco-support-bounces@mit.edu" target="_blank" class="">ecco-support-bounces@mit.edu</a><span class="Apple-converted-space"> </span>on behalf of<span class="Apple-converted-space"> </span><a href="mailto:Dimitris.Menemenlis@jpl.nasa.gov" target="_blank" class="">Dimitris.Menemenlis@jpl.nasa.<wbr class="">gov</a>> wrote:<br class=""><br class=""> Hi Martha, I don't have answer for your question. I return discussion to ECCO support, and also include Daria in cc, in case someone else has advice for you.<br class=""><br class=""> On May 1, 2017, at 6:39 PM, Martha Buckley <<a href="mailto:marthabuckley@gmail.com" target="_blank" class="">marthabuckley@gmail.com</a>> wrote:<br class=""><br class=""><br class=""><br class=""> Thanks for the replies. It seems that all the ECCO products use the 0.8C density equivalent. <br class=""><br class=""><br class=""> The reason I am asking is that I am calculating predictability timescales for temperature integrated over the maximum climatological (e.g., wintertime) MLD for a variety of different ocean data products. I am trying to be as comparable as possible between<br class=""> products in the definition of MLD, but the MLD definitions for the products varies wildly (e.g. 0.4 kg/m^3 for ECDA and 0.8 C for ECCO v4). The question is whether the most comparable approach is to apply the same criteria to all models/data (using for example<br class=""> monthly output) or whether there is a better approach.<span class="Apple-converted-space"> </span><br class=""> Thanks!<br class=""> Martha<br class=""><br class=""><br class=""><br class=""> ---------- Forwarded message ----------<br class=""> From: Wang, Ou (329B) <<a href="mailto:Ou.Wang@jpl.nasa.gov" target="_blank" class="">Ou.Wang@jpl.nasa.gov</a>><br class=""> Date: Mon, May 1, 2017 at 1:43 PM<br class=""> Subject: Re: [ecco-support] ECCO Mixed Layer Depths criteria<br class=""> To: "ECCO support list, wider membership" <<a href="mailto:ecco-support@mit.edu" target="_blank" class="">ecco-support@mit.edu</a>><br class=""><br class=""><br class=""> The same default parameters were used for ECCOv4:<br class=""> hMixCriteria = -.8 _d 0<br class=""> dRhoSmall = 1. _d -6<br class=""> hMixSmooth = 0. _d 0<br class=""><br class=""><br class=""> The criteria of the density equivalent of 0.8C is applied to the model profiles at the model time step.<span class="Apple-converted-space"> </span><br class=""><br class=""><br class=""> Ou Wang<br class=""><br class=""><br class=""> From: <<a href="mailto:ecco-support-bounces@mit.edu" target="_blank" class="">ecco-support-bounces@mit.edu</a>> on behalf of "Menemenlis, Dimitris (329C)" <<a href="mailto:Dimitris.Menemenlis@jpl.nasa.gov" target="_blank" class="">Dimitris.Menemenlis@jpl.nasa.<wbr class="">gov</a>><br class=""> Reply-To: "ECCO support list, wider membership" <<a href="mailto:ecco-support@mit.edu" target="_blank" class="">ecco-support@mit.edu</a>><br class=""> Date: Monday, May 1, 2017 at 10:19 AM<br class=""> To: "ECCO support list, wider membership" <<a href="mailto:ecco-support@mit.edu" target="_blank" class="">ecco-support@mit.edu</a>><br class=""> Subject: Re: [ecco-support] ECCO Mixed Layer Depths criteria<br class=""><br class=""><br class=""><br class=""> Don’t know what criterion was used for ECCOv4<br class=""> For cube92, the default parameters were used:<br class=""><br class=""><br class=""> C hMixCriteria:: criteria for mixed-layer diagnostic<br class=""> C dRhoSmall :: parameter for mixed-layer diagnostic<br class=""> C hMixSmooth :: Smoothing parameter for mixed-layer diag (default=0=no smoothing)<br class=""><br class=""> hMixCriteria = -.8 _d 0<br class=""> dRhoSmall = 1. _d -6<br class=""> hMixSmooth = 0. _d 0<br class=""><br class=""> C-- First method :<br class=""> C where the potential density (ref.lev=surface) is larger than<br class=""> C surface density plus Delta_Rho = hMixCriteria * Alpha(surf)<br class=""> C = density of water which is -hMixCriteria colder than surface water<br class=""> C (see Kara, Rochford, and Hurlburt JGR 2000 for default criterion)<br class=""><br class=""><br class=""><br class=""> See model/src/calc_oce_mxlayer for other options.<br class=""><br class=""><br class=""> Another available definition (when using KPP scheme)<br class=""> is KPPhbl, which is the “mixing” layer depth.<br class=""> The daily maximum of KPPhbl is a very good proxy for<br class=""> mixed layer depth most places, except seasonally-sea-ice-covered and polar-winter regions.<br class=""> (Kara’s criterion is also problematic for polar regions, which tend to be salt-stratified.)<br class=""><br class=""><br class=""><br class=""> On May 1, 2017, at 9:50 AM, Martha Buckley <<a href="mailto:marthabuckley@gmail.com" target="_blank" class="">marthabuckley@gmail.com</a>> wrote:<br class=""><br class=""> One follow up question regarding mixed layer depths in ECCO. <br class=""><br class=""><br class=""> What is the criteria for the variable "MXLDEPTH" in the ecco v4 product. I think I recall being told it is the density equivalent of 0.8C. Two questions:<br class=""> (1) Is this the criteria that was used? It the criteria applied to model profiles at the model timestep?<br class=""> (2) If so, isn't this criteria a bit large? 0.125 kg/m^3 or 0.5C seem common MLD criteria when applying criteria to monthly-averaged data. However, criteria that are smaller (~0.03 kg/m^3 or ~0.2 C) are typically applied when using raw observational<br class=""> profiles. <br class=""><br class=""><br class=""> I know that the MLD criteria is somewhat arbitrary, but I want an idea of the rational for the chosen criteria.<span class="Apple-converted-space"> </span><br class=""> Thanks!<br class=""> Martha<br class=""><br class=""><br class=""> On Fri, Apr 28, 2017 at 6:04 PM, HONG ZHANG<span class="Apple-converted-space"> </span><br class=""> <<a href="mailto:hong.zhang@ucla.edu" target="_blank" class="">hong.zhang@ucla.edu</a>> wrote:<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""> -------- Forwarded Message --------<br class=""> Subject: Re: [ecco-support] ECCO2 Mixed Layer Depths ConcernsDate: Fri, 28 Apr 2017 12:48:59 -0700From: HONG ZHANG<span class="Apple-converted-space"> </span><br class=""> <<a href="mailto:hong.zhang@ucla.edu" target="_blank" class="">hong.zhang@ucla.edu</a>> <mailto:<a href="mailto:hong.zhang@ucla.edu" target="_blank" class="">hong.zhang@ucla.edu</a>><wbr class="">To: Dimitris Menemenlis<span class="Apple-converted-space"> </span><br class=""> <<a href="mailto:dmenemenlis@gmail.com" target="_blank" class="">dmenemenlis@gmail.com</a>> <mailto:<a href="mailto:dmenemenlis@gmail.com" target="_blank" class="">dmenemenlis@gmail.com</a>><wbr class="">, Moisan, John R. (WFF-610W)<span class="Apple-converted-space"> </span><br class=""> <<a href="mailto:john.r.moisan@nasa.gov" target="_blank" class="">john.r.moisan@nasa.gov</a>> <mailto:<a href="mailto:john.r.moisan@nasa.gov" target="_blank" class="">john.r.moisan@nasa.gov</a><wbr class="">>CC: Daria Halkides<span class="Apple-converted-space"> </span><br class=""> <<a href="mailto:halkides@esr.org" target="_blank" class="">halkides@esr.org</a>> <mailto:<a href="mailto:halkides@esr.org" target="_blank" class="">halkides@esr.org</a>><br class=""><br class=""> On 4/28/17 9:37 AM, Dimitris Menemenlis wrote:<br class=""><br class=""><br class=""> Daria or Hong, would you have looked at MXLDEPTH depth in the “cube92” solution?<br class=""> <a href="ftp://ecco2.jpl.nasa.gov/data1/cube/cube92/lat_lon/quart_90S_90N/MXLDEPTH.nc/" target="_blank" class="">ftp://ecco2.jpl.nasa.gov/<wbr class="">data1/cube/cube92/lat_lon/<wbr class="">quart_90S_90N/MXLDEPTH.nc/</a><br class=""><br class=""><br class=""> John, why not use the ECCO-V4, which is data-constrained?<br class=""> <a href="http://www.ecco-group.org/products.htm" target="_blank" class="">http://www.ecco-group.org/<wbr class="">products.htm</a><br class=""> ECCO v4 is on a 1-degree grid, which is closer to your desired 2x2 degree grid.<br class=""><br class=""> Dimitris Menemenlis<span class="Apple-converted-space"> </span><br class=""><br class=""><br class=""> On Apr 26, 2017, at 9:41 AM, Moisan, John R. (WFF-610W) <<a href="mailto:john.r.moisan@nasa.gov" target="_blank" class="">john.r.moisan@nasa.gov</a>> wrote:<br class=""><br class=""> Hello ECCO2’ers,<br class=""><br class=""> I would like to make use of the Mixed-Layer Depth model outputs from the global model runs 1992-present (Nov. 2016) to drive some simple 1D biogeochemical models for CDOM. I have binned to model outputs from the jpg site into monthly 2x2 deg. grids. One thing<br class=""> of concern that I have is the trends in the global MLDs (see attached figure), which are more extreme in the SO. Is there an issue with maintaining correct ocean boundary layers, esp. in the SO for this model? Has this been observed yet?<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""> Hi John,<br class=""> You're right, MLD from cube92 run has bias in late stage<br class=""> because it's a free run for the later period.<br class=""> As Dimitris said, you can look at other ECCO product provided at the website.<br class=""> Please see attached fig showing MLD from ECCOv4 and ECCO llc270 for an example,<br class=""> where red line is for ECCOv4 (1992-2012) and blue line for ECCO llc270 (2001-2015) in upper plot<br class=""> also MLD for cube92 is also shown in lower plot for reference.<br class=""><br class=""> cheers<br class=""> Hong<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""> ___________________________<wbr class="">____________________<br class=""> ecco-support mailing list<br class=""> <a href="mailto:ecco-support@mit.edu" target="_blank" class="">ecco-support@mit.edu</a><br class=""> <a href="http://mailman.mit.edu/mailman/listinfo/ecco-support" target="_blank" class="">http://mailman.mit.edu/<wbr class="">mailman/listinfo/ecco-support</a><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""> --<span class="Apple-converted-space"> </span><br class=""> Martha W. Buckley<br class=""> <a href="mailto:marthabuckley@gmail.com" target="_blank" class="">marthabuckley@gmail.com</a><br class=""> <a href="http://sites.google.com/site/marthabuckley/home" target="_blank" class="">http://sites.google.com/<wbr class="">site/marthabuckley/home</a><br class=""><br class=""><br class=""><br class=""><br class=""> ___________________________<wbr class="">____________________<br class=""> ecco-support mailing list<br class=""> <a href="mailto:ecco-support@mit.edu" target="_blank" class="">ecco-support@mit.edu</a><br class=""> <a href="http://mailman.mit.edu/mailman/listinfo/ecco-support" target="_blank" class="">http://mailman.mit.edu/<wbr class="">mailman/listinfo/ecco-support</a><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""> ___________________________<wbr class="">____________________<br class=""> ecco-support mailing list<br class=""> <a href="mailto:ecco-support@mit.edu" target="_blank" class="">ecco-support@mit.edu</a><br class=""> <a href="http://mailman.mit.edu/mailman/listinfo/ecco-support" target="_blank" class="">http://mailman.mit.edu/<wbr class="">mailman/listinfo/ecco-support</a><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""> --<span class="Apple-converted-space"> </span><br class=""> Martha W. Buckley<br class=""> <a href="mailto:marthabuckley@gmail.com" target="_blank" class="">marthabuckley@gmail.com</a><br class=""> <a href="http://sites.google.com/site/marthabuckley/home" target="_blank" class="">http://sites.google.com/<wbr class="">site/marthabuckley/home</a><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class="">______________________________<wbr class="">_________________<br class="">ecco-support mailing list<br class=""><a href="mailto:ecco-support@mit.edu" target="_blank" class="">ecco-support@mit.edu</a><br class=""><a href="http://mailman.mit.edu/mailman/listinfo/ecco-support" target="_blank" class="">http://mailman.mit.edu/<wbr class="">mailman/listinfo/ecco-support</a><br class=""></div></div></blockquote></div><br class=""></div></div></div></div></blockquote></div><br class=""><br clear="all" class=""><div class=""><br class=""></div>--<span class="Apple-converted-space"> </span><br class=""><div class="gmail_signature" data-smartmail="gmail_signature">Martha W. Buckley<br class=""><a href="mailto:marthabuckley@gmail.com" target="_blank" class="">marthabuckley@gmail.com</a><br class=""><a href="http://sites.google.com/site/marthabuckley/home" target="_blank" class="">http://sites.google.com/site/marthabuckley/home</a></div></div></div></div></blockquote></div><br class=""></div></body></html>