Publications

Export 110 results:
Sort by: Author [ Title  (Asc)] Type Year
A B C D E F G H I J K L M N O P Q R S T U V W X Y Z 
I
Lubin, D, Chen B, Bromwich DH, Somerville RCJ, Lee WH, Hines KM.  1998.  The impact of Antarctic cloud radiative properties on a GCM climate simulation. Journal of Climate. 11:447-462.   10.1175/1520-0442(1998)011<0447:tioacr>2.0.co;2   AbstractWebsite

A sensitivity study to evaluate the impact upon regional and hemispheric climate caused by changing the optical properties of clouds over the Antarctic continent is conducted with the NCAR Community Model version 2 (CCM2). Sensitivity runs are performed in which radiation interacts with ice clouds with particle sizes of 10 and 40 mu m rather than with the standard 10-mu m water clouds. The experiments are carried out for perpetual January conditions with the diurnal cycle considered. The effects of these cloud changes on the Antarctic radiation budget are examined by considering cloud forcing at the top of the atmosphere and net radiation at the surface. Changes of the cloud radiative properties to those of 10-mu m ice clouds over Antarctica have significant Impacts on regional climate: temperature increases throughout the Antarctic troposphere by 1 degrees-2 degrees C and total cloud fraction over Antarctica is smaller than that of the control at low levels but is larger than that of the control in the mid- to upper troposphere. As a result of Antarctic warming and changes in the north-south temperature gradient, the drainage flows at the surface as well as the meridional mass circulation are weakened. Similarly, the circumpolar trough weakens significantly by 4-8 hPa and moves northward by about 4 degrees-5 degrees latitude. This regional mass field adjustment halves the strength of the simulated surface westerly winds. As a result of indirect thermodynamic and dynamic effects, significant changes are observed in the zonal mean circulation and eddies in the middle latitudes. In fact, the simulated impacts of the Antarctic cloud radiative alteration are not confined to the Southern Hemisphere. The meridional mean mass flux, zonal wind, and latent heat release exhibit statistically significant changes in the Tropics and even extratropics of the Northern Hemisphere. The simulation with radiative properties of 40-mu m ice clouds produces colder surface temperatures over Antarctica by up to 3 degrees C compared to the control. Otherwise, the results of the 40-mu m ice cloud simulation are similar to those of the 10-mu m ice cloud simulation.

Yang, Y, Russell LM, Xu L, Lou SJ, Lamjiri MA, Somerville RCJ, Miller AJ, Cayan DR, DeFlorio MJ, Ghan SJ, Liu Y, Singh B, Wang HL, Yoon JH, Rasch PJ.  2016.  Impacts of ENSO events on cloud radiative effects in preindustrial conditions: Changes in cloud fraction and their dependence on interactive aerosol emissions and concentrations. Journal of Geophysical Research-Atmospheres. 121:6321-6335.   10.1002/2015jd024503   AbstractWebsite

We use three 150 year preindustrial simulations of the Community Earth System Model to quantify the impacts of El Nino-Southern Oscillation (ENSO) events on shortwave and longwave cloud radiative effects (CRESW and CRELW). Compared to recent observations from the Clouds and the Earth's Radiant Energy System data set, the model simulation successfully reproduces larger variations of CRESW and CRELW over the tropics. The ENSO cycle is found to dominate interannual variations of cloud radiative effects. Simulated cooling (warming) effects from CRESW (CRELW) are strongest over the tropical western and central Pacific Ocean during warm ENSO events, with the largest difference between 20 and 60 W m(-2), with weaker effects of 10-40 W m(-2) over Indonesian regions and the subtropical Pacific Ocean. Sensitivity tests show that variations of cloud radiative effects are mainly driven by ENSO-related changes in cloud fraction. The variations in midlevel and high cloud fractions each account for approximately 20-50% of the interannual variations of CRESW over the tropics and almost all of the variations of CRELW between 60 degrees S and 60 degrees N. The variation of low cloud fraction contributes to most of the variations of CRESW over the midlatitude oceans. Variations in natural aerosol concentrations explained 10-30% of the variations of both CRESW and CRELW over the tropical Pacific, Indonesian regions, and the tropical Indian Ocean. Changes in natural aerosol emissions and concentrations enhance 3-5% and 1-3% of the variations of cloud radiative effects averaged over the tropics.

