ALBERT

Description: The Geophysical Fluid Dynamics Laboratory (GFDL) has developed a coupled general circulation model (CM3) for the atmosphere, oceans, land, and sea ice. The goal of CM3 is to address emerging issues in climate change, including aerosol–cloud interactions, chemistry–climate interactions, and coupling between the troposphere and stratosphere. The model is also designed to serve as the physical system component of earth system models and models for decadal prediction in the near-term future—for example, through improved simulations in tropical land precipitation relative to earlier-generation GFDL models. This paper describes the dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component (AM3) of this model. Relative to GFDL AM2, AM3 includes new treatments of deep and shallow cumulus convection, cloud droplet activation by aerosols, subgrid variability of stratiform vertical velocities for droplet activation, and atmospheric chemistry driven by emissions with advective, convective, and turbulent transport. AM3 employs a cubed-sphere implementation of a finite-volume dynamical core and is coupled to LM3, a new land model with ecosystem dynamics and hydrology. Its horizontal resolution is approximately 200 km, and its vertical resolution ranges approximately from 70 m near the earth’s surface to 1 to 1.5 km near the tropopause and 3 to 4 km in much of the stratosphere. Most basic circulation features in AM3 are simulated as realistically, or more so, as in AM2. In particular, dry biases have been reduced over South America. In coupled mode, the simulation of Arctic sea ice concentration has improved. AM3 aerosol optical depths, scattering properties, and surface clear-sky downward shortwave radiation are more realistic than in AM2. The simulation of marine stratocumulus decks remains problematic, as in AM2. The most intense 0.2% of precipitation rates occur less frequently in AM3 than observed. The last two decades of the twentieth century warm in CM3 by 0.32°C relative to 1881–1920. The Climate Research Unit (CRU) and Goddard Institute for Space Studies analyses of observations show warming of 0.56° and 0.52°C, respectively, over this period. CM3 includes anthropogenic cooling by aerosol–cloud interactions, and its warming by the late twentieth century is somewhat less realistic than in CM2.1, which warmed 0.66°C but did not include aerosol–cloud interactions. The improved simulation of the direct aerosol effect (apparent in surface clear-sky downward radiation) in CM3 evidently acts in concert with its simulation of cloud–aerosol interactions to limit greenhouse gas warming.

Description: The Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), CloudSat radar, and the Moderate Resolution Imaging Spectroradiometer (MODIS) cloud data on the A-Train constellation complemented with the European Centre for Medium-Range Forecasts (ECMWF) analyses are used to investigate the cloud and boundary layer structure across a 10° wide cross section starting at 5°S near the international date line and extending to 35°N near the California coast from March 2008 to February 2009. The mean large-scale inversion height and low-level cloud tops, which correspond very closely to each other, are very shallow (∼500 m) over cold SSTs and high static stability near California and deepen southwestward (to a maximum of ∼1.5–2.0 km) along the cross section as SSTs rise. Deep convection near the ITCZ occurs at a surface temperature close to 298 K. While the boundary layer relative humidity (RH) is nearly constant where a boundary layer is well defined, it drops sharply near cloud top in stratocumulus regions, corresponding with strong thermal inversions and water vapor decrease, such that the maximum (−∂RH/∂z) marks the boundary layer cloud top very well. The magnitude correlates well with low cloud frequency during March–May (MAM), June–August (JJA), and September–November (SON) (r 2 = 0.85, 0.88, and 0.86, respectively). Also, CALIPSO and MODIS isolated low cloud frequency generally agree quite well, but CloudSat senses only slightly more than one-third of the low clouds as observed by the other sensors, as many clouds are shallower than 1 km and thus cannot be discerned with CloudSat due to contamination from the strong signal from surface clutter. Mean tropospheric ω between 300 and 700 hPa is examined from the ECMWF Year of Tropical Convection (YOTC) analysis dataset, and during JJA and SON, strong rising motion in the middle troposphere is confined to a range of 2-m surface temperatures between 297 and 300 K, consistent with previous studies that show a narrow range of SSTs over which deep ascent occurs. During December–February (DJF), large-scale ascending motion extends to colder SSTs and high boundary layer stability. A slightly different boundary layer stability metric is derived, the difference of moist static energy (MSE) at the middle point of the inversion (or at 700 hPa if no inversion exists) and the surface, referred to as ΔMSE. The utility of ΔMSE is its prediction of isolated uniform low cloud frequency, with very high r 2 values of 0.93 and 0.88, respectively, for the MODIS and joint lidar plus radar product during JJA but significantly lower values during DJF (0.46 and 0.40), with much scatter. To quantify the importance of free tropospheric dynamics in modulating the ΔMSE–low cloud relationships, the frequency as a function of ΔMSE of rising motion profiles (ω 〈 −0.05 Pa s−1) is added to the observed low cloud frequency for a maximum hypothetical low cloud frequency. Doing this greatly reduces the interseasonal differences and holds promise for using ΔMSE for parameterization schemes and examining low cloud feedbacks.

