The automation and evaluation of nested clade phylogeographic analysis

Nested clade phylogeographic analysis (NCPA) is a popular method for reconstructing the demographic history of spatially distributed populations from genetic data. Although some parts of the analysis are automated, there is no unique and widely followed algorithm for doing this in its entirety, beginning with the data, and ending with the inferences drawn from the data. This article describes a method that automates NCPA, thereby providing a framework for replicating analyses in an objective way. To do so, a number of decisions need to be made so that the automated implementation is representative of previous analyses. We review how the NCPA procedure has evolved since its inception and conclude that there is scope for some variability in the manual application of NCPA. We apply the automated software to three published datasets previously analyzed manually and replicate many details of the manual analyses, suggesting that the current algorithm is representative of how a typical user will perform NCPA. We simulate a large number of replicate datasets for geographically distributed, but entirely random-mating, populations. These are then analyzed using the automated NCPA algorithm. Results indicate that NCPA tends to give a high frequency of false positives. In our simulations we observe that 14% of the clades give a conclusive inference that a demographic event has occurred, and that 75% of the datasets have at least one clade that gives such an inference. This is mainly due to the generation of multiple statistics per clade, of which only one is required to be significant to apply the inference key. We survey the inferences that have been made in recent publications and show that the most commonly inferred processes (restricted gene flow with isolation by distance and contiguous range expansion) are those that are commonly inferred in our simulations. However, published datasets typically yield a richer set of inferences with NCPA than obtained in our random-mating simulations, and further testing of NCPA with models of structured populations is necessary to examine its accuracy.

Proteomic analysis of plasma membrane vesicles isolated from the rat renal cortex

Polarized epithelial cells are responsible for the vectorial transport of solutes and have a key role in maintaining body fluid and electrolyte homeostasis. Such cells contain structurally and functionally distinct plasma membrane domains. Brush border and basolateral membranes of renal and intestinal epithelial cells can be separated using a number of different separation techniques, which allow their different transport functions and receptor expressions to be studied. In this communication, we report a proteomic analysis of these two membrane segments, apical and basolateral, obtained from the rat renal cortex isolated by two different methods: differential centrifugation and free-flow electrophoresis. The study was aimed at assessing the nature of the major proteins isolated by these two separation techniques. Two analytical strategies were used: separation by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) at the protein level or by cation-exchange high-performance liquid chromatography (HPLC) after proteolysis (i.e., at the peptide level). Proteolytic peptides derived from the proteins present in gel pieces or from HPLC fractions after proteolysis were sequenced by on-line liquid chromatography-tandem mass spectrometry (LC-MS/MS). Several hundred proteins were identified in each membrane section. In addition to proteins known to be located at the apical and basolateral membranes, several novel proteins were also identified. In particular, a number of proteins with putative roles in signal transduction were identified in both membranes. To our knowledge, this is the first reported study to try and characterize the membrane proteome of polarized epithelial cells and to provide a data set of the most abundant proteins present in renal proximal tubule cell membranes.

Mechanisms and function of flower and inflorescence reversion

