On the interaction between a radiatively damped planetary wave and the zonally averaged circulation in the middle atmosphere

The interaction between a planetary wave damped by cooling to space and the zonally averaged circulation in the middle atmosphere is examined for a steady-state situation in middle latitudes. Quasi-geostrophic scaling of Type 2 is assumed (i.e. the space scales are planetary and the superrotation is small).A set of mean equations is derived for this scaling which is complementary to the set of perturbation equations previously studied. The mean equations show that a planetary wave induces a mean meridional circulation which is balanced by an eddy momentum forcing function and a mean diabatic heating which is balanced by an eddy heat flux forcing function. The vertical eddy fluxes enter the forcing at the same order as the horizontal eddy fluxes.An analytical wave solution is found for the case of an atmosphere in constant superrotation. The eddy fluxes and forcing functions are evaluated for this special case. It is found that they are very sensitive to the values of the radiative damping coefficient and the superrotation. Since the damping coefficient depends on the ozone concentration and the intensity of the solar ultraviolet flux, the results suggest that changes in these quantities can strongly modify the wave-mean flow interaction in the middle atmosphere. Possible implications for climate change are discussed.     

A theory of enhanced saturation of the gravity wave spectrum due to increases in atmospheric stability

In this paper we consider a vertical wavenumber spectrum of vertically propagating gravity waves impinging on a rapid increase in atmospheric stability. If the high-wavenumber range is saturated below the increase, as is usually observed, then the compression of vertical scales as the waves enter a region of higher stability results in that range becoming supersaturated, that is, the spectral amplitude becomes larger than the saturation limit. The supersaturated wave energy must then dissipate in a vertical distance of the order of a wavelength, resulting in an enhanced turbulent energy dissipation rate. If the wave spectrum is azimuthally anisotropic, the dissipation also results in an enhanced vertical divergence of the vertical flux of horizontal momentum and enhanced wave drag in the same region. Estimates of the enhanced dissipation rates and radar reflectivities appear to be consistent with the enhancements observed near the high-latitude summer mesopause. Estimates of the enhanced mean flow acceleration appear to be consistent with the wave drag that is needed near the tropopause and the high-latitude summer mesopause in large-scale models of the atmosphere. Thus, this process may play a significant role in determining the global effects of gravity waves on the large-scale circulation.     

A consistent derivation of the wave-energy equation from basic hydrodynamic principles

Based on a decomposition of the velocity into mean flow, turbulent and wave components, momentum and hereafter a wave-energy equation is derived. It contains a turbulent energy dissipation term which is closed by applying a wave-related mixing length model and linear wave theory solutions. This closure produces a non-linear turbulent wave-energy dissipation including the wave energy in a 5/2 power law. The theory is able to predict correctly the shape of deep-water wave spectra according to Phillips' similarity law.Responsible Editor: Hans Burchard     

Climatological Ertel’s potential-vorticity flux and mean meridional circulation in the extratropical troposphere – lower stratosphere

We investigated to what extent the isentropic, non-geostrophic formulation of zonally averaged circulation derived for stratospheric conditions is applicable to climatological transport in the extratropical troposphere and lower stratosphere. The study is based on 10 years of daily data of ECMWF analysis and on the ECHAM3 climate model of the German Climate Computing Centre. The main result is a scalar isentropic mixing coefficient, K_yy, and a mean meridional transport circulation consistently derived from the same data base. For both data sources, isentropic mean meridional circulation is derived from horizontal mass flow rate for 4 representative months. Alternatively, a mean meridional circulation is calculated from total diabatic heating rates of the ECHAM3 model. It is shown that only the latter is in good agreement with the ECMWF mean meridional circulation. Isentropic analysis also comprises the seasonal cycle of the climatological meridional gradient and flux of Ertels potential vorticity (PV). Application of Tungs flux-gradient relation yields that for all seasons K_yy is positive in height-latitude regions where statistical significance is reached. Large K_yy values, marking regions of more efficient mixing, have been found in the subtropical vertical band of weak westerly wind and in mid-latitudes in regions of upward-propagating baroclinic wave activity in the middle and upper troposphere. Based on the ECMWF data and results of baroclinic-wave behaviour, strong indications are presented that positive zonally averaged PV flux polewards of the jet core in the NH is strengthened by stationary waves and nonlinear effects. Reduced eddy transport is apparent in winter and spring slightly below the subtropical tropopause jet. The seasonal cycle of K_yy from ECHAM3 data is to a great extent in agreement with the result based on ECMWF analysis. In the model, reduced interannual variability enlarges the height-latitude range where sign of K_yy is significant.     