Iacobellis, SF, Somerville RCJ.  2000.  Implications of microphysics for cloud-radiation parameterizations: Lessons from TOGA COARE. Journal of the Atmospheric Sciences. 57:161-183.   10.1175/1520-0469(2000)057<0161:iomfcr>2.0.co;2   AbstractWebsite

A single-column model (SCM) and observational data collected during TOGA COARE were used to investigate the sensitivity of model-produced cloud properties and radiative fluxes to the representation of cloud microphysics in the cloud-radiation parameterizations. Four 78-day SCM numerical experiments were conducted for the atmospheric column overlying the COARE Intensive Flux Array. Each SCM experiment used a different cloud-radiation parameterization with a different representation of cloud microphysics. All the SCM experiments successfully reproduced most of the observed temporal variability in precipitation, cloud fraction, shortwave and longwave cloud forcing, and downwelling surface shortwave flux. The magnitude and temporal variability of the downward surface longwave flux was overestimated by all the SCM experiments. This bins is probably due to clouds forming too low in the model atmosphere. Time-averaged model results were used to examine the sensitivity of model performance to the differences between the four cloud-radiation parameterization packages. The SCM versions that calculated cloud amount as a function of cloud liquid water, instead of using a relative humidity-based cloud scheme, produced smaller amounts of both low and deep convective clouds. Additionally, larger high (cirrus) cloud emissivities were obtained with interactive cloud liquid water schemes than with the relative humidity-based scheme. Surprisingly. calculating cloud optical properties as a function of cloud liquid water amount, instead of parameterizing them based on temperature, humidity, and pressure, resulted in relatively little change in radiative fluxes. However. model radiative fluxes were sensitive to the specification of the effective cloud droplet radius. Optically thicker low clouds and optically thinner high clouds were produced when an interactive effective cloud droplet radius scheme was used instead of specifying a constant value. Comparison of model results to both surface and satellite observations revealed that model experiments that calculated cloud properties as a function of cloud liquid water produced more realistic cloud amounts and radiative fluxes. The most realistic vertical distribution of clouds was obtained from the SCM experiment that included the most complete representation of cloud microphysics. Due to the limitations of SCMs. the above conclusions are model dependent and need to be tested in a general circulation model.

Berque, J, Lubin D, Somerville RCJ.  2004.  Infrared radiative properties of the Antarctic plateau from AVHRR data. Part I: Effect of the snow surface. Journal of Applied Meteorology. 43:350-362.   10.1175/1520-0450(2004)043<0350:irpota>2.0.co;2   AbstractWebsite

The effective scene temperature, or "brightness temperature," measured in channel 3 (3.5-3.9 m m) of the Advanced Very High Resolution Radiometer (AVHRR) is shown to be sensitive, in principle, to the effective particle size of snow grains on the Antarctic plateau, over the range of snow grain sizes reported in field studies. In conjunction with a discrete ordinate method radiative transfer model that couples the polar atmosphere with a scattering and absorbing snowpack, the thermal infrared channels of the AVHRR instrument can, therefore, be used to estimate effective grain size at the snow surface over Antarctica. This is subject to uncertainties related to the modeled top-of-atmosphere bidirectional reflectance distribution function resulting from the possible presence of sastrugi and to lack of complete knowledge of snow crystal shapes and habits as they influence the scattering phase function. However, when applied to NOAA-11 and NOAA-12 AVHRR data from 1992, the snow grain effective radii of order 50 mum are retrieved, consistent with field observations, with no apparent discontinuity between two spacecraft having different viewing geometries. Retrieved snow grain effective radii are 10-20-mum larger when the snow grains are modeled as hexagonal solid columns rather than as spheres with a Henyey-Greenstein phase function. Despite the above-mentioned uncertainties, the retrievals are consistent enough that one should be able to monitor climatically significant changes in surface snow grain size due to major precipitation events. It is also shown that a realistic representation of the surface snow grain size is critical when retrieving the optical depth and effective particle radius of clouds for the optically thin clouds most frequently encountered over the Antarctic plateau.

