Vertical propagation of information in a middle atmosphere data assimilation system by gravity-wave drag feedbacks

The mesospheric response to the 2002 Antarctic Stratospheric Sudden Warming (SSW) is analysed using the Canadian Middle Atmosphere Model Data Assimilation System (CMAM-DAS), where it represents a vertical propagation of information from the observations into the data-free mesosphere. The CMAM-DAS simulates a cooling in the lowest part of the mesosphere which is accomplished by resolved motions, but which is extended to the mid- to upper mesosphere by the response of the model's non-orographic gravity-wave drag parameterization to the change in zonal winds. The basic mechanism is that elucidated by Holton consisting of a net eastward wave-drag anomaly in the mesosphere during the SSW, although in this case there is a net upwelling in the polar mesosphere. Since the zonal-mean mesospheric response is shown to be predictable, this demonstrates that variations in the mesospheric state can be slaved to the lower atmosphere through gravity-wave drag.

Salmon consumption during pregnancy alters fatty acid composition and secretory IgA concentration in human breast milk

Fish oil supplementation during pregnancy alters breast milk composition, but there is little information about the impact of oily fish consumption. We determined whether increased salmon consumption during pregnancy alters breast milk fatty acid composition and immune factors. Women (n = 123) who rarely ate oily fish were randomly assigned to consume their habitual diet or to consume 2 portions of farmed salmon per week from 20 wk of pregnancy until delivery. The salmon provided 3.45 g long-chain (LC) (n-3) PUFA/wk. Breast milk fatty acid composition and immune factors [soluble CD14, transforming growth factor-b (TGFb)1, TGFb2, and secretory IgA] were analyzed at 1, 5, 14, and 28 d postpartum (PP). Breast milk from the salmon group had higher proportions of EPA (80%), docosapentaenoic acid (30%), and DHA (90%) on d 5 PP compared with controls (P < 0.01). The LC (n-6) PUFA:LC (n-3) PUFA ratio was lower for the salmon group on all days of PP sampling (P < 0.004), although individual (n-6) PUFA proportions, including arachidonic acid, did not differ. All breast milk immune factors decreased between d 1 and 28 PP (P < 0.001). Breast milk secretory IgA (sIgA) was lower in the salmon group (d 1–28 PP; P = 0.006). Salmon consumption during pregnancy, at the current recommended intakes, increases the LC (n-3) PUFA concentration of breast milk in early lactation, thus improving the supply of these important fatty acids to the breast-fed neonate. The consequence of the lower breast milk concentration of sIgA in the salmon group is not clear.

Effects of urban ventilation patterns on the carbon monoxide concentration in a high-rise mega city

Carbon monoxide (CO) concentration data from 1999–2006, monitored at 5 different pollution stations in a high-rise mega city (Hong Kong), were collected and investigated. The spatio-temporal characteristics of urban CO concentration profiles were obtained. A new approach was put forward to examine the relationship between urban CO concentration and different wind flow patterns. Rather than relying on the meteorological data from a single weather station, usually adopted in previous studies, four weather stations on the boundary of Hong Kong territory were used in the present study so as to identify 16 different wind flow patterns, among which a typical urban heat island circulation (UHIC) can be distinguished. Higher concentrations were observed to be associated with the flow pattern of an inflow from Lau Fau Shan (LFS) station which is located in the northwest of Hong Kong. This suggests that the ability of dilution for north-to-west wind is relatively weak due to the pollutants carried from outside Hong Kong. The effectiveness of wind speed on the alleviation of urban concentration is dependent on the initial concentration of the approaching wind. The increase of wind speed of north-to-west wind from 0 m/s to 6 m/s has little effect on the reduction of urban CO concentration, especially on the non-roadside stations. By contrast, for the southerly marine wind, pollution concentration decreases sharply with an increase in the wind speed. It was also found that urban heat island circulation (UHIC) is conducive of the accumulation of pollutants, especially at night. There exists a positive correlation between CO concentration and UHI intensity. This correlation is much stronger at night compared to during the day. Keywords: urban pollution monitoring, urban ventilation pattern, urban heat island circulation, mega city

Asymptotic gravity wave drag expressions for non-hydrostatic rotating flow over a ridge

