Using seasonal hindcasts to understand the origin of the equatorial cold tongue bias in CGCMs and its impact on ENSO

The cold equatorial SST bias in the tropical Pacific that is persistent in many coupled OAGCMs severely impacts the fidelity of the simulated climate and variability in this key region, such as the ENSO phenomenon. The classical bias analysis in these models usually concentrates on multi-decadal to centennial time series needed to obtain statistically robust features. Yet, this strategy cannot fully explain how the models errors were generated in the first place. Here, we use seasonal re-forecasts (hindcasts) to track back the origin of this cold bias. As such hindcasts are initialized close to observations, the transient drift leading to the cold bias can be analyzed to distinguish pre-existing errors from errors responding to initial ones. A time sequence of processes involved in the advent of the final mean state errors can then be proposed. We apply this strategy to the ENSEMBLES-FP6 project multi-model hindcasts of the last decades. Four of the five AOGCMs develop a persistent equatorial cold tongue bias within a few months. The associated systematic errors are first assessed separately for the warm and cold ENSO phases. We find that the models are able to reproduce either El Niño or La Niña close to observations, but not both. ENSO composites then show that the spurious equatorial cooling is maximum for El Niño years for the February and August start dates. For these events and at this time of the year, zonal wind errors in the equatorial Pacific are present from the beginning of the simulation and are hypothesized to be at the origin of the equatorial cold bias, generating too strong upwelling conditions. The systematic underestimation of the mixed layer depth in several models can also amplify the growth of the SST bias. The seminal role of these zonal wind errors is further demonstrated by carrying out ocean-only experiments forced by the AOCGCMs daily 10-meter wind. In a case study, we show that for several models, this forcing is sufficient to reproduce the main SST error patterns seen after 1 month in the AOCGCM hindcasts.

A 27 day persistence model of near-Earth solar wind conditions: a long lead-time forecast and a benchmark for dynamical models

Geomagnetic activity has long been known to exhibit approximately 27 day periodicity, resulting from solar wind structures repeating each solar rotation. Thus a very simple near-Earth solar wind forecast is 27 day persistence, wherein the near-Earth solar wind conditions today are assumed to be identical to those 27 days previously. Effective use of such a persistence model as a forecast tool, however, requires the performance and uncertainty to be fully characterized. The first half of this study determines which solar wind parameters can be reliably forecast by persistence and how the forecast skill varies with the solar cycle. The second half of the study shows how persistence can provide a useful benchmark for more sophisticated forecast schemes, namely physics-based numerical models. Point-by-point assessment methods, such as correlation and mean-square error, find persistence skill comparable to numerical models during solar minimum, despite the 27 day lead time of persistence forecasts, versus 2–5 days for numerical schemes. At solar maximum, however, the dynamic nature of the corona means 27 day persistence is no longer a good approximation and skill scores suggest persistence is out-performed by numerical models for almost all solar wind parameters. But point-by-point assessment techniques are not always a reliable indicator of usefulness as a forecast tool. An event-based assessment method, which focusses key solar wind structures, finds persistence to be the most valuable forecast throughout the solar cycle. This reiterates the fact that the means of assessing the “best” forecast model must be specifically tailored to its intended use.

Decadal predictions of the cooling and freshening of the North Atlantic in the 1960s and the role of ocean circulation

In the 1960s North Atlantic sea surface temperatures (SST) cooled rapidly. The magnitude of the cooling was largest in the North Atlantic subpolar gyre (SPG), and was coincident with a rapid freshening of the SPG. Here we analyze hindcasts of the 1960s North Atlantic cooling made with the UK Met Office’s decadal prediction system (DePreSys), which is initialised using observations. It is shown that DePreSys captures—with a lead time of several years—the observed cooling and freshening of the North Atlantic SPG. DePreSys also captures changes in SST over the wider North Atlantic and surface climate impacts over the wider region, such as changes in atmospheric circulation in winter and sea ice extent. We show that initialisation of an anomalously weak Atlantic Meridional Overturning Circulation (AMOC), and hence weak northward heat transport, is crucial for DePreSys to predict the magnitude of the observed cooling. Such an anomalously weak AMOC is not captured when ocean observations are not assimilated (i.e. it is not a forced response in this model). The freshening of the SPG is also dominated by ocean salt transport changes in DePreSys; in particular, the simulation of advective freshwater anomalies analogous to the Great Salinity Anomaly were key. Therefore, DePreSys suggests that ocean dynamics played an important role in the cooling of the North Atlantic in the 1960s, and that this event was predictable.