Xu, L, Pierce DW, Russell LM, Miller AJ, Somerville RCJ, Twohy CH, Ghan SJ, Singh B, Yoon J-H, Rasch PJ.  2015.  Interannual to decadal climate variability of sea salt aerosols in the coupled climate model CESM1.0. Journal of Geophysical Research: Atmospheres. :2014JD022888.   10.1002/2014JD022888   AbstractWebsite

This study examines multi-year climate variability associated with sea salt aerosols and their contribution to the variability of shortwave cloud forcing (SWCF) using a 150-year simulation for pre-industrial conditions of the Community Earth System Model version 1.0 (CESM1). The results suggest that changes in sea salt and related cloud and radiative properties on interannual timescales are dominated by the ENSO cycle. Sea salt variability on longer (interdecadal) timescales is associated with low-frequency variability in the Pacific Ocean similar to the interdecadal Pacific Oscillation (IPO), but does not show a statistically significant spectral peak. A multivariate regression suggests that sea salt aerosol variability may contribute to SWCF variability in the tropical Pacific, explaining up to 20-30% of the variance in that region. Elsewhere, there is only a small sea salt aerosol influence on SWCF through modifying cloud droplet number and liquid water path that contributes to the change of cloud effective radius and cloud optical depth (and hence cloud albedo), producing a multi-year aerosol-cloud-wind interaction.

Xie, SC, Xu KM, Cederwall RT, Bechtold P, Delgenio AD, Klein SA, Cripe DG, Ghan SJ, Gregory D, Iacobellis SF, Krueger SK, Lohmann U, Petch JC, Randall DA, Rotstayn LD, Somerville RCJ, Sud YC, Von Salzen K, Walker GK, Wolf A, Yio JJ, Zhang GJ, Zhang MG.  2002.  Intercomparison and evaluation of cumulus parametrizations under summertime midlatitude continental conditions. Quarterly Journal of the Royal Meteorological Society. 128:1095-1135.   10.1256/003590002320373229   AbstractWebsite

This study reports the Single-Column Model (SCM) part of the Atmospheric Radiation Measurement (ARM)/the Global Energy and Water Cycle Experiment (GEWEX) Cloud System Study (GCSS) joint SCM and Cloud-Resolving Model (CRM) Case 3 intercomparison study, with a focus on evaluation Of Cumulus parametrizations used in SCMs. Fifteen SCMs are evaluated under summertime midlatitude continental conditions using data collected at the ARM Southern Great Plains site during the summer 1997 Intensive Observing Period. Results from ten CRMs are also used to diagnose problems in the SCMs. It is shown that most SCMs can generally capture well the convective events that were well-developed within the SCM domain, while most of them have difficulties in simulating the occurrence of those convective events that only occurred within a small part of the domain. All models significantly underestimate the surface stratiform precipitation. A third of them produce large errors in surface precipitation and thermodynamic structures. Deficiencies in convective triggering mechanisms are thought to be one of the major reasons. Using a triggering mechanism that is based on the vertical integral of parcel buoyant energy without additional appropriate constraints results in overactive convection, which in turn leads to large systematic warm/dry biases in the troposphere. It is also shown that a non-penetrative convection scheme can underestimate the depth of instability for midlatitude convection, which leads to large systematic cold/moist biases in the troposphere. SCMs agree well quantitatively with CRMs in the updraught mass fluxes, while most models significantly underestimate the downdraught mass fluxes. Neglect of mesoscale updraught and downdraught mass fluxes in the SCMs contributes considerably to the discrepancies between the SCMs and the CRMs. In addition, uncertainties in the diagnosed mass fluxes in the CRMs and deficiencies with cumulus parametrizations are not negligible. Similar results are obtained in the sensitivity tests when different forcing approaches are used. Finally. sensitivity tests from an SCM indicate that its simulations can be greatly improved when its triggering mechanism and closure assumption are improved.