Asymptotic expressions are derived for the mountain wave drag in flow with constant wind and static stability over a ridge when both rotation and non-hydrostatic effects are important. These expressions, which are much more manageable than the corresponding exact drag expressions (when these do exist) are found to provide accurate approximations to the drag, even when non-hydrostatic and rotation effects are strong, despite having been developed for cases where these effects are weak. The derived expressions are compared with approximations to the drag found previously, and their asymptotic behaviour in various limits is studied.

A linear model of gravity wave drag for hydrostatic sheared flow over elliptical mountains

An analytical model of orographic gravity wave drag due to sheared flow past elliptical mountains is developed. The model extends the domain of applicability of the well-known Phillips model to wind profiles that vary relatively slowly in the vertical, so that they may be treated using a WKB approximation. The model illustrates how linear processes associated with wind profile shear and curvature affect the drag force exerted by the airflow on mountains, and how it is crucial to extend the WKB approximation to second order in the small perturbation parameter for these effects to be taken into account. For the simplest wind profiles, the normalized drag depends only on the Richardson number, Ri, of the flow at the surface and on the aspect ratio, γ, of the mountain. For a linear wind profile, the drag decreases as Ri decreases, and this variation is faster when the wind is across the mountain than when it is along the mountain. For a wind that rotates with height maintaining its magnitude, the drag generally increases as Ri decreases, by an amount depending on γ and on the incidence angle. The results from WKB theory are compared with exact linear results and also with results from a non-hydrostatic nonlinear numerical model, showing in general encouraging agreement, down to values of Ri of order one.

Linear criteria for gravity-wave breaking in resonant stratified flow over a ridge

Using linear theory, it is shown that, in resonant flow over a 2D mountain ridge, such as exists when a layer of uniform wind is topped by an environmental critical level, the conditions for internal gravity-wave breaking are different from those determined in previous studies for non-resonant flows. For Richardson numbers in the shear layer not exceeding 2.25, two zones of flow overturning exist, respectively below and downstream and above and upstream of the expected locations. Flow overturning occurs for values of the dimensionless height of the ridge smaller than those required for a uniform wind profile. These results may have implications for the physical understanding of high-drag states.

Resonant gravity-wave drag enhancement in linear stratified flow over mountains

High-drag states produced in stratified flow over a 2D ridge and an axisymmetric mountain are investigated using a linear, hydrostatic, analytical model. A wind profile is assumed where the background velocity is constant up to a height z1 and then decreases linearly, and the internal gravity-wave solutions are calculated exactly. In flow over a 2D ridge, the normalized surface drag is given by a closed-form analytical expression, while in flow over an axisymmetric mountain it is given by an expression involving a simple 1D integral. The drag is found to depend on two dimensionless parameters: a dimensionless height formed with z_1, and the Richardson number, Ri, in the shear layer. The drag oscillates as z_1 increases, with a period of half the hydrostatic vertical wavelength of the gravity waves. The amplitude of this modulation increases as Ri decreases. This behaviour is due to wave reflection at z_1. Drag maxima correspond to constructive interference of the upward- and downward-propagating waves in the region z < z_1, while drag minima correspond to destructive interference. The reflection coefficient at the interface z = z_1 increases as Ri decreases. The critical level, z_c, plays no role in the drag amplification. A preliminary numerical treatment of nonlinear effects is presented, where z_c appears to become more relevant, and flow over a 2D ridge qualitatively changes its character. But these effects, and their connection with linear theory, still need to be better understood.

The effect of wind shear and curvature on the gravity wave drag produced by a ridge

The analytical model proposed by Teixeira, Miranda, and Valente is modified to calculate the gravity wave drag exerted by a stratified flow over a 2D mountain ridge. The drag is found to be more strongly affected by the vertical variation of the background velocity than for an axisymmetric mountain. In the hydrostatic approximation, the corrections to the drag due to this effect do not depend on the detailed shape of the ridge as long as this is exactly 2D. Besides the drag, all the perturbed quantities of the flow at the surface, including the pressure, may be calculated analytically.

Sources of variation in the ampicillin-resistant Escherichia coli concentration in the feces of organic broiler chickens

