\n",
"
"
],
"text/plain": [
"<xarray.DataArray 'explained_variance_ratio' (mode: 4)> Size: 32B\n",
"array([0.82575411, 0.03647744, 0.02654952, 0.01653281])\n",
"Coordinates:\n",
" * mode (mode) int64 32B 1 2 3 4\n",
"Attributes: (12/16)\n",
" model: EOF analysis\n",
" software: xeofs\n",
" version: 2.3.2\n",
" date: 2024-06-07 15:13:21\n",
" n_modes: 4\n",
" center: True\n",
" ... ...\n",
" feature_name: feature\n",
" random_state: None\n",
" verbose: False\n",
" compute: True\n",
" solver: auto\n",
" solver_kwargs: {}\n",
"
"
],
"text/plain": [
"<xarray.DataArray 'explained_variance_ratio' (mode: 4)> Size: 32B\n",
"array([0.12226184, 0.04733037, 0.034539 , 0.03201537])\n",
"Coordinates:\n",
" * mode (mode) int64 32B 1 2 3 4\n",
"Attributes: (12/16)\n",
" model: EOF analysis\n",
" software: xeofs\n",
" version: 2.3.2\n",
" date: 2024-06-07 15:14:36\n",
" n_modes: 4\n",
" center: True\n",
" ... ...\n",
" feature_name: feature\n",
" random_state: None\n",
" verbose: False\n",
" compute: True\n",
" solver: auto\n",
" solver_kwargs: {}\n",
"
"
],
"text/plain": [
"<xarray.DataArray 'explained_variance_ratio' (mode: 4)> Size: 32B\n",
"array([0.32955431, 0.08712222, 0.05315635, 0.04717559])\n",
"Coordinates:\n",
" * mode (mode) int64 32B 1 2 3 4\n",
"Attributes: (12/16)\n",
" model: EOF analysis\n",
" software: xeofs\n",
" version: 2.3.2\n",
" date: 2024-06-07 15:14:39\n",
" n_modes: 4\n",
" center: True\n",
" ... ...\n",
" feature_name: feature\n",
" random_state: None\n",
" verbose: False\n",
" compute: True\n",
" solver: auto\n",
" solver_kwargs: {}![\"ONI\"](\"images/oni.png\")
"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Pacific Decadal Oscillation (PDO)\n",
"\n",
"Here we restrict our domain to the North Pacific. You can [read more about PDO here](https://www.ncei.noaa.gov/access/monitoring/pdo/)."
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {},
"outputs": [],
"source": [
"np_ssta_dt = ssta_dt.where((ssta_dt.lat < 70) & (ssta_dt.lat > 20) & (ssta_dt.lon > 120) & (ssta_dt.lon < 260), drop=True)"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [],
"source": [
"np_model = xe.models.EOF(n_modes=4, use_coslat=True)\n",
"np_model.fit(np_ssta_dt, dim='time')\n",
"np_eofs = np_model.components()\n",
"np_pcs = np_model.scores()\n",
"np_expvar = np_model.explained_variance_ratio()"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": [
"\n",
"
"
],
"text/plain": [
"<xarray.DataArray 'explained_variance_ratio' (mode: 4)> Size: 32B\n",
"array([0.17464065, 0.11074887, 0.072165 , 0.06702556])\n",
"Coordinates:\n",
" * mode (mode) int64 32B 1 2 3 4\n",
"Attributes: (12/16)\n",
" model: EOF analysis\n",
" software: xeofs\n",
" version: 2.3.2\n",
" date: 2024-06-07 15:14:57\n",
" n_modes: 4\n",
" center: True\n",
" ... ...\n",
" feature_name: feature\n",
" random_state: None\n",
" verbose: False\n",
" compute: True\n",
" solver: auto\n",
" solver_kwargs: {}