Flower and inflorescence reversion involve a switch from floral development back to vegetative development, thus rendering flowering a phase in an ongoing growth pattern rather than a terminal act of the meristem. Although it can be considered an unusual event, reversion raises questions about the nature and function of flowering. It is linked to environmental conditions and is most often a response to conditions opposite to those that induce flowering. Research on molecular genetic mechanisms underlying plant development over the last 15 years has pinpointed some of the key genes involved in the transition to flowering and flower development. Such investigations have also uncovered mutations which reduce floral maintenance or alter the balance between vegetative and floral features of the plant. How this information contributes to an understanding of floral reversion is assessed here. One issue that arises is whether floral commitment (defined as the ability to continue flowering when inductive conditions no longer exist) is a developmental switch affecting the whole plant or is a mechanism which assigns autonomy to individual meristems. A related question is whether floral or vegetative development is the underlying default pathway of the plant. This review begins by considering how studies of flowering in Arabidopsis thaliana have aided understanding of mechanisms of floral maintenance. Arabidopsis has not been found to revert to leaf production in any of the conditions or genetic backgrounds analysed to date. A clear-cut reversion to leaf production has, however, been described in Impatiens balsamina. It is proposed that a single gene controls whether Impatiens reverts or can maintain flowering when inductive conditions are removed, and it is inferred that this gene functions to control the synthesis or transport of a leaf-generated signal. But it is also argued that the susceptibility of Impatiens to reversion is a consequence of the meristem-based mechanisms controlling development of the flower in this species. Thus, in Impatiens, a leaf-derived signal is critical for completion of flowering and can be considered to be the basis of a plant-wide floral commitment that is achieved without accompanying meristem autonomy. The evidence, derived from in vitro and other studies, that similar mechanisms operate in other species is assessed. It is concluded that most species (including Arabidopsis) are less prone to reversion because signals from the leaf are less ephemeral, and the pathways driving flower development have a high level of redundancy that generates meristem autonomy even when leaf-derived signals are weak. This gives stability to the flowering process, even where its initiation is dependent on environmental cues. On this interpretation, Impatiens reversion appears as an anomaly resulting from an unusual combination of leaf signalling and meristem regulation. Nevertheless, it is shown that the ability to revert can serve a function in the life history strategy (perenniality) or reproductive habit (pseudovivipary) of many plants. In these instances reversion has been assimilated into regular plant development and plays a crucial role there.

Proteomic analysis of plasma membrane vesicles isolated from the rat renal cortex

Polarized epithelial cells are responsible for the vectorial transport of solutes and have a key role in maintaining body fluid and electrolyte homeostasis. Such cells contain structurally and functionally distinct plasma membrane domains. Brush border and basolateral membranes of renal and intestinal epithelial cells can be separated using a number of different separation techniques, which allow their different transport functions and receptor expressions to be studied. In this communication, we report a proteomic analysis of these two membrane segments, apical and basolateral, obtained from the rat renal cortex isolated by two different methods: differential centrifugation and free-flow electrophoresis. The study was aimed at assessing the nature of the major proteins isolated by these two separation techniques. Two analytical strategies were used: separation by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) at the protein level or by cation-exchange high-performance liquid chromatography (HPLC) after proteolysis (i.e., at the peptide level). Proteolytic peptides derived from the proteins present in gel pieces or from HPLC fractions after proteolysis were sequenced by on-line liquid chromatography-tandem mass spectrometry (LC-MS/MS). Several hundred proteins were identified in each membrane section. In addition to proteins known to be located at the apical and basolateral membranes, several novel proteins were also identified. In particular, a number of proteins with putative roles in signal transduction were identified in both membranes. To our knowledge, this is the first reported study to try and characterize the membrane proteome of polarized epithelial cells and to provide a data set of the most abundant proteins present in renal proximal tubule cell membranes.

Factors controlling the distribution of ozone in the West African lower troposphere during the AMMA (African Monsoon Multidisciplinary Analysis) wet season campaign

Ozone and its precursors were measured on board the Facility for Airborne Atmospheric Measurements (FAAM) BAe 146 Atmospheric Research Aircraft during the monsoon season 2006 as part of the African Monsoon Multidisciplinary Analysis (AMMA) campaign. One of the main features observed in the west African boundary layer is the increase of the ozone mixing ratios from 25 ppbv over the forested area (south of 12 degrees N) up to 40 ppbv over the Sahelian area. We employ a two-dimensional ( latitudinal versus vertical) meteorological model coupled with an O-3-NOx-VOC chemistry scheme to simulate the distribution of trace gases over West Africa during the monsoon season and to analyse the processes involved in the establishment of such a gradient. Including an additional source of NO over the Sahelian region to account for NO emitted by soils we simulate a mean NOx concentration of 0.7 ppbv at 16 degrees N versus 0.3 ppbv over the vegetated region further south in reasonable agreement with the observations. As a consequence, ozone is photochemically produced with a rate of 0.25 ppbv h(-1) over the vegetated region whilst it reaches up to 0.75 ppbv h(-1) at 16 degrees N. We find that the modelled gradient is due to a combination of enhanced deposition to vegetation, which decreases the ozone levels by up to 11 pbbv, and the aforementioned enhanced photochemical production north of 12 degrees N. The peroxy radicals required for this enhanced production in the north come from the oxidation of background CO and CH4 as well as from VOCs. Sensitivity studies reveal that both the background CH4 and partially oxidised VOCs, produced from the oxidation of isoprene emitted from the vegetation in the south, contribute around 5-6 ppbv to the ozone gradient. These results suggest that the northward transport of trace gases by the monsoon flux, especially during nighttime, can have a significant, though secondary, role in determining the ozone gradient in the boundary layer. Convection, anthropogenic emissions and NO produced from lightning do not contribute to the establishment of the discussed ozone gradient.

