.. This is a comment in RestructuredText format (two periods and a space).
.. Note that all "statements" and "paragraphs" need to be separated by a blank
line. This means the source code can be hard-wrapped to 80 columns for ease
of reading. Multi-line comments or commands like this need to be indented by
exactly three spaces.
.. Underline with '='s to set top-level heading:
https://docutils.sourceforge.io/docs/user/rst/quickref.html#section-structure
Tropical Pacific Sea Level Diagnostic Documentation
===================================================
Last update: 11/16/2020
Sea level rise is closely related to climate variability and change. It has
important socio-economic impacts on many coastal cities and island nations.
The largest sea level variability occurs in the tropical Pacific, with a
magnitude of ~200mm on interannual timescales during El Niño. This sea level
variability is superimposed on long-term sea level trends and closely related
to global temperature evolution and other climate indices. A detailed
understanding of climate model’s ability in simulating sea level variability
and change in the tropical Pacific is crucial for future projections of
climate and sea level.
.. Underline with '-'s to make a second-level heading.
Version & Contact info
----------------------
.. '-' starts items in a bulleted list:
https://docutils.sourceforge.io/docs/user/rst/quickref.html#bullet-lists
- Version/revision information: version 1 (11/16/2020)
- PI (Jianjun Yin, University of Arizona, yin@arizona.edu)
- Developer/point of contact (Chia-Wei Hsu, University of Arizona, chiaweih@arizona.edu)
.. Underline with '^'s to make a third-level heading.
Open source copyright agreement
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
The MDTF framework is distributed under the LGPLv3 license (see LICENSE.txt).
Unless you've distirbuted your script elsewhere, you don't need to change this.
Functionality
-------------
The main script generates the tropical Pacific dynamic sea level
and wind stress curl scatter plots at different time scales
due to their strong dependency from Ekman pumping/suction
and barotropic response over the ocean.
Python function used
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
- spherical_area.cal_area : generate area array based on the lon lat of data
- dynamical_balance2.curl_var_3d : calculate wind stress curl in obs (for Dataset with time dim)
- dynamical_balance2.curl_var : calculate wind stress curl in obs (for Dataset without time dim)
- dynamical_balance2.curl_tau_3d : calculate wind stress curl in model (for Dataset with time dim)
- dynamical_balance2.curl_tau : calculate wind stress curl in model (for Dataset without time dim)
- xr_ufunc.da_linregress : linregress for Dataset with time dim
Required programming language and libraries
-------------------------------------------
The programming language is python version 3 or up. The third-party libraries
include "matplotlib", "xarray", "cartopy","cftime","numpy". The conda environment
can be set to _MDTF_python3_base.
Required model output variables
-------------------------------
With monthly frequency from the model output. This diagnostic needs
input model variables
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
- tauuo (surface wind stress in the x direction from native model output resolution/grid)
- tauvo (surface wind stress in the y direction from native model output resolution/grid)
- zos (dynamic sea level height in the model from native model output resolution/grid)
The script is written based on the CESM2-OMIP1 download provided by CMIP6-OMIP
hosted by WCRP.
The dimension of all variable is 3-D with (time,nlat,nlon) in dimension and 2-D
array for lat and lon as coordinate.
Required observational data
-------------------------------
This diagnostic needs
input observational variables
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
- adt (absolute dynamic topography from CMEMS)
preprocessing from daily to monthly mean is needed (use 'io_cmems_adt.py')
- tx (surface wind stress in the x direction from WASwind)
no preprocessing needed
- ty (surface wind stress in the y direction from WASwind)
no preprocessing needed
data access :
**********************
- adt :
Ftp server is the fastest way to manage download
`http://marine.copernicus.eu/services-portfolio/access-to-products/ `__
search for product ID - "SEALEVEL_GLO_PHY_L4_REP_OBSERVATIONS_008_047"
Need to download the daily data with adt (absolute dynamic topography) available
- tx,ty :
`https://www.riam.kyushu-u.ac.jp/oed/tokinaga/waswind.html `__
The dimension of all variable is 3-D with 2-D in space and time
References
----------
.. Note this syntax, which sets the "anchor" for the hyperlink: two periods, one
space, one underscore, the reference tag, and a colon, then a blank line.
.. _ref-Hsu:
1. C.-W. Hsu et al. (2020): A Mechanistic Analysis of Tropical Pacific Dynamic
Sea Level in GFDL-OM4 under OMIP-I and OMIP-II Forcings. *GMD*, under review.
2. S. M. Griffies et al. (2016): OMIP contribution to CMIP6: experimental and
diagnostic protocol for the physical component of the Ocean Model Intercomparison
Project. *GMD*, `https://doi.org/10.5194/gmd-9-3231-2016 `__
3. S. Kobayashi et al., (2015): The JRA-55 Reanalysis: General Specifications and Basic Characteristics.
*Journal of the Meteorological Society of Japan. Ser. II*,
`https://doi.org/10.2151/jmsj.2015-001 `__
4. W. G. Large and S. G. Yeager, (2009): The global climatology of an interannually varying air–sea flux data set.
*Climate Dynamics*, `https://doi.org/10.1007/s00382-008-0441-3 `__
More about this diagnostic
--------------------------
The sea level over the tropical Pacific is a key indicator reflecting vertically
integrated heat distribution over the ocean. We find persisting mean state dynamic
sea level (DSL) bias along 9◦N even with updated wind forcing in JRA55-do relative to CORE.
The mean state bias is related to biases in wind stress forcing and geostrophic currents
in the 4◦N to 9◦N latitudinal band. The simulation forced by JRA55-do significantly reduces
the bias in DSL trend over the northern tropical Pacific relative to CORE. In the CORE forcing,
the anomalous westerly wind trend in the eastern tropical Pacific causes an underestimated
DSL trend across the entire Pacific basin along 10◦N. The simulation forced by JRA55-do
significantly reduces the bias in DSL trend over the northern tropical Pacific relative to CORE.
We also identify a bias in the 10 easterly wind trend along 20◦N in both JRA55-do and CORE,
thus motivating future improvement. In JRA55-do, an accurate Rossby wave initiated in the eastern
tropical Pacific at seasonal time scale corrects a biased seasonal variability of the northern
equatorial counter-current in the CORE simulation. Both CORE and JRA55-do generate realistic
DSL variation during El Niño. We find an asymmetry in the DSL pattern on two sides of the equator
is strongly related to wind stress curl that follows the sea level pressure evolution during El Niño.