Multiple stationary solutions of an irradiated slab

A mathematical model describing the heat budget of an irradiated medium is introduced. The one-dimensional form of the equations and boundary conditions are presented and analysed. Heat transport at one face of the slab occurs by absorption (and reflection) of an incoming beam of short-wave radiation with a fraction of this radiation penetrating into the body of the slab, a diffusive heat flux in the slab and a prescribed incoming heat flux term. The other face of the slab is immersed in its own melt and is considered to be a free surface. Here, temperature continuity is prescribed and evolution of the surface is determined by a Stefan condition. These boundary conditions are flexible enough to describe a range of situations such as a laser shining on an opaque medium, or the natural environment of polar sea ice or lake ice. A two-stream radiation model is used which replaces the simple Beer’s law of radiation attenuation frequently used for semi-infinite domains. The stationary solutions of the governing equations are sought and it is found that there exists two possible stationary solutions for a given set of boundary conditions and a range of parameter choices. It is found that the existence of two stationary solutions is a direct result of the model of radiation absorption, due to its effect on the albedo of the medium. A linear stability analysis and numerical calculations indicate that where two stationary solutions exist, the solution corresponding to a larger thickness is always stable and the solution corresponding to a smaller thickness is unstable. Numerical simulations reveal that when there are two solutions, if the slab is thinner than the smaller stationary thickness it will melt completely, whereas if the slab is thicker than the smaller stationary thickness it will evolve toward the larger stationary thickness. These results indicate that other mechanisms (e.g. wave-induced agglomeration of crystals) are necessary to grow a slab from zero initial thickness in the parameter regime that yields two stationary solutions.

Monte Carlo field-theoretic simulations for melts of symmetric diblock copolymer

Monte Carlo field-theoretic simulations (MCFTS) are performed on melts of symmetric diblock copolymer for invariant polymerization indexes extending down to experimentally relevant values of N̅ ∼ 10^4. The simulations are performed with a fluctuating composition field, W_−(r), and a pressure field, W_+(r), that follows the saddle-point approximation. Our study focuses on the disordered-state structure function, S(k), and the order−disorder transition (ODT). Although shortwavelength fluctuations cause an ultraviolet (UV) divergence in three dimensions, this is readily compensated for with the use of an effective Flory−Huggins interaction parameter, χ_e. The resulting S(k) matches the predictions of renormalized one-loop (ROL) calculations over the full range of χ_eN and N̅ examined in our study, and agrees well with Fredrickson−Helfand (F−H) theory near the ODT. Consistent with the F−H theory, the ODT is discontinuous for finite N̅ and the shift in (χ_eN)_ODT follows the predicted N̅^−1/3 scaling over our range of N̅.

Simulations of the impact of orbital forcing and ocean on the Asian summer monsoon during the Holocene

The importance of orbital forcing and ocean impact on the Asian summer monsoon in the Holocene is investigated by comparing simulations with a fully coupled ocean–atmosphere general circulation model (FOAM) and with the atmospheric component of this model (FSSTAM) forced with prescribed modern sea-surface temperatures (SSTs). The results show: (1) the ocean amplifies the orbitally-induced increase in African monsoon precipitation, makes somewhat increase in southern India and damps the increase over the southeastern China. (2) The ocean could change the spatial distribution and local intensity of the orbitally-induced latitudinal atmospheric oscillation over the southeastern China and the subtropical western Pacific Ocean. (3) The orbital forcing mostly enhances the Asian summer precipitation in the FOAM and FSSTAM simulations. However, the ocean reduces the orbitally-induced summer precipitation and postpones the time of summer monsoon onset over the Asian monsoon region. (4) The orbital forcing considerably enhances the intensity of upper divergence, which is amplified by ocean further, over the eastern hemisphere. But the divergence is weaker in the FOAM simulations than in the FSSTAM simulations when the orbital forcing is fixed. (5) The orbital forcing can enhance the amplitude of precipitation variability over the subtropical Africa, the southeastern China and northwestern China, inversely, reduce it over central India and North China in the FOAM and FSSTAM simulations. The ocean obviously reduces the amplitude of precipitation variability over most of the Asian monsoon regions in the fixed orbital forcing simulations. (6) The areas characterized by increased summer precipitation in the long-term mean are mostly characterized by increased amplitude of short-term variability, whereas regions characterized by decreased precipitation are primarily characterized by decreased amplitude of short-term variability. However, the influences of orbital forcing or dynamical ocean on regional climate depend on the model.

Results of PMIP2 coupled simulations of the mid-Holocene and Last Glacial Maximum, Part 2: feedbacks with emphasis on the location of the ITCZ and mid- and high latitudes heat budget