Life cycle analysis of UK coal fired power plants

This paper examines the life cycle GHG emissions from existing UK pulverized coal power plants. The life cycle of the electricity Generation plant includes construction, operation and decommissioning. The operation phase is extended to upstream and downstream processes. Upstream processes include the mining and transport of coal including methane leakage and the production and transport of limestone and ammonia, which are necessary for flue gas clean up. Downstream processes, on the other hand, include waste disposal and the recovery of land used for surface mining. The methodology used is material based process analysis that allows calculation of the total emissions for each process involved. A simple model for predicting the energy and material requirements of the power plant is developed. Preliminary calculations reveal that for a typical UK coal fired plant, the life cycle emissions amount to 990 g CO2-e/kWh of electricity generated, which compares well with previous UK studies. The majority of these emissions result from direct fuel combustion (882 g/kWh 89%) with methane leakage from mining operations accounting for 60% of indirect emissions. In total, mining operations (including methane leakage) account for 67.4% of indirect emissions, while limestone and other material production and transport account for 31.5%. The methodology developed is also applied to a typical IGCC power plant. It is found that IGCC life cycle emissions are 15% less than those from PC power plants. Furthermore, upon investigating the influence of power plant parameters on life cycle emissions, it is determined that, while the effect of changing the load factor is negligible, increasing efficiency from 35% to 38% can reduce emissions by 7.6%. The current study is funded by the UK National Environment Research Council (NERC) and is undertaken as part of the UK Carbon Capture and Storage Consortium (UKCCSC). Future work will investigate the life cycle emissions from other power generation technologies with and without carbon capture and storage. The current paper reveals that it might be possible that, when CCS is employed. the emissions during generation decrease to a level where the emissions from upstream processes (i.e. coal production and transport) become dominant, and so, the life cycle efficiency of the CCS system can be significantly reduced. The location of coal, coal composition and mining method are important in determining the overall impacts. In addition to studying the net emissions from CCS systems, future work will also investigate the feasibility and technoeconomics of these systems as a means of carbon abatement.

High pressure induced water uptake characteristics of Thai glutinous rice

Glutinous rice (or sticky rice) has to be soaked in water over an extended period of time before cooking. Soaking provides some of the water needed for starch gelatinisation to occur during cooking. The extent of water uptake during soaking is known to be influenced by temperature. This paper explores the use of very high pressures up to 600 MPa to accelerate water uptake kinetics during soaking. Changes occurring in length, diameter and moisture content were determined as a function of soaking time, pressure and temperature. The results show that length and diameter are positively correlated with all three parameters. However, the expansion ratios are not very high: the maximum length expansion ratio observed was 1.2, while the maximum diameter expansion ratio was 1. 1. Given these low values, it was possible to model water uptake kinetics by using the well-known Fickian model applied to a finite cylinder, assuming uniform average dimensions and effective diffusion coefficient. The results showed that the overall rates of water uptake and the equilibrium moisture content increased with pressure and temperature. The effective diffusion coefficient, on the other hand, did not follow the same trend. Temperature influenced the effective diffusion coefficient below 300 MPa, but had a marginal effect at higher pressures. Moreover, the effective diffusion coefficient increased with temperature between 20 and 50 degrees C, but dropped at higher temperatures. This drop can be attributed to the gelatinisation of starch, which restricts the transport of water. Regardless, it is possible to increase the quantity of water absorbed by rice and the rate at which it is absorbed, by using high pressures and temperatures. (c) 2004 Elsevier Ltd. All rights reserved.

Export of dissolved organic carbon from an upland peatland during storm events: implications for flux estimates

