Joseph Smagorinsky

Bio: Joseph Smagorinsky is an academic researcher from Princeton University. The author has contributed to research in topics: Primitive equations & Troposphere. The author has an hindex of 14, co-authored 17 publications receiving 14142 citations. Previous affiliations of Joseph Smagorinsky include Environmental Science Services Administration & Institute for Advanced Study.

General circulation experiments with the primitive equations

Abstract: An extended period numerical integration of a baroclinic primitive equation model has been made for the simulation and the study of the dynamics of the atmosphere's general circulation. The solution corresponding to external gravitational propagation is filtered by requiring the vertically integrated divergence to vanish identically. The vertical structure permits as dependent variables the horizontal wind at two internal levels and a single temperature, with the static stability entering as a parameter. The incoming radiation is a function of latitude only corresponding to the annual mean, and the outgoing radiation is taken to be a function of the local temperature. With the requirement for thermal equilibrium, the domain mean temperature is specified as a parameter. The role of condensation is taken into account only as it effectively reduces the static stability. All other external sources and sinks of heat are assumed to balance each other locally, and are thus omitted. The kinematics are th...

Simulated climatology of a general circulation model with a hydrologic cycle

Abstract: The results of a numerical time integration of a hemispheric general circulation model of the atmosphere with moist processes and a uniform earth's surface has already been published by Manabe, Smagorinsky, and Strickler. In this study, the integration is repeated after halving the midlatitude grid size from approximately 500 to 250 km. This increase in the resolution of the horizontal finite differences markedly improves the features of the model atmosphere. For example, the system of fronts and the associated cyclone families in the high resolution atmosphere is much more realistic than that in the low resolution atmosphere. Furthermore, the general magnitude and the spectral distribution of eddy kinetic energy are in better agreement with the actual atmosphere as a result of the improvement in resolution. In order to explain these improvements, an extensive analysis of the energetics of both the low and high resolution atmospheric models is carried out. It is shown that these improvements are ...

Numerical Results from a Nine-Level General Circulation Model of the ATMOSPHERE1

Abstract: The “primitive equations of motion” are adopted for this study. The nine levels of the model are distributed so as to resolve surface boundary layer fluxes as well as radiative transfer by ozone, carbon dioxide, and water vapor. The lower boundary is a kinematically uniform land surface without any heat capacity. The stabilizing effect of moist convection is implicitly incorporated into the model by requiring an adjustment of the lapse rate whenever it exceeds the moist adiabatic value. The numerical integrations are performed for the mean annual conditions over a hemisphere starting with an isothermal atmosphere at rest. The spatial distribution of gaseous absorbers is assumed to have the annual mean value of the actual atmosphere and to be constant with time. A quasi-equilibrium is attained about which a cyclic energy variation occurs with an irregular period of about 2 weeks. The dominant wave number of the meridional component of the wind is 5 to 6 in the troposphere but is reduced to about 3...

The dynamical influence of large‐scale heat sources and sinks on the quasi‐stationary mean motions of the atmosphere

Abstract: The classical convective theory for the influence of large-scale zonally asymmetric heating and cooling on the quasi-stationary mean motions is re-examined and shown to be inconsistent with observation. Scale considerations indicate that the classical mechanism cannot be applied to heating on a scale larger than that responsible for sea-breezes. In analogy with the theory of orographically produced stationary disturbances, heat sources and sinks of planetary dimensions produce quasi-static and quasi-geostrophic perturbations. A thermally active model is constructed in which all motions except those of large scale are filtered out a priori. The small perturbations produced by an idealized large-scale distribution of heat sources and sinks acting on a horizontally uniform baroclinic zonal current confined between rigid horizontal planes are then studied. The influence of surface friction is taken into account. Calculations are also performed on a model with no friction and on the assumption that entropy is advected horizontally. The results generally show good agreement with those for the more general model. When the stratosphere is taken into account, the disturbances in the upper troposphere are reduced in magnitude. The conditions for resonance in the absence of friction are determined. The passage through the resonant state is shown to modify significantly the character of the motions. The observed normals are discussed at some length in light of the theoretical computations. It is found that the magnitude and spatial distribution of the model disturbances display substantial correspondence with observation when allowance is made for orographic effects. Finally a method is suggested for quantitatively deducing the three-dimensional field of non-adiabatic heating and cooling from the observed normal maps by dynamical considerations.

A dynamic subgrid‐scale eddy viscosity model

Abstract: One major drawback of the eddy viscosity subgrid‐scale stress models used in large‐eddy simulations is their inability to represent correctly with a single universal constant different turbulent fields in rotating or sheared flows, near solid walls, or in transitional regimes. In the present work a new eddy viscosity model is presented which alleviates many of these drawbacks. The model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model is based on an algebraic identity between the subgrid‐scale stresses at two different filtered levels and the resolved turbulent stresses. The subgrid‐scale stresses obtained using the proposed model vanish in laminar flow and at a solid boundary, and have the correct asymptotic behavior in the near‐wall region of a turbulent boundary layer. The results of large‐eddy simulations of transitional and turbulent channel flow that use the proposed model are in good agreement with the direct simulation data.

