Publications

Export 25 results:
Sort by: Author Title Type [ Year  (Desc)]
2018
Amaya, DJ, Siler N, Xie SP, Miller AJ.  2018.  The interplay of internal and forced modes of Hadley Cell expansion: lessons from the global warming hiatus. Climate Dynamics. 51:305-319.   10.1007/s00382-017-3921-5   AbstractWebsite

The poleward branches of the Hadley Cells and the edge of the tropics show a robust poleward shift during the satellite era, leading to concerns over the possible encroachment of the globe's subtropical dry zones into currently temperate climates. The extent to which this trend is caused by anthropogenic forcing versus internal variability remains the subject of considerable debate. In this study, we use a Joint EOF method to identify two distinct modes of tropical width variability: (1) an anthropogenically-forced mode, which we identify using a 20-member simulation of the historical climate, and (2) an internal mode, which we identify using a 1000-year pre-industrial control simulation. The forced mode is found to be closely related to the top of the atmosphere radiative imbalance and exhibits a long-term trend since 1860, while the internal mode is essentially indistinguishable from the El Nio Southern Oscillation. Together these two modes explain an average of 70% of the interannual variability seen in model "edge indices" over the historical period. Since 1980, the superposition of forced and internal modes has resulted in a period of accelerated Hadley Cell expansion and decelerated global warming (i.e., the "hiatus"). A comparison of the change in these modes since 1980 indicates that by 2013 the signal has emerged above the noise of internal variability in the Southern Hemisphere, but not in the Northern Hemisphere, with the latter also exhibiting strong zonal asymmetry, particularly in the North Atlantic. Our results highlight the important interplay of internal and forced modes of tropical width change and improve our understanding of the interannual variability and long-term trend seen in observations.

2017
Miller, AJ, Collins M, Gualdi S, Jensen TG, Misra V, Pezzi LP, Pierce DW, Putrasahan D, Seo H, Tseng YH.  2017.  Coupled ocean-atmosphere modeling and predictions. Journal of Marine Research. 75:361-402. AbstractWebsite

Key aspects of the current state of the ability of global and regional climate models to represent dynamical processes and precipitation variations are summarized. Interannual, decadal, and global-warming timescales, wherein the influence of the oceans is relevant and the potential for predictability is highest, are emphasized. Oceanic influences on climate occur throughout the ocean and extend over land to affect many types of climate variations, including monsoons, the El Nino Southern Oscillation, decadal oscillations, and the response to greenhouse gas emissions. The fundamental ideas of coupling between the ocean-atmosphere-land system are explained for these modes in both global and regional contexts. Global coupled climate models are needed to represent and understand the complicated processes involved and allow us to make predictions over land and sea. Regional coupled climate models are needed to enhance our interpretation of the fine-scale response. The mechanisms by which large-scale, low-frequency variations can influence shorter timescale variations and drive regional-scale effects are also discussed. In this light of these processes, the prospects for practical climate predictability are also presented.

Gan, BL, Wu LX, Jia F, Li SJ, Cai WJ, Nakamura H, Alexander MA, Miller AJ.  2017.  On the response of the Aleutian Low to greenhouse warming. Journal of Climate. 30:3907-3925.   10.1175/jcli-d-15-0789.1   AbstractWebsite

Past and future changes in the Aleutian low are investigated by using observation-based sea level pressure (SLP) datasets and CMIP5 models. It is found that the Aleutian low intensity, measured by the North Pacific Index (NPI), has significantly strengthened during the twentieth century, with the observed centennial trend double the modeled counterpart for the multimodel average of historical simulations, suggesting compound signals of anthropogenic warming and natural variability. As climate warms under the strongest future warming scenario, the climatological-mean Aleutian low will continue to intensify and expand northward, as manifested in the significant decrease (-1.3 hPa) of the multimodel-averaged NPI, which is 1.6 times its unforced internal variability, and the increase in the central area of low pressure (SLP < 999.0 hPa), which expands about 7 times that in the twentieth century. A suite of idealized experiments further demonstrates that the deepening of the Aleutian low can be driven by an El Nino-like warming of the tropical Pacific sea surface temperature (SST), with a reduction in the climatological-mean zonal SST gradient, which overshadows the dampening effect of a weakened wintertime land-ocean thermal contrast on the Aleutian low change in a warmer climate. While the projected deepening of Aleutian low on multimodel average is robust, individual model portrayals vary primarily in magnitude. Intermodel difference in surface warming amplitude over the Asian continent, which is found to explain about 31% of the variance of the NPI changes across models, has a greater contribution than that in the spatial pattern of tropical Pacific SST warming (which explains about 23%) to model uncertainty in the projection of Aleutian low intensity.

