Glazer, R. H., & Misra, V. (2018). Ice versus liquid water saturation in simulations of the Indian summer monsoon.
Climate Dynamics, .
Goff, J. A., & Arbic, B. K. (2010). Global prediction of abyssal hill roughness statistics for use in ocean models from digital maps of paleo-spreading rate, paleo-ridge orientation, and sediment thickness.
Ocean Modelling, 32(1-2), 36–43.
Gouillon, F. (2010).
Internal Wave Propagation and Numerically Induced Diapycnal Mixing in Oceanic General Circulation Models. Ph.D. thesis, Florida State University, Tallahassee, FL.
Abstract: Numerical ocean models have become powerful tools for providing a realistic view of the ocean state and for describing ocean processes that are difficult to observe. Recent improvements in model performance focus on simulating realistic ocean interior mixing rates, as ocean mixing is the main physical process that creates water masses and maintains their properties. Below the mixed layer, diapycnal mixing primarily arises from the breaking of internal waves, whose energy is largely supplied by winds and tides. This is particularly true in abyssal regions, where the barotropic tide interacts with rough topography and where high levels of diapycnal mixing have been recorded (e.g., the Hawaiian Archipelago). Many studies have discussed the representation of internal wave generation, propagation, and evolution in ocean numerical models. Expanding on these studies, this work seeks to better understand the representation of internal wave dynamics, energetics, and their associated mixing in several different classes of widely used ocean models (e.g., the HYbrid Coordinate Ocean Model, HYCOM; the Regional Ocean Modeling System, ROMS; and the MIT general circulation model, MITgcm). First, a multi-model study investigates the representation of internal waves for a wide spectrum of numerical choices, such as the horizontal and vertical resolution, the vertical coordinate, and the choice of the numerical advection scheme. Idealized configurations are compared to their corresponding analytical solutions. Some preliminary results of realistic baroclinic tidal simulations are shown for the Gulf of Mexico. Second, the spurious diapycnal mixing that exists in models with fixed vertical coordinates (i.e., geopotential and terrain following) is documented and quantified. This purely numerical error arises because, in fixed-coordinate models, the numerical framework cannot properly maintain the adiabatic properties of an advected water parcel. This unrealistic mixing of water masses can be a source of major error in both regional and global ocean models. We use the tracer flux method to compute the spurious diapycnal diffusivities for both a lockexchange scenario and a propagating internal wave field using all three models. Results for the lock exchange experiments are compared to the results of a recent study. Our results, obtained by using three different model classes, provide a comprehensive analysis of the impact of model resolution choice and numerical framework on the magnitude of the spurious diapycnal mixing and the representation of internal waves.
Guimond, S. (2010).
Tropical Cyclone Inner-Core Dynamics: A Latent Heat Retrieval and Its Effects on Intensity and Structure Change; and the Impacts of Effective Diffusion on the Axisymmetrization Process. Ph.D. thesis, Florida State University, Tallahassee, FL.
Abstract: Despite the fact that latent heating in cloud systems drives many atmospheric circulations, including tropical cyclones, little is known of its magnitude and structure due in large part to inadequate observations. In this work, a reasonably high-resolution (2 km), four-dimensional airborne Doppler radar retrieval of the latent heat of condensation is presented for rapidly intensifying Hurricane Guillermo (1997). Several advancements in the retrieval algorithm are shown including: (1) analyzing the scheme within the dynamically consistent framework of a numerical model, (2) identifying algorithm sensitivities through the use of ancillary data sources and (3) developing a precipitation budget storage term parameterization. The determination of the saturation state is shown to be an important part of the algorithm for updrafts of ~ 5 m s-1 or less. The uncertainties in the magnitude of the retrieved heating are dominated by errors in the vertical velocity. Using a combination of error propagation and Monte Carlo uncertainty techniques, biases were found to be small, and randomly distributed errors in the heating magnitude were ~16 % for updrafts greater than 5 m s-1 and ~156 % for updrafts of 1 m s- 1. The impact of the retrievals is assessed by inserting the heating into realistic numerical simulations at 2 km resolution and comparing the generated wind structure to the Doppler radar observations of Guillermo. Results show that using the latent heat retrievals outperforms a simulation that relies on a state-of-the-art microphysics scheme (Reisner and Jeffery 2009), in terms of wind speed root-mean-square errors, explained variance and eye/eyewall structure. The incorrect transport of water vapor (a function of the sub-grid model and the numerical approximations to advection) and the restrictions on the magnitude of heat release that ensure the present model's stability are suggested as sources of error in the simulation without the retrievals. Motivated by the latent heat retrievals, the dynamics of vortex axisymmetrization from the perspective of thermal anomalies is investigated using an idealized, non-linear atmospheric model (HIGRAD). Attempts at reproducing the results of previous work (Nolan and Grasso 2003; NG03) revealed a discrepancy with the impacts of purely asymmetric forcing. While NG03 found that purely asymmetric heating led to a negligible, largely negative impact on the vortex intensification, in the present study the impacts of asymmetries are found to have an important, largely positive role. Absolute angular momentum budgets revealed that the essential difference between the present work and that of NG03 was the existence of a significant, axisymmetric secondary circulation in the basic-state vortex used in the HIGRAD simulations. This secondary circulation was larger than that present in NG03's simulations. The spin-up of the vortex caused by the asymmetric thermal anomalies was dominated by the axisymmetric fluxes of angular momentum at all times, indicating fundamentally different evolution of asymmetries in the presence of radial flow. Radial momentum budgets were performed to elucidate the mechanisms responsible for the formation of the physically significant secondary circulation. Results show that explicit (sub-grid) diffusion in the model was producing a gradient wind imbalance, which drives a radial inflow and associated secondary circulation in an attempt to re-gain balance. In addition, the production of vorticity anomalies from the asymmetric heating was found to be sensitive to the eddy diffusivity, with large differences between HIGRAD and the widely used WRF model for the exact same value of this uncertain parameter.
Hong, S. - Y., Park, H., Cheong, H. - B., Kim, J. - E. E., Koo, M. - S., Jang, J., et al. (2013). The Global/Regional Integrated Model system (GRIMs).
Asia-Pacific J Atmos Sci, 49(2), 219–243.
Kvaleberg, E. (2004).
Generation of Cold Core Filaments and Eddies Through Baroclinic Instability on a Continental Shelf. Ph.D. thesis, Florida State University, Tallahassee, FL.
Abstract: The formation of cold core filaments on an idealized continental shelf is investigated using a numerical model to simulate the ocean's response to surface cooling. A horizontal density gradient forms because of uneven buoyancy loss due to the sloping bottom, and this gradient induces an alongshelf current in thermal wind balance, that in time becomes unstable. As the instabilities grow, filaments, and later eddies, are generated so that dense water near the coast is mixed offshore. Scaling arguments of the filament wavelength indicate that the current is baroclinically unstable, and an analytical model of the frontal expansion with time is in very good agreement with the simulations. This study was inspired by satellite observations of sea surface temperature on the West Florida Shelf during the winter months, in which it is clearly seen that cold core filaments extend from a thermal front. Numerical experiments are therefore designed to allow for reliable comparisons with conditions in this region.
Le Sommer, J., Chassignet, E. P., & Wallcraft, A. J. (2018). Ocean Circulation Modeling for Operational Oceanography: Current Status and Future Challenges. In and J. Verron J. Tintoré A. Pascual E. P. Chassignet (Ed.),
New Frontiers in Operational Oceanography (pp. 289–305). Tallahassee, FL: GODAE OceanView.
Abstract: This chapter focuses on ocean circulation models used in operational oceanography, physical oceanography and climate science. Ocean circulation models area particular branch of ocean numerical modeling that focuses on the representation of ocean physical properties over spatial scales ranging from the global scale to less than a kilometer and time scales ranging from hours to decades. As such, they are an essential build-ing block for operational oceanography systems and their design receives a lot of attention from operational and research centers.
Magaldi, M. G., Özgökmen, T. M., Griffa, A., Chassignet, E. P., Iskandarani, M., & Peters, H. (2008). Turbulent flow regimes behind a coastal cape in a stratified and rotating environment.
Ocean Modelling, 25(1-2), 65–82.
Morey, S. L. (2002). The spring transition from horizontal to vertical thermal stratification on a midlatitude continental shelf.
J. Geophys. Res., 107(C8).
Morey, S. L., & Dukhovskoy, D. S. (2013). A downscaling method for simulating deep current interactions with topography – Application to the Sigsbee Escarpment.
Ocean Modelling, 69, 50–63.