On the climate response of the low-latitude Pacific Ocean to changes in the global freshwater cycle

Under global warming, the predicted intensification of the global freshwater cycle will modify the net freshwater flux at the ocean surface. Since the freshwater flux maintains ocean salinity structures, changes to the density-driven ocean circulation are likely. A modified ocean circulation could further alter the climate, potentially allowing rapid changes, as seen in the past. The relevant feedback mechanisms and timescales are poorly understood in detail, however, especially at low latitudes where the effects of salinity are relatively subtle. In an attempt to resolve some of these outstanding issues, we present an investigation of the climate response of the low-latitude Pacific region to changes in freshwater forcing. Initiated from the present-day thermohaline structure, a control run of a coupled ocean-atmosphere general circulation model is compared with a perturbation run in which the net freshwater flux is prescribed to be zero over the ocean. Such an extreme experiment helps to elucidate the general adjustment mechanisms and their timescales. The atmospheric greenhouse gas concentrations are held constant, and we restrict our attention to the adjustment of the upper 1,000 m of the Pacific Ocean between 40°N and 40°S, over 100 years. In the perturbation run, changes to the surface buoyancy, near-surface vertical mixing and mixed-layer depth are established within 1 year. Subsequently, relative to the control run, the surface of the low-latitude Pacific Ocean in the perturbation run warms by an average of 0.6°C, and the interior cools by up to 1.1°C, after a few decades. This vertical re-arrangement of the ocean heat content is shown to be achieved by a gradual shutdown of the heat flux due to isopycnal (i.e. along surfaces of constant density) mixing, the vertical component of which is downwards at low latitudes. This heat transfer depends crucially upon the existence of density-compensating temperature and salinity gradients on isopycnal surfaces. The timescale of the thermal changes in the perturbation run is therefore set by the timescale for the decay of isopycnal salinity gradients in response to the eliminated freshwater forcing, which we demonstrate to be around 10-20 years. Such isopycnal heat flux changes may play a role in the response of the low-latitude climate to a future accelerated freshwater cycle. Specifically, the mechanism appears to represent a weak negative sea surface temperature feedback, which we speculate might partially shield from view the anthropogenically-forced global warming signal at low latitudes. Furthermore, since the surface freshwater flux is shown to play a role in determining the ocean's thermal structure, it follows that evaporation and/or precipitation biases in general circulation models are likely to cause sea surface temperature biases.

The role of nonhydrostatic dynamics in controlling development of a surface ocean front

Numerical studies of surface ocean fronts forced by inhomogeneous buoyancy loss show nonhydrostatic convective plumes coexisting with baroclinic eddies. The character of the vertical overturning depends sensitively on the treatment of the vertical momentum equation in the model. It is less well known how the frontal evolution over scales of O(10 km) is affected by these dynamics. Here, we compare highly resolved numerical experiments using nonhydrostatic and hydrostatic models and the convective-adjustment parametrization. The impact of nonhydrostatic processes on average cross-frontal transfer is weak compared to the effect of the O(1 km) scale baroclinic motions. For water-mass distribution and formation rate nonhydrostatic dynamics have similar influence to the baroclinic eddies although adequate resolution of the gradients in forcing fluxes is more important. The overall implication is that including nonhydrostatic surface frontal dynamics in ocean general circulation models will have only a minor effect on scales of O(1 km) and greater.

Vertical current flow through extensive layer clouds

The global atmospheric electrical circuit sustains a vertical current density between the ionosphere and the Earth's surface, the existence of which is well-established from measurements made in fair-weather conditions. In overcast, but non-thunderstorm, non-precipitating conditions, the current travels through the cloud present, despite cloud layers having low electrical conductivity. For extensive layer clouds, this leads to space charge at the upper and lower cloud boundaries. Using a combination of atmospheric electricity and solar radiation measurements at three UK sites, vertical current measurements have been categorised into clear, broken, and overcast cloud conditions. This approach shows that the vertical “fair weather” current is maintained despite the presence of cloud. In fully overcast conditions with thick cloud, the vertical current is reduced compared to thin cloud overcast conditions, associated with the cloud's resistance contributions. Contribution of cloud to the columnar resistance depends both on cloud thickness, and the cloud's height.

Testing the influence of small crystals on ice size spectra using Doppler lidar observations