2016
Lou, SJ, Russell LM, Yang Y, Xu L, Lamjiri MA, DeFlorio MJ, Miller AJ, Ghan SJ, Liu Y, Singh B.  2016.  Impacts of the East Asian Monsoon on springtime dust concentrations over China. Journal of Geophysical Research-Atmospheres. 121:8137-8152.   10.1002/2016jd024758   AbstractWebsite

We use 150year preindustrial simulations of the Community Earth System Model to quantify the impacts of the East Asian Monsoon strength on interannual variations of springtime dust concentrations over China. The simulated interannual variations in March-April-May (MAM) dust column concentrations range between 20-40% and 10-60% over eastern and western China, respectively. The dust concentrations over eastern China correlate negatively with the East Asian Monsoon (EAM) index, which represents the strength of monsoon, with a regionally averaged correlation coefficient of -0.64. Relative to the strongest EAM years, MAM dust concentrations in the weakest EAM years are higher over China, with regional relative differences of 55.6%, 29.6%, and 13.9% in the run with emissions calculated interactively and of 33.8%, 10.3%, and 8.2% over eastern, central, and western China, respectively, in the run with prescribed emissions. Both interactive run and prescribed emission run show the similar pattern of climate change between the weakest and strongest EAM years. Strong anomalous northwesterly and westerly winds over the Gobi and Taklamakan deserts during the weakest EAM years result in larger transport fluxes, and thereby increase the dust concentrations over China. These differences in dust concentrations between the weakest and strongest EAM years (weakest-strongest) lead to the change in the net radiative forcing by up to -8 and -3Wm(-2) at the surface, compared to -2.4 and +1.2Wm(-2) at the top of the atmosphere over eastern and western China, respectively.

Newman, M, Alexander MA, Ault TR, Cobb KM, Deser C, Di Lorenzo E, Mantua NJ, Miller AJ, Minobe S, Nakamura H, Schneider N, Vimont DJ, Phillips AS, Scott JD, Smith CA.  2016.  The Pacific Decadal Oscillation, Revisited. Journal of Climate. 29:4399-4427.   10.1175/jcli-d-15-0508.1   AbstractWebsite

The Pacific decadal oscillation (PDO), the dominant year-round pattern of monthly North Pacific sea surface temperature (SST) variability, is an important target of ongoing research within the meteorological and climate dynamics communities and is central to the work of many geologists, ecologists, natural resource managers, and social scientists. Research over the last 15 years has led to an emerging consensus: the PDO is not a single phenomenon, but is instead the result of a combination of different physical processes, including both remote tropical forcing and local North Pacific atmosphere-ocean interactions, which operate on different time scales to drive similar PDO-like SST anomaly patterns. How these processes combine to generate the observed PDO evolution, including apparent regime shifts, is shown using simple autoregressive models of increasing spatial complexity. Simulations of recent climate in coupled GCMs are able to capture many aspects of the PDO, but do so based on a balance of processes often more independent of the tropics than is observed. Finally, it is suggested that the assessment of PDO-related regional climate impacts, reconstruction of PDO-related variability into the past with proxy records, and diagnosis of Pacific variability within coupled GCMs should all account for the effects of these different processes, which only partly represent the direct forcing of the atmosphere by North Pacific Ocean SSTs.

2015
Amaya, DJ, Xie SP, Miller AJ, McPhaden MJ.  2015.  Seasonality of tropical Pacific decadal trends associated with the 21st century global warming hiatus. Journal of Geophysical Research-Oceans. 120:6782-6798.   10.1002/2015jc010906   AbstractWebsite

Equatorial Pacific changes during the transition from a nonhiatus period (pre-1999) to the present global warming hiatus period (post-1999) are identified using a combination of reanalysis and observed data sets. Results show increased surface wind forcing has excited significant changes in wind-driven circulation. Over the last two decades, the core of the Equatorial Undercurrent intensified at a rate of 6.9 cm s(-1) decade(-1). Similarly, equatorial upwelling associated with the shallow meridional overturning circulation increased at a rate of 2.0 x 10(-4) cm s(-1) decade(-1) in the central Pacific. Further, a seasonal dependence is identified in the sea surface temperature trends and in subsurface dynamics. Seasonal variations are evident in reversals of equatorial surface flow trends, changes in subsurface circulation, and seasonal deepening/shoaling of the thermocline. Anomalous westward surface flow drives cold-water zonal advection from November to February, leading to surface cooling from December through May. Conversely, eastward surface current anomalies in June-July drive warm-water zonal advection producing surface warming from July to November. An improved dynamical understanding of how the tropical Pacific Ocean responds during transitions into hiatus events, including its seasonal structure, may help to improve future predictability of decadal climate variations.