Zhang, C, Wang M, Morrison H, Somerville RCJ, Zhang K, Liu X, Li J-LF.  2014.  Investigating ice nucleation in cirrus clouds with an aerosol-enabled Multiscale Modeling Framework. Journal of Advances in Modeling Earth Systems. 6:998-1015.   10.1002/2014MS000343   Abstract

In this study, an aerosol-dependent ice nucleation scheme has been implemented in an aerosol-enabled Multiscale Modeling Framework (PNNL MMF) to study ice formation in upper troposphere cirrus clouds through both homogeneous and heterogeneous nucleation. The MMF model represents cloud scale processes by embedding a cloud-resolving model (CRM) within each vertical column of a GCM grid. By explicitly linking ice nucleation to aerosol number concentration, CRM-scale temperature, relative humidity and vertical velocity, the new MMF model simulates the persistent high ice supersaturation and low ice number concentration (10–100/L) at cirrus temperatures. The new model simulates the observed shift of the ice supersaturation PDF toward higher values at low temperatures following the homogeneous nucleation threshold. The MMF model predicts a higher frequency of midlatitude supersaturation in the Southern Hemisphere and winter hemisphere, which is consistent with previous satellite and in situ observations. It is shown that compared to a conventional GCM, the MMF is a more powerful model to simulate parameters that evolve over short time scales such as supersaturation. Sensitivity tests suggest that the simulated global distribution of ice clouds is sensitive to the ice nucleation scheme and the distribution of sulfate and dust aerosols. Simulations are also performed to test empirical parameters related to auto-conversion of ice crystals to snow. Results show that with a value of 250 μm for the critical diameter, Dcs, that distinguishes ice crystals from snow, the model can produce good agreement with the satellite-retrieved products in terms of cloud ice water path and ice water content, while the total ice water is not sensitive to the specification of Dcs value.

Zhao, Z, Kooperman GJ, Pritchard MS, Russell LM, Somerville RCJ.  2014.  Investigating impacts of forest fires in Alaska and western Canada on regional weather over the northeastern United States using CAM5 global simulations to constrain transport to a WRF-Chem regional domain. Journal of Geophysical Research: Atmospheres. 119:2013JD020973.   10.1002/2013JD020973   AbstractWebsite

An aerosol-enabled globally driven regional modeling system has been developed by coupling the National Center for Atmospheric Research's Community Atmosphere Model version 5 (CAM5) with the Weather Research and Forecasting model with chemistry (WRF-Chem). In this modeling system, aerosol-enabled CAM5, a state-of-the-art global climate model is downscaled to provide coherent meteorological and chemical boundary conditions for regional WRF-Chem simulations. Aerosol particle emissions originating outside the WRF-Chem domain can be a potentially important nonlocal aerosol source. As a test case, the potential impacts of nonlocal forest fire aerosols on regional precipitation and radiation were investigated over the northeastern United States during the summer of 2004. During this period, forest fires in Alaska and western Canada lofted aerosol particles into the midtroposphere, which were advected across the United States. WRF-Chem simulations that included nonlocal biomass burning aerosols had domain-mean aerosol optical depths that were nearly three times higher than those without, which reduced peak downwelling domain-mean shortwave radiation at the surface by ~25 W m−2. In this classic twin experiment design, adding nonlocal fire plume led to near-surface cooling and changes in cloud vertical distribution, while variations in domain-mean cloud liquid water path were negligible. The higher aerosol concentrations in the simulation with the fire plume resulted in a ~10% reduction in domain-mean precipitation coincident with an ~8% decrease in domain-mean CAPE. A suite of simulations was also conducted to explore sensitivities of meteorological feedbacks to the ratio of black carbon to total plume aerosols, as well as to overall plume concentrations. Results from this ensemble revealed that plume-induced near-surface cooling and CAPE reduction occur in a wide range of conditions. The response of moist convection was very complex because of strong thermodynamic internal variability.

Somerville, RCJ.  2014.  Is learning about climate change like having a colonoscopy? Earth's Future. 2:119-121.: Wiley Periodicals, Inc.   10.1002/2013EF000169   Abstract

Many people avoid having valuable medical tests from fear of the results
People resist learning about climate change, fearing unpleasant consequences
Research suggests addressing these concerns early, aids in communication