Causes and underlying dynamic processes of the mid-winter suppression in the North Pacific storm track

Baroclinic wave activity in the North Pacific exhibit peaks in late fall and early spring, and a local minimum in midwinter, when by linear baroclinic instability theory it should attain its maximum. This counterintuitive phenomenon, or"midwinter suppression"(MWM) as called, is investigated with a functional analysis apparatus, multiscale window transform(MWT), and the MWT-based theory of canonical transfer and localized multi-scale energetics analysis, together with a feature tracking technique, using the data from the European Centre for Medium-Range Weather Forecasts ReAnalysis(ERA-40). It is found that the MWM results from a variety of different physical processes, including baroclinic canonical transfer, diabatic effect, energy flux divergence, and frictional dissipation. On one hand, baroclinic canonical transfer and diabatic effect achieve their respective maxima in late fall. More transient available potential energy is produced and then converted to transient kinetic energy, resulting in a stronger storm track in late fall than in midwinter. On the other hand, in early spring, although baroclinic instability and buoyancy conversion are weak, energy flux convergences are substantially strengthened, leading to a net energy inflow into the storm track. Meanwhile, frictional dissipation is greatly reduced in spring; as a result, less transient energy is dissipated in early spring than in midwinter. It is further found that the weakening of baroclinic canonical transfer in midwinter(compared to late fall) is due to the far distance between the storm and the jet stream(located at its southernmost point), which suppresses the interaction between them. Regarding the increase in energy flux convergence in early spring, it appears to originate from the increase(enhancement) in the number(strength) of storms from the upstream into the Pacific.     

The Tsunami of 26 December, 2004: Numerical Modeling and Energy Considerations

A numerical model for the global tsunamis computation constructed by Kowalik et al. (2005), is applied to the tsunami of 26 December, 2004 in the World Ocean from 80°S to 69°N with spatial resolution of one minute. Because the computational domain includes close to 200 million grid points, a parallel version of the code was developed and run on a Cray X1 supercomputer. An energy flux function is used to investigate energy transfer from the tsunami source to the Atlantic and Pacific Oceans. Although the first energy input into the Pacific Ocean was the primary (direct) wave, reflections from the Sri Lankan and eastern shores of Maldives were a larger source. The tsunami traveled from Indonesia, around New Zealand, and into the Pacific Ocean by various routes. The direct path through the deep ocean to North America carried miniscule energy, while the stronger signal traveled a considerably longer distance via South Pacific ridges as these bathymetric features amplified the energy flux vectors. Travel times for these amplified energy fluxes are much longer than the arrival of the first wave. These large fluxes are organized in the wave-like form when propagating between Australia and Antarctica. The sources for the larger fluxes are multiple reflections from the Seychelles, Maldives and a slower direct signal from the Bay of Bengal. The energy flux into the Atlantic Ocean shows a different pattern since the energy is pumped into this domain through the directional properties of the source function. The energy flow into the Pacific Ocean is approximately 75% of the total flow to the Atlantic Ocean. In many locations along the Pacific and Atlantic coasts, the first arriving signal, or forerunner, has lower amplitude than the main signal which often is much delayed. Understanding this temporal distribution is important for an application to tsunami warning and prediction.     

Instability of a geostrophic front and its energetics