2013
Moon, JH, Song YT, Bromirski PD, Miller AJ.  2013.  Multidecadal regional sea level shifts in the Pacific over 1958-2008. Journal of Geophysical Research-Oceans. 118:7024-7035.   10.1002/2013jc009297   AbstractWebsite

Altimeter data have significantly improved our understanding of regional sea level variability and trends, but their relatively short records do not allow either evaluation of the ocean state prior to 1993 or multidecadal low-frequency signals in the ocean. Here we characterize and quantify the multidecadal regional sea level rise (rSLR) and related ocean heat content in the Pacific from a non-Boussinesq ocean circulation model in comparison with data sets from altimeters, two sea level reconstructions, and in situ ocean profiles from 1958 to 2008. We show that the rSLR trends have undergone two shifts, during the mid-1970s and in the early 1990s, with an east-west dipole pattern in the tropical Pacific. In each of these phases, rSLR accelerated on one side of the Pacific, but decelerated on the other side. The multidecadal sea level shifts can be explained by the dynamical (steric) upper-ocean responses to the surface wind forcing associated with the Pacific Decadal Oscillation (PDO), with negligible contributions from internal (depth-integrated) ocean mass changes. Additional model experimentation further confirms that the Pacific wind stress trend over the recent two decades has played an important role in strengthening the rSLR in the western Pacific while suppressing the rSLR in the eastern Pacific. The climate-forced large-scale rSLR variability is likely to impose a long-term and uneven impact on coastal communities.

2011
Bromirski, PD, Miller AJ, Flick RE, Auad G.  2011.  Dynamical suppression of sea level rise along the Pacific coast of North America: Indications for imminent acceleration. Journal of Geophysical Research-Oceans. 116   10.1029/2010jc006759   AbstractWebsite

Long-term changes in global mean sea level (MSL) rise have important practical implications for shoreline and beach erosion, coastal wetlands inundation, storm surge flooding, and coastal development. Altimetry since 1993 indicates that global MSL rise has increased about 50% above the 20th century rise rate, from 2 to 3 mm yr(-1). At the same time, both tide gauge measurements and altimetry indicate virtually no increase along the Pacific coast of North America during the satellite epoch. Here we show that the dynamical steric response of North Pacific eastern boundary ocean circulation to a dramatic change in wind stress curl, tau(xy), which occurred after the mid-1970s regime shift, can account for the suppression of regional sea level rise along this coast since 1980. Alarmingly, mean tau(xy) over the North Pacific recently reached levels not observed since before the mid-1970s regime shift. This change in wind stress patterns may be foreshadowing a Pacific Decadal Oscillation regime shift, causing an associated persistent change in basin-scale tau(xy) that may result in a concomitant resumption of sea level rise along the U.S. West Coast to global or even higher rates.

Yeh, SW, Kang YJ, Noh Y, Miller AJ.  2011.  The North Pacific climate transitions of the winters of 1976/77 and 1988/89. Journal of Climate. 24:1170-1183.   10.1175/2010jcli3325.1   AbstractWebsite

This paper examines characteristic changes in North Pacific sea surface temperature (SST) variability during the boreal winter (December-February) for two subperiods (1956-88 and 1977-2009) during which the 1976/77 and the 1988/89 climate transitions occurred. It is found that the Pacific decadal oscillation (PDO)-like SST variability plays a dominant role in the 1976/77 climate transition, while both the North Pacific Gyre Oscillation (NPGO)-like and PDO-like SST variability contribute to the 1988/89 climate transition. Furthermore, the leading mode changes from PDO-like SST variability during the period 1956-88 to NPGO-like SST variability during the period 1977-2009, indicative of an enhancement of NPGO-like SST variability since 1988. Changes in sea level pressure across the 1976/77 climate transition project strongly onto the Aleutian low pressure system. But sea level pressure changes across the 1988/89 climate transition project primarily onto the North Pacific Oscillation, which is associated with remote changes in the Arctic Oscillation over the polar region as well. This contributes to enhancing the NPGO-like SST variability after 1988. The authors also analyze the output from an ensemble of Tropical Ocean and Global Atmosphere (TOGA) experiments in which the observed SSTs are inserted only at grid points in the tropics between 20 degrees S and 20 degrees N. The results indicate that the changes in the North Pacific atmosphere in the 1976/77 climate transition are mostly due to the tropics, whereas those in the 1988/89 climate transition are not.

2009
Capotondi, A, Combes V, Alexander MA, Di Lorenzo E, Miller AJ.  2009.  Low-frequency variability in the Gulf of Alaska from coarse and eddy-permitting ocean models. Journal of Geophysical Research-Oceans. 114   10.1029/2008jc004983   AbstractWebsite

