Gregory J. Hakim

Bio: Gregory J. Hakim is an academic researcher from University of Washington. The author has contributed to research in topics: Data assimilation & Climate model. The author has an hindex of 34, co-authored 100 publications receiving 3246 citations. Previous affiliations of Gregory J. Hakim include State University of New York System & National Center for Atmospheric Research.

Consistent multi-decadal variability in global temperature reconstructions and simulations over the Common Era.

Abstract: Multi-decadal surface temperature changes may be forced by natural as well as anthropogenic factors, or arise unforced from the climate system. Distinguishing these factors is essential for estimating sensitivity to multiple climatic forcings and the amplitude of the unforced variability. Here we present 2,000-year-long global mean temperature reconstructions using seven different statistical methods that draw from a global collection of temperature-sensitive paleoclimate records. Our reconstructions display synchronous multi-decadal temperature fluctuations, which are coherent with one another and with fully forced CMIP5 millennial model simulations across the Common Era. The most significant attribution of pre-industrial (1300-1800 CE) variability at multi-decadal timescales is to volcanic aerosol forcing. Reconstructions and simulations qualitatively agree on the amplitude of the unforced global mean multi-decadal temperature variability, thereby increasing confidence in future projections of climate change on these timescales. The largest warming trends at timescales of 20 years and longer occur during the second half of the 20th century, highlighting the unusual character of the warming in recent decades.

Boundary Conditions for Limited-Area Ensemble Kalman Filters

Abstract: One aspect of implementing a limited-area ensemble Kalman filter (EnKF) involves the specification of a suitable ensemble of lateral boundary conditions. Two classes of methods to populate a boundary condition ensemble are proposed. In the first class, the ensemble of boundary conditions is provided by an EnKF on a larger domain and is approximately a random draw from the probability distribution function for the forecast (or analysis) on the limited-area domain boundary. The second class perturbs around a deterministic estimate of the state using assumed spatial and temporal covariance relationships. Methods in the second class are relatively flexible and easy to implement. Experiments that test the utility of these methods are performed for both an idealized low-dimensional model and limited-area simulations using the Weather Research and Forecasting (WRF) model; all experiments employ simulated observations under the perfect model assumption. The performance of the ensemble boundary condition ...

Ensemble-Based Sensitivity Analysis

Abstract: The sensitivity of forecasts to observations is evaluated using an ensemble approach with data drawn from a pseudo-operational ensemble Kalman filter. For Gaussian statistics and a forecast metric defined as a scalar function of the forecast variables, the effect of observations on the forecast metric is quantified by changes in the metric mean and variance. For a single observation, expressions for these changes involve a product of scalar quantities, which can be rapidly evaluated for large numbers of observations. This technique is applied to determining climatological forecast sensitivity and predicting the impact of observations on sea level pressure and precipitation forecast metrics. The climatological 24-h forecast sensitivity of the average pressure over western Washington State shows a region of maximum sensitivity to the west of the region, which tilts gently westward with height. The accuracy of ensemble sensitivity predictions is tested by withholding a single buoy pressure observation from this region and comparing this perturbed forecast with the control case where the buoy is assimilated. For 30 cases, there is excellent agreement between these forecast differences and the ensemble predictions, as measured by the forecast metric. This agreement decreases for increasing numbers of observations. Nevertheless, by using statistical confidence tests to address sampling error, the impact of thousands of observations on forecast-metric variance is shown to be well estimated by a subset of the O(100) most significant observations.

Comparing Adjoint- and Ensemble-Sensitivity Analysis with Applications to Observation Targeting

Abstract: The sensitivity of numerical weather forecasts to small changes in initial conditions is estimated using ensemble samples of analysis and forecast errors. Ensemble sensitivity is defined here by linear regression of analysis errors onto a given forecast metric. It is shown that ensemble sensitivity is proportional to the projection of the analysis-error covariance onto the adjoint-sensitivity field. Furthermore, the ensemble-sensitivity approach proposed here involves a small calculation that is easy to implement. Ensemble- and adjoint-based sensitivity fields are compared for a representative wintertime flow pattern near the west coast of North America for a 90-member ensemble of independent initial conditions derived from an ensemble Kalman filter. The forecast metric is taken for simplicity to be the 24-h forecast of sea level pressure at a single point in western Washington State. Results show that adjoint and ensemble sensitivities are very different in terms of location, scale, and magnitud...

The last millennium climate reanalysis project: Framework and first results