A set of coupled ocean-atmosphere(-vegetation) simulations using state of the art climate models is now available for the Last Glacial Maximum (LGM) and the Mid-Holocene (MH) through the second phase of the Paleoclimate Modeling Intercomparison Project (PMIP2). Here we quantify the latitudinal shift of the location of the Intertropical Convergence Zone (ITCZ) in the tropical regions during boreal summer and the change in precipitation in the northern part of the ITCZ. For both periods the shift is more pronounced over the continents and East Asia. The maritime continent is the region where the largest spread is found between models. We also clearly establish that the larger the increase in the meridional temperature gradient in the tropical Atlantic during summer at the MH, the larger the change in precipitation over West Africa. The vegetation feedback is however not as large as found in previous studies, probably due to model differences in the control simulation. Finally, we show that the feedback from snow and sea-ice at mid and high latitudes contributes for half of the cooling in the Northern Hemisphere for the LGM, with the remaining being achieved by the reduced CO2 and water vapour in the atmosphere. For the MH the snow and albedo feedbacks strengthen the spring cooling and enhance the boreal summer warming, whereas water vapour reinforces the late summer warming. These feedbacks are modest in the Southern Hemisphere. For the LGM most of the surface cooling is due to CO2 and water vapour.

Results of PMIP2 coupled simulations of the mid-Holocene and Last Glacial Maximum, Part 1: experiments and large-scale features

A set of coupled ocean-atmosphere simulations using state of the art climate models is now available for the Last Glacial Maximum and the Mid-Holocene through the second phase of the Paleoclimate Modeling Intercomparison Project (PMIP2). This study presents the large-scale features of the simulated climates and compares the new model results to those of the atmospheric models from the first phase of the PMIP, for which sea surface temperature was prescribed or computed using simple slab ocean formulations. We consider the large-scale features of the climate change, pointing out some of the major differences between the different sets of experiments. We show in particular that systematic differences between PMIP1 and PMIP2 simulations are due to the interactive ocean, such as the amplification of the African monsoon at the Mid-Holocene or the change in precipitation in mid-latitudes at the LGM. Also the PMIP2 simulations are in general in better agreement with data than PMIP1 simulations.

Technical note: a simple procedure for removing temporal discontinuities in ERA-Interim upper stratospheric temperatures for use in nudged chemistry-climate model simulations

This note describes a simple procedure for removing unphysical temporal discontinuities in ERA-Interim upper stratospheric global mean temperatures in March 1985 and August 1998 that have arisen due to changes in satellite radiance data used in the assimilation. The derived temperature adjustments (offsets) are suitable for use in stratosphere-resolving chemistry-climate models that are nudged (relaxed) to ERA-Interim winds and temperatures. Simulations using a nudged version of the Canadian Middle Atmosphere Model (CMAM) show that the inclusion of the temperature adjustments produces temperature time series that are devoid of the large jumps in 1985 and 1998. Due to its strong temperature dependence, the simulated upper stratospheric ozone is also shown to vary smoothly in time, unlike in a nudged simulation without the adjustments where abrupt changes in ozone occur at the times of the temperature jumps. While the adjustments to the ERA-Interim temperatures remove significant artefacts in the nudged CMAM simulation, spurious transient effects that arise due to water vapour and persist for about 5 yr after the 1979 switch to ERA-Interim data are identified, underlining the need for caution when analysing trends in runs nudged to reanalyses.

The Landau free energy of hard ellipses obtained from microscopic simulations

Systems of two-dimensional hard ellipses of varying aspect ratios and packing fractions are studied by Monte Carlo simulations in the generalised canonical ensemble. From this microscopic model, we extract a coarse-grained macroscopic Landau-de Gennes free energy as a function of packing fraction and orientational order parameter. We separate the free energy into the ideal orientational entropy of non-interacting two-dimensional spins and an excess free energy associated with excluded volume interactions. We further explore the isotropic-nematic phase transition using our empirical expression for the free energy and find that the nature of the phase transition is continuous for the aspect ratios we studied.

The water balance of northern Africa during the mid-Holocene: an evaluation of the 6ka BP PMIP simulations

Runoff fields over northern Africa (10–25°N, 20°W–30°E) derived from 17 atmospheric general circulation models driven by identical 6 ka BP orbital forcing, sea surface temperatures, and CO2 concentration have been analyzed using a hydrological routing scheme (HYDRA) to simulate changes in lake area. The AGCM-simulated runoff produced six-fold differences in simulated lake area between models, although even the largest simulated changes considerably underestimate the observed changes in lake area during the mid-Holocene. The inter-model differences in simulated lake area are largely due to differences in simulated runoff (the squared correlation coefficient, R2, is 0.84). Most of these differences can be attributed to differences in the simulated precipitation (R2=0.83). The higher correlation between runoff and simulated lake area (R2=0.92) implies that simulated differences in evaporation have a contributory effect. When runoff is calculated using an offline land-surface scheme (BIOME3), the correlation between runoff and simulated lake area is (R2=0.94). Finally, the spatial distribution of simulated precipitation can exert an important control on the overall response.