Abstract

The instability of a current with a geostrophic surface density front is investigated by means of a reduced gravity model having a velocity profile with nearly uniform potential vorticity. It is shown that currents are unstable when the mean potential vorticity decreases toward the surface front at the critical point of the frontal trapped waves investigated by Paldor (1983). This instability is identical with that demonstrated by Killworth (1983) in the longwave limit.

The cross-stream component of mass flux and the rates of energy conversions among the five energy forms defined by Orlanski (1968) are also calculated. The main results are as follows, (a) The mass flux toward the surface front is positive near the front and negative around the critical point. The positive mass flux near the front does not vanish at the position of the undisturbed surface front, so that the mean position of the front moves outward and the region of the strong current spreads. (b) The potential energy of the mean flow integrated over the fluid is released through the work done by the force of the pressure gradient of the mean flow on the fluid, and is converted into the kinetic energy of the mean flow. (c) In the critical layer, the mean flow is rapidly accelerated with the growth of the unstable wave. This acceleration is caused by the rapid phase shift of the unstable wave in the critical layer.     

Characterization of seismic hazard and structural response by energy flux

Seismic safety of structures depends on the structure's ability to absorb the seismic energy that is transmitted from ground to structure. One parameter that can be used to characterize seismic energy is the energy flux. Energy flux is defined as the amount of energy transmitted per unit time through a cross-section of a medium, and is equal to kinetic energy multiplied by the propagation velocity of seismic waves. The peak or the integral of energy flux can be used to characterize ground motions. By definition, energy flux automatically accounts for site amplification. Energy flux in a structure can be studied by formulating the problem as a wave propagation problem. For buildings founded on layered soil media and subjected to vertically incident plane shear waves, energy flux equations are derived by modeling the building as an extension of the layered soil medium, and considering each story as another layer. The propagation of energy flux in the layers is described in terms of the upgoing and downgoing energy flux in each layer, and the energy reflection and transmission coefficients at each interface. The formulation results in a pair of simple finite-difference equations for each layer, which can be solved recursively starting from the bedrock. The upgoing and downgoing energy flux in the layers allows calculation of the energy demand and energy dissipation in each layer. The methodology is applicable to linear, as well as nonlinear structures.     

Energy balance of simple elastodynamic sources

Complete relations are derived for energy and energy flux of elastic waves generated by an isotropic and double-couple source in a perfectly elastic, homogeneous, isotropic, and unbounded medium. In the energy balance of elastodynamic sources near-field waves play an essential role, transforming static energy into wave energy, andvice versa. For explosive and dislocation sources, the sources surface radiates a positive wave energy that is partially distributed to the medium transforming into static energy. For implosive and antidislocation sources, the source surface generates elastic waves, but it does not necessarily imply that it also radiates a positive wave energy. The energy transported by waves can originate in gradual transformation of the static-to-wave energy during propagation of waves through a stressed medium.On leave from Geophysical Institute, Czech Academy of Sciences, Boní II/1401, 41 31, Praha 4 Czech Republic     

Interaction of multiple courses of wave‐induced fluid flow in layered porous media

