Geophysical Fluid Dynamics Laboratory

About: Geophysical Fluid Dynamics Laboratory is a facility organization based out in Princeton, New Jersey, United States. It is known for research contribution in the topics: Climate model & Climate change. The organization has 525 authors who have published 2432 publications receiving 264545 citations. The organization is also known as: GFDL.

The KPP Boundary Layer Scheme for the Ocean: Revisiting Its Formulation and Benchmarking One-Dimensional Simulations Relative to LES

Abstract: We evaluate the Community ocean Vertical Mixing (CVMix) project version of the K-profile parameterization (KPP). For this purpose, one-dimensional KPP simulations are compared across a suite of oceanographically relevant regimes against large eddy simulations (LES). The LES is forced with horizontally uniform boundary fluxes and has horizontally uniform initial conditions, allowing its horizontal average to be compared to one-dimensional KPP tests. We find the standard configuration of KPP [Danabasoglu et al., 2006] consistent with LES across many forcing regimes, supporting the physical basis of KPP. Our evaluation motivates recommendations for ”best practices” for using KPP within ocean circulation models, and identifies areas where further research is warranted. Further, our test suite can be used as a baseline for evaluation of a broad suite of boundary layer models. The original treatment of KPP recommends the matching of interior diffusivities and their gradients to the KPP predicted values computed in the ocean surface boundary layer (OSBL). However, we find that difficulties in representing derivatives of rapidly changing diffusivities near the base of the OSBL can lead to loss of simulation fidelity. We propose two alternative approaches: (1) match to the internal predicted diffusivity along, (2) set the KPP diffusivity to zero at the OSBL base. Although computationally simpler, the second alternative is sensitive to implementation details and we offer methods to prevent the emergence of numerical high frequency noise. We find the KPP entrainment buoyancy flux to be sensitive to vertical grid resolution and details of how to diagnose the KPP boundary layer depth. We modify the KPP turbulent shear velocity parameterization to reduce resolution dependence. Additionally, our results show that the KPP parameterized non-local tracer flux is incomplete due to the assumption that it solely redistributes the surface tracer flux. However, examination of the LES vertical turbulent scalar flux budgets show that non-local fluxes can exist in the absence of surface tracer fluxes. This result motivates further studies of the non-local flux parameterization. Draft from March 9, 2018

How ice shelf morphology controls basal melting

Abstract: [1] The response of ice shelf basal melting to climate is a function of ocean temperature, circulation, and mixing in the open ocean and the coupling of this external forcing to the sub–ice shelf circulation. Because slope strongly influences the properties of buoyancy-driven flow near the ice shelf base, ice shelf morphology plays a critical role in linking external, subsurface heat sources to the ice. In this paper, the slope-driven dynamic control of local and area-integrated melting rates is examined under a wide range of ocean temperatures and ice shelf shapes, with an emphasis on smaller, steeper ice shelves. A 3-D numerical ocean model is used to simulate the circulation underneath five idealized ice shelves, forced with subsurface ocean temperatures ranging from −2.0°C to 1.5°C. In the sub–ice shelf mixed layer, three spatially distinct dynamic regimes are present. Entrainment of heat occurs predominately under deeper sections of the ice shelf; local and area-integrated melting rates are most sensitive to changes in slope in this “initiation” region. Some entrained heat is advected upslope and used to melt ice in the “maintenance” region; however, flow convergence in the “outflow” region limits heat loss in flatter portions of the ice shelf. Heat flux to the ice exhibits (1) a spatially nonuniform, superlinear dependence on slope and (2) a shape- and temperature-dependent, internally controlled efficiency. Because the efficiency of heat flux through the mixed layer decreases with increasing ocean temperature, numerical simulations diverge from a simple quadratic scaling law.

Impact of meteorology and emissions on methane trends, 1990-2004

Abstract: [1] Over the past century, atmospheric methane (CH4) rose dramatically before leveling off in the late 1990s. The processes controlling this trend are poorly understood, limiting confidence in projections of future CH4. The MOZART-2 global tropospheric chemistry model qualitatively captures the observed CH4 trend (increasing in the early 1990s and then leveling off) with constant emissions. From 1991–1995 to 2000–2004, the CH4 lifetime versus tropospheric OH decreases by 1.6%, reflecting increases in OH and temperature. The rise in OH stems from an increase in lightning NOx as parameterized in the model. A simulation including annually varying anthropogenic and wetland CH4 emissions, as well as the changes in meteorology, best reproduces the observed CH4 distribution, trend, and seasonal cycles. Projections of future CH4 abundances should consider climate-driven changes in CH4 sources and sinks.

Sensitivity of nitrate aerosols to ammonia emissions and to nitrate chemistry: implications for present and future nitrate optical depth

Abstract: . We update and evaluate the treatment of nitrate aerosols in the Geophysical Fluid Dynamics Laboratory (GFDL) atmospheric model (AM3). Accounting for the radiative effects of nitrate aerosols generally improves the simulated aerosol optical depth, although nitrate concentrations at the surface are biased high. This bias can be reduced by increasing the deposition of nitrate to account for the near-surface volatilization of ammonium nitrate or by neglecting the heterogeneous production of nitric acid to account for the inhibition of N2O5 reactive uptake at high nitrate concentrations. Globally, uncertainties in these processes can impact the simulated nitrate optical depth by up to 25 %, much more than the impact of uncertainties in the seasonality of ammonia emissions (6 %) or in the uptake of nitric acid on dust (13 %). Our best estimate for fine nitrate optical depth at 550 nm in 2010 is 0.006 (0.005–0.008). In wintertime, nitrate aerosols are simulated to account for over 30 % of the aerosol optical depth over western Europe and North America. Simulated nitrate optical depth increases by less than 30 % (0.0061–0.010) in response to projected changes in anthropogenic emissions from 2010 to 2050 (e.g., −40 % for SO2 and +38 % for ammonia). This increase is primarily driven by greater concentrations of nitrate in the free troposphere, while surface nitrate concentrations decrease in the midlatitudes following lower concentrations of nitric acid. With the projected increase of ammonia emissions, we show that better constraints on the vertical distribution of ammonia (e.g., convective transport and biomass burning injection) and on the sources and sinks of nitric acid (e.g., heterogeneous reaction on dust) are needed to improve estimates of future nitrate optical depth.

Impacts of a Pinatubo-size volcanic eruption on ENSO

Abstract: Observations and model simulations of the climate responses to strong explosive low-latitude volcanic eruptions suggest a significant increase in the likelihood of El Nino during the eruption and post-eruption years, though model results have been inconclusive and have varied in magnitude and even sign. In this study, we test how this spread of responses depends on the initial phase of El Nino-Southern Oscillation (ENSO) in the eruption year, and on the eruption's seasonal timing. We employ the GFDL CM2.1 global coupled general circulation model to investigate the impact of the Pinatubo 1991 eruption, assuming that in 1991 ENSO would otherwise be in Central or Eastern Pacific El Nino, La Nina, or neutral phases. We obtain statistically significant El Nino responses in a year after the eruption for all cases except La Nina, which shows no response in the eastern equatorial Pacific. The eruption has a weaker impact on Eastern Pacific El Ninos than on Central Pacific El Ninos. We find that the ocean dynamical thermostat, and (to a lesser extent) wind changes due to land-ocean temperature gradients, are the main feedbacks affecting El Nino development after the eruption. The El Nino responses to eruptions occurring in summer are more pronounced than for winter and spring eruptions. That the climate response depends on eruption season and initial ENSO phase may help to reconcile apparent inconsistencies among previous studies.