L
Engquist, BE, Osher S, Somerville RCJ.  1985.  Large-Scale Computations in Fluid Mechanics. Lectures in Applied Mathematics. :779.: American Mathematical Society AbstractWebsite
n/a
Somerville, RCJ, Jouzel J.  2007.  Le groupe intergouvernemental d'experts sur l'evolution du climat: le consensus a l'échelle planétaire. Comprendre le changement climatique. ( André JC, Fellous JL, Gautier C, Eds.).:27-44., Paris: O. Jacob Abstract
n/a
Roads, JO, Somerville RCJ.  1984.  Linear Predictability: Effects of Stationary Forcing. Aip Conference Proceedings. :557-570. AbstractWebsite
n/a
M
Razafimpanilo, H, Frouin R, Iacobellis SF, Somerville RCJ.  1995.  Methodology for estimating burned area from AVHRR reflectance data. Remote Sensing of Environment. 54:273-289.   10.1016/0034-4257(95)00154-9   AbstractWebsite

Two methods are described to determine burned area from Advanced Very High Resolution Radiometer (AVHRR) data. The first method, or the ''linear method,'' employs Channel 2 reflectance, R(2), and is based on the nearly linear relationship between the fraction of pixel burned, P, and R(2). The second method, or the ''nonlinear method,'' employs the Normalized Difference Vegetation Index (NDVI) derived from Channels 1 and 2 reflectances, and is based on the nonlinear relationship P=f(NDVI), a polynomial of order 2 in NDVI. The coefficients of the polynomial are parameterized as a function of the NDVI of the background before the fire event. Radiative transfer simulations indicate that the linear method, unlike the nonlinear method, must be applied to top-of-atmosphere reflectances that have been corrected for atmospheric influence. Sensitivity studies suggest that the methods are subject to some limitations. To avoid discontinuity problems, the original background (just before the fire) must be characterized by a Channel 2 reflectance above 0.07 and by a positive NDVI. To separate the useful signal from atmospheric effects, the fire scar must occupy at least 20% and 12% of the pixel area in the case of savanna and green vegetation (e.g., forest), respectively When applied to uniform pixels, the mean relative error on the fraction of area burned is about 20% for the linear method and 10% for the nonlinear method. The linear method gives better results for nonuniform pixels, but neither method can be used when the pixel contains low reflectance backgrounds (e.g., water).

Xu, KM, Zhang MH, Eitzen MA, Ghan SJ, Klein SA, Wu XQ, Xie SC, Branson M, Delgenio AD, Iacobellis SF, Khairoutdinov M, Lin WY, Lohmann U, Randall DA, Somerville RCJ, Sud YC, Walker GK, Wolf A, Yio JJ, Zhang JH.  2005.  Modeling springtime shallow frontal clouds with cloud-resolving and single-column models. Journal of Geophysical Research-Atmospheres. 110   10.1029/2004jd005153   AbstractWebsite

This modeling study compares the performance of eight single-column models (SCMs) and four cloud-resolving models (CRMs) in simulating shallow frontal cloud systems observed during a short period of the March 2000 Atmospheric Radiation Measurement (ARM) intensive operational period. Except for the passage of a cold front at the beginning of this period, frontal cloud systems are under the influence of an upper tropospheric ridge and are driven by a persistent frontogenesis over the Southern Great Plains and moisture transport from the northwestern part of the Gulf of Mexico. This study emphasizes quantitative comparisons among the model simulations and with the ARM data, focusing on a 27-hour period when only shallow frontal clouds were observed. All CRMs and SCMs simulate clouds in the observed shallow cloud layer. Most SCMs also produce clouds in the middle and upper troposphere, while none of the CRMs produce any clouds there. One possible cause for this is the decoupling between cloud condensate and cloud fraction in nearly all SCM parameterizations. Another possible cause is the weak upper tropospheric subsidence that has been averaged over both descending and ascending regions. Significantly different cloud amounts and cloud microphysical properties are found in the model simulations. All CRMs and most SCMs underestimate shallow clouds in the lowest 125 hPa near the surface, but most SCMs overestimate the cloud amount above this layer. These results are related to the detailed formulations of cloud microphysical processes and fractional cloud parameterizations in the SCMs, and possibly to the dynamical framework and two-dimensional configuration of the CRMs. Although two of the CRMs with anelastic dynamical frameworks simulate the shallow frontal clouds much better than the SCMs, the CRMs do not necessarily perform much better than the SCMs for the entire period when deep and shallow frontal clouds are present.