An eddy-permitting ocean model of the northeast Pacific is used to examine the ocean adjustment to changing wind forcing in the Gulf of Alaska (GOA) at interannual-to-decadal timescales. It is found that the adjustment of the ocean model in the presence of mesoscale eddies is similar to that obtained with coarse-resolution models. Local Ekman pumping plays a key role in forcing pycnocline depth variability and, to a lesser degree, sea surface height (SSH) variability in the center of the Alaska gyre and in some areas of the eastern and northern GOA. Westward Rossby wave propagation is evident in the SSH field along some latitudes but is less noticeable in the pycnocline depth field. Differences between SSH and pycnocline depth are also found when considering their relationship with the local forcing and leading modes of climate variability in the northeast Pacific. In the central GOA pycnocline depth variations are more clearly related to changes in the local Ekman pumping than SSH. While SSH is marginally correlated with both Pacific Decadal Oscillation (PDO) and North Pacific Gyre Oscillation (NPGO) indices, the pycnocline depth evolution is primarily related to NPGO variability. The intensity of the mesoscale eddy field increases with increasing circulation strength. The eddy field is generally more energetic after the 1976-1977 climate regime shift, when the gyre circulation intensified. In the western basin, where eddies primarily originate from intrinsic instabilities of the flow, variations in eddy kinetic energy are statistically significant correlated with the PDO index, indicating that eddy statistics may be inferred, to some degree, from the characteristics of the large-scale flow.

2008
Haidvogel, DB, Arango H, Budgell WP, Cornuelle BD, Curchitser E, Di Lorenzo E, Fennel K, Geyer WR, Hermann AJ, Lanerolle L, Levin J, McWilliams JC, Miller AJ, Moore AM, Powell TM, Shchepetkin AF, Sherwood CR, Signell RP, Warner JC, Wilkin J.  2008.  Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System. Journal of Computational Physics. 227:3595-3624.   10.1016/j.jcp.2007.06.016   AbstractWebsite

Systematic improvements in algorithmic design of regional ocean circulation models have led to significant enhancement in simulation ability across a wide range of space/time scales and marine system types. As an example, we briefly review the Regional Ocean Modeling System, a member of a general class of three-dimensional, free-surface, terrain-following numerical models. Noteworthy characteristics of the ROMS computational kernel include: consistent temporal averaging of the barotropic mode to guarantee both exact conservation and constancy preservation properties for tracers; redefined barotropic pressure-gradient terms to account for local variations in the density field; vertical interpolation performed using conservative parabolic splines; and higher-order, quasi-monotone advection algorithms. Examples of quantitative skill assessment are shown for a tidally driven estuary, an ice-covered high-latitude sea, a wind- and buoyancy-forced continental shelf, and a mid-latitude ocean basin. The combination of moderate-order spatial approximations, enhanced conservation properties, and quasi-monotone advection produces both more robust and accurate, and less diffusive, solutions than those produced in earlier terrain-following ocean models. Together with advanced methods of data assimilation and novel observing system technologies, these capabilities constitute the necessary ingredients for multi-purpose regional ocean prediction systems. (c) 2007 Elsevier Inc. All rights reserved.

Seo, H, Jochum M, Murtugudde R, Miller AJ, Roads JO.  2008.  Precipitation from African easterly waves in a coupled model of the tropical Atlantic. Journal of Climate. 21:1417-1431.   10.1175/2007jcli1906.1   AbstractWebsite

A regional coupled climate model is configured for the tropical Atlantic to explore the role of synopticscale African easterly waves (AEWs) on the simulation of mean precipitation in the marine intertropical convergence zone (ITCZ). Sensitivity tests with varying atmospheric resolution in the coupled model show that these easterly waves are well represented with comparable amplitudes on both fine and coarse grids of the atmospheric model. Significant differences in the model simulations are found in the precipitation fields, however, where heavy rainfall events occur in the region of strong cyclonic shear of the easterly waves only on the higher-resolution grid. This is because the low-level convergence due to the waves is much larger and more realistic in the fine-resolution simulation, which enables heavier precipitation events that skew the rainfall distributions toward longer tails. The variability in rainfall on these time scales accounts for more than 60%-70% of the total variability. As a result, the simulation of mean rainfall in the ITCZ and its seasonal migration improves in the higher-resolution case. This suggests that capturing these transient waves and the resultant strong low-level convergence is one of the key ingredients for improving the simulation of precipitation in global coupled climate models.

2006
Auad, G, Miller A, Di Lorenzo E.  2006.  Long-term forecast of oceanic conditions off California and their biological implications. Journal of Geophysical Research-Oceans. 111   10.1029/2005jc003219   AbstractWebsite

