Curl free pressure gradients over orography in a solution of the fully compressible Euler equations with implicit treatment of acoustic and gravity waves

Steep orography can cause noisy solutions and instability in models of the atmosphere. A new technique for modelling flow over orography is introduced which guarantees curl free gradients on arbitrary grids, implying that the pressure gradient term is not a spurious source of vorticity. This mimetic property leads to better hydrostatic balance and better energy conservation on test cases using terrain following grids. Curl-free gradients are achieved by using the co-variant components of velocity over orography rather than the usual horizontal and vertical components. In addition, gravity and acoustic waves are treated implicitly without the need for mean and perturbation variables or a hydrostatic reference profile. This enables a straightforward description of the implicit treatment of gravity waves. Results are presented of a resting atmosphere over orography and the curl-free pressure gradient formulation is advantageous. Results of gravity waves over orography are insensitive to the placement of terrain-following layers. The model with implicit gravity waves is stable in strongly stratified conditions, with N∆t up to at least 10 (where N is the Brunt-V ̈ais ̈al ̈a frequency). A warm bubble rising over orography is simulated and the curl free pressure gradient formulation gives much more accurate results for this test case than a model without this mimetic property.

Impact of non-hydrostatic effects and trapped lee waves on mountain wave drag in directionally sheared flow

The orographic gravity wave drag produced in flow over an axisymmetric mountain when both vertical wind shear and non-hydrostatic effects are important was calculated using a semi-analytical two-layer linear model, including unidirectional or directional constant wind shear in a layer near the surface, above which the wind is constant. The drag behaviour is determined by partial wave reflection at the shear discontinuity, wave absorption at critical levels (both of which exist in hydrostatic flow), and total wave reflection at levels where the waves become evanescent (an intrinsically non-hydrostatic effect), which produces resonant trapped lee wave modes. As a result of constructive or destructive wave interference, the drag oscillates with the thickness of the constant-shear layer and the Richardson number within it (Ri), generally decreasing at low Ri and when the flow is strongly non-hydrostatic. Critical level absorption, which increases with the angle spanned by the wind velocity in the constant-shear layer, shields the surface from reflected waves, keeping the drag closer to its hydrostatic limit. While, for the parameter range considered here, the drag seldom exceeds this limit, a substantial drag fraction may be produced by trapped lee waves, particularly when the flow is strongly non-hydrostatic, the lower layer is thick and Ri is relatively high. In directionally sheared flows with Ri = O(1), the drag may be misaligned with the surface wind in a direction opposite to the shear, a behaviour which is totally due to non-trapped waves. The trapped lee wave drag, whose reaction force on the atmosphere is felt at low levels, may therefore have a distinctly different direction from the drag associated with vertically propagating waves, which acts on the atmosphere at higher levels.

Analysis of wave phenomena in a morphogenetic mechanochemical model and an application to post-fertilization waves on eggs

Wave solutions to a mechanochemical model for cytoskeletal activity are studied and the results applied to the waves of chemical and mechanical activity that sweep over an egg shortly after fertilization. The model takes into account the calcium-controlled presence of actively contractile units in the cytoplasm, and consists of a viscoelastic force equilibrium equation and a conservation equation for calcium. Using piecewise linear caricatures, we obtain analytic solutions for travelling waves on a strip and demonstrate uiat the full nonlinear system behaves as predicted by the analytic solutions. The equations are solved on a sphere and the numerical results are similar to the analytic solutions. We indicate how the speed of the waves can be used as a diagnostic tool with which the chemical reactivity of the egg surface can be measured.

Ion acceleration at both the interior and exterior Alfvén waves associated with the magnetopause reconnection site: Signatures in cusp precipitation