Different theoretical and laboratory studies on the propagation of elastic waves in layered hydrocarbon reservoir have shown characteristic velocity dispersion and attenuation of seismic waves. The wave‐induced fluid flow between mesoscopic‐scale heterogeneities (larger than the pore size but smaller than the predominant wavelengths) is the most important cause of attenuation for frequencies below 1 kHz. Most studies on mesoscopic wave‐induced fluid flow in the seismic frequency band are based on the representative elementary volume, which does not consider interaction of fluid flow due to the symmetrical structure of representative elementary volume. However, in strongly heterogeneous media with unsymmetrical structures, different courses of wave‐induced fluid flow may lead to the interaction of the fluid flux in the seismic band; this has not yet been explored. This paper analyses the interaction of different courses of wave‐induced fluid flow in layered porous media. We apply a one‐dimensional finite‐element numerical creep test based on Biot's theory of consolidation to obtain the fluid flux in the frequency domain. The characteristic frequency of the fluid flux and the strain rate tensor are introduced to characterise the interaction of different courses of fluid flux. We also compare the behaviours of characteristic frequencies and the strain rate tensor on two scales: the local scale and the global scale. It is shown that, at the local scale, the interaction between different courses of fluid flux is a dynamic process, and the weak fluid flux and corresponding characteristic frequencies contain detailed information about the interaction of the fluid flux. At the global scale, the averaged strain rate tensor can facilitate the identification of the interaction degree of the fluid flux for the porous medium with a random distribution of mesoscopic heterogeneities, and the characteristic frequency of the fluid flux is potentially related to that of the peak attenuation. The results are helpful for the prediction of the distribution of oil–gas patches based on the statistical properties of phase velocities and attenuation in layered porous media with random disorder.     

Comparison of geostrophic and nonlinear balanced winds from LIMS data and implications for derived dynamical quantities

Nimbus 7 LIMS geopotential height data are utilized to infer the rotational wind distribution in the Northern Hemisphere stratosphere and lower mesosphere during a period of substantial wave-mean flow interaction in January, 1979. Rotational winds are derived from the application of a successive relaxation numerical procedure which incorporates the spherical polar coordinate iterative algorithm ofPaegle andTomlinson (1975) for the nondivergent nonlinear balance equation. Optimum convergence of the numerical solutions is found to occur when under-relaxation is utilized. The LIMS height analyses were also latitudinally smoothed and constrained to obey the ellipticity criterion for spherical coordinates. The balanced winds are compared with geostrophically derived values and within situ radiosonde reports for 100 mb to 10 mb over Berlin.From a localized perspective, the Berlin-LIMS comparison indicates that radiosonde and balanced wind vectors exhibit somewhat closer agreement in direction than is associated with the geostrophic estimates. However, substantial quantitative differences between radiosonde, balanced, and geostrophic wind speeds are also evident, suggesting that caution should be exercised in the local application of derived winds, as for example in the quantitative interpretation of trajectories derived from satellite height analyses during periods of enhanced stratospheric wave activity.On a longitudinally averaged basis, balanced zonal-mean wind speeds are typically 20% weaker than geostrophic values in polar latitudes, and as much as 50% weaker in tropical and midlatitude regions. Meridional balanced wind velocities, at a given longitude, are generally within ±10% of geostrophic values. Although these alterations in horizontal wind components result in only modest differences between balanced and geostrophic meridional eddy heat fluxes, a more substantial change appears in the meridional eddy momentum flux analysis. The corresponding patterns of Eliassen-Palm flux divergence are found to be somewhat more (less) intense for the balanced wind case in the stratosphere (lower mesosphere) in polar latitudes.     

非绝热加热对大气局地扰动位能的影响和机理

本文运用统计相关和奇异值分解方法,系统地考察了非绝热加热对大气局地扰动位能的影响特征和机理问题.分析结果表明,热带地区海表温度异常和扰动位能的耦合相关特征与厄尔尼诺和南方涛动变率的关系密切,扰动位能在热带外地区的耦合模态空间型呈现出与北太平洋-北美大气遥相关型极为相似的分布特征.通过进一步考察大气中视热源与扰动位能的耦合特征,发现在低纬地区仍主要反映了厄尔尼诺和南方涛动的影响机制,在热带外地区,视热源异常的耦合模态则表现出与北半球环状模类似的特点,这表明热带外地区大气的局地能量有效性与热带非绝热加热的遥强迫以及局地性热源加热强迫有着极为密切的联系.     

Oblique wave diffraction by a flexible floating structure in the presence of a submerged flexible structure