[1] The impact of global warming due to an increased content of atmospheric CO(2) is studied by forcing a numerical eddy-resolving ocean model with wind stresses, heat fluxes, and open boundary conditions obtained from a state-of-the-art coupled model. Specifically, we have compared the 1986 - 1996 and 2040 - 2050 decades to describe and analyze the changes attained by several oceanographic variables in the California Current System. A richer atmosphere in CO(2) leads to increased sea surface and near-surface temperatures in the model domain and to an increased stratification along the coast that, however, is not strong enough to overcome the effect of increased upwelling favorable winds. A mild oceanic cooling is forecast below the 70-m depth, in agreement with recent studies of global warming trends. Near-surface vertical velocities increase by about 30% in April, but our simulations also forecast anomalous offshore transports in most of the coastal areas. The eddy kinetic energy decreases, on an annual mean, along the California Current main path with maximum negative anomalies in winter. The upward displacement of the 26.5 isopycnal surface, especially in the northern half of our study area, leads to an increase in the concentration of nutrients in the subsurface. The agreement between some results of this forecasting study and recent published findings would suggest that the current global warming trend would persist in the study area with similar changes to those observed over the last half century.

Mestas-Nunez, AM, Miller AJ.  2006.  Interdecadal variability and climate change in the eastern tropical Pacific: A review. Progress in Oceanography. 69:267-284.   10.1016/j.pocean.2006.03.011   AbstractWebsite

In this paper, we review interdecadal climatic variability in the eastern tropical Pacific Ocean. This variability dominates the climatic fluctuations in the North Pacific on scales between ENSO and the centennial trend and is commonly referred to as the Pacific Decadal Oscillation or PDO. We include a historical overview and a summary of observational work that describes the surface, tropospheric and subsurface signatures of this variability. Descriptions of interdecadal variability are incomplete at best, mostly due to limitations in the observational record. We emphasize that the well-known "ENSO-like" sea surface temperature (SST) pattern describing the PDO may not be an accurate representation. In the eastern tropical Pacific, the SST maxima are displaced north and south of the equator with larger amplitudes in the northern branch near the coast of North America, which has significant implications for the troposphere-driven circulations. Several mechanisms have been proposed to explain the PDO. We review these mechanisms and models, which capture our present level of understanding of the problem. We conclude by reporting there is little evidence of both multidecadal variability and the centennial trend in the eastern tropical Pacific. This paper is part of a comprehensive review of the oceanography of the eastern tropical Pacific. (c) 2006 Elsevier Ltd. All rights reserved.

2005
Miller, AJ, Di Lorenzo E, Neilson DJ, Kim HJ, Capotondi A, Alexander MA, Bograd SJ, Schwing FB, Mendelssohn R, Hedstrom K, Musgrave DL.  2005.  Interdecadal changes in mesoscale eddy variance in the Gulf of Alaska circulation: Possible implications for the Steller sea lion ecline. Atmosphere-Ocean. 43:231-240.   10.3137/ao.430303   AbstractWebsite

A distinct change in the ocean circulation of the Gulf of Alaska after the 1976-77 climate shift is studied in an eddy-permitting primitive equation model forced by observed wind stresses from. 1951-99. When the Aleutian Low strengthens after 1976-77, strong changes occur in the mean velocity, of the Alaskan Stream and in its associated mesoscale eddy field. In contrast, the Alaska Current and the eddy flows in the eastern Gulf remain relatively unchanged after the shift. Since mesoscale eddies provide a possible mechanism for transporting nutrient-rich open-ocean waters to the productive shelf region, the flow of energy through the food web may have been altered by this physical oceanographic change. This climate-driven mechanism, which has a characteristic east-west spatial asymmetry, may potentially help to explain changes in forage fish quality, in diet diversity of Steller sea lions whose populations have declined precipitously since the mid-1970s in the western Gulf while remaining stable in the eastern Gulf.

Di Lorenzo, E, Miller AJ, Schneider N, McWilliams JC.  2005.  The warming of the California current system: Dynamics and ecosystem implications. Journal of Physical Oceanography. 35:336-362.   10.1175/jpo-2690.1   AbstractWebsite