We present evidence for the acceleration of magnetospheric ions by reflection off two Alfvén waves, launched by the reconnection site into the inflow regions on both sides of the reconnecting magnetopause. The “exterior” wave stands in the inflow from the magnetosheath and is the magnetopause, in the sense that the majority of the field rotation occurs there. The other, “interior” wave stands in the inflow region on the magnetospheric side of the boundary. The population reflected by the interior wave is the more highly energized of the two and appears at low altitudes on open field lines, immediately equatorward of the cusp precipitation. In addition, we identify the population of magnetosheath ions transmitted across the exterior Alfvén wave, as well as a population of magnetospheric ions which are accelerated, after transmission through the interior wave, by reflection off the exterior wave. The ion populations near the X line are modeled and, with allowance for time-of-flight effects, are also derived from observations in the dayside auroral ionosphere. Agreement between observed and theoretical spectra is very good and the theory also explains the observed total fluxes and average energies of the precipitations poleward of the open/closed field line boundary. The results offer a physical interpretation of all the various classifications of precipitation into the dayside ionosphere (central plasma sheet, dayside boundary plasma sheet, void, low-latitude boundary layer, cusp, and mantle) and allow the conditions in both the magnetosphere and the magnetosheath adjacent to the X line to be studied.

Simple M-factor algorithm for improved estimation of the basic maximum usable frequency of radio waves reflected from the ionospheric F-region

The equations of Milsom are evaluated, giving the ground range and group delay of radio waves propagated via the horizontally stratified model ionosphere proposed by Bradley and Dudeney. Expressions for the ground range which allow for the effects of the underlying E- and F1-regions are used to evaluate the basic maximum usable frequency or M-factors for single F-layer hops. An algorithm for the rapid calculation of the M-factor at a given range is developed, and shown to be accurate to within 5%. The results reveal that the M(3000)F2-factor scaled from vertical-incidence ionograms using the standard URSI procedure can be up to 7.5% in error. A simple addition to the algorithm effects a correction to ionogram values to make these accurate to 0.5%.

Extracting planetary waves from geomagnetic time series using Empirical Mode Decomposition

Empirical Mode Decomposition is presented as an alternative to traditional analysis methods to decompose geomagnetic time series into spectral components. Important comments on the algorithm and its variations will be given. Using this technique, planetary wave modes of 5-, 10-, and 16-day mean periods can be extracted from magnetic field components of three different stations in Germany. In a second step, the amplitude modulation functions of these wave modes can be shown to contain significant contribution from solar cycle variation through correlation with smoothed sunspot numbers. Additionally, the data indicate connections with geomagnetic jerk occurrences, supported by a second set of data providing reconstructed near-Earth magnetic field for 150 years. Usually attributed to internal dynamo processes within the Earth's outer core, the question of who is impacting whom will be briefly discussed here.

A route to systematic error in forecasts of Rossby waves

Recent work has shown that both the amplitude of upper-level Rossby waves and the tropopause sharpness decrease with forecast lead time for several days in some operational weather forecast systems. In this contribution, the evolution of error growth in a case study of this forecast error type is diagnosed through analysis of operational forecasts and hindcast simulations. Potential vorticity (PV) on the 320-K isentropic surface is used to diagnose Rossby waves. The Rossby-wave forecast error in the operational ECMWF high-resolution forecast is shown to be associated with errors in the forecast of a warm conveyor belt (WCB) through trajectory analysis and an error metric for WCB outflows. The WCB forecast error is characterised by an overestimation of WCB amplitude, a location of the WCB outflow regions that is too far to the southeast, and a resulting underestimation of the magnitude of the negative PV anomaly in the outflow. Essentially the same forecast error development also occurred in all members of the ECMWF Ensemble Prediction System and the Met Office MOGREPS-15 suggesting that in this case model error made an important contribution to the development of forecast error in addition to initial condition error. Exploiting this forecast error robustness, a comparison was performed between the realised flow evolution, proxied by a sequence of short-range simulations, and a contemporaneous forecast. Both the proxy to the realised flow and the contemporaneous forecast a were produced with the Met Office Unified Model enhanced with tracers of diabatic processes modifying potential temperature and PV. Clear differences were found in the way potential temperature and PV are modified in the WCB between proxy and forecast. These results demonstrate that differences in potential temperature and PV modification in the WCB can be responsible for forecast errors in Rossby waves.

