Optimization of a sea ice model using basinwide observations of Arctic sea ice thickness, extent, and velocity

A stand-alone sea ice model is tuned and validated using satellite-derived, basinwide observations of sea ice thickness, extent, and velocity from the years 1993 to 2001. This is the first time that basin-scale measurements of sea ice thickness have been used for this purpose. The model is based on the CICE sea ice model code developed at the Los Alamos National Laboratory, with some minor modifications, and forcing consists of 40-yr ECMWF Re-Analysis (ERA-40) and Polar Exchange at the Sea Surface (POLES) data. Three parameters are varied in the tuning process: Ca, the air–ice drag coefficient; P*, the ice strength parameter; and α, the broadband albedo of cold bare ice, with the aim being to determine the subset of this three-dimensional parameter space that gives the best simultaneous agreement with observations with this forcing set. It is found that observations of sea ice extent and velocity alone are not sufficient to unambiguously tune the model, and that sea ice thickness measurements are necessary to locate a unique subset of parameter space in which simultaneous agreement is achieved with all three observational datasets.

Influence of mass diffusion in sea ice dynamical models

The mixing of floes of different thickness caused by repeated deformation of the ice cover is modeled as diffusion, and the mass balance equation for sea ice accounting for mass diffusion is developed. The effect of deformational diffusion on the ice thickness balance is shown to reach 1% of the divergence effect, which describes ridging and lead formation. This means that with the same accuracy the mass balance equation can be written in terms of mean velocity rather than mean mass-weighted velocity, which one should correctly use for a multicomponent fluid such as sea ice with components identified by floe thickness. Mixing (diffusion) of sea ice also occurs because of turbulent variations in wind and ocean drags that are unresolved in models. Estimates of the importance of turbulent mass diffusion on the dynamic redistribution of ice thickness are determined using empirical data for the turbulent diffusivity. For long-time-scale prediction (≫5 days), where unresolved atmospheric motion may have a length scale on the order of the Arctic basin and the time scale is larger than the synoptic time scale of atmospheric events, turbulent mass diffusion can exceed 10% of the divergence effect. However, for short-time-scale prediction, for example, 5 days, the unresolved scales are on the order of 100 km, and turbulent diffusion is about 0.1% of the divergence effect. Because inertial effects are small in the dynamics of the sea ice pack, diffusive momentum transfer can be disregarded.

Efficient fully nonlinear data assimilation for geophysical fluid dynamics

A potential problem with Ensemble Kalman Filter is the implicit Gaussian assumption at analysis times. Here we explore the performance of a recently proposed fully nonlinear particle filter on a high-dimensional but simplified ocean model, in which the Gaussian assumption is not made. The model simulates the evolution of the vorticity field in time, described by the barotropic vorticity equation, in a highly nonlinear flow regime. While common knowledge is that particle filters are inefficient and need large numbers of model runs to avoid degeneracy, the newly developed particle filter needs only of the order of 10-100 particles on large scale problems. The crucial new ingredient is that the proposal density cannot only be used to ensure all particles end up in high-probability regions of state space as defined by the observations, but also to ensure that most of the particles have similar weights. Using identical twin experiments we found that the ensemble mean follows the truth reliably, and the difference from the truth is captured by the ensemble spread. A rank histogram is used to show that the truth run is indistinguishable from any of the particles, showing statistical consistency of the method.

Electron acceleration and parallel electric fields due to kinetic Alfvén waves in plasma with similar thermal and Alfvén speeds

We investigate electron acceleration due to shear Alfven waves in a collissionless plasma for plasma parameters typical of 4–5RE radial distance from the Earth along auroral field lines. Recent observational work has motivated this study, which explores the plasma regime where the thermal velocity of the electrons is similar to the Alfven speed of the plasma, encouraging Landau resonance for electrons in the wave fields. We use a self-consistent kinetic simulation model to follow the evolution of the electrons as they interact with a short-duration wave pulse, which allows us to determine the parallel electric field of the shear Alfven wave due to both electron inertia and electron pressure effects. The simulation demonstrates that electrons can be accelerated to keV energies in a modest amplitude sub-second period wave. We compare the parallel electric field obtained from the simulation with those provided by fluid approximations.

Power transfer in velocity-vortex acceleration

A theoretical model is presented of an electron acceleration-as-oscillator method derived from the work of Joseph Larmor unified with J. Clerk Maxwell’s theory of vorticity for the displacement of radiation into free-space at an antenna interface.

Tracking (poly)phenol components from raspberries in ileal fluid