Long-term changes in the observed temperature and salinity along the southern California coast are studied using a four-dimensional space-time analysis of the 52-yr (1949-2000) California Cooperative Oceanic Fisheries Investigations (CalCOFI) hydrography combined with a sensitivity analysis of an eddy-permitting primitive equation ocean model under various forcing scenarios. An overall warming trend of 1.3 degrees C in the ocean surface, a deepening in the depth of the mean thermocline (18 m), and increased stratification between 1950 and 1999 are found to be primarily forced by large-scale decadal fluctuations in surface heat fluxes combined with horizontal advection by the mean currents. After 1998 the surface heat fluxes suggest the beginning of a period of cooling, consistent with colder observed ocean temperatures. Salinity changes are decoupled from temperature and appear to be controlled locally in the coastal ocean by horizontal advection by anomalous currents. A cooling trend of -0.5 degrees C in SST is driven in the ocean model by the 50-yr NCEP wind reanalysis, which contains a positive trend in upwelling-favorable winds along the southern California coast. A net warming trend of +1 degrees C in SST occurs, however, when the effects of observed surface heat fluxes are included as forcing functions in the model. Within 50-100 km of the coast, the ocean model simulations show that increased stratification/deepening of the thermocline associated with the warming reduces the efficiency of coastal upwelling in advecting subsurface waters to the ocean surface, counteracting any effects of the increased strength of the upwelling winds. Such a reduction in upwelling efficiency leads in the model to a freshening of surface coastal waters. Because salinity and nutrients at the coast have similar distributions this must reflect a reduction of the nutrient supply at the coast, which is manifestly important in explaining the observed decline in zooplankton concentration. The increased winds also drive an intensification of the mean currents of the southern California Current System (SCCS). Model mesoscale eddy variance significantly increases in recent decades in response to both the stronger upwelling winds and the warmer upper-ocean temperatures, suggesting that the stability properties of the SCCS have also changed.

Capotondi, A, Alexander MA, Deser C, Miller AJ.  2005.  Low-frequency pycnocline variability in the northeast Pacific. Journal of Physical Oceanography. 35:1403-1420.   10.1175/jpo2757.1   AbstractWebsite

The output from an ocean general circulation model (OGCM) driven by observed surface forcing is used in conjunction with simpler dynamical models to examine the physical mechanisms responsible for interannual to interdecadal pycnocline variability in the northeast Pacific Ocean during 1958-97, a period that includes the 1976-77 climate shift. After 1977 the pycnocline deepened in a broad band along the coast and shoaled in the central part of the Gulf of Alaska. The changes in pycnocline depth diagnosed from the model are in agreement with the pycnocline depth changes observed at two ocean stations in different areas of the Gulf of Alaska. A simple Ekman pumping model with linear damping explains a large fraction of pycnocline variability in the OGCM. The fit of the simple model to the OGCM is maximized in the central part of the Gulf of Alaska, where the pycnocline variability produced by the simple model can account for similar to 70%-90% of the pycnocline depth variance in the OGCM. Evidence of westward-propagating Rossby waves is found in the OGCM, but they are not the dominant signal. On the contrary, large-scale pycnocline depth anomalies have primarily a standing character, thus explaining the success of the local Ekman pumping model. The agreement between the Ekman pumping model and OGCM deteriorates in a large band along the coast, where propagating disturbances within the pycnocline, due to either mean flow advection or boundary waves, appear to play an important role in pycnocline variability. Coastal propagation of pycnocline depth anomalies is especially relevant in the western part of the Gulf of Alaska, where local Ekman pumping-induced changes are anticorrelated with the OGCM pycnocline depth variations. The pycnocline depth changes associated with the 1976-77 climate regime shift do not seem to be consistent with Sverdrup dynamics, raising questions about the nature of the adjustment of the Alaska Gyre to low-frequency wind stress variability.

2004
Auad, G, Miller AJ, Roads JO.  2004.  Pacific Ocean forecasts. Journal of Marine Systems. 45:75-90.   10.1016/j.jmarsys.2003.11.010   AbstractWebsite

A primitive equation Pacific Ocean model forced by wind stresses and heat fluxes is used to obtain uncoupled forecasts of sea surface temperature (SST), heat storage (upper 400 m), and surface currents. The forecasts are displayed in real-time on the web http://ecpc.ucsd.edu) and are compared against observations obtained from the Reynolds (SST) and Joint Environmental Data Analysis Center, JEDAC http://jedac.ucsd.edu), (0- to 400-m temperature) data sets. The resulting forecast skill, for both total and anomalous fields, are reasonably good given the simplicity of our methodology and the fact that feedback processes between ocean and atmosphere are absent. SST forecasts are equal and even superior to anomaly persistence forecasts in some regions during some seasons. Given this skill, which depends both on model performance and on quality and sampling density of the observations, we are beginning to develop various applications for these experimental forecasts. (C) 2004 Elsevier B.V. All rights reserved.

2001
Auad, G, Miller AJ, Roads JO, Cayan D.  2001.  Pacific Ocean wind stress and surface heat flux anomalies from NCEP reanalysis and observations: Cross-statistics and ocean model responses. Journal of Geophysical Research-Oceans. 106:22249-22265.   10.1029/2000jc000264   AbstractWebsite