N
Hathaway, DH, Somerville RCJ.  1986.  Nonlinear Interactions between Convection, Rotation and Flows with Vertical Shear. Journal of Fluid Mechanics. 164:91-&.   10.1017/s0022112086002483   AbstractWebsite

A three-dimensional and time-dependent numerical model is used to study the nonlinear interactions between thermal convective motions, rotation, and imposed flows with vertical shear. All cases have Rayleigh numbers of 104 and Prandtl numbers of 1.0. Rotating cases have Taylor numbers of 104.For the non-rotating cases, the effects of the shear on the convection produce longitudinal rolls aligned with the shear flow and a downgradient flux of momentum. The interaction between the convection and the shear flow decreases the shear in the interior of the fluid layer while adding kinetic energy to the convective motions. For unit Prandtl number the dimensionless flux of momentum is equal to the dimensionless flux of heat.For rotating cases with vertical rotation vectors, the shear flow favours rolls aligned with the shear and produces a downgradient flux of momentum. However, the Coriolis force turns the flow induced by the convection to produce a more complicated shear that changes direction with height. As in the non-rotating cases, the convective motions become more energetic by extracting energy from the mean flow. For Richardson numbers larger than about − 1.0, the dominant source of eddy kinetic energy is the shear flow rather than buoyancy.For rotating cases with tilted rotation vectors the results depend upon the direction of the shear. For weak shear, convective rolls aligned with the rotation vector are favoured. When the shear flow is directed to the east along the top, the rolls become broader and the convection weaker. For large shear in this direction, the convective motions are quenched by the competition between the shear flow and the tilted rotation vector. When the shear flow is directed to the west along the top, strong shear produces rolls aligned with the shear. The heat and momentum fluxes become large and can exceed those found in the absence of a tilted rotation vector. Countergradient fluxes of momentum can also be produced.

Somerville, RCJ.  1967.  A Nonlinear Spectral Model of Convection in a Fluid Unevenly Heated from Below. Journal of the Atmospheric Sciences. 24:665-676.: American Meteorological Society   10.1175/1520-0469(1967)024<0665:ansmoc>2.0.co;2   AbstractWebsite

A two-dimensional form of the Boussinesq equations is integrated numerically for the case of a rectangular channel with a temperature gradient maintained along the bottom. The side walls are insulating, the top wall has a constant temperature, and the velocity obeys free boundary conditions on all four walls. The fields of stream function and temperature departure are represented by truncated double Fourier series, and integration of the initial-value problem for the spectral amplitudes results in steady states which agree qualitatively with those of previous experimental and theoretical investigations.Calculations are presented at two levels of truncation (wave numbers 2 and 3) for a wide range of Prandtl numbers and a moderate range of horizontal Rayleigh numbers and top temperatures. For sufficiently large gravitational stability, a single asymmetric convection cell develops. Its intensity and asymmetry increase markedly with increasing horizontal Rayleigh number, decrease with increasing top temperature, and respond very slightly to changes in Prandtl number. As the top temperature is decreased below the temperature of the warm side of the bottom, however, the possibility is indicated that the single cell may be modified by a Bénard-like multi-cellular structure.

Iacobellis, SF, Frouin R, Razafimpanilo H, Somerville RCJ, Piper SC.  1994.  North African savanna fires and atmospheric carbon dioxide. Journal of Geophysical Research-Atmospheres. 99:8321-8334.   10.1029/93jd03339   AbstractWebsite