Currently, there are limited published data for the population dynamics of antimicrobial-resistant commensal bacteria. This study was designed to evaluate both the proportions of the Escherichia coli populations that are resistant to ampicillin at the level of the individual chicken on commercial broiler farms and the feasibility of obtaining repeated measures of fecal E. coli concentrations. Short-term temporal variation in the concentration of fecal E. coli was investigated, and a preliminary assessment was made of potential factors involved in the shedding of high numbers of ampicillin-resistant E. coli by growing birds in the absence of the use of antimicrobial drugs. Multilevel linear regression modeling revealed that the largest component of random variation in log-transformed fecal E. coli concentrations was seen between sampling occasions for individual birds. The incorporation of fixed effects into the model demonstrated that the older, heavier birds in the study were significantly more likely (P = 0.0003) to shed higher numbers of ampicillin-resistant E. coli. This association between increasing weight and high shedding was not seen for the total fecal E. coli population (P = 0.71). This implies that, in the absence of the administration of antimicrobial drugs, the proportion of fecal E. coli that was resistant to ampicillin increased as the birds grew. This study has shown that it is possible to collect quantitative microbiological data on broiler farms and that such data could make valuable contributions to risk assessments concerning the transfer of resistant bacteria between animal and human populations.

Recycling and resensitization of the neurokinin 1 receptor: influence of agonist concentration and Rab GTPases

Substance P (SP) induces endocytosis and recycling of the neurokinin 1 receptor (NK1R) in endothelial cells and spinal neurons at sites of inflammation and pain, and it is thus important to understand the mechanism and function of receptor trafficking. We investigated how the SP concentration affects NK1R trafficking and determined the role of Rab GTPases in trafficking. NK1R trafficking was markedly influenced by the SP concentration. High SP (10 nM) induced translocation of the NK1R and beta-arrestin 1 to perinuclear sorting endosomes containing Rab5a, where NK1R remained for >60 min. Low SP (1 nM) induced translocation of the NK1R to early endosomes located immediately beneath the plasma membrane that also contained Rab5a and beta-arrestin 1, followed by rapid recycling of the NK1R. Overexpression of Rab5a promoted NK1R translocation to perinuclear sorting endosomes, whereas the GTP binding-deficient mutant Rab5aS34N caused retention of the NK1R in superficial early endosomes. NK1R translocated from superficial early endosomes to recycling endosomes containing Rab4a and Rab11a, and Rab11aS25N inhibited NK1R recycling. Rapid NK1R recycling coincided with resensitization of SP-induced Ca2+ mobilization and with the return of surface SP binding sites. Resensitization was minimally affected by inhibition of vacuolar H(+)-ATPase and phosphatases but was markedly suppressed by disruption of Rab4a and Rab11a. Thus, whereas beta-arrestins mediate NK1R endocytosis, Rab5a regulates translocation between early and sorting endosomes, and Rab4a and Rab11a regulate trafficking through recycling endosomes. We have thus identified a new function of Rab5a as a control protein for directing concentration-dependent trafficking of the NK1R into different intracellular compartments and obtained evidence that Rab4a and Rab11a contribute to G-protein-coupled receptor recycling from early endosomes.

Is missing orographic gravity wave drag near 60°S the cause of the stratospheric zonal wind biases in chemistry–climate models?

Nearly all chemistry–climate models (CCMs) have a systematic bias of a delayed springtime breakdown of the Southern Hemisphere (SH) stratospheric polar vortex, implying insufficient stratospheric wave drag. In this study the Canadian Middle Atmosphere Model (CMAM) and the CMAM Data Assimilation System (CMAM-DAS) are used to investigate the cause of this bias. Zonal wind analysis increments from CMAMDAS reveal systematic negative values in the stratosphere near 608S in winter and early spring. These are interpreted as indicating a bias in the model physics, namely, missing gravity wave drag (GWD). The negative analysis increments remain at a nearly constant height during winter and descend as the vortex weakens, much like orographic GWD. This region is also where current orographic GWD parameterizations have a gap in wave drag, which is suggested to be unrealistic because of missing effects in those parameterizations. These findings motivate a pair of free-runningCMAMsimulations to assess the impact of extra orographicGWDat 608S. The control simulation exhibits the cold-pole bias and delayed vortex breakdown seen in the CCMs. In the simulation with extra GWD, the cold-pole bias is significantly reduced and the vortex breaks down earlier. Changes in resolved wave drag in the stratosphere also occur in response to the extra GWD, which reduce stratospheric SH polar-cap temperature biases in late spring and early summer. Reducing the dynamical biases, however, results in degraded Antarctic column ozone. This suggests that CCMs that obtain realistic column ozone in the presence of an overly strong and persistent vortex may be doing so through compensating errors.