Most of the dissolved organic carbon (DOC) exported from catchments is transported during storm events. Accurate assessments of DOC fluxes are essential to understand long-term trends in the transport of DOC from terrestrial to aquatic systems, and also the loss of carbon from peatlands to determine changes in the source/sink status of peatland carbon stores. However, many long-term monitoring programmes collect water samples at a frequency (e.g. weekly/monthly) less than the time period of a typical storm event (typically <1–2 days). As widespread observations in catchments dominated by organo-mineral soils have shown that both concentration and flux of DOC increases during storm events, lower frequency monitoring could result in substantial underestimation of DOC flux as the most dynamic periods of transport are missed. However, our intensive monitoring study in a UK upland peatland catchment showed a contrasting response to these previous studies. Our results showed that (i) DOC concentrations decreased during autumn storm events and showed a poor relationship with flow during other seasons; and that (ii) this decrease in concentrations during autumn storms caused DOC flux estimates based on weekly monitoring data to be over-estimated, rather than under-estimated, because of over rather than under estimation of the flow-weighted mean concentration used in flux calculations. However, as DOC flux is ultimately controlled by discharge volume, and therefore rainfall, and the magnitude of change in discharge was greater than the magnitude of decline in concentrations, DOC flux increased during individual storm events. The implications for long-term DOC trends are therefore contradictory, as increased rainfall could increase flux but cause an overall decrease in DOC concentrations from peatland streams. Care needs to be taken when interpreting long-term trends in DOC flux rather than concentration; as flux is calculated from discharge estimates, and discharge is controlled by rainfall, DOC flux and rainfall/discharge will always be well correlated.

Localization to the nucleolus is a common feature of coronavirus nucleoproteins and the protein may disrupt host cell division

The subcellular localization of transmissible gastroenteritis virus (TGEV) and mouse hepatitis virus (MHV) (group I and group II coronaviruses, respectively) nucleoproteins (N proteins) were examined by confocal microscopy. The proteins were shown to localize either to the cytoplasm alone or to the cytoplasm and a structure in the nucleus. This feature was confirmed to be the nucleolus by using specific antibodies to nucleolin, a major component of the nucleolus, and by confocal microscopy to image sections through a cell expressing N protein. These findings are consistent with our previous report for infectious bronchitis virus (group III coronavirus) (J. A. Hiscox et al., J. Virol. 75:506-512, 2001), indicating that nucleolar localization of the N protein is a common feature of the coronavirus family and is possibly of functional significance. Nucleolar localization signals were identified in the domain III region of the N protein from all three coronavirus groups, and this suggested that transport of N protein to the nucleus might be an active process. In addition, our results suggest that the N protein might function to disrupt cell division. Thus, we observed that approximately 30% of cells transfected with the N protein appeared to be undergoing cell division. The most likely explanation for this is that the N protein induced a cell cycle delay or arrest, most likely in the G2/M phase. In a fraction of transfected cells expressing coronavirus N proteins, we observed multinucleate cells and dividing cells with nucleoli (which are only present during interphase). These findings are consistent with the possible inhibition of cytokinesis in these cells.

Evaluating the ability of a numerical weather prediction model to forecast tracer concentrations during ETEX 2

In this paper the meteorological processes responsible for transporting tracer during the second ETEX (European Tracer EXperiment) release are determined using the UK Met Office Unified Model (UM). The UM predicted distribution of tracer is also compared with observations from the ETEX campaign. The dominant meteorological process is a warm conveyor belt which transports large amounts of tracer away from the surface up to a height of 4 km over a 36 h period. Convection is also an important process, transporting tracer to heights of up to 8 km. Potential sources of error when using an operational numerical weather prediction model to forecast air quality are also investigated. These potential sources of error include model dynamics, model resolution and model physics. In the UM a semi-Lagrangian monotonic advection scheme is used with cubic polynomial interpolation. This can predict unrealistic negative values of tracer which are subsequently set to zero, and hence results in an overprediction of tracer concentrations. In order to conserve mass in the UM tracer simulations it was necessary to include a flux corrected transport method. Model resolution can also affect the accuracy of predicted tracer distributions. Low resolution simulations (50 km grid length) were unable to resolve a change in wind direction observed during ETEX 2, this led to an error in the transport direction and hence an error in tracer distribution. High resolution simulations (12 km grid length) captured the change in wind direction and hence produced a tracer distribution that compared better with the observations. The representation of convective mixing was found to have a large effect on the vertical transport of tracer. Turning off the convective mixing parameterisation in the UM significantly reduced the vertical transport of tracer. Finally, air quality forecasts were found to be sensitive to the timing of synoptic scale features. Errors in the position of the cold front relative to the tracer release location of only 1 h resulted in changes in the predicted tracer concentrations that were of the same order of magnitude as the absolute tracer concentrations.