Description: The authors exploit three years of data from the CloudSat and Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellites to document for the first time the seasonally varying vertical structure of cloudiness throughout Antarctica and the high-latitude Southern Ocean. The results provide a baseline reference of Southern Hemisphere high-latitude cloudiness for future observational and modeling studies, and they highlight several previously undocumented aspects and key features of Antarctic cloudiness. The key features of high-latitude Southern Hemisphere cloudiness documented here include 1) a pronounced seasonal cycle in cloudiness over the high-latitude Southern Hemisphere, with higher cloud incidences generally found during the winter season over both the Southern Ocean and Antarctica; 2) two distinct maxima in vertical profiles of cloud incidence over the Southern Ocean, one centered near the surface and another centered in the upper troposphere; 3) a nearly discontinuous drop-off in cloudiness near 8 km over much of the continent that peaks during autumn, winter, and spring; 4) large east–west gradients in upper-level cloudiness in the vicinity of the Antarctic Peninsula that peak during the austral spring season; and 5) evidence that cloudiness in the polar stratosphere is marked not by a secondary maximum at stratospheric levels but by a nearly monotonic decrease with height from the tropopause. Key results are interpreted in the context of the seasonally varying profiles of vertical motion and static stability and compared with results of previous studies.

Description: This paper documents time mean simulation characteristics from the ocean and sea ice components in a new coupled climate model developed at the NOAA Geophysical Fluid Dynamics Laboratory (GFDL). The GFDL Climate Model version 3 (CM3) is formulated with effectively the same ocean and sea ice components as the earlier CM2.1 yet with extensive developments made to the atmosphere and land model components. Both CM2.1 and CM3 show stable mean climate indices, such as large-scale circulation and sea surface temperatures (SSTs). There are notable improvements in the CM3 climate simulation relative to CM2.1, including a modified SST bias pattern and reduced biases in the Arctic sea ice cover. The authors anticipate SST differences between CM2.1 and CM3 in lower latitudes through analysis of the atmospheric fluxes at the ocean surface in corresponding Atmospheric Model Intercomparison Project (AMIP) simulations. In contrast, SST changes in the high latitudes are dominated by ocean and sea ice effects absent in AMIP simulations. The ocean interior simulation in CM3 is generally warmer than in CM2.1, which adversely impacts the interior biases.

Description: Continuous estimates of the oceanic meridional heat transport in the Atlantic are derived from the Rapid Climate Change–Meridional Overturning Circulation (MOC) and Heatflux Array (RAPID–MOCHA) observing system deployed along 26.5°N, for the period from April 2004 to October 2007. The basinwide meridional heat transport (MHT) is derived by combining temperature transports (relative to a common reference) from 1) the Gulf Stream in the Straits of Florida; 2) the western boundary region offshore of Abaco, Bahamas; 3) the Ekman layer [derived from Quick Scatterometer (QuikSCAT) wind stresses]; and 4) the interior ocean monitored by “endpoint” dynamic height moorings. The interior eddy heat transport arising from spatial covariance of the velocity and temperature fields is estimated independently from repeat hydrographic and expendable bathythermograph (XBT) sections and can also be approximated by the array. The results for the 3.5 yr of data thus far available show a mean MHT of 1.33 ± 0.40 PW for 10-day-averaged estimates, on which time scale a basinwide mass balance can be reasonably assumed. The associated MOC strength and variability is 18.5 ± 4.9 Sv (1 Sv ≡ 106 m3 s−1). The continuous heat transport estimates range from a minimum of 0.2 to a maximum of 2.5 PW, with approximately half of the variance caused by Ekman transport changes and half caused by changes in the geostrophic circulation. The data suggest a seasonal cycle of the MHT with a maximum in summer (July–September) and minimum in late winter (March–April), with an annual range of 0.6 PW. A breakdown of the MHT into “overturning” and “gyre” components shows that the overturning component carries 88% of the total heat transport. The overall uncertainty of the annual mean MHT for the 3.5-yr record is 0.14 PW or about 10% of the mean value.