Sensitivity of simulated climate to conservation of momentum in gravity wave drag parameterization

The Canadian Middle Atmosphere Model is used to examine the sensitivity of simulated climate to conservation of momentum in gravity wave drag parameterization. Momentum conservation requires that the parameterized gravity wave momentum flux at the top of the model be zero and corresponds to the physical boundary condition of no momentum flux at the top of the atmosphere. Allowing momentum flux to escape the model domain violates momentum conservation. Here the impact of momentum conservation in two sets of model simulations is investigated. In the first set, the simulation of present-day climate for two model-lid height configurations, 0.001 and 10 hPa, which are identical below 10 hPa, is considered. The impact of momentum conservation on the climate with the model lid at 0.001 hPa is minimal, which is expected because of the small amount of gravity wave momentum flux reaching 0.001 hPa. When the lid is lowered to 10 hPa and momentum is conserved, there is only a modest impact on the climate in the Northern Hemisphere; however, the Southern Hemisphere climate is more adversely affected by the deflection of resolved waves near the model lid. When momentum is not conserved in the 10-hPa model the climate is further degraded in both hemispheres, particularly in winter at high latitudes, and the impact of momentum conservation extends all the way to the surface. In the second set of simulations, the impact of momentum conservation and model-lid height on the modeled response to ozone depletion in the Southern Hemisphere is considered, and it is found that the response can display significant sensitivity to both factors. In particular, both the lower-stratospheric polar temperature and surface responses are significantly altered when the lid is lowered, with the effect being most severe when momentum is not conserved. The implications with regard to the current round of Intergovernmental Panel on Climate Change model projections are discussed.

Angular momentum conservation and gravity wave drag parameterization: implications for climate models

The robustness of the parameterized gravity wave response to an imposed radiative perturbation in the middle atmosphere is examined. When momentum is conserved and for reasonable gravity wave drag parameters, the response to a polar cooling induces polar downwelling above the region of the imposed cooling, with consequent adiabatic warming. This response is robust to changes in the gravity wave source spectrum, background flow, gravity wave breaking criterion, and model lid height. When momentum is not conserved, either in the formulation or in the implementation of the gravity wave drag parameterization, the response becomes sensitive to the above-mentioned factors—in particular to the model lid height. The spurious response resulting from nonconservation is found to be nonnegligible in terms of the total gravity wave drag–induced downwelling.

Estimation of optimal gravity wave parameters for climate models using data assimilation

There is a current need to constrain the parameters of gravity wave drag (GWD) schemes in climate models using observational information instead of tuning them subjectively. In this work, an inverse technique is developed using data assimilation principles to estimate gravity wave parameters. Because mostGWDschemes assume instantaneous vertical propagation of gravity waves within a column, observations in a single column can be used to formulate a one-dimensional assimilation problem to estimate the unknown parameters. We define a cost function that measures the differences between the unresolved drag inferred from observations (referred to here as the ‘observed’ GWD) and the GWD calculated with a parametrisation scheme. The geometry of the cost function presents some difficulties, including multiple minima and ill-conditioning because of the non-independence of the gravity wave parameters. To overcome these difficulties we propose a genetic algorithm to minimize the cost function, which provides a robust parameter estimation over a broad range of prescribed ‘true’ parameters. When real experiments using an independent estimate of the ‘observed’ GWD are performed, physically unrealistic values of the parameters can result due to the non-independence of the parameters. However, by constraining one of the parameters to lie within a physically realistic range, this degeneracy is broken and the other parameters are also found to lie within physically realistic ranges. This argues for the essential physical self-consistency of the gravity wave scheme. A much better fit to the observed GWD at high latitudes is obtained when the parameters are allowed to vary with latitude. However, a close fit can be obtained either in the upper or the lower part of the profiles, but not in both at the same time. This result is a consequence of assuming an isotropic launch spectrum. The changes of sign in theGWDfound in the tropical lower stratosphere, which are associated with part of the quasi-biennial oscillation forcing, cannot be captured by the parametrisation with optimal parameters.