The Kain–Fritsch Convective Parameterization: An Update

Abstract: Numerous modifications to the Kain‐Fritsch convective parameterization have been implemented over the last decade. These modifications are described, and the motivating factors for the changes are discussed. Most changes were inspired by feedback from users of the scheme (primarily numerical modelers) and interpreters of the model output (mainly operational forecasters). The specific formulation of the modifications evolved from an effort to produce desired effects in numerical weather prediction while also rendering the scheme more faithful to observations and cloud-resolving modeling studies.

A tensorial approach to computational continuum mechanics using object-oriented techniques

Abstract: In this article the principles of the field operation and manipulation (FOAM) C++ class library for continuum mechanics are outlined. Our intention is to make it as easy as possible to develop reliable and efficient computational continuum-mechanics codes: this is achieved by making the top-level syntax of the code as close as possible to conventional mathematical notation for tensors and partial differential equations. Object-orientation techniques enable the creation of data types that closely mimic those of continuum mechanics, and the operator overloading possible in C++ allows normal mathematical symbols to be used for the basic operations. As an example, the implementation of various types of turbulence modeling in a FOAM computational-fluid-dynamics code is discussed, and calculations performed on a standard test case, that of flow around a square prism, are presented. To demonstrate the flexibility of the FOAM library, codes for solving structures and magnetohydrodynamics are also presented with appropriate test case results given. © 1998 American Institute of Physics.

Semi-Supervised Learning

Abstract: In the field of machine learning, semi-supervised learning (SSL) occupies the middle ground, between supervised learning (in which all training examples are labeled) and unsupervised learning (in which no label data are given). Interest in SSL has increased in recent years, particularly because of application domains in which unlabeled data are plentiful, such as images, text, and bioinformatics. This first comprehensive overview of SSL presents state-of-the-art algorithms, a taxonomy of the field, selected applications, benchmark experiments, and perspectives on ongoing and future research. Semi-Supervised Learning first presents the key assumptions and ideas underlying the field: smoothness, cluster or low-density separation, manifold structure, and transduction. The core of the book is the presentation of SSL methods, organized according to algorithmic strategies. After an examination of generative models, the book describes algorithms that implement the low-density separation assumption, graph-based methods, and algorithms that perform two-step learning. The book then discusses SSL applications and offers guidelines for SSL practitioners by analyzing the results of extensive benchmark experiments. Finally, the book looks at interesting directions for SSL research. The book closes with a discussion of the relationship between semi-supervised learning and transduction. Adaptive Computation and Machine Learning series

Oceanic vertical mixing: A review and a model with a nonlocal boundary layer parameterization

Abstract: If model parameterizations of unresolved physics, such as the variety of upper ocean mixing processes, are to hold over the large range of time and space scales of importance to climate, they must be strongly physically based. Observations, theories, and models of oceanic vertical mixing are surveyed. Two distinct regimes are identified: ocean mixing in the boundary layer near the surface under a variety of surface forcing conditions (stabilizing, destabilizing, and wind driven), and mixing in the ocean interior due to internal waves, shear instability, and double diffusion (arising from the different molecular diffusion rates of heat and salt). Mixing schemes commonly applied to the upper ocean are shown not to contain some potentially important boundary layer physics. Therefore a new parameterization of oceanic boundary layer mixing is developed to accommodate some of this physics. It includes a scheme for determining the boundary layer depth h, where the turbulent contribution to the vertical shear of a bulk Richardson number is parameterized. Expressions for diffusivity and nonlocal transport throughout the boundary layer are given. The diffusivity is formulated to agree with similarity theory of turbulence in the surface layer and is subject to the conditions that both it and its vertical gradient match the interior values at h. This nonlocal “K profile parameterization” (KPP) is then verified and compared to alternatives, including its atmospheric counterparts. Its most important feature is shown to be the capability of the boundary layer to penetrate well into a stable thermocline in both convective and wind-driven situations. The diffusivities of the aforementioned three interior mixing processes are modeled as constants, functions of a gradient Richardson number (a measure of the relative importance of stratification to destabilizing shear), and functions of the double-diffusion density ratio, Rρ. Oceanic simulations of convective penetration, wind deepening, and diurnal cycling are used to determine appropriate values for various model parameters as weak functions of vertical resolution. Annual cycle simulations at ocean weather station Papa for 1961 and 1969–1974 are used to test the complete suite of parameterizations. Model and observed temperatures at all depths are shown to agree very well into September, after which systematic advective cooling in the ocean produces expected differences. It is argued that this cooling and a steady salt advection into the model are needed to balance the net annual surface heating and freshwater input. With these advections, good multiyear simulations of temperature and salinity can be achieved. These results and KPP simulations of the diurnal cycle at the Long-Term Upper Ocean Study (LOTUS) site are compared with the results of other models. It is demonstrated that the KPP model exchanges properties between the mixed layer and thermocline in a manner consistent with observations, and at least as well or better than alternatives.