Analysis of the vertical velocity of ice crystals observed with a 1.5micron Doppler lidar from a continuous sample of stratiform ice clouds over 17 months show that the distribution of Doppler velocity varies strongly with temperature, with mean velocities of 0.2m/s at -40C, increasing to 0.6m/s at -10C due to particle growth and broadening of the size spectrum. We examine the likely influence of crystals smaller than 60microns by forward modelling their effect on the area-weighted fall speed, and comparing the results to the lidar observations. The comparison strongly suggests that the concentration of small crystals in most clouds is much lower than measured in-situ by some cloud droplet probes. We argue that the discrepancy is likely due to shattering of large crystals on the probe inlet, and that numerous small particles should not be included in numerical weather and climate model parameterizations.

The GlobMODEL demonstrator: assimilation of new satellite products in an operational meteorological center

As part of its Data User Element programme, the European Space Agency funded the GlobMODEL project which aimed at investigating the scientific, technical, and organizational issues associated with the use and exploitation of remotely-sensed observations, particularly from new sounders. A pilot study was performed as a "demonstrator" of the GlobMODEL idea, based on the use of new data, with a strong European heritage, not yet assimilated operationally. Two parallel assimilation experiments were performed, using either total column ozone or ozone profiles retrieved at the Royal Netherlands Meteorological Institute (KNMI) from the Ozone Monitoring Instrument (OMI). In both cases, the impact of assimilating OMI data in addition to the total ozone columns from the SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY) on the European Centre for Medium Range Weather Forecasts (ECMWF) ozone analyses was assessed by means of independent measurements. We found that the impact of OMI total columns is mainly limited to the region between 20 and 80 hPa, and is particularly important at high latitudes in the Southern hemisphere where the stratospheric ozone transport and chemical depletion are generally difficult to model with accuracy. Furthermore, the assimilation experiments carried out in this work suggest that OMI DOAS (Differential Optical Absorption Spectroscopy) total ozone columns are on average larger than SCIAMACHY total columns by up to 3 DU, while OMI total columns derived from OMI ozone profiles are on average about 8 DU larger than SCIAMACHY total columns. At the same time, the demonstrator brought to light a number of issues related to the assimilation of atmospheric composition profiles, such as the shortcomings arising when the vertical resolution of the instrument is not properly accounted for in the assimilation. The GlobMODEL demonstrator accelerated scientific and operational utilization of new observations and its results - prompted ECMWF to start the operational assimilation of OMI total column ozone data.

The effect of mechanical stirring on horizontal convection

An important experimental result, as yet poorly understood, is that mechanical stirring can significantly enhance the strength of horizontal convection. A contentious issue is whether this necessarily implies that the mechanical stirring replaces the buoyancy forcing as the main source of energy driving the observed overturning circulation, as has been suggested for the Atlantic meridional overturning circulation (AMOC). In this paper, rigorous energetics considerations and idealized numerical experiments reveal that the rate at which the surface buoyancy forcing supplies energy to the fluid, as measured by the production rate of available potential energy G(APE), does not solely depend upon the buoyancy forcing, as is often implicitly assumed, but also upon the vertical stratification, such that the deeper the thermocline depth, the larger G(APE). This suggests that mechanical stirring enhances horizontal convection because it causes more energy to be extracted from the buoyancy forcing. It does so by enhancing turbulent mixing, which allows surface heating to reach greater depths, which increases the thermocline depth and hence G(APE). This paper therefore proposes a new hypothesis, namely that mechanically stirred horizontal convection and the AMOC are best described as mechanically controlled heat engines.

Latitude and longitude vertical disparities

The literature on vertical disparity is complicated by the fact that several different definitions of the term “vertical disparity” are in common use, often without a clear statement about which is intended or a widespread appreciation of the properties of the different definitions. Here, we examine two definitions of retinal vertical disparity: elevation-latitude and elevation-longitude disparities. Near the fixation point, these definitions become equivalent, but in general, they have quite different dependences on object distance and binocular eye posture, which have not previously been spelt out. We present analytical approximations for each type of vertical disparity, valid for more general conditions than previous derivations in the literature: we do not restrict ourselves to objects near the fixation point or near the plane of regard, and we allow for non-zero torsion, cyclovergence, and vertical misalignments of the eyes. We use these expressions to derive estimates of the latitude and longitude vertical disparities expected at each point in the visual field, averaged over all natural viewing. Finally, we present analytical expressions showing how binocular eye position—gaze direction, convergence, torsion, cyclovergence, and vertical misalignment—can be derived from the vertical disparity field and its derivatives at the fovea.

Seasonal trends in the stable isotopic composition of snow and meltwater runoff in a subarctic catchment at Okstindan, Norway