The (poly)phenols in ileal fluid after ingestion of raspberries were analysed by targeted and non-targeted LC-MSn approaches. Targeted approaches identified major anthocyanin and ellagitannin components at varying recoveries and with considerable inter-individual variation. Non-targeted LC-MSn analysis using an Orbitrap mass spectrometer gave exact mass MS data which was sifted using a software program to select peaks that changed significantly after supplementation. This method confirmed the recovery of the targeted components but also identified novel raspberry-specific metabolites. Some components (including ellagitannin and previously unidentified proanthocyanidin derivatives) may have arisen from raspberry seeds that survived intact in ileal samples. Other components include potential breakdown products of anthocyanins, unidentified components and phenolic metabolites formed in either the gut epithelia or after absorption into the circulatory system and efflux back into the gut lumen. The possible physiological roles of the ileal metabolites in the large bowel are discussed.

Solution calorimetry as a tool for investigating drug interaction with intestinal fluid

Solution calorimetry offers a reproducible technique for measuring the enthalpy of solution (ΔsolH) of a solute dissolving into a solvent. The ΔsolH of two solutes, propranolol HCl and mannitol were determined in simulated intestinal fluid (SIF) solutions designed to model the fed and fasted states within the gut, and in Hanks’ balanced salt solution (HBSS) of varying pH. The bile salt and lipid within the SIF solutions formed mixed micelles. Both solutes exhibited endothermic reactions in all solvents. The ΔsolH for propranolol HCl in the SIF solutions differed from those in the HBSS and was lower in the fed state than the fasted state SIF solution, revealing an interaction between propranolol and the micellar phase in both SIF solutions. In contrast, for mannitol the ΔsolH was constant in all solutions indicating minimal interaction between mannitol and the micellar phases of the SIF solutions. In this study, solution calorimetry proved to be a simple method for measuring the enthalpy associated with the dissolution of model drugs in complex biological media such as SIF solutions. In addition, the derived power–time curves allowed the time taken for the powdered solutes to form solutions to be estimated.

Rheological and electrical properties of NaY zeolite electrorheological fluid

The relations between the rheological and electrical properties of NaY zeolite electrorheological fluid and its solid phase are studied. It is found that then exist complex relations between its electrical and theological properties. The temperature spectra of dielectric properties of the fluid under high AC electric field are strongly field strength dependent. The relation between the DC conductivity of the fluid and the exciting electric field is experimentally presented as log sigma =A+BE1/2, when A is a strong function, but B, a very weak function of temperature. The shear stress of the fluid under a fixed electric field and temperature decreases with shear rate. A relaxation time for the adsorbed charges is estimated to be about 0.3 to 6.6 s in the temperature range from 280 to 380 K. The relaxation time qualitatively corresponds to the shear rate at which the shear stress begins to drop. The time dependent leaking current of the ER fluids under DC electric field is also measured. The conductivity increase is mainly caused by the structure evolution of particles. The experimental results can he explained with the calculations of Davis (J. Appl. Phys. 81(1997) pp.1985-1991) and Martin (J. Chem. Phys. 110(1999) pp.4854-4866). It is predicted that the NaY zeolite ER fluid strength would get degraded slowly.

Lattice Boltzmann simulation of viscous fluid systems with elastic boundaries

A lattice Boltzmann model able to simulate viscous fluid systems with elastic and movable boundaries is proposed. By introducing the virtual distribution function at the boundary, the Galilean invariance is recovered for the full system. As examples of application, the how in elastic vessels is simulated with the pressure-radius relationship similar to that of the pulmonary blood vessels. The numerical results for steady how are in good agreement with the analytical prediction, while the simulation results for pulsative how agree with the experimental observation of the aortic flows qualitatively. The approach has potential application in the study of the complex fluid systems such as the suspension system as well as the arterial blood flow.

Structure and viscoelasticity of an electrorheological fluid in oscillatory shear: computer simulation investigation

The three-dimensional molecular dynamics simulation method has been used to study the dynamic responses of an electrorheological (ER) fluid in oscillatory shear. The structure and related viscoelastic behaviour of the fluid are found to be sensitive to the amplitude of the strain. With the increase of the strain amplitude, the structure formed by the particles changes from isolated columns to sheet-like structures which may be perpendicular or parallel to the oscillating direction. Along with the structure evolution, the field-induced moduli decrease significantly with an increase in strain amplitude. The viscoelastic behaviour of the structures obtained in the cases of different strain amplitudes was examined in the linear response regime and an evident structure dependence of the moduli was found. The reason for this lies in the anisotropy of the arrangement of the particles in these structures. Short-range interactions between the particles cannot be neglected in determining the viscoelastic behaviour of ER fluids at small strain amplitude, especially for parallel sheets. The simulation results were compared with available experimental data and good agreement was reached for most of them.