Description: Hurricane Gustav (2008) made landfall in southern Louisiana on 1 September 2008 with its eye never closer than 75 km to New Orleans, but its waves and storm surge threatened to flood the city. Easterly tropical-storm-strength winds impacted the region east of the Mississippi River for 12–15 h, allowing for early surge to develop up to 3.5 m there and enter the river and the city’s navigation canals. During landfall, winds shifted from easterly to southerly, resulting in late surge development and propagation over more than 70 km of marshes on the river’s west bank, over more than 40 km of Caernarvon marsh on the east bank, and into Lake Pontchartrain to the north. Wind waves with estimated significant heights of 15 m developed in the deep Gulf of Mexico but were reduced in size once they reached the continental shelf. The barrier islands further dissipated the waves, and locally generated seas existed behind these effective breaking zones. The hardening and innovative deployment of gauges since Hurricane Katrina (2005) resulted in a wealth of measured data for Gustav. A total of 39 wind wave time histories, 362 water level time histories, and 82 high water marks were available to describe the event. Computational models—including a structured-mesh deepwater wave model (WAM) and a nearshore steady-state wave (STWAVE) model, as well as an unstructured-mesh “simulating waves nearshore” (SWAN) wave model and an advanced circulation (ADCIRC) model—resolve the region with unprecedented levels of detail, with an unstructured mesh spacing of 100–200 m in the wave-breaking zones and 20–50 m in the small-scale channels. Data-assimilated winds were applied using NOAA’s Hurricane Research Division Wind Analysis System (H*Wind) and Interactive Objective Kinematic Analysis (IOKA) procedures. Wave and surge computations from these models are validated comprehensively at the measurement locations ranging from the deep Gulf of Mexico and along the coast to the rivers and floodplains of southern Louisiana and are described and quantified within the context of the evolution of the storm.

Description: Although the upper-layer dynamics of the Loop Current and eddies in the Gulf of Mexico are well studied, the understanding of how they are coupled to the deep flows is limited. In this work, results from a numerical model are analyzed to classify the expansion, shedding, retraction, and deep-coupling cycle (the Loop Current cycle) according to the vertical mass flux across the base of the Loop. Stage A is the “Loop reforming” period, with downward flux and deep divergence under the Loop Current. Stage B is the “incipient shedding,” with strong upward flux and deep convergence. Stage C is the “eddy migration,” with waning upward flux and deep throughflow from the western Gulf into the Yucatan Channel. Because of the strong deep coupling between the eastern and western Gulf, the Loop’s expansion is poorly correlated with deep flows through the Yucatan Channel. Stage A is longest and the mean vertical flux under the Loop Current is downward. Therefore, because the net circulation around the abyssal basin is zero, the abyssal gyre in the western Gulf is cyclonic. The gyre’s strength is strongest when the Loop Current is reforming and weakest after an eddy is shed. The result suggests that the Loop Current cycle can force a low-frequency [time scales ∼ shedding periods; O(months)] abyssal oscillation in the Gulf of Mexico.

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 40 (2010): 2743–2756, doi:10.1175/2010JPO4339.1.