Seasonal variations in the stable isotopic composition of snow and meltwater were investigated in a sub-arctic, mountainous, but non-glacial, catchment at Okstindan in northern Norway based on analyses of delta(18)O and deltaD. Samples were collected during four field periods (August 1998; April 1999; June 1999 and August 1999) at three sites lying on an altitudinal transect (740-970 m a.s.l.). Snowpack data display an increase in the mean values of delta(18)O (increasing from a mean value of - 13.51 to - 11.49% between April and August), as well as a decrease in variability through the melt period. Comparison with a regional meteoric water line indicates that the slope of the delta(18)O - deltaD line for the snowpacks decreases over the same period, dropping from 7.49 to approximately 6.2. This change points to the role of evaporation in snowpack ablation and is confirmed by the vertical profile of deuterium excess. Snowpack seepage data, although limited, also suggest reduced values of deltaD, as might be associated with local evaporation during meltwater generation. In general, meltwaters were depleted in delta(18)O relative to the source snowpack at the peak of the melt (June), but later in the year (August) the difference between the two was not statistically significant. The diurnal pattern of isotopic composition indicates that the most depleted meltwaters coincide with the peak in temperature and, hence, meltwater production.

The three-dimensional vortical nature of atmospheric and oceanic turbulent flows

Using a novel numerical method at unprecedented resolution, we demonstrate that structures of small to intermediate scale in rotating, stratified flows are intrinsically three-dimensional. Such flows are characterized by vortices (spinning volumes of fluid), regions of large vorticity gradients, and filamentary structures at all scales. It is found that such structures have predominantly three-dimensional dynamics below a horizontal scale LLR, where LR is the so-called Rossby radius of deformation, equal to the characteristic vertical scale of the fluid H divided by the ratio of the rotational and buoyancy frequencies f/N. The breakdown of two-dimensional dynamics at these scales is attributed to the so-called "tall-column instability" [D. G. Dritschel and M. de la Torre Juárez, J. Fluid. Mech. 328, 129 (1996)], which is active on columnar vortices that are tall after scaling by f/N, or, equivalently, that are narrow compared with LR. Moreover, this instability eventually leads to a simple relationship between typical vertical and horizontal scales: for each vertical wave number (apart from the vertically averaged, barotropic component of the flow) the average horizontal wave number is equal to f/N times the vertical wave number. The practical implication is that three-dimensional modeling is essential to capture the behavior of rotating, stratified fluids. Two-dimensional models are not valid for scales below LR. ©1999 American Institute of Physics.

A singular vector perspective of 4DVAR: The spatial structure and evolution of baroclinic weather systems

The extent to which the four-dimensional variational data assimilation (4DVAR) is able to use information about the time evolution of the atmosphere to infer the vertical spatial structure of baroclinic weather systems is investigated. The singular value decomposition (SVD) of the 4DVAR observability matrix is introduced as a novel technique to examine the spatial structure of analysis increments. Specific results are illustrated using 4DVAR analyses and SVD within an idealized 2D Eady model setting. Three different aspects are investigated. The first aspect considers correcting errors that result in normal-mode growth or decay. The results show that 4DVAR performs well at correcting growing errors but not decaying errors. Although it is possible for 4DVAR to correct decaying errors, the assimilation of observations can be detrimental to a forecast because 4DVAR is likely to add growing errors instead of correcting decaying errors. The second aspect shows that the singular values of the observability matrix are a useful tool to identify the optimal spatial and temporal locations for the observations. The results show that the ability to extract the time-evolution information can be maximized by placing the observations far apart in time. The third aspect considers correcting errors that result in nonmodal rapid growth. 4DVAR is able to use the model dynamics to infer some of the vertical structure. However, the specification of the case-dependent background error variances plays a crucial role.

The diurnal cycle of convection and atmospheric tides in an aquaplanet

The diurnal cycle of tropical convection and its relationship to the atmospheric tides is investigated using an aquaplanet GCM. The diurnal and semidiurnal harmonics of precipitation are both found to contribute significantly to the total diurnal variability of precipitation in the model, which is broadly consistent with observations of the diurnal cycle of convection over the open ocean. The semidiurnal tide is found to be the dominant forcing for the semidiurnal harmonic of precipitation. In contrast the diurnal tide plays only a small role in forcing the diurnal harmonic of precipitation, which is dominated by the variations in shortwave and longwave heating. In both the diurnal and semidiurnal harmonics, the feedback onto the convection by the humidity tendencies due to the convection is found to be important in determining the phase of the harmonics. Further experiments show that the diurnal cycle of precipitation is sensitive to the choice of closure in the convection scheme. While the surface pressure signal of the simulated atmospheric tides in the model agree well with both theory and observations in their magnitude and phase, sensitivity experiments suggest that the role of the stratospheric ozone in forcing the semidiurnal tide is much reduced compared to theoretical predictions. Furthermore, the influence of the cloud radiative effects seems small. It is suggested that the radiative heating profile in the troposphere, associated primarily with the water vapor distribution, is more important than previously thought for driving the semidiurnal tide. However, this result may be sensitive to the vertical resolution and extent of the model.