The development of a data-driven application benchmarking approach to performance modelling

Performance modelling is a useful tool in the lifeycle of high performance scientific software, such as weather and climate models, especially as a means of ensuring efficient use of available computing resources. In particular, sufficiently accurate performance prediction could reduce the effort and experimental computer time required when porting and optimising a climate model to a new machine. In this paper, traditional techniques are used to predict the computation time of a simple shallow water model which is illustrative of the computation (and communication) involved in climate models. These models are compared with real execution data gathered on AMD Opteron-based systems, including several phases of the U.K. academic community HPC resource, HECToR. Some success is had in relating source code to achieved performance for the K10 series of Opterons, but the method is found to be inadequate for the next-generation Interlagos processor. The experience leads to the investigation of a data-driven application benchmarking approach to performance modelling. Results for an early version of the approach are presented using the shallow model as an example.

An evaluation of decadal probability forecasts from state-of-the-art climate models

While state-of-the-art models of Earth's climate system have improved tremendously over the last 20 years, nontrivial structural flaws still hinder their ability to forecast the decadal dynamics of the Earth system realistically. Contrasting the skill of these models not only with each other but also with empirical models can reveal the space and time scales on which simulation models exploit their physical basis effectively and quantify their ability to add information to operational forecasts. The skill of decadal probabilistic hindcasts for annual global-mean and regional-mean temperatures from the EU Ensemble-Based Predictions of Climate Changes and Their Impacts (ENSEMBLES) project is contrasted with several empirical models. Both the ENSEMBLES models and a “dynamic climatology” empirical model show probabilistic skill above that of a static climatology for global-mean temperature. The dynamic climatology model, however, often outperforms the ENSEMBLES models. The fact that empirical models display skill similar to that of today's state-of-the-art simulation models suggests that empirical forecasts can improve decadal forecasts for climate services, just as in weather, medium-range, and seasonal forecasting. It is suggested that the direct comparison of simulation models with empirical models becomes a regular component of large model forecast evaluations. Doing so would clarify the extent to which state-of-the-art simulation models provide information beyond that available from simpler empirical models and clarify current limitations in using simulation forecasting for decision support. Ultimately, the skill of simulation models based on physical principles is expected to surpass that of empirical models in a changing climate; their direct comparison provides information on progress toward that goal, which is not available in model–model intercomparisons.

Modelling primary producer interaction and composition: an example of a UK lowland river

Nutrient enrichment and drought conditions are major threats to lowland rivers causing ecosystem degradation and composition changes in plant communities. The controls on primary producer composition in chalk rivers are investigated using a new model and existing data from the River Frome (UK) to explore abiotic and biotic interactions. The growth and interaction of four primary producer functional groups (suspended algae, macrophytes, epiphytes, sediment biofilm) were successfully linked with flow, nutrients (N, P), light and water temperature such that the modelled biomass dynamics of the four groups matched that of the observed. Simulated growth of suspended algae was limited mainly by the residence time of the river rather than in-stream phosphorus concentrations. The simulated growth of the fixed vegetation (macrophytes, epiphytes, sediment biofilm) was overwhelmingly controlled by incoming solar radiation and light attenuation in the water column. Nutrients and grazing have little control when compared to the other physical controls in the simulations. A number of environmental threshold values were identified in the model simulations for the different producer types. The simulation results highlighted the importance of the pelagic–benthic interactions within the River Frome and indicated that process interaction defined the behaviour of the primary producers, rather than a single, dominant driver. The model simulations pose interesting questions to be considered in the next iteration of field- and laboratory based studies.

Validation of a priori CME arrival predictions made using real-time heliospheric imager observations

Between December 2010 and March 2013, volunteers for the Solar Stormwatch (SSW) Citizen Science project have identified and analyzed coronal mass ejections (CMEs) in the near real-time Solar Terrestrial Relations Observatory Heliospheric Imager observations, in order to make “Fearless Forecasts” of CME arrival times and speeds at Earth. Of the 60 predictions of Earth-directed CMEs, 20 resulted in an identifiable Interplanetary CME (ICME) at Earth within 1.5–6 days, with an average error in predicted transit time of 22 h, and average transit time of 82.3 h. The average error in predicting arrival speed is 151 km s−1, with an average arrival speed of 425km s−1. In the same time period, there were 44 CMEs for which there are no corresponding SSW predictions, and there were 600 days on which there was neither a CME predicted nor observed. A number of metrics show that the SSW predictions do have useful forecast skill; however, there is still much room for improvement. We investigate potential improvements by using SSW inputs in three models of ICME propagation: two of constant acceleration and one of aerodynamic drag. We find that taking account of interplanetary acceleration can improve the average errors of transit time to 19 h and arrival speed to 77 km s−1.