The effect of north African savanna fires on atmospheric CO2 is investigated using a tracer transport model. The model uses winds from operational numerical weather prediction analyses and provides CO2 Concentrations as a function of space and time. After a spin-up period of several years, biomass-burning sources are added, and model experiments are run for an additional year, utilizing various estimates of CO2 sources. The various model experiments show that biomass burning in the north African savannas significantly affects CO2 concentrations in South America. The effect is more pronounced during the period from January through March, when biomass burning in South America is almost nonexistent. During this period, atmospheric CO2 concentrations in parts of South America typically may increase by 0.5 to 0.75 ppm at 970 mbar, the average pressure of the lowest model layer. These figures are above the probable uncertainty level, as model runs with biomass-burning sources estimated from independent studies using distinct data sets and techniques indicate. From May through September, when severe biomass burning occurs in South America, the effect of north African savanna fires over South America has become generally small at 970 mbar, but north of the equator it may be of the same magnitude or larger than the effect of South American fires. The CO2 concentration increase in the extreme northern and southern portions of South America, however, is mostly due to southern African fires, whose effect may be 2-3 times larger than the effect of South American fires at 970 mbar. Even in the central part of the continent, where local biomass-burning emissions are maximum, southern African fires contribute to at least 15% of the CO2 concentration increase at 970 mbar. At higher levels in the atmosphere, less CO2 emitted by north African savanna fires reaches South America, and at 100 mbar no significant amount of CO2 is transported across the Atlantic Ocean. The vertical structure of the CO2 concentration increase due to biomass burning differs substantially, depending on whether sources are local or remote. A prominent maximum Of CO2 concentration increase in the lower layers characterizes the effect of local sources, whereas a more homogenous profile of CO2 concentration increase characterizes the effect of remote sources. The results demonstrate the strong remote effects of African biomass burning which, owing to the general circulation of the atmosphere, are felt as far away as South America.

Hathaway, DH, Somerville RCJ.  1985.  Numerical simulation in three space dimensions of time-dependent thermal convection in a rotating fluid. Lectures in Applied Mathematics. 22:309-319. Abstract

Three-dimensional time-dependent convection in a plane layer of fluid, uniformly heated from below and subject to vertical shear and to rotation about an axis tilted from the vertical, was simulated by the numerical solution of the Boussinesq equations, including all Coriolis terms. Rotation about a vertical axis produces smaller convection cells with diminished heat fluxes and considerable vorticity. When the rotation axis is tilted from the vertical to represent tropical latitudes, the convection cells become elongated in a N-S direction. Imposed flows with constant vertical shear produce convective rolls aligned with the mean flow. When the rotation vector is tilted from the vertical, the competing effects due to rotation and shear can stabilize the convective motions.

Somerville, RCJ, Galchen T.  1979.  A Numerical Simulation of Convection with Mean Vertical Motion. Journal of the Atmospheric Sciences. 36:805-815.   10.1175/1520-0469(1979)036<0805:nsocwm>2.0.co;2   AbstractWebsite

The flow in a convectively unstable layer of fluid may be strongly influenced by large-scale ascent or descent. We consider cellular convection between horizontal surfaces on which vertical velocity is maintained at a constant value. Using an efficient numerical model to simulate the evolution of the convection in three space dimensions and time, we investigate the effect of the imposed vertical velocity on the flow.For moderately supercritical values of the Rayleigh number and for Prandtl numbers near unity, convection is known to occur in the form of steady rolls if the specified mean vertical motion is zero, i.e., in the case of the conventional Bénard problem for a Boussinesq fluid. Our model also produces rolls under these circumstances. For sufficiently large values of the imposed vertical velocity, however, the numerically simulated rolls are replaced by polygonal cells in which the direction of flow depends on whether ascent or descent is prescribed at the boundaries, in accordance with recent theoretical and laboratory results of R. Krishnamurti. We have also investigated the dependence of the convection on the Rayleigh and Prandtl numbers within limited ranges of these parameters, and we discuss several aspects of agreement and disagreement among analytical theory, laboratory experiment and numerical simulation.

Somerville, RCJ, Lipps FB.  1973.  A Numerical Study in Three Space Dimensions of Bénard Convection in a Rotating Fluid. Journal of the Atmospheric Sciences. 30:590-596.: American Meteorological Society   10.1175/1520-0469(1973)030<0590:ansits>2.0.co;2   AbstractWebsite