Wind stresses and surface heat fluxes over the Pacific Ocean from the National Center for Environmental Prediction (NCEP) reanalysis and the comprehensive Ocean-Atmosphere Data Set (COADS) (blended with FSU tropical wind stresses) are compared over a common time interval (1958-1997) in their statistics anal in the responses that they induce in sea surface temperature (SST) and heat storage when used to force an ocean model. Wind stress anomalies from the two data sets are well correlated in the midlatitude extratropics, especially in the highly sampled North Pacific. In the tropics and subtropics, low correlations were found between the two wind stress data sets. The amplitudes of the stress variations of the two data sets are similar in midlatitudes, but in the tropics NCEP wind stresses are weaker than the LOADS/FSU stresses, especially on interannual timescales. Surface heat flux anomalies from the two data sets are well correlated on interannual and shorter timescales in the North Pacific Ocean poleward of 20 degreesN, but they are poorly correlated elsewhere and on decadal timescales. In the extratropics the amplitudes of the heat flux variations of the two data sets are comparable, but in the tropics the NCEP heat fluxes are weaker than those of CORDS. Ocean model hindcasts driven by bath data sets are also compared: The midlatitude SST hindcasts were superior when using the NCEP flux anomalies while tropical SST hindcasts were equally skillful for the two hindcasts when considering all climatic timescales. The spatial and temporal sampling rates of the LOADS observations and their consequent impacts on constraining the NCEP reanalysis appear to be the main factors controlling the results found here.

2000
Miller, AJ, Schneider N.  2000.  Interdecadal climate regime dynamics in the North Pacific Ocean: theories, observations and ecosystem impacts. Progress in Oceanography. 47:355-379.   10.1016/s0079-6611(00)00044-6   AbstractWebsite

Basin-scale variations in oceanic physical variables are thought to organize patterns of biological response across the Pacific Ocean over decadal time scales. Different physical mechanisms can be responsible for the diverse basin-scale patterns of sea-surface temperature (SST), mixed-layer depth, thermocline depth, and horizontal currents, although they are linked in various ways. In light of various theories and observations, we interpret observed basinwide patterns of decadal-scale variations in upper-ocean temperatures. Evidence so far indicates that large-scale perturbations of the Aleutian Low generate temperature anomalies in the central and eastern North Pacific through the combined action of net surface heat Aux, turbulent mixing and Ekman advection. The surface-forced temperature anomalies in the central North Pacific subduct and propagate southwestwards in the ocean thermocline to the subtropics but apparently do not reach the equator. The large-scale Ekman pumping resulting from changes of the Aleutian Low forces western-intensified thermocline depth anomalies that are approximately consistent with Sverdrup theory. These thermocline changes are associated with SST anomalies in the Kuroshio/Oyashio Extension that are of the same sign as those in the central North Pacific, but lagged by several years. The physics of the possible feedback from the SST anomalies to the Aleutian Low, which might close a coupled ocean-atmosphere mode of decadal variability, is poorly understood and is an area of active research. The possible responses of North Pacific Ocean ecosystems to these complicated physical patterns is summarized. (C) 2000 Elsevier Science Ltd. All rights reserved.

1999
Miller, AJ, Cornuelle BD.  1999.  Forecasts from fits of frontal fluctuations. Dynamics of Atmospheres and Oceans. 29:305-333.   10.1016/s0377-0265(99)00009-3   AbstractWebsite

A primitive equation ocean model is fit with strong constraints to non-synoptic hydrographic surveys in an unstable frontal current region, the Iceland-Faeroe Front. The model is first initialized from a time-independent objective analysis of non-synoptic data (spanning 2 to 6 days). A truncated set of eddy-scale basis functions is used to represent the initial error in temperature, salinity, and velocity. A series of model integrations, each perturbed with one basis function for one dependent variable in one layer, is used to determine the sensitivity to the objective-analysis initial state of the match to the non-synoptic hydrographic data. A new initial condition is then determined from a generalized inverse of the sensitivity matrix and the process is repeated to account for non-linearity. The method is first tested in 'identical twin' experiments to demonstrate the adequacy of the basis functions in representing initial condition error and the convergence of the method to the true solution. The approach is then applied to observations gathered in August 1993 in the Iceland-Faeroe Front. Model fits are successful in improving the match to the true data, leading to dynamically consistent evolution scenarios. However, the forecast skill (here defined as the variance of the model-data differences) of the model runs from the optimized initial condition is not superior to less sophisticated methods of initialization, probably due to inadequate initialization data. The limited verification data in the presence of strong frontal slopes may not be sufficient to establish Forecast skill, so that it must be judged subjectively or evaluated by other quantitative measures. (C) 1999 Elsevier Science B.V. All rights reserved.

1996
Nese, JM, Miller AJ, Dutton JA.  1996.  The nature of predictability enhancement in a low-order ocean-atmosphere model. Journal of Climate. 9:2167-2172.   10.1175/1520-0442(1996)009<2167:tnopei>2.0.co;2   AbstractWebsite