Coupled versus uncoupled hindcast simulations of the Madden-Julian oscillation in the Year of Tropical Convection

This study investigates the impact of a full interactive ocean on daily initialised 15 day hindcasts of the Madden-Julian oscillation (MJO), measured against a Met Office Unified Model (MetUM) atmosphere control simulation (AGCM) during a 3 month period of the Year of Tropical Convection (YOTC). Results indicated that the coupled configuration (CGCM) extends MJO predictability over that of the AGCM, by up to 3-5 days. Propagation is improved in the CGCM, which we partly attribute to a more realistic phase relationship between sea surface temperature (SST) and convection. In addition, the CGCM demonstrates skill in representing downwelling oceanic Kelvin and Rossby waves which warm SSTs along their trajectory, with the potential to feed back on the atmosphere. These results imply that an ocean model capable of simulating internal ocean waves may be required to capture the full effect of air-sea coupling for the MJO.

Creep and recovery of magnetorheological fluids: experiments and simulations

A direct comparative study on the creep-recovery behavior of conventional MR fluids is carried out using magnetorheometry and particle-level simulations. Two particle concentrations are investigated (ϕ=0.05 and 0.30) at two different magnetic field strengths (53 kA•m-1 and 173 kA•m-1) in order to match the yield stresses developed in both systems for easier comparison. Simulations are mostly started with random initial structures with some additional tests of using preassembled single chains in the low concentration case. Experimental and simulation data are in good qualitative agreement. The results demonstrate three regions in the creep curves: i) In the initial viscoelastic region, the chain-like (at ϕ=0.05) or percolated three-dimensional network (at ϕ=0.30) structures fill up the gap and the average cluster size remains constant; ii) Above a critical strain of 10 %, in the retardation region, these structures begin to break and rearrange under shear. At large enough imposed stress values, they transform into thin sheet-like or thick lamellar structures, depending on the particle concentration; iii) Finally in the case of larger strain values either the viscosity diverges (at low stress values) or reaches a constant low value (at high stress values), showing a clear bifurcation behavior. For stresses below the bifurcation point the MR fluid is capable to recover the strain by a certain fraction. However, no recovery is observed for large stress values.

A novel ensemble method for retrieving properties of warm cloud in 3-D using ground-based scanning radar and zenith radiances

We present a novel method for retrieving high-resolution, three-dimensional (3-D) nonprecipitating cloud fields in both overcast and broken-cloud situations. The method uses scanning cloud radar and multiwavelength zenith radiances to obtain gridded 3-D liquid water content (LWC) and effective radius (re) and 2-D column mean droplet number concentration (Nd). By using an adaption of the ensemble Kalman filter, radiances are used to constrain the optical properties of the clouds using a forward model that employs full 3-D radiative transfer while also providing full error statistics given the uncertainty in the observations. To evaluate the new method, we first perform retrievals using synthetic measurements from a challenging cumulus cloud field produced by a large-eddy simulation snapshot. Uncertainty due to measurement error in overhead clouds is estimated at 20% in LWC and 6% in re, but the true error can be greater due to uncertainties in the assumed droplet size distribution and radiative transfer. Over the entire domain, LWC and re are retrieved with average error 0.05–0.08 g m-3 and ~2 μm, respectively, depending on the number of radiance channels used. The method is then evaluated using real data from the Atmospheric Radiation Measurement program Mobile Facility at the Azores. Two case studies are considered, one stratocumulus and one cumulus. Where available, the liquid water path retrieved directly above the observation site was found to be in good agreement with independent values obtained from microwave radiometer measurements, with an error of 20 g m-2.

Computer simulations of the structure of colloidal ferrofluids

The structure of a ferrofluid under the influence of an external magnetic field is expected to become anisotropic due to the alignment of the dipoles into the direction of the external field, and subsequently to the formation of particle chains due to the attractive head to tail orientations of the ferrofluid particles. Knowledge about the structure of a colloidal ferrofluid can be inferred from scattering data via the measurement of structure factors. We have used molecular-dynamics simulations to investigate the structure of both monodispersed and polydispersed ferrofluids. The results for the isotropic structure factor for monodispersed samples are similar to previous data by Camp and Patey that were obtained using an alternative Monte Carlo simulation technique, but in a different parameter region. Here we look in addition at bidispersed samples and compute the anisotropic structure factor by projecting the q vector onto the XY and XZ planes separately, when the magnetic field was applied along the z axis. We observe that the XY- plane structure factor as well as the pair distribution functions are quite different from those obtained for the XZ plane. Further, the two- dimensional structure factor patterns are investigated for both monodispersed and bidispersed samples under different conditions. In addition, we look at the scaling exponents of structure factors. Our results should be of value to interpret scattering data on ferrofluids obtained under the influence of an external field.