Optimization of the physical refining of sunflower oil concerning the final contents of trans-fatty acids

The final contents of total and individual trans-fatty acids of sunflower oil, produced during the deacidification step of physical refining were obtained using a computational simulation program that considered cis-trans isomerization reaction features for oleic, linoleic, and linolenic acids attached to the glycerol part of triacylglycerols. The impact of process variables, such as temperature and liquid flow rate, and of equipment configuration parameters, such as liquid height, diameter, and number of stages, that influence the retention time of the oil in the equipment was analyzed using the response-surface methodology (RSM). The computational simulation and the RSM results were used in two different optimization methods, aiming to minimize final levels of total and individual trans-fatty acids (trans-FA), while keeping neutral oil loss and final oil acidity at low values. The main goal of this work was to indicate that computational simulation, based on a careful modeling of the reaction system, combined with optimization could be an important tool for indicating better processing conditions in industrial physical refining plants of vegetable oils, concerning trans-FA formation.

Mass transport effects in the borohydride oxidation reaction-Influence of the residence time on the reaction onset and faradaic efficiency

The borohydride oxidation reaction (BOR) was studied on Pt and Au electrodes by cyclic voltammetry in dilute alkaline borohydride solutions (0.1 M NaOH + 10(-3) mol L(-1) NaBH(4)). More specifically, the electrodes were considered as either Vulcan XC72-supported Pt or Au (noted as Pt/C and Au/C, respectively) active layers or smooth Pt or Au surfaces, the latter possibly being covered by a layer of (non-metalized) Vulcan XC72 carbon powder. The BOR onset potential and the number of electrons (n(e-)) exchanged per BH(4)(-) anion (faradaic efficiency) were investigated for these electrodes, to determine whether the residence time of reaction intermediates (at the electrode surface or inside the porous layer) does influence the overall reaction pathway/completion. For the carbon-supported platinum, n(e-) strongly depends on the thickness of the active layer. While thin (ca. 0.5 mu m-thick) Pt/C active layers yield n(e-) < 4, thick layers (approximately 3 mu m) yield n(e-)approximate to 8, which can be ascribed to the sufficient residence time of the molecules formed within the active layer (H(2), by heterogeneous hydrolysis, or BOR intermediates) enabling further (near-complete) oxidation. This puts into evidence that not only the nature of the electrocatalyst is important to reach high BOR efficiency, but also the structure/thickness of the active layer. The same trend applies for Au/C active layers and for smooth Pt or Au surfaces covered with a layer of (inactive) Vulcan XC72. In addition, the BOR onset usually shifts negative when the reaction intermediates are trapped, which suggests that some of the intermediates are more easily oxidized than BH(4)(-) itself; based on literature data, BH(3)OH(-) species is a likely candidate. (C) 2011 Elsevier B.V. All rights reserved.

Load balancing solution and evaluation of F5 content switch equipment

The Thesis focused on hardware based Load balancing solution of web traffic through a load balancer F5 content switch. In this project, the implemented scenario for distributing HTTPtraffic load is based on different CPU usages (processing speed) of multiple member servers.Two widely used load balancing algorithms Round Robin (RR) and Ratio model (weighted Round Robin) are implemented through F5 load balancer. For evaluating the performance of F5 content switch, some experimental tests has been taken on implemented scenarios using RR and Ratio model load balancing algorithms. The performance is examined in terms of throughput (bits/sec) and Response time of member servers in a load balancing pool. From these experiments we have observed that Ratio Model load balancing algorithm is most suitable in the environment of load balancing servers with different CPU usages as it allows assigning the weight according to CPU usage both in static and dynamic load balancing of servers.

Accelerated aging of thick glass second surface silvered reflectors under sandstorm conditions