Predictability of frontal waves and cyclones

The statistical properties and skill in predictions of objectively identified and tracked cyclonic features (frontal waves and cyclones) are examined in MOGREPS-15, the global 15-day version of the Met Office Global and Regional Ensemble Prediction System (MOGREPS). The number density of cyclonic features is found to decline with increasing lead-time, with analysis fields containing weak features which are not sustained past the first day of the forecast. This loss of cyclonic features is associated with a decline in area averaged enstrophy with increasing lead time. Both feature number density and area averaged enstrophy saturate by around 7 days into the forecast. It is found that the feature number density and area averaged enstrophy of forecasts produced using model versions that include stochastic energy backscatter saturate at higher values than forecasts produced without stochastic physics. The ability of MOGREPS-15 to predict the locations of cyclonic features of different strengths is evaluated at different spatial scales by examining the Brier Skill (relative to the analysis climatology) of strike probability forecasts: the probability that a cyclonic feature center is located within a specified radius. The radius at which skill is maximised increases with lead time from 650km at 12h to 950km at 7 days. The skill is greatest for the most intense features. Forecast skill remains above zero at these scales out to 14 days for the most intense cyclonic features, but only out to 8 days when all features are included irrespective of intensity.

Improving the health forecasting alert system for cold weather and heat-waves in England: a proof-of-concept using temperature-mortality relationships

Objectives In this study a prototype of a new health forecasting alert system is developed, which is aligned to the approach used in the Met Office’s (MO) National Severe Weather Warning Service (NSWWS). This is in order to improve information available to responders in the health and social care system by linking temperatures more directly to risks of mortality, and developing a system more coherent with other weather alerts. The prototype is compared to the current system in the Cold Weather and Heatwave plans via a case-study approach to verify its potential advantages and shortcomings. Method The prototype health forecasting alert system introduces an “impact vs likelihood matrix” for the health impacts of hot and cold temperatures which is similar to those used operationally for other weather hazards as part of the NSWWS. The impact axis of this matrix is based on existing epidemiological evidence, which shows an increasing relative risk of death at extremes of outdoor temperature beyond a threshold which can be identified epidemiologically. The likelihood axis is based on a probability measure associated with the temperature forecast. The new method is tested for two case studies (one during summer 2013, one during winter 2013), and compared to the performance of the current alert system. Conclusions The prototype shows some clear improvements over the current alert system. It allows for a much greater degree of flexibility, provides more detailed regional information about the health risks associated with periods of extreme temperatures, and is more coherent with other weather alerts which may make it easier for front line responders to use. It will require validation and engagement with stakeholders before it can be considered for use.

Activation of p21-activated protein kinase alpha (alpha PAK) by hyperosmotic shock in neonatal ventricular myocytes.

The p21-activated protein kinases (PAKs) may participate in signalling from Cdc42/Rac1 to the stress-regulated MAPKs (SAPKs/JNKs and p38-/HOG-1-related-MAPKs). We characterized the expression and regulation of alpha PAK in cultured ventricular myocytes. alpha PAK was specifically immunoprecipitated from myocyte extracts. High basal alpha PAK activity was detected in unstimulated myocytes. Its activity was increased rapidly (<30 s) by hyperosmotic shock in the presence of okadaic acid, and was maximal by 3 min (187 +/- 7% relative to unstimulated cells). Endothelin-1 and interleukin-1beta, which also activate SAPKs/JNKs, did not increase alpha PAK activity and presumably act through different PAK isoforms or other mechanisms.

Stimulation of multiple mitogen-activated protein kinase sub-families by oxidative stress and phosphorylation of the small heat shock protein, HSP25/27, in neonatal ventricular myocytes.

