A multithickness sea ice model accounting for sliding friction

A multithickness sea ice model explicitly accounting for the ridging and sliding friction contributions to sea ice stress is developed. Both ridging and sliding contributions depend on the deformation type through functions adopted from the Ukita and Moritz kinematic model of floe interaction. In contrast to most previous work, the ice strength of a uniform ice sheet of constant ice thickness is taken to be proportional to the ice thickness raised to the 3/2 power, as is revealed in discrete element simulations by Hopkins. The new multithickness sea ice model for sea ice stress has been implemented into the Los Alamos “CICE” sea ice model code and is shown to improve agreement between model predictions and observed spatial distribution of sea ice thickness in the Arctic.

Assimilation of sea-ice concentration in a global climate model — physical and statistical aspects

We investigate the initialization of Northern-hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates significantly reduces assimilation error both in identical-twin experiments and when assimilating sea-ice observations, reducing the concentration error by a factor of four to six, and the thickness error by a factor of two. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that the strong dependence of thermodynamic ice growth on ice concentration necessitates an adjustment of mean ice thickness in the analysis update. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that proportional mean-thickness updates are superior to the other two methods considered and enable us to assimilate sea ice in a global climate model using simple Newtonian relaxation.

Assimilation of sea-ice concentration in a global climate model — physical and statistical aspects

We investigate the initialisation of Northern Hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates leads to good assimilation performance for sea-ice concentration and thickness, both in identical-twin experiments and when assimilating sea-ice observations. The simulation of other Arctic surface fields in the coupled model is, however, not significantly improved by the assimilation. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that an adjustment of mean ice thickness in the analysis update is essential to arrive at plausible state estimates. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that assimilation with proportional mean-thickness updates outperforms the other two methods considered. The method described here is very simple to implement, and gives results that are sufficiently good to be used for initialising sea ice in a global climate model for seasonal to decadal predictions.

Fold-change detection in a whole-pathway model of Escherichia coli chemotaxis

There has been recent interest in sensory systems that are able to display a response which is proportional to a fold change in stimulus concentration, a feature referred to as fold-change detection (FCD). Here, we demonstrate FCD in a recent whole-pathway mathematical model of Escherichia coli chemotaxis. FCD is shown to hold for each protein in the signalling cascade and to be robust to kinetic rate and protein concentration variation. Using a sensitivity analysis, we find that only variations in the number of receptors within a signalling team lead to the model not exhibiting FCD. We also discuss the ability of a cell with multiple receptor types to display FCD and explain how a particular receptor configuration may be used to elucidate the two experimentally determined regimes of FCD behaviour. All findings are discussed in respect of the experimental literature.

Global forecasting of thermal health hazards: the skill of probabilistic predictions of the Universal Thermal Climate Index (UTCI)

Although over a hundred thermal indices can be used for assessing thermal health hazards, many ignore the human heat budget, physiology and clothing. The Universal Thermal Climate Index (UTCI) addresses these shortcomings by using an advanced thermo-physiological model. This paper assesses the potential of using the UTCI for forecasting thermal health hazards. Traditionally, such hazard forecasting has had two further limitations: it has been narrowly focused on a particular region or nation and has relied on the use of single ‘deterministic’ forecasts. Here, the UTCI is computed on a global scale,which is essential for international health-hazard warnings and disaster preparedness, and it is provided as a probabilistic forecast. It is shown that probabilistic UTCI forecasts are superior in skill to deterministic forecasts and that despite global variations, the UTCI forecast is skilful for lead times up to 10 days. The paper also demonstrates the utility of probabilistic UTCI forecasts on the example of the 2010 heat wave in Russia.

The impact of two coupled cirrus microphysics-radiation parameterizations on the temperature and specific humidity biases in the tropical tropopause layer in a climate model

The impact of two different coupled cirrus microphysics-radiation parameterizations on the zonally averaged temperature and humidity biases in the tropical tropopause layer (TTL) of a Met Office climate model configuration is assessed. One parameterization is based on a linear coupling between a model prognostic variable, the ice mass mixing ratio, qi, and the integral optical properties. The second is based on the integral optical properties being parameterized as functions of qi and temperature, Tc, where the mass coefficients (i.e. scattering and extinction) are parameterized as nonlinear functions of the ratio between qi and Tc. The cirrus microphysics parameterization is based on a moment estimation parameterization of the particle size distribution (PSD), which relates the mass moment (i.e. second moment if mass is proportional to size raised to the power of 2 ) of the PSD to all other PSD moments through the magnitude of the second moment and Tc. This same microphysics PSD parameterization is applied to calculate the integral optical properties used in both radiation parameterizations and, thus, ensures PSD and mass consistency between the cirrus microphysics and radiation schemes. In this paper, the temperature-non-dependent and temperature-dependent parameterizations are shown to increase and decrease the zonally averaged temperature biases in the TTL by about 1 K, respectively. The temperature-dependent radiation parameterization is further demonstrated to have a positive impact on the specific humidity biases in the TTL, as well as decreasing the shortwave and longwave biases in the cloudy radiative effect. The temperature-dependent radiation parameterization is shown to be more consistent with TTL and global radiation observations.

Exact vortex solutions in an extended Skyrme-Faddeev model