Concentrated solar power systems are expected to be sited in desert locations where the direct normal irradiation is above 1800 kWh/m2.year. These systems include large solar collector assemblies, which account for a significant share of the investment cost. Solarreflectors are the main components of these solar collector assemblies and dust/sand storms may affect their reflectance properties, either by soiling or by surface abrasion. While soiling can be reverted by cleaning, surface abrasion is a non reversible degradation.The aim of this project was to study the accelerated aging of second surface silvered thickglass solar reflectors under simulated sandstorm conditions and develop a multi-parametric model which relates the specular reflectance loss to dust/sand storm parameters: wind velocity, dust concentration and time of exposure. This project focused on the degradation caused by surface abrasion.Sandstorm conditions were simulated in a prototype environmental test chamber. Material samples (6cm x 6cm) were exposed to Arizona coarse test dust. The dust stream impactedthese material samples at a perpendicular angle. Both wind velocity and dust concentrationwere maintained at a stable level for each accelerated aging test. The total exposure time in the test chamber was limited to 1 hour. Each accelerated aging test was interrupted every 4 minutes to measure the specular reflectance of the material sample after cleaning.The accelerated aging test campaign had to be aborted prematurely due to a contamination of the dust concentration sensor. A robust multi-parametric degradation model could thus not be derived. The experimental data showed that the specular reflectance loss decreasedeither linearly or exponentially with exposure time, so that a degradation rate could be defined as a single modeling parameter. A correlation should be derived to relate this degradation rate to control parameters such as wind velocity and dust/sand concentration.The sandstorm chamber design would have to be updated before performing further accelerated aging test campaigns. The design upgrade should improve both the reliability of the test equipment and the repeatability of accelerated aging tests. An outdoor exposure test campaign should be launched in deserts to learn more about the intensity, frequencyand duration of dust/sand storms. This campaign would also serve to correlate the results of outdoor exposure tests with accelerated exposure tests in order to develop a robust service lifetime prediction model for different types of solar reflector materials.

Evaluation of a high temperature solar thermal seasonal borehole storage

A solar thermal system with seasonal borehole storage for heating of a residential area in Anneberg, Sweden, approximately 10 km north of Stockholm, has been in operation since late 2002. Originally, the project was part of the EU THERMIE project “Large-scale Solar Heating Systems for Housing Developments” (REB/0061/97) and was the first solar heating plant in Europe with borehole storage in rock not utilizing a heat pump. Earlier evaluations of the system show lower performance than the preliminary simulation study, with residents complaining of a high use of electricity for domestic hot water (DHW) preparation and auxiliary heating. One explanation mentioned in the earlier evaluations is that the borehole storage had not yet reached “steady state” temperatures at the time of evaluation. Many years have passed since then and this paper presents results from a new evaluation. The main aim of this work is to evaluate the current performance of the system based on several key figures, as well as on system function based on available measurement data. The analysis show that though the borehole storage now has reached a quasi-steady state and operates as intended, the auxiliary electricity consumption is much higher than the original design values largely due to high losses in the distribution network, higher heat loads as well as lower solar gains.

Radiation properties of coil-coated steel in building envelope surfaces and the influence on building thermal performance

Recent studies have shown that the optical properties of building exterior surfaces are important in terms of energy use and thermal comfort. While the majority of the studies are related to exterior surfaces, the radiation properties of interior surfaces are less thoroughly investigated. Development in the coil-coating industries has now made it possible to allocate different optical properties for both exterior and interior surfaces of steel-clad buildings. The aim of this thesis is to investigate the influence of surface radiation properties with the focus on the thermal emittance of the interior surfaces, the modeling approaches and their consequences in the context of the building energy performance and indoor thermal environment. The study consists of both numerical and experimental investigations. The experimental investigations include parallel field measurements on three similar test cabins with different interior and exterior surface radiation properties in Borlänge, Sweden, and two ice rink arenas with normal and low emissive ceiling in Luleå, Sweden. The numerical methods include comparative simulations by the use of dynamic heat flux models, Building Energy Simulation (BES), Computational Fluid Dynamics (CFD) and a coupled model for BES and CFD. Several parametric studies and thermal performance analyses were carried out in combination with the different numerical methods. The parallel field measurements on the test cabins include the air, surface and radiation temperatures and energy use during passive and active (heating and cooling) measurements. Both measurement and comparative simulation results indicate an improvement in the indoor thermal environment when the interior surfaces have low emittance. In the ice rink arenas, surface and radiation temperature measurements indicate a considerable reduction in the ceiling-to-ice radiation by the use of low emittance surfaces, in agreement with a ceiling-toice radiation model using schematic dynamic heat flux calculations. The measurements in the test cabins indicate that the use of low emittance surfaces can increase the vertical indoor air temperature gradients depending on the time of day and outdoor conditions. This is in agreement with the transient CFD simulations having the boundary condition assigned on the exterior surfaces. The sensitivity analyses have been performed under different outdoor conditions and surface thermal radiation properties. The spatially resolved simulations indicate an increase in the air and surface temperature gradients by the use of low emittance coatings. This can allow for lower air temperature at the occupied zone during the summer. The combined effect of interior and exterior reflective coatings in terms of energy use has been investigated by the use of building energy simulation for different climates and internal heat loads. The results indicate possible energy savings by the smart choice of optical properties on interior and exterior surfaces of the building. Overall, it is concluded that the interior reflective coatings can contribute to building energy savings and improvement of the indoor thermal environment. This can be numerically investigated by the choice of appropriate models with respect to the level of detail and computational load. This thesis includes comparative simulations at different levels of detail.