We investigated the activation of three subfamilies of mitogen-activated protein kinases (MAPKs), namely the stress-activated protein kinases/c-Jun N-terminal kinases (SAPKs/JNKs), the extracellularly responsive kinases (ERKs) and p38-MAPK, by oxidative stress as exemplified by H2O2 in primary cultures of neonatal rat ventricular myocytes. The 46 and 54 kDa species of SAPKs/JNKs were activated 5- and 10-fold, respectively, by 0.1 mM H2O2 (the maximally effective concentration). Maximal activation occurred at 15-30 min, but was still detectable after 2 h. Both ERK1 and ERK2 were activated 16-fold by 0.1 mM H2O2 with a similar time course to the SAPKs/JNKs, and this was comparable with their activation by 1 microM PMA, the most powerful activator of ERKs that we have so far identified in these cells. The activation of ERKs by H2O2 was inhibited by PD98059, which inhibits the activation of MAPK (or ERK) kinases, and by the protein kinase C (PKC) inhibitor, GF109203X. ERK activation was also inhibited by down-regulation of PMA-sensitive PKC isoforms. p38-MAPK was activated by 0.1 mM H2O2 as shown by an increase in its phosphorylation. However, maximal phosphorylation (activation) was more rapid (<5 min) than for the SAPKs/JNKs or the ERKs. We studied the downstream consequences of p38-MAPK activation by examining activation of MAPK-activated protein kinase 2 (MAPKAPK2) and phosphorylation of the MAPKAPK2 substrate, the small heat shock protein HSP25/27. As with p38-MAPK, MAPKAPK2 was rapidly activated (maximal within 5 min) by 0.1 mM H2O2. This activation was abolished by 10 microM SB203580, a selective inhibitor of certain p38-MAPK isoforms. The phosphorylation of HSP25/27 rapidly followed activation of MAPKAPK2 and was also inhibited by SB203580. Phosphorylation of HSP25/27 was associated with a decrease in its aggregation state. These data indicate that oxidative stress is a powerful activator of all three MAPK subfamilies in neonatal rat ventricular myocytes. Activation of all three MAPKs has been associated with the development of the hypertrophic phenotype. However, stimulation of p38-MAPK and the consequent phosphorylation of HSP25/27 may also be important in cardioprotection.

Effects of oxidative stress on the cardiac myocyte proteome: modifications to peroxiredoxins and small heat shock proteins

Endogenous oxidative stress is a likely cause of cardiac myocyte death in vivo. We examined the early (0-2 h) changes in the proteome of isolated cardiac myocytes from neonatal rats exposed to H2O2 (0.1 mM), focussing on proteins with apparent molecular masses of between 20 and 30 kDa. Proteins were separated by two-dimensional gel electrophoresis (2DGE), located by silver-staining and identified by mass spectrometry. Incorporation of [35S]methionine or 32Pi was also studied. For selected proteins, transcript abundance was examined by reverse transcriptase-polymerase chain reaction. Of the 38 protein spots in the region, 23 were identified. Two families showed changes in 2DGE migration or abundance with H2O2 treatment: the peroxiredoxins and two small heat shock protein (Hsp) family members: heat shock 27 kDa protein 1 (Hsp25) and alphaB-crystallin. Peroxiredoxins shifted to lower pI values and this was probably attributable to 'over-oxidation' of active site Cys-residues. Hsp25 also shifted to lower pI values but this was attributable to phosphorylation. alphaB-crystallin migration was unchanged but its abundance decreased. Transcripts encoding peroxiredoxins 2 and 5 increased significantly. In addition, 10 further proteins were identified. For two (glutathione S-transferase pi, translationally-controlled tumour protein), we could not find any previous references indicating their occurrence in cardiac myocytes. We conclude that exposure of cardiac myocytes to oxidative stress causes post-translational modification in two protein families involved in cytoprotection. These changes may be potentially useful diagnostically. In the short term, oxidative stress causes few detectable changes in global protein abundance as assessed by silver-staining.