Evaluation of a Lagrangian box model using field measurements from EASE (Eastern Atlantic Summer Experiment) 1996

The Cambridge Tropospheric Trajectory model of Chemistry and Transport (CiTTyCAT), a Lagrangian chemistry model, has been evaluated using atmospheric chemical measurements collected during the East Atlantic Summer Experiment 1996 (EASE '96). This field campaign was part of the UK Natural Environment Research Council's (NERC) Atmospheric Chemistry Studies in the Oceanic Environment (ACSOE) programme, conducted at Mace Head, Republic of Ireland, during July and August 1996. The model includes a description of gas-phase tropospheric chemistry, and simple parameterisations for surface deposition, mixing from the free troposphere and emissions. The model generally compares well with the measurements and is used to study the production and loss of O3 under a variety of conditions. The mean difference between the hourly O3 concentrations calculated by the model and those measured is 0.6 ppbv with a standard deviation of 8.7 ppbv. Three specific air-flow regimes were identified during the campaign – westerly, anticyclonic (easterly) and south westerly. The westerly flow is typical of background conditions for Mace Head. However, on some occasions there was evidence of long-range transport of pollutants from North America. In periods of anticyclonic flow, air parcels had collected emissions of NOx and VOCs immediately before arriving at Mace Head, leading to O3 production. The level of calculated O3 depends critically on the precise details of the trajectory, and hence on the emissions into the air parcel. In several periods of south westerly flow, low concentrations of O3 were measured which were consistent with deposition and photochemical destruction inside the tropical marine boundary layer.

Report of the 2010 assessment of the scientific assessment panel: future ozone and its impact on surface UV

SCIENTIFIC SUMMARY Globally averaged total column ozone has declined over recent decades due to the release of ozone-depleting substances (ODSs) into the atmosphere. Now, as a result of the Montreal Protocol, ozone is expected to recover from the effects of ODSs as ODS abundances decline in the coming decades. However, a number of factors in addition to ODSs have led to and will continue to lead to changes in ozone. Discriminating between the causes of past and projected ozone changes is necessary, not only to identify the progress in ozone recovery from ODSs, but also to evaluate the effectiveness of climate and ozone protection policy options. Factors Affecting Future Ozone and Surface Ultraviolet Radiation • At least for the next few decades, the decline of ODSs is expected to be the major factor affecting the anticipated increase in global total column ozone. However, several factors other than ODS will affect the future evolution of ozone in the stratosphere. These include changes in (i) stratospheric circulation and temperature due to changes in long-lived greenhouse gas (GHG) abundances, (ii) stratospheric aerosol loading, and (iii) source gases of highly reactive stratospheric hydrogen and nitrogen compounds. Factors that amplify the effects of ODSs on ozone (e.g., stratospheric aerosols) will likely decline in importance as ODSs are gradually eliminated from the atmosphere. • Increases in GHG emissions can both positively and negatively affect ozone. Carbon dioxide (CO2)-induced stratospheric cooling elevates middle and upper stratospheric ozone and decreases the time taken for ozone to return to 1980 levels, while projected GHG-induced increases in tropical upwelling decrease ozone in the tropical lower stratosphere and increase ozone in the extratropics. Increases in nitrous oxide (N2O) and methane (CH4) concentrations also directly impact ozone chemistry but the effects are different in different regions. • The Brewer-Dobson circulation (BDC) is projected to strengthen over the 21st century and thereby affect ozone amounts. Climate models consistently predict an acceleration of the BDC or, more specifically, of the upwelling mass flux in the tropical lower stratosphere of around 2% per decade as a consequence of GHG abundance increases. A stronger BDC would decrease the abundance of tropical lower stratospheric ozone, increase poleward transport of ozone, and could reduce the atmospheric lifetimes of long-lived ODSs and other trace gases. While simulations showing faster ascent in the tropical lower stratosphere to date are a robust feature of chemistry-climate models (CCMs), this has not been confirmed by observations and the responsible mechanisms remain unclear. • Substantial ozone losses could occur if stratospheric aerosol loading were to increase in the next few decades, while halogen levels are high. Stratospheric aerosol increases may be caused by sulfur contained in volcanic plumes entering the stratosphere or from human activities. The latter might include attempts to geoengineer the climate system by enhancing the stratospheric aerosol layer. The ozone losses mostly result from enhanced heterogeneous chemistry on stratospheric aerosols. Enhanced aerosol heating within the stratosphere also leads to changes in temperature and circulation that affect ozone. • Surface ultraviolet (UV) levels will not be affected solely by ozone changes but also by the effects of climate change and by air quality change in the troposphere. These tropospheric effects include changes in clouds, tropospheric aerosols, surface reflectivity, and tropospheric sulfur dioxide (SO2) and nitrogen dioxide (NO2). The uncertainties in projections of these factors are large. Projected increases in tropospheric ozone are more certain and may lead to reductions in surface erythemal (“sunburning”) irradiance of up to 10% by 2100. Changes in clouds may lead to decreases or increases in surface erythemal irradiance of up to 15% depending on latitude. Expected Future Changes in Ozone Full ozone recovery from the effects of ODSs and return of ozone to historical levels are not synonymous. In this chapter a key target date is chosen to be 1980, in part to retain the connection to previous Ozone Assessments. Noting, however, that decreases in ozone may have occurred in some regions of the atmosphere prior to 1980, 1960 return dates are also reported. The projections reported on in this chapter are taken from a recent compilation of CCM simulations. The ozone projections, which also form the basis for the UV projections, are limited in their representativeness of possible futures since they mostly come from CCM simulations based on a single GHG emissions scenario (scenario A1B of Emissions Scenarios. A Special Report of Working Group III of the Intergovernmental Panel on Climate Change, Cambridge University Press, 2000) and a single ODS emissions scenario (adjusted A1 of the previous (2006) Ozone Assessment). Throughout this century, the vertical, latitudinal, and seasonal structure of the ozone distribution will be different from what it was in 1980. For this reason, ozone changes in different regions of the atmosphere are considered separately. • The projections of changes in ozone and surface clear-sky UV are broadly consistent with those reported on in the 2006 Assessment. • The capability of making projections and attribution of future ozone changes has been improved since the 2006 Assessment. Use of CCM simulations from an increased number of models extending through the entire period of ozone depletion and recovery from ODSs (1960–2100) as well as sensitivity simulations have allowed more robust projections of long-term changes in the stratosphere and of the relative contributions of ODSs and GHGs to those changes. • Global annually averaged total column ozone is projected to return to 1980 levels before the middle of the century and earlier than when stratospheric halogen loading returns to 1980 levels. CCM projections suggest that this early return is primarily a result of GHG-induced cooling of the upper stratosphere because the effects of circulation changes on tropical and extratropical ozone largely cancel. Global (90°S–90°N) annually averaged total column ozone will likely return to 1980 levels between 2025 and 2040, well before the return of stratospheric halogens to 1980 levels between 2045 and 2060. • Simulated changes in tropical total column ozone from 1960 to 2100 are generally small. The evolution of tropical total column ozone in models depends on the balance between upper stratospheric increases and lower stratospheric decreases. The upper stratospheric increases result from declining ODSs and a slowing of ozone destruction resulting from GHG-induced cooling. Ozone decreases in the lower stratosphere mainly result from an increase in tropical upwelling. From 1960 until around 2000, a general decline is simulated, followed by a gradual increase to values typical of 1980 by midcentury. Thereafter, although total column ozone amounts decline slightly again toward the end of the century, by 2080 they are no longer expected to be affected by ODSs. Confidence in tropical ozone projections is compromised by the fact that simulated decreases in column ozone to date are not supported by observations, suggesting that significant uncertainties remain. • Midlatitude total column ozone is simulated to evolve differently in the two hemispheres. Over northern midlatitudes, annually averaged total column ozone is projected to return to 1980 values between 2015 and 2030, while for southern midlatitudes the return to 1980 values is projected to occur between 2030 and 2040. The more rapid return to 1980 values in northern midlatitudes is linked to a more pronounced strengthening of the poleward transport of ozone due to the effects of increased GHG levels, and effects of Antarctic ozone depletion on southern midlatitudes. By 2100, midlatitude total column ozone is projected to be above 1980 values in both hemispheres. • October-mean Antarctic total column ozone is projected to return to 1980 levels after midcentury, later than in any other region, and yet earlier than when stratospheric halogen loading is projected to return to 1980 levels. The slightly earlier return of ozone to 1980 levels (2045–2060) results primarily from upper stratospheric cooling and resultant increases in ozone. The return of polar halogen loading to 1980 levels (2050–2070) in CCMs is earlier than in empirical models that exclude the effects of GHG-induced changes in circulation. Our confidence in the drivers of changes in Antarctic ozone is higher than for other regions because (i) ODSs exert a strong influence on Antarctic ozone, (ii) the effects of changes in GHG abundances are comparatively small, and (iii) projections of ODS emissions are more certain than those for GHGs. Small Antarctic ozone holes (areas of ozone <220 Dobson units, DU) could persist to the end of the 21st century. • March-mean Arctic total column ozone is projected to return to 1980 levels two to three decades before polar halogen loading returns to 1980 levels, and to exceed 1980 levels thereafter. While CCM simulations project a return to 1980 levels between 2020 and 2035, most models tend not to capture observed low temperatures and thus underestimate present-day Arctic ozone loss such that it is possible that this return date is biased early. Since the strengthening of the Brewer-Dobson circulation through the 21st century leads to increases in springtime Arctic column ozone, by 2100 Arctic ozone is projected to lie well above 1960 levels. Uncertainties in Projections • Conclusions dependent on future GHG levels are less certain than those dependent on future ODS levels since ODS emissions are controlled by the Montreal Protocol. For the six GHG scenarios considered by a few CCMs, the simulated differences in stratospheric column ozone over the second half of the 21st century are largest in the northern midlatitudes and the Arctic, with maximum differences of 20–40 DU between the six scenarios in 2100. • There remain sources of uncertainty in the CCM simulations. These include the use of prescribed ODS mixing ratios instead of emission fluxes as lower boundary conditions, the range of sea surface temperatures and sea ice concentrations, missing tropospheric chemistry, model parameterizations, and model climate sensitivity. • Geoengineering schemes for mitigating climate change by continuous injections of sulfur-containing compounds into the stratosphere, if implemented, would substantially affect stratospheric ozone, particularly in polar regions. Ozone losses observed following large volcanic eruptions support this prediction. However, sporadic volcanic eruptions provide limited analogs to the effects of continuous sulfur emissions. Preliminary model simulations reveal large uncertainties in assessing the effects of continuous sulfur injections. Expected Future Changes in Surface UV. While a number of factors, in addition to ozone, affect surface UV irradiance, the focus in this chapter is on the effects of changes in stratospheric ozone on surface UV. For this reason, clear-sky surface UV irradiance is calculated from ozone projections from CCMs. • Projected increases in midlatitude ozone abundances during the 21st century, in the absence of changes in other factors, in particular clouds, tropospheric aerosols, and air pollutants, will result in decreases in surface UV irradiance. Clear-sky erythemal irradiance is projected to return to 1980 levels on average in 2025 for the northern midlatitudes, and in 2035 for the southern midlatitudes, and to fall well below 1980 values by the second half of the century. However, actual changes in surface UV will be affected by a number of factors other than ozone. • In the absence of changes in other factors, changes in tropical surface UV will be small because changes in tropical total column ozone are projected to be small. By the middle of the 21st century, the model projections suggest surface UV to be slightly higher than in the 1960s, very close to values in 1980, and slightly lower than in 2000. The projected decrease in tropical total column ozone through the latter half of the century will likely result in clear-sky surface UV remaining above 1960 levels. Average UV irradiance is already high in the tropics due to naturally occurring low total ozone columns and high solar elevations. • The magnitude of UV changes in the polar regions is larger than elsewhere because ozone changes in polar regions are larger. For the next decades, surface clear-sky UV irradiance, particularly in the Antarctic, will continue to be higher than in 1980. Future increases in ozone and decreases in clear-sky UV will occur at slower rates than those associated with the ozone decreases and UV increases that occurred before 2000. In Antarctica, surface clear-sky UV is projected to return to 1980 levels between 2040 and 2060, while in the Arctic this is projected to occur between 2020 and 2030. By 2100, October surface clear-sky erythemal irradiance in Antarctica is likely to be between 5% below to 25% above 1960 levels, with considerable uncertainty. This is consistent with multi-model-mean October Antarctic total column ozone not returning to 1960 levels by 2100. In contrast, by 2100, surface clear-sky UV in the Arctic is projected to be 0–10% below 1960 levels.

Sensitivity analysis of a cyanobacterial growth and movement model under two different flow regimes

Bloom-forming and toxin-producing cyanobacteria remain a persistent nuisance across the world. Modelling cyanobacterial behaviour in freshwaters is an important tool for understanding their population dynamics and predicting the location and timing of the bloom events in lakes, reservoirs and rivers. A new deterministic–mathematical model was developed, which simulates the growth and movement of cyanobacterial blooms in river systems. The model focuses on the mathematical description of the bloom formation, vertical migration and lateral transport of colonies within river environments by taking into account the major factors that affect the cyanobacterial bloom formation in rivers including light, nutrients and temperature. A parameter sensitivity analysis using a one-at-a-time approach was carried out. There were two objectives of the sensitivity analysis presented in this paper: to identify the key parameters controlling the growth and movement patterns of cyanobacteria and to provide a means for model validation. The result of the analysis suggested that maximum growth rate and day length period were the most significant parameters in determining the population growth and colony depth, respectively.

Moisture transport in mid-latitude cyclones

We discuss how synoptic-scale variability controls the transport of atmospheric water vapour by mid-latitude cyclones. Idealised simulations are used to investigate quantitatively what factors determine the magnitude of cyclone moisture transport. It is demonstrated that large-scale ascent on the warmconveyor belt and shallow cumulus convection are equally important for ventilating moisture from the boundary layer into the free troposphere, and that ventilated moisture can be transported large distances eastwards and polewards by the cyclone, before being returned to the surface as precipitation. The initial relative humidity is shown to have little affect on the ability of the cyclone to transport moisture, whilst the absolute temperature and meridional temperature gradient provide much stronger controls. Scaling arguments are presented to quantify the dependence of moisture transport on large-scale and boundary-layer parameters. It is shown that ventilation by shallow convection and warm-conveyor belt advection vary in the same way with changes to large-scale parameters. However, shallow convective ventilation has a much stronger dependence on boundary-layer parameters than warm-conveyor belt ventilation.

The climate response to CO2 of the Hadley Centre coupled AOGCM with and without flux adjustment

Coupled atmosphere‐ocean general circulation models have a tendency to drift away from a realistic climatology. The modelled climate response to an increase of CO2 concentration may be incorrect if the simulation of the current climate has significant errors, so in many models, including ours, the drift is counteracted by applying prescribed fluxes of heat and fresh water at the ocean‐atmosphere interface in addition to the calculated surface exchanges. Since the additional fluxes do not have a physical basis, the use of this technique of “flux adjustment” itself introduces some uncertainty in the simulated response to increased CO2. We find that the global‐average temperature response of our model to CO2 increasing at 1% per year is about 30% less without flux adjustment than with flux adjustment. The geographical patterns of the response are similar, indicating that flux adjustment is not causing any gross distortion. The reduced size of the response is due to more effective vertical transport of heat into the ocean, and a somewhat smaller climate sensitivity. Although the response in both cases lies within the generally accepted range for the climate sensitivity, systematic uncertainties of this size are clearly undesirable, and the best strategy for future development is to improve the climate model in order to reduce the need for flux adjustment.