We construct exact vortex solutions in 3+1 dimensions to a theory which is an extension, due to Gies, of the Skyrme-Faddeev model, and that is believed to describe some aspects of the low energy limit of the pure SU(2) Yang-Mills theory. Despite the efforts in the last decades those are the first exact analytical solutions to be constructed for such type of theory. The exact vortices appear in a very particular sector of the theory characterized by special values of the coupling constants, and by a constraint that leads to an infinite number of conserved charges. The theory is scale invariant in that sector, and the solutions satisfy Bogomolny type equations. The energy of the static vortex is proportional to its topological charge, and waves can travel with the speed of light along them, adding to the energy a term proportional to a U(1) No ether charge they create. We believe such vortices may play a role in the strong coupling regime of the pure SU(2) Yang-Mills theory.

The shoving model for the glass-former LiCl center dot 6H(2)O: A molecular dynamics simulation study

Molecular dynamics (MD) simulations of LiCl center dot 6H(2)O Showed that the diffusion coefficient D, and also I lie structural relaxation time , follow a power law at high temperatures, D(-1) proportional to (T - T(0))(-mu), with the same experimental parameters for viscosity (T(0) = 207 K, mu = 2.08). Decoupling between D and occurs at T(x) similar to 1.1 T(0). High frequency acoustic excitations for the LiCl center dot 6H(2)O model were obtained by the calculation of time correlation functions of mass current fluctuations. The temperature dependence of the instantaneous shear modulus, G,(T), was considered in the shoving model for supercooled liquids [J.C. Dyre, T. Christensen, N.B. Olsen, J. Non-Cryst. Solids 352 (2006) 4635] resulting in a linear relationship log (D(-1)) vs. G root T. (C) 2009 Elsevier B.V. All rights reserved.

Short-term temporal changes of soil carbon losses after tillage described by a first-order decay model

Tillage stimulates soil carbon (C) losses by increasing aeration, changing temperature and moisture conditions, and thus favoring microbial decomposition. In addition, soil aggregate disruption by tillage exposes once protected organic matter to decomposition. We propose a model to explain carbon dioxide (CO2) emission after tillage as a function of the no-till emission plus a correction due to the tillage disturbance. The model assumes that C in the readily decomposable organic matter follows a first-order reaction kinetics equation as: dC(sail)(t)/dt = -kC(soil)(t) and that soil C-CO2 emission is proportional to the C decay rate in soil, where C-soil(t) is the available labile soil C (g m(-2)) at any time (t). Emissions are modeled in terms soil C available to decomposition in the tilled and non-tilled plots, and a relationship is derived between no-till (F-NT) and tilled (F-Gamma) fluxes, which is: F-T = a1F(NT)e(-a2t), where t is time after tillage. Predicted and observed fluxes showed good agreement based on determination coefficient (R-2), index of agreement and model efficiency, with R-2 as high as 0.97. The two parameters included in the model are related to the difference between the decay constant (k factor) of tilled and no-till plots (a(2)) and also to the amount of labile carbon added to the readily decomposable soil organic matter due to tillage (a,). These two parameters were estimated in the model ranging from 1.27 and 2.60 (a(1)) and - 1.52 x 10(-2) and 2.2 x 10(-2) day(-1) (a(2)). The advantage is that temporal variability of tillage-induced emissions can be described by only one analytical function that includes the no-till emission plus an exponential term modulated by tillage and environmentally dependent parameters. (C) 2008 Elsevier B.V. All rights reserved.

Using of model updating considering damped frequency response function

The search for better performance in the structural systems has been taken to more refined models, involving the analysis of a growing number of details, which should be correctly formulated aiming at defining a representative model of the real system. Representative models demand a great detailing of the project and search for new techniques of evaluation and analysis. Model updating is one of this technologies, it can be used to improve the predictive capabilities of computer-based models. This paper presents a FRF-based finite element model updating procedure whose the updating variables are physical parameters of the model. It includes the damping effects in the updating procedure assuming proportional and none proportional damping mechanism. The updating parameters are defined at an element level or macro regions of the model. So, the parameters are adjusted locally, facilitating the physical interpretation of the adjusting of the model. Different tests for simulated and experimental data are discussed aiming at defining the characteristics and potentialities of the methodology.

Development of a rabbit's urethral sphincter deficiency animal model for anatomical-functional evaluation

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A simplified Fermi Accelerator Model under quadratic frictional force

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Can Drag Force Suppress Fermi Acceleration in a Bouncer Model?

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Non-uniform drag force on the Fermi accelerator model

Some dynamical properties of a particle suffering the action of a generic drag force are obtained for a dissipative Fermi Acceleration model. The dissipation is introduced via a viscous drag force, like a gas, and is assumed to be proportional to a power of the velocity: F alpha -nu(gamma). The dynamics is described by a two-dimensional nonlinear area-contracting mapping obtained via the solution of Newton's second law of motion. We prove analytically that the decay of high energy is given by a continued fraction which recovers the following expressions: (i) linear for gamma = 1; (ii) exponential for gamma = 2; and (iii) second-degree polynomial type for gamma = 1.5. Our results are discussed for both the complete version and the simplified version. The procedure used in the present paper can be extended to many different kinds of system, including a class of billiards problems.

A two band model for superconductivity: probing interband pair formation

We propose a two band model for superconductivity. It turns out that the simplest nontrivial case considers solely interband scattering, and both bands can be modeled as symmetric (around the Fermi level) and flat, thus each band is completely characterized by its half-band width Wn (n=1,2). A useful dimensionless parameter is d, proportional to W2 - W1. The case delta = 0 retrieves the conventional BCS model. We probe the specific heat, the ratio gap over critical temperature, the thermodynamic critical field and tunneling conductance as functions of d and temperature (from zero to Tc). We compare our results with experimental results for MgB2 and good quantitative agreement is obtained, indicating the relevance of interband coupling. Work in progress also considers the inclusion of band hybridization and general interband as well as intra-band scattering mechanisms.