The cumulus-capped boundary layer. I: Modelling transports in the cloud layer

Scalar-flux budgets have been obtained from large-eddy simulations (LESs) of the cumulus-capped boundary layer. Parametrizations of the terms in the budgets are discussed, and two parametrizations for the transport term in the cloud layer are proposed. It is shown that these lead to two models for scalar transports by shallow cumulus convection. One is equivalent to the subsidence detrainment form of convective tendencies obtained from mass-flux parametrizations of cumulus convection. The second is a flux-gradient relationship that is similar in form to the non-local parametrizations of turbulent transports in the dry-convective boundary layer. Using the fluxes of liquid-water potential temperature and total water content from the LES, it is shown that both models are reasonable diagnostic relations between fluxes and the vertical gradients of the mean fields. The LESs used in this study are for steady-state convection and it is possible to treat the fluxes of conserved thermodynamic variables as independent, and ignore the effects of condensation. It is argued that a parametrization of cumulus transports in a model of the cumulus-capped boundary layer should also include an explicit representation of condensation. A simple parametrization of the liquid-water flux in terms of conserved variables is also derived.

The influence of hilly terrain on canopy-atmosphere carbon dioxide exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO2 (N-ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO2 sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO2 exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO2. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N-ee = F-c, the vertical turbulent flux of CO2. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F-c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F-c-at canopy top differs by more than 20% compared to the flat-terrain case.

Effects of fluctuating daily surface fluxes on the time-mean oceanic circulation

The effect of fluctuating daily surface fluxes on the time-mean oceanic circulation is studied using an empirical flux model. The model produces fluctuating fluxes resulting from atmospheric variability and includes oceanic feedbacks on the fluxes. Numerical experiments were carried out by driving an ocean general circulation model with three different versions of the empirical model. It is found that fluctuating daily fluxes lead to an increase in the meridional overturning circulation (MOC) of the Atlantic of about 1 Sv and a decrease in the Antarctic circumpolar current (ACC) of about 32 Sv. The changes are approximately 7% of the MOC and 16% of the ACC obtained without fluctuating daily fluxes. The fluctuating fluxes change the intensity and the depth of vertical mixing. This, in turn, changes the density field and thus the circulation. Fluctuating buoyancy fluxes change the vertical mixing in a non-linear way: they tend to increase the convective mixing in mostly stable regions and to decrease the convective mixing in mostly unstable regions. The ACC changes are related to the enhanced mixing in the subtropical and the mid-latitude Southern Ocean and reduced mixing in the high-latitude Southern Ocean. The enhanced mixing is related to an increase in the frequency and the depth of convective events. As these events bring more dense water downward, the mixing changes lead to a reduction in meridional gradient of the depth-integrated density in the Southern Ocean and hence the strength of the ACC. The MOC changes are related to more subtle density changes. It is found that the vertical mixing in a latitudinal strip in the northern North Atlantic is more strongly enhanced due to fluctuating fluxes than the mixing in a latitudinal strip in the South Atlantic. This leads to an increase in the density difference between the two strips, which can be responsible for the increase in the Atlantic MOC.

Atmospheric electricity coupling between earthquake regions and the ionosphere

We propose a mechanism to explain suggested links between seismic activity and ionospheric changes detected overhead. Specifically, we explain changes in the natural extremely low-frequency (ELF) radio noise recently observed in the topside ionosphere aboard the DEMETER satellite at night, before major earthquakes. Our mechanism utilises increased electrical conductivity of surface layer air before a major earthquake, which reduces the surface-ionosphere electrical resistance. This increases the vertical fair weather current, and (to maintain continuity of electron flow) lowers the ionosphere. Magnitudes of crucial parameters are estimated and found to be consistent with observations. Natural variability in ionospheric and atmospheric electrical properties is evaluated, and may be overcome using a hybrid detection approach. Suggested experiments to investigate the mechanism involve measuring the cut-off frequency of ELF “tweeks”, the amplitude and phase of very low frequency radio waves in the Earth–ionosphere waveguide, or medium frequency radar, incoherent scatter or rocket studies of the lower ionospheric electron density.