Brazilian indexing and inertial inflation: evidence from time-varying estimates of an inflation transfer function

This paper is a theoretica1 and empirica1 study of the re1ationship between indexing po1icy and feedback mechanisms in the inflationary adjustment process in Brazil. The focus of our study is on two policy issues: (1) did the Brazilian system of indexing of interest rates, the exchange rate, and wages make inflation so dependent on its own past values that it created a significant feedback process and inertia in the behaviour of inflation in and (2) was the feedback effect of past inf1ation upon itself so strong that dominated the effect of monetary/fiscal variables upon current inflation? This paper develops a simple model designed to capture several "stylized facts" of Brazi1ian indexing po1icy. Separate ru1es of "backward indexing" for interest rates, the exchange rate, and wages, reflecting the evolution of po1icy changes in Brazil, are incorporated in a two-sector model of industrial and agricultural prices. A transfer function derived irom this mode1 shows inflation depending on three factors: (1) past values of inflation, (2) monetary and fiscal variables, and (3) supply- .shock variables. The indexing rules for interest rates, the exchange rate, and wages place restrictions on the coefficients of the transfer function. Variations in the policy-determined parameters of the indexing rules imply changes in the coefficients of the transfer function for inflation. One implication of this model, in contrast to previous results derived in analytically simpler models of indexing, is that a higher degree of indexing does not make current inflation more responsive to current monetary shocks. The empirical section of this paper studies the central hypotheses of this model through estimation of the inflation transfer function with time-varying parameters. The results show a systematic non-random variation of the transfer function coefficients closely synchronized with changes in the observed values of the wage-indexing parameters. Non-parametric tests show the variation of the transfer function coefficients to be statistically significant at the time of the changes in wage indexing rules in Brazil. As the degree of indexing increased, the inflation feadback coefficients increased, while the effect of external price and agricultura shocs progressively increased and monetary effects progressively decreased.

The Effect of Residence Time on the Tensile Properties of Superelastic and Thermal Activated Ni-Ti Orthodontic Wires

Since the 1980s, different devices based on superelastic alloys have been developed to fulfill orthodontic applications. Particularly in the last decades several researches have been carried out to evaluate the mechanical behavior of Ni-Ti alloys, including their tensile, torsion and fatigue properties. However, studies regarding the dependence of elastic properties on residence time of Ni-Ti wires in the oral cavity are scarce. Such approach is essential since metallic alloys are submitted to mechanical stresses during orthodontic treatment as well as pH and temperature fluctuations. The goal of the present contribution is to provide elastic stress-strain results to guide the orthodontic choice between martensitic thermal activated and austenitic superelastic Ni-Ti alloys. From the point of view of an orthodontist, the selection of appropriate materials and the correct maintenance of the orthodontic apparatus are essential needs during clinical treatment. The present work evaluated the elastic behavior of Ni-Ti alloy wires with diameters varying from 0.014 to 0.020 inches, submitted to hysteresis tensile tests with 8% strain. Tensile tests were performed after periods of use of 1, 2 and 3 months in the oral cavity of patients submitted to orthodontic treatment. The results from the hysteresis tests allowed to exam the strain range covered by isostress lines upon loading and unloading, as well as the residual strain after unloading for both superelastic and thermal activated Ni-Ti wires. Superelastic Ni-Ti wires exhibited higher load isostress values compared to thermal activated wires. It was found that such differences in the load isostress values can increase with increasing residence time.