Influence of the size distribution of particles on the viscous property of an electrorheological fluid

The influence of the size distribution of particles on the viscous property of an electrorheological fluid has been investigated by the molecular dynamic simulation method. The shear stress of the fluid is found to decrease with the increase of the variance sigma(2) of the Gaussian distribution of the particle size, and then reach a steady value when sigma is larger than 0.5. This phenomenon is attributed to the influence of the particle size distribution on the dynamic structural evolution in the fluid as well as the strength of the different chain-like structures formed by the particles.

The importance of vertical velocity variability for estimates of the indirect aerosol effects

The activation of aerosols to form cloud droplets is dependent upon vertical velocities whose local variability is not typically resolved at the GCM grid scale. Consequently, it is necessary to represent the subgrid-scale variability of vertical velocity in the calculation of cloud droplet number concentration. This study uses the UK Chemistry and Aerosols community model (UKCA) within the Hadley Centre Global Environmental Model (HadGEM3), coupled for the first time to an explicit aerosol activation parameterisation, and hence known as UKCA-Activate. We explore the range of uncertainty in estimates of the indirect aerosol effects attributable to the choice of parameterisation of the subgrid-scale variability of vertical velocity in HadGEM-UKCA. Results of simulations demonstrate that the use of a characteristic vertical velocity cannot replicate results derived with a distribution of vertical velocities, and is to be discouraged in GCMs. This study focuses on the effect of the variance (σw2) of a Gaussian pdf (probability density function) of vertical velocity. Fixed values of σw (spanning the range measured in situ by nine flight campaigns found in the literature) and a configuration in which σw depends on turbulent kinetic energy are tested. Results from the mid-range fixed σw and TKE-based configurations both compare well with observed vertical velocity distributions and cloud droplet number concentrations. The radiative flux perturbation due to the total effects of anthropogenic aerosol is estimated at −1.9 W m−2 with σw = 0.1 m s−1, −2.1 W m−2 with σw derived from TKE, −2.25 W m−2 with σw = 0.4 m s−1, and −2.3 W m−2 with σw = 0.7 m s−1. The breadth of this range is 0.4 W m−2, which is comparable to a substantial fraction of the total diversity of current aerosol forcing estimates. Reducing the uncertainty in the parameterisation of σw would therefore be an important step towards reducing the uncertainty in estimates of the indirect aerosol effects. Detailed examination of regional radiative flux perturbations reveals that aerosol microphysics can be responsible for some climate-relevant radiative effects, highlighting the importance of including microphysical aerosol processes in GCMs.

Comment on “B_y fluctuations in the magnetosheath and azimuthal flow velocity transients in the dayside ionosphere” by Newell and Sibeck

Newell and Sibeck [1993] (hereafter N&S) list some objections to our interpretation of dayside auroral transients and associated azimuthal flow bursts in terms of pulsed reconnection [e.g. Lockwood et al., 1989; 1993a]. They present what they term an “apparently overlooked” alternative explanation in terms of steady reconnection and fluctuations in the magnitude of the By component of the magnetosheath field. The objections of N&S can all be answered by reference to our previous publications and their alternative explanation was only “overlooked” in so far as it fails to explain the observations. Here we discuss just some of the reasons why the objections of N&S are invalid, and then give reasons why the events are not simply due to magnetosheath |By| changes.

Dayside ionospheric convection changes in response to long-period interplanetary Magnetic field oscillations: Determination of the ionospheric phase velocity

Ground magnetic field perturbations recorded by the CANOPUS magnetometer network in the 7 to 13 MLT sector are used to examine how reconfigurations of the dayside polar ionospheric flow take place in response to north-south changes of the IMF. During the 6-hour interval in question IMF Bz oscillates between ±7 nT with about a 1-hour period. Corresponding variations in the ground magnetic disturbance are observed which we infer are due to changes in ionospheric flow. Cross correlation of the data obtained from two ground stations at 73.5° magnetic latitude, but separated by ∼2 hours in MLT, shows that changes in the flow are initiated in the prenoon sector (∼10 MLT) and then spread outward toward dawn and dusk with a phase speed of ∼5 km s−1 over the longitude range ∼8 to 12 MLT, slowing to ∼2 km s−1 outside this range. Cross correlating the data from these ground stations with IMP 8 IMF Bz records produces a MLT variation in the ground response delay relative to the IMF which is compatible with these deduced phase speeds. We interpret these observations in terms of the ionospheric response to the onset, expansion and decay of magnetic reconnection at the dayside magnetopause.