The primitive, nonlinear, Boussinesq equations of motion, continuity and thermodynamic energy are integrated numerically in three space dimensions and time to study convection driven by unstable vertical density gradients and subject to Coriolis forces. Parameter values are chosen to permit quantitative comparison with data from laboratory experiments for rotating Bénard convection in water. The model realistically simulates the structure of the convection cells, their horizontal scale, and the mean vertical heat transport. The experimentally observed phenomenon of a non-monotone dependence of heat transport on rotation rate is reproduced and shown to be a consequence of the rotational constraint on the wavelength of the cells.

Galchen, T, Somerville RCJ.  1975.  Numerical-Solution of Navier-Stokes Equations with Topography. Journal of Computational Physics. 17:276-310.   10.1016/0021-9991(75)90054-6   AbstractWebsite

A finite difference scheme for solving the equations of fluid motion in a generalized coordinate system has been constructed. The scheme conserves mass and all the first integral moments of the motion. The scheme also advectively “almost conserves” second moments, in that the magnitude of implicit numerical smoothing is typically about an order smaller than explicit viscosity and diffusion. Calculations with the model support the theoretical conjecture that the difference scheme is stable whenever the analogous Cartesian scheme is stable. The scheme has been used to calculate dry atmospheric convection due to differential heating between top and bottom of mountainous terrain. The general small-scale characteristics of mountain up-slope winds have been simulated. In addition, the results have demonstrated the crucial role played by the eddy diffusivities and the environmental stability, in determining both the quantitative and the qualitative features of the circulation.

O
Galchen, T, Somerville RCJ.  1975.  On the Use of a Coordinate Transformation for Solution of Navier-Stokes Equations. Journal of Computational Physics. 17:209-228.   10.1016/0021-9991(75)90037-6   AbstractWebsite

The equations of fluid motion have been formulated in a generalized noncartesian, non-orthogonal coordinate system. A particular coordinate transformation, which transforms a domain with an irregular lower boundary into a cube, has been constructed. The transformed system, unlike the original one, has flat boundaries and homogeneous boundary conditions. Where the topography is flat, the original and transformed systems are identical, and extra terms do not appear. A finite difference scheme for solving the transformed equations has been constructed and will be described in a subsequent issue of this journal.

Pritchard, MS, Moncrieff MW, Somerville RCJ.  2011.  Orogenic Propagating Precipitation Systems over the United States in a Global Climate Model with Embedded Explicit Convection. Journal of the Atmospheric Sciences. 68:1821-1840.   10.1175/2011jas3699.1   AbstractWebsite

In the lee of major mountain chains worldwide, diurnal physics of organized propagating convection project onto seasonal and climate time scales of the hydrologic cycle, but this phenomenon is not represented in conventional global climate models (GCMs). Analysis of an experimental version of the superparameterized (SP) Community Atmosphere Model (CAM) demonstrates that propagating orogenic nocturnal convection in the central U.S. warm season is, however, representable in GCMs that use the embedded explicit convection model approach [i.e., multiscale modeling frameworks (MMFs)]. SP-CAM admits propagating organized convective systems in the lee of the Rockies during synoptic conditions similar to those that generate mesoscale convective systems in nature. The simulated convective systems exhibit spatial scales, phase speeds, and propagation speeds comparable to radar observations, and the genesis mechanism in the model agrees qualitatively with established conceptual models. Convective heating and condensate structures are examined on both resolved scales in SP-CAM, and coherently propagating cloud "metastructures" are shown to transcend individual cloud-resolving model arrays. In reconciling how this new mode of diurnal convective variability is admitted in SP-CAM despite the severe idealizations in the cloud-resolving model configuration, an updated discussion is presented of what physics may transcend the re-engineered scale interface in MMFs. The authors suggest that the improved diurnal propagation physics in SP-CAM are mediated by large-scale first-baroclinic gravity wave interactions with a prognostic organization life cycle, emphasizing the physical importance of preserving "memory" at the inner resolved scale.

P
Somerville, RCJ.  1977.  Pattern Recognition Techniques for Weather Forecast Verification. Contributions to Atmospheric Physics [Beitraege zur Physik der Atmosphaere.], Wiesbaden, Germany. 50:403-410. Abstract
n/a