Expansion formulae associated with the interaction of oblique surface gravity waves with a floating flexible plate in the presence of a submerged horizontal flexible structure are derived using Green’s integral theorem in water of finite and infinite water depths. The associated Green’s functions are derived using the fundamental solution associated with the reduced wave equation. The integral forms of the Green’s functions and the velocity potentials are advantageous over the eigenfunction expansion method in situation when the roots of the dispersion relation coalesce. As an application of the expansion formulae, diffraction of oblique waves by a finite floating elastic plate in the presence of a submerged horizontal flexible membrane is investigated in water of finite depth. The accuracy of the numerical computation is demonstrated by analysing the convergence of the complex amplitude of the reflected waves and the energy relation. Effect of the submerged membrane on the diffraction of surface waves is studied by analysing the reflection and transmission coefficients for various parametric values. Further, the derivation of long wave equation under shallow water approximation is derived in a direct manner in the appendix. The concept and methodology can be easily extended to deal with acoustic wave interaction with flexible structures and related problems of mathematical physics and engineering.     

Maximum seismic displacement of inelastic systems based on energy concept

The energy balance and energy input of Single-Degree-of-Freedom (SDF) systems under earthquake motion is studied for elastic and inelastic systems. The maximum displacement of an inelastic system is related to that of an elastic system having the same initial stiffness and mass by considering the earthquake energy input per cycle of oscillation. With an assumption that the cyclic energy input is equal for both elastic and inelastic system for intermediate- and long-period systems, a simplified relation is suggested. Newmark's equal energy rule is shown to be the upper bound of the simplified relation; the lower bound of maximum response displacement is also derived in this paper. The numerical analysis results were mostly shown to fall between the proposed upper and lower bounds. A separate approximate relation is proposed for short-period systems. The reason for divergence from the suggested relations is discussed for short-period systems. Copyright © 1999 John Wiley & Sons Ltd.     

Diurnal cycle of surface energy fluxes in high mountain terrain: High-resolution fully coupled atmosphere-hydrology modelling and impact of lateral flow

Water and energy fluxes are inextricably interlinked within the interface of the land surface and the atmosphere. In the regional earth system models, the lower boundary parameterization of land surface neglects lateral hydrological processes, which may inadequately depict the surface water and energy fluxes variations, thus affecting the simulated atmospheric system through land-atmosphere feedbacks. Therefore, the main objective of this study is to evaluate the hydrologically enhanced regional climate modelling in order to represent the diurnal cycle of surface energy fluxes in high spatial and temporal resolution. In this study, the Weather Research and Forecasting model (WRF) and coupled WRF Hydrological modelling system (WRF-Hydro) are applied in a high alpine catchment in Northeastern Tibetan Plateau, the headwater area of the Heihe River. By evaluating and intercomparing model results by both models, the role of lateral flow processes on the surface energy fluxes dynamics is investigated. The model evaluations suggest that both WRF and coupled WRF-Hydro reasonably represent the diurnal variations of the near-surface meteorological fields, surface energy fluxes and hourly partitioning of available energy. By incorporating additional lateral flow processes, the coupled WRF-Hydro simulates higher surface soil moisture over the mountainous area, resulting in increased latent heat flux and decreased sensible heat flux of around 20–50 W/m² in their diurnal peak values during summertime, although the net radiation and ground heat fluxes remain almost unchanged. The simulation results show that the diurnal cycle of surface energy fluxes follows the local terrain and vegetation features. This highlights the importance of consideration of lateral flow processes over areas with heterogeneous terrain and land surfaces.     

真实基流中非线性Rossby波演变特征(二):能量、结构演变及初始场影响

针对非线性的准地转正压位涡方程,利用自行设计的差分格式和高斯函数拟合得到的真实基流分布,数值研究了线性和非线性Rossby波流场结构和总能量的演变以及初值对总能量演变的影响.发现在非线性的真实基流中,线性和非线性Rossby波的相对总能量出现振荡型增长或衰减,非线性波动的振荡周期明显小于线性波动,非线性项不仅抑制能量的...     

Transport and the seasonal variation of ozone