Description: Analysis of modern and historical observations demonstrates that the temperature of the intermediate-depth (150–900 m) Atlantic water (AW) of the Arctic Ocean has increased in recent decades. The AW warming has been uneven in time; a local 1°C maximum was observed in the mid-1990s, followed by an intervening minimum and an additional warming that culminated in 2007 with temperatures higher than in the 1990s by 0.24°C. Relative to climatology from all data prior to 1999, the most extreme 2007 temperature anomalies of up to 1°C and higher were observed in the Eurasian and Makarov Basins. The AW warming was associated with a substantial (up to 75–90 m) shoaling of the upper AW boundary in the central Arctic Ocean and weakening of the Eurasian Basin upper-ocean stratification. Taken together, these observations suggest that the changes in the Eurasian Basin facilitated greater upward transfer of AW heat to the ocean surface layer. Available limited observations and results from a 1D ocean column model support this surmised upward spread of AW heat through the Eurasian Basin halocline. Experiments with a 3D coupled ice–ocean model in turn suggest a loss of 28–35 cm of ice thickness after 50 yr in response to the 0.5 W m−2 increase in AW ocean heat flux suggested by the 1D model. This amount of thinning is comparable to the 29 cm of ice thickness loss due to local atmospheric thermodynamic forcing estimated from observations of fast-ice thickness decline. The implication is that AW warming helped precondition the polar ice cap for the extreme ice loss observed in recent years.

Description: This study was supported by JAMSTEC (IP and VI), NOAA (IP, VI, and ID), NSF (IP,VA,VI, ID, JT, andMS),NASA(IP andVI), BMBF (ID), and UK NERC (SB) grants.

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 40 (2010): 2605–2623, doi:10.1175/2010JPO4132.1.

Description: Steady scale-invariant solutions of a kinetic equation describing the statistics of oceanic internal gravity waves based on wave turbulence theory are investigated. It is shown in the nonrotating scale-invariant limit that the collision integral in the kinetic equation diverges for almost all spectral power-law exponents. These divergences come from resonant interactions with the smallest horizontal wavenumbers and/or the largest horizontal wavenumbers with extreme scale separations. A small domain is identified in which the scale-invariant collision integral converges and numerically find a convergent power-law solution. This numerical solution is close to the Garrett–Munk spectrum. Power-law exponents that potentially permit a balance between the infrared and ultraviolet divergences are investigated. The balanced exponents are generalizations of an exact solution of the scale-invariant kinetic equation, the Pelinovsky–Raevsky spectrum. A small but finite Coriolis parameter representing the effects of rotation is introduced into the kinetic equation to determine solutions over the divergent part of the domain using rigorous asymptotic arguments. This gives rise to the induced diffusion regime. The derivation of the kinetic equation is based on an assumption of weak nonlinearity. Dominance of the nonlocal interactions puts the self-consistency of the kinetic equation at risk. However, these weakly nonlinear stationary states are consistent with much of the observational evidence.

Description: This research is supported by NSF CMG Grants 0417724, 0417732 and 0417466. YL is also supported by NSF DMS Grant 0807871 and ONR Award N00014-09-1-0515.

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 40 (2010): 2679–2695, doi:10.1175/2010JPO4395.1.

Description: Observations of stratification and currents between June 2007 and March 2009 reveal a strong overflow between 400- and 570-m depth from the Panay Strait into the Sulu Sea. The overflow water is derived from approximately 400 m deep in the South China Sea. Temporal mean velocity is greater than 0.75 m s−1 at 50 m above the 570-m Panay Sill. Empirical orthogonal function analysis of a mooring time series shows that the flow is dominated by the bottom overflow current with little seasonal variance. The overflow does not descend below 1250 m in the Sulu Sea but rather settles above high-salinity deep water derived from the Sulawesi Sea. The mean observed overflow transport at the sill is 0.32 × 106 m3 s−1. The observed transport was used to calculate a bulk diapycnal diffusivity of 4.4 × 10−4 m2 s−1 within the Sulu Sea slab (575–1250 m) ventilated from Panay Strait. Analysis of Froude number variation across the sill shows that the flow is hydraulically controlled. A suitable hydraulic control model shows overflow transport equivalent to the observed overflow. Thorpe-scale estimates show turbulent dissipation rates up to 5 × 10−7 W kg−1 just downstream of the supercritical to subcritical flow transition, suggesting a hydraulic jump downstream of the sill.

Description: This work was supported by the Office of Naval Research Grant N00014-09-1-0582 to Lamont-Doherty Earth Observatory of Columbia University; Grants ONR-13759000 and N00014-09-1-0582 to the Woods Hole Oceanographic Institution; Grant ONR-N00014-06-1-0690 to Scripps Institute of Oceanography; and a National Defense Science and Engineering Graduate Fellowship.