A low-order moist general circulation model of the coupled ocean-atmosphere system is reexamined to determine the source of short-term predictability enhancement that occurs when an oceanic circulation is activated. The predictability enhancement is found to originate predominantly in thermodynamic processes involving changes in the mean hydrologic cycle of the model, which arise because the mean sea surface temperature is altered by the oceanic circulation. Thus, time-dependent sea surface temperature anomalies forced by anomalous geostrophic currents in the altered mean conditions do not contribute to the dominant ocean-atmosphere feedback mechanism that causes the predictability enhancement in the model.

Miller, AJ.  1996.  Recent advances in California Current modeling: Decadal and interannual thermocline variations. California Cooperative Oceanic Fisheries Investigations Reports. 37:69-79. AbstractWebsite

Some recent advances in large-scale modeling of the California Current and its interaction with basin-scale circulation and forcing are summarized. The discussion concentrates on a decadal-scale change and interannual-scale variations identified in the thermocline off the California Coast. The western-intensified decadal-scale change is part of a basinwide change in the North Pacific thermocline from the early 1970s to the early 1980s, which has been observed and modeled. The decadal change is driven by a basin-scale change in wind stress curl (Ekman pumping) and is associated with a deepening of the thermocline off California but no significant change in the strength of the California Current. Interannual variations of the thermocline off the California Coast, which tend to be associated with ENSO, have also been observed and modeled. High-resolution models often exhibit a coastal-trapped Kelvin-like wave arriving from the equatorial zone, but even a coarse-resolution model can capture aspects of the midlatitude wind-forced thermocline signals that propagate westward on ENSO time scales.

Miller, AJ, Lermusiaux PFJ, Poulain PM.  1996.  A topographic-Rossby mode resonance over the Iceland-Faeroe Ridge. Journal of Physical Oceanography. 26:2735-2747.   10.1175/1520-0485(1996)026<2735:atmrot>2.0.co;2   AbstractWebsite

An array of current meter moorings along 12 degrees W on the southern side of the Iceland-Faeroe Ridge reveals a narrowband barotropic oscillation with period 1.8 days in spectra of velocity. The signal is coherent over at least 55-km scales and propagates phase with shallow water on the right (toward the northwest). Velocity ellipses tend to be elongated (crossing contours of f/H) and rotate anticyclonically. Solutions of the rigid-lid barotropic shallow-water equations predict the occurrence of a topographic-Rossby normal mode on the south side of the ridge with spatial Scales exceeding 250 km and with intrinsic period near 1.84 days. This fundamental mode of the south side of the ridge has predicted spatial structure, phase propagation, and velocity ellipses consistent with the observed oscillation. The frictional amplitude e-folding decay time for this normal mode is estimated from the observations to be 13 days. The observed ocean currents are significantly coherent with zonal wind stress fluctuations (but not with wind stress curl) in the relevant period band, which indicates the oscillation is wind forced. This appears to be the first clear evidence of a stochastically forced resonant barotropic topographic-Rossby normal mode in the ocean.

1994
Miller, AJ, Cayan DR, Barnett TP, Graham NE, Oberhuber JM.  1994.  Interdecadal variability of the Pacific Ocean: model response to observed heat flux and wind stress anomalies. Climate Dynamics. 9:287-302.   10.1007/bf00204744   AbstractWebsite

Variability of the Pacific Ocean is examined in numerical simulations with an ocean general circulation model forced by observed anomalies of surface heat flux, wind stress and turbulent kinetic energy (TKE) over the period 1970-88. The model captures the 1976-77 winter time climate shift in sea surface temperature, as well as its monthly, seasonal and longer term variability as evidenced in regional time series and empirical orthogonal function analyses. Examination of the surface mixed-layer heat budget reveals that the 1976-77 shift was caused by a unique concurrance of sustained heat flux input anomalies and very strong horizontal advection anomalies during a multi-month period preceding the shift in both the central Pacific region (where cooling occurred) and the California coastal region (where warming occurred). In the central Pacific, the warm conditions preceding and the cold conditions following the shift tend to be maintained by anomalous vertical mixing due to increases in the atmospheric momentum flux (TKE input) into the mixed layer (which deepens in the model after the shift) from the early 1970s to the late 1970s and 1980s. Since the ocean model does not contain feedback to the atmosphere and it succeeds in capturing the major features of the 1976-77 shift, it appears that the midlatitude part of the shift was driven by the atmosphere, although effects of midlatitude ocean-atmosphere feedback are still possible. The surface mixed-layer heat budget also reveals that, in the central Pacific, the effects of heat flux input and vertical mixing anomalies are comparable in amplitude while horizontal advection anomalies are roughly half that size. In the California coastal region, in contrast, where wind variability is much weaker than in the central Pacific, horizontal advection and vertical mixing effects on the mixed-layer heat budget are only one-quarter the size of typical heat flux input anomalies.