The transport mechanisms responsible for the seasonal behavior of total ozone are deduced from the comparison of model results to stratospheric data. The seasonal transport is dominated by a combination of the diabatic circulation and transient planetary wave activity acting on a diffusively and photochemically determined background state. The seasonal variation is not correctly modeled as a diffusive process. The buildup of total ozone at high latitudes during winter is dependent upon transient planetary wave activity of sufficient strength to cause the breakdown of the polar vortex. While midwinter warmings are responsible for enhanced ozone transport to high latitudes, the final warming marking the transition from zonal mean westerlies to zonal mean easterlies is the most important event leading to the spring maximum. The final warming is not followed by reacceleration of the mean flow; so that the ozone transport associated with this event is more pronounced than that associated with midwinter warmings.     

Tracking Earth’s Energy: From El Ni?o to Global Warming

The state of knowledge and outstanding issues with respect to the global mean energy budget of planet Earth are described, along with the ability to track changes over time. Best estimates of the main energy components involved in radiative transfer and energy flows through the climate system do not satisfy physical constraints for conservation of energy without adjustments. The main issues relate to the downwelling longwave (LW) radiation and the hydrological cycle, and thus the surface evaporative cooling. It is argued that the discrepancy is 18% of the surface latent energy flux, but only 4% of the downwelling LW flux and, for various reasons, it is most likely that the latter is astray in some calculations, including many models, although there is also scope for precipitation estimates to be revised. Beginning in 2000, the top-of-atmosphere radiation measurements provide stable estimates of the net global radiative imbalance changes over a decade, but after 2004 there is “missing energy” as the observing system of the changes in ocean heat content, melting of land ice, and so on is unable to account for where it has gone. Based upon a number of climate model experiments for the twenty-first century where there are stases in global surface temperature and upper ocean heat content in spite of an identifiable global energy imbalance, we infer that the main sink of the missing energy is likely the deep ocean below 275?m depth.     

Hemispheric simulation of the Asian summer monsoon

A three-level, -plane, filtered model is used to simulate the Northern Hemisphere summer monsoon. A time-averaged initial state, devoid of sub-planetary scale waves, is integrated through 30 days on a 5° latitude-longitude grid. Day 25 through day 30 integrations are then repeated on a 2.5° grid. The planetary-scale waves are forced by time-independent, spatially varying diabatic heating. Energy is extracted via internal and surface frictional processes. Orography is excluded to simplify synoptic-scale energy sources.During integration the model energy first increases, but stabilizes near day 10. Subsequent flow patterns closely resemble the hemisphere summer monsoon. Climatological features remain quasi-stationary. At 200 mb high pressure dominates the land area, large-scale troughs are found over the Atlantic and Pacific Oceans, the easterly jet forms south of Asia, and subtropical jets develop in the westerlies. At 800 mb subtropical highs dominate the oceans and the monsoon trough develops over the Asian land mass. The planetary scales at all levels develop a realistic cellular structure from the passage of transient synoptic-scale features, e.g., a baroclinic cyclone track develops near 55°N and westward propagating waves form in the easterlies.Barotropic redistribution of kinetic energy is examined over a low-latitude zonal strip using a Fourier wave-space. In contrast to higher latitudes where the zonal flow and both longer and shorter waves are fed by barotropic energy redistribution from the baroclinically unstable wavelengths, the low-latitude waves have a planetary-scale kinetic energy source. Wave numbers 1 and 2 maintain both the zonal flow and all shorter scales via barotropic transfers. Transient and standing wave processes are examined individually and in combination.Wave energy accumulates at wave numbers 7 and 8 at 200 mb and at wave number 11 in the lower troposphere. The 800-mb waves are thermally indirect and in the mean they give energy to the zonal flow. These characteristics agree with atmospheric observation. The energy source for these waves is the three wave barotropic transfer. The implications of examining barotropic processes in a Fourier wave-space, vice the more common approach of separating the flow into a mean plus a deviation are discussed.