Probabilistic parameterisation of the surface mass balance--elevation feedback in regional climate model simulations of the Greenland ice sheet

We present a new parameterisation that relates surface mass balance (SMB: the sum of surface accumulation and surface ablation) to changes in surface elevation of the Greenland ice sheet (GrIS) for the MAR (Modèle Atmosphérique Régional: Fettweis, 2007) regional climate model. The motivation is to dynamically adjust SMB as the GrIS evolves, allowing us to force ice sheet models with SMB simulated by MAR while incorporating the SMB–elevation feedback, without the substantial technical challenges of coupling ice sheet and climate models. This also allows us to assess the effect of elevation feedback uncertainty on the GrIS contribution to sea level, using multiple global climate and ice sheet models, without the need for additional, expensive MAR simulations. We estimate this relationship separately below and above the equilibrium line altitude (ELA, separating negative and positive SMB) and for regions north and south of 77� N, from a set of MAR simulations in which we alter the ice sheet surface elevation. These give four “SMB lapse rates”, gradients that relate SMB changes to elevation changes. We assess uncertainties within a Bayesian framework, estimating probability distributions for each gradient from which we present best estimates and credibility intervals (CI) that bound 95% of the probability. Below the ELA our gradient estimates are mostly positive, because SMB usually increases with elevation: 0.56 (95% CI: −0.22 to 1.33) kgm−3 a−1 for the north, and 1.91 (1.03 to 2.61) kgm−3 a−1 for the south. Above the ELA, the gradients are much smaller in magnitude: 0.09 (−0.03 to 0.23) kgm−3 a−1 in the north, and 0.07 (−0.07 to 0.59) kgm−3 a−1 in the south, because SMB can either increase or decrease in response to increased elevation. Our statistically founded approach allows us to make probabilistic assessments for the effect of elevation feedback uncertainty on sea level projections (Edwards et al., 2014).

Aerosol microphysics simulations of the Mt.~Pinatubo eruption with the UM-UKCA composition-climate model

We use a stratosphere–troposphere composition–climate model with interactive sulfur chemistry and aerosol microphysics, to investigate the effect of the 1991 Mount Pinatubo eruption on stratospheric aerosol properties. Satellite measurements indicate that shortly after the eruption, between 14 and 23 Tg of SO2 (7 to 11.5 Tg of sulfur) was present in the tropical stratosphere. Best estimates of the peak global stratospheric aerosol burden are in the range 19 to 26 Tg, or 3.7 to 6.7 Tg of sulfur assuming a composition of between 59 and 77 % H2SO4. In light of this large uncertainty range, we performed two main simulations with 10 and 20 Tg of SO2 injected into the tropical lower stratosphere. Simulated stratospheric aerosol properties through the 1991 to 1995 period are compared against a range of available satellite and in situ measurements. Stratospheric aerosol optical depth (sAOD) and effective radius from both simulations show good qualitative agreement with the observations, with the timing of peak sAOD and decay timescale matching well with the observations in the tropics and mid-latitudes. However, injecting 20 Tg gives a factor of 2 too high stratospheric aerosol mass burden compared to the satellite data, with consequent strong high biases in simulated sAOD and surface area density, with the 10 Tg injection in much better agreement. Our model cannot explain the large fraction of the injected sulfur that the satellite-derived SO2 and aerosol burdens indicate was removed within the first few months after the eruption. We suggest that either there is an additional alternative loss pathway for the SO2 not included in our model (e.g. via accommodation into ash or ice in the volcanic cloud) or that a larger proportion of the injected sulfur was removed via cross-tropopause transport than in our simulations. We also critically evaluate the simulated evolution of the particle size distribution, comparing in detail to balloon-borne optical particle counter (OPC) measurements from Laramie, Wyoming, USA (41° N). Overall, the model captures remarkably well the complex variations in particle concentration profiles across the different OPC size channels. However, for the 19 to 27 km injection height-range used here, both runs have a modest high bias in the lowermost stratosphere for the finest particles (radii less than 250 nm), and the decay timescale is longer in the model for these particles, with a much later return to background conditions. Also, whereas the 10 Tg run compared best to the satellite measurements, a significant low bias is apparent in the coarser size channels in the volcanically perturbed lower stratosphere. Overall, our results suggest that, with appropriate calibration, aerosol microphysics models are capable of capturing the observed variation in particle size distribution in the stratosphere across both volcanically perturbed and quiescent conditions. Furthermore, additional sensitivity simulations suggest that predictions with the models are robust to uncertainties in sub-grid particle formation and nucleation rates in the stratosphere.