Abstract: An “offline” approach to DA is used, where static ensemble samples are drawn from existing CMIP climate-model simulations to serve as the prior estimate of climate variables. We use linear, univariate forward models (“proxy system models (PSMs)”) that map climate variables to proxy measurements by fitting proxy data to 2 m air temperature from gridded instrumental temperature data; the linear PSMs are then used to predict proxy values from the prior estimate. Results for the LMR are compared against six gridded instrumental temperature data sets and 25% of the proxy records are withheld from assimilation for independent verification. Results show broad agreement with previous reconstructions of Northern Hemisphere mean 2 m air temperature, with millennial-scale cooling, a multicentennial warm period around 1000 C.E., and a cold period coincident with the Little Ice Age (circa 1450–1800 C.E.). Verification against gridded instrumental data sets during 1880–2000 C.E. reveals greatest skill in the tropics and lowest skill over Northern Hemisphere land areas. Verification against independent proxy records indicates substantial improvement relative to the model (prior) data without proxy assimilation. As an illustrative example, we present multivariate reconstructed fields for a singular event, the 1808/1809 “mystery” volcanic eruption, which reveal global cooling that is strongly enhanced locally due to the presence of the Pacific-North America wave pattern in the 500 hPa geopotential height field.

Numerical study of convection observed during the Winter Monsoon Experiment using a mesoscale two-dimensional model [presentation]

Abstract: Abstract A two-dimensional version of the Pennsylvania State University mesoscale model has been applied to Winter Monsoon Experiment data in order to simulate the diurnally occurring convection observed over the South China Sea. The domain includes a representation of part of Borneo as well as the sea so that the model can simulate the initiation of convection. Also included in the model are parameterizations of mesoscale ice phase and moisture processes and longwave and shortwave radiation with a diurnal cycle. This allows use of the model to test the relative importance of various heating mechanisms to the stratiform cloud deck, which typically occupies several hundred kilometers of the domain. Frank and Cohen's cumulus parameterization scheme is employed to represent vital unresolved vertical transports in the convective area. The major conclusions are: Ice phase processes are important in determining the level of maximum large-scale heating and vertical motion because there is a strong anvil componen...

The Ensemble Kalman Filter: Theoretical formulation and practical implementation

Abstract: The purpose of this paper is to provide a comprehensive presentation and interpretation of the Ensemble Kalman Filter (EnKF) and its numerical implementation. The EnKF has a large user group, and numerous publications have discussed applications and theoretical aspects of it. This paper reviews the important results from these studies and also presents new ideas and alternative interpretations which further explain the success of the EnKF. In addition to providing the theoretical framework needed for using the EnKF, there is also a focus on the algorithmic formulation and optimal numerical implementation. A program listing is given for some of the key subroutines. The paper also touches upon specific issues such as the use of nonlinear measurements, in situ profiles of temperature and salinity, and data which are available with high frequency in time. An ensemble based optimal interpolation (EnOI) scheme is presented as a cost-effective approach which may serve as an alternative to the EnKF in some applications. A fairly extensive discussion is devoted to the use of time correlated model errors and the estimation of model bias.

The Self-Organizing Map

Abstract: An overview of the self-organizing map algorithm, on which the papers in this issue are based, is presented in this article.

Geophysical fluid dynamics.

Abstract: The potential for sea ice-albedo feedback to give rise to nonlinear climate change in the Arctic Ocean – defined as a nonlinear relationship between polar and global temperature change or, equivalently, a time-varying polar amplification – is explored in IPCC AR4 climate models. Five models supplying SRES A1B ensembles for the 21 st century are examined and very linear relationships are found between polar and global temperatures (indicating linear Arctic Ocean climate change), and between polar temperature and albedo (the potential source of nonlinearity). Two of the climate models have Arctic Ocean simulations that become annually sea ice-free under the stronger CO 2 increase to quadrupling forcing. Both of these runs show increases in polar amplification at polar temperatures above-5 o C and one exhibits heat budget changes that are consistent with the small ice cap instability of simple energy balance models. Both models show linear warming up to a polar temperature of-5 o C, well above the disappearance of their September ice covers at about-9 o C. Below-5 o C, surface albedo decreases smoothly as reductions move, progressively, to earlier parts of the sunlit period. Atmospheric heat transport exerts a strong cooling effect during the transition to annually ice-free conditions. Specialized experiments with atmosphere and coupled models show that the main damping mechanism for sea ice region surface temperature is reduced upward heat flux through the adjacent ice-free oceans resulting in reduced atmospheric heat transport into the region.