A 1D microphysical cloud model for Earth,and Earth-like exoplanets: Liquid water and water ice clouds in the convective troposphere

One significant difference between the atmospheres of stars and exoplanets is the presence of condensed particles (clouds or hazes) in the atmosphere of the latter. In current 1D models clouds and hazes are treated in an approximate way by raising the surface albedo, or adopting measured Earth cloud properties. The former method introduces errors to the modeled spectra of the exoplanet, as clouds shield the lower atmosphere and thus modify the spectral features. The latter method works only for an exact Earth-analog, but it is challenging to extend to other planets.The main goal of this paper is to develop a self-consistent microphysical cloud model for 1D atmospheric codes, which can reproduce some observed properties of Earth, such as the average albedo, surface temperature, and global energy budget. The cloud model is designed to be computationally efficient, simple to implement, and applicable for a wide range of atmospheric parameters for planets in the habitable zone.We use a 1D, cloud-free, radiative–convective, and photochemical equilibrium code originally developed by Kasting, Pavlov, Segura, and collaborators as basis for our cloudy atmosphere model. The cloud model is based on models used by the meteorology community for Earth’s clouds. The free parameters of the model are the relative humidity and number density of condensation nuclei, and the precipitation efficiency. In a 1D model, the cloud coverage cannot be self-consistently determined, thus we treat it as a free parameter.We apply this model to Earth (aerosol number density 100 cm^?3, relative humidity 77%, liquid cloud fraction 40%, and ice cloud fraction 25%) and find that a precipitation efficiency of 0.8 is needed to reproduce the albedo, average surface temperature and global energy budget of Earth. We perform simulations to determine how the albedo and the climate of a planet is influenced by the free parameters of the cloud model. We find that the planetary climate is most sensitive to changes in the liquid water cloud fraction and precipitation efficiency.The advantage of our cloud model is that the cloud height and the droplet sizes are self-consistently calculated, both of which influence the climate and albedo of exoplanets.     

A comparative study of the modelling of cement hydration and cement-rock laboratory experiments

The use of cement and concrete as fracture grouting or as tunnel seals in a geological disposal facility for radioactive wastes creates potential issues concerning chemical reactivity. From a long-term safety perspective, it is desirable to be able model these interactions and changes quantitatively. The ‘Long-term Cement Studies’ (LCS) project was formulated with an emphasis on in situ field experiments with more realistic boundary conditions and longer time scales compared with former experiments. As part of the project programme, a modelling inter-comparison has been conducted, involving the modelling of two experiments describing cement hydration on one hand and cement-rock reaction on the other, with teams representing the NDA (UK), Posiva (Finland), and JAEA (Japan).This modelling exercise showed that the dominant reaction pathways in the two experiments are fairly well understood and are consistent between the different modelling teams, although significant differences existed amongst the precise parameterisation (e.g. reactive surface areas, dependences of rate upon pH, types of secondary minerals), and in some instances, processes (e.g. partition of alkali elements between solids and liquid during cement hydration; kinetic models of cement hydration). It was not conclusive if certain processes such as surface complexation (preferred by some modellers, but not by others) played a role in the cement-rock experiment or not. These processes appear to be more relevant at early times in the experiment and the evolution at longer timescales was not affected. The observed permeability profile with time could not be matched. The fact that no secondary minerals could be observed and that the precipitated mass calculated during the simulations is minor might suggest that the permeability reduction does not have a chemical origin, although a small amount of precipitates at pore throats could have a large impact on permeability.The modelling exercises showed that there is an interest in keeping the numerical models as simple as possible and trying to obtain a reasonable fit with a minimum of processes, minerals and parameters. However, up-scaling processes and model parameterisation to the timescales appropriate to repository safety assessment are of considerable concern. Future modelling exercises of this type should focus on a suitable natural or industrial analogue that might aid assessing mineral-fluid reactions at these longer timescales.     

Development of ST:REAM: a reach‐based stream power balance approach for predicting alluvial river channel adjustment

River channel sediment dynamics are important in integrated catchment management because changes in channel morphology resulting from sediment transfer have important implications for many river functions. However, application of existing approaches that account for catchment‐scale sediment dynamics has been limited, largely due to the difficulty in obtaining data necessary to support them. It is within this context that this study develops a new, reach‐based, stream power balance approach for predicting river channel adjustment. The new approach, named ST:REAM (sediment transport: reach equilibrium assessment method), is based upon calculations of unit bed area stream power (ω) derived from remotely sensed slope, width and discharge datasets. ST:REAM applies a zonation algorithm to values of ω that are spaced every 50 m along the catchment network in order to divide the branches of the network up into relatively homogenous reaches. ST:REAM then compares each reach's ω value with the ω of its upstream neighbour in order to predict whether or not the reach is likely to be either erosion dominated or deposition dominated. The paper describes the application of ST:REAM to the River Taff in South Wales, UK. This test study demonstrated that ST:REAM can be rapidly applied using remotely sensed data that are available across many river catchments and that ST:REAM correctly predicted the status of 87.5% of sites within the Taff catchment that field observations had defined as being either erosion or deposition dominated. However, there are currently a number of factors that limit the usefulness of ST:REAM, including inconsistent performance and the need for additional, resource intensive, data to be collected to both calibrate the model and aid interpretation of its results. Copyright © 2014 John Wiley & Sons, Ltd.     

Displaced non-Keplerian orbits using impulsive thrust

This paper investigates new families of displaced, highly non-Keplerian orbits in the two-body problem and artificial equilibria in the circular restricted three-body problem. The families of orbits presented extend prior work by using periodic impulses to generate displaced orbits rather than continuous thrust. The new displaced orbits comprise a sequence of individual Keplerian arcs whose intersection is continuous in position, with discontinuities in velocity removed using impulses. For frequent impulses the new families of orbits approximate continuous thrust non-Keplerian orbits found in previous studies. To generate approximations to artificial equilibria in the circular restricted three-body problem, periodic impulses are used to generate a sequence of connected three-body arcs which begin and terminate at a fixed position in the rotating frame of reference. Again, these families of orbits reduce to the families of artificial equilibria found using continuous thrust.     

A two-layer αω-dynamo model with dynamic feedback on the ω-effect

In an attempt to produce a simple representation of an interface dynamo, I examine a dynamo model composed of two one-dimensional (radially averaged) pseudo-spherical layers, one in the convection zone and possessing an α-effect, and the other in the tachocline and possessing an ω-effect. The two layers communicate by means of an analogue of Newton's law of cooling, and a dynamical back-reaction of the magnetic field on ω is provided. Extensive bifurcation diagrams are calculated for three separate values of η, the ratio of magnetic diffusivities of the two layers. I find recognizable similarities to, but also dramatic differences from, the comparable one-layer model examined by Roald &38; Thomas. In particular, the solar-like dynamo mode found previously is no longer stable in the two-layer version; in its place there is a sequence of periodic, quasi-periodic and chaotic modes probably created in a homoclinic bifurcation. These differences are important enough to provide support for the view that the solar dynamo cannot be meaningfully modelled in one dimension.     

The effect of sideband ratio on line intensity for Herschel/HIFI

The Heterodyne Instrument for the Far Infrared (HIFI) on board the Herschel Space Observatory is composed of a set of fourteen double sideband mixers. We discuss the general problem of the sideband ratio (SBR) determination and the impact of an imbalanced sideband ratio on the line calibration in double sideband heterodyne receivers. The HIFI SBR is determined from a combination of data taken during pre-launch gas cell tests and in-flight. The results and some of the calibration artefacts discovered in the gas cell test data are presented here along with some examples of how these effects appear in science data taken in orbit.     

High Resolution Radio Observations of Luminous IR Galaxies

High resolution radio observations provide a unique tool for investigation ofultraluminous infrared galaxies, due to the combination of high angularresolution and source transparency available at centimeter wavelengths. Thetypical angular size of the IR emitting region in the most luminous systems isof order 0.1 arcsec, while dust absorption appears to obscure this region fromdirect view at wavelengths shorter than 100 microns. Radio observationswith VLBI and MERLIN are revealing compact radio structures indicative of bothstarburst and AGN activity in these systems. Additional observations arerequired to clarify relationships between radio structures and the energygeneration mechanisms, and to realize the full potential of radio interferometryfor ULIRG studies.     

Artificial halo orbits for low-thrust propulsion spacecraft

We consider periodic halo orbits about artificial equilibrium points (AEP) near to the Lagrange points L ₁ and L ₂ in the circular restricted three body problem, where the third body is a low-thrust propulsion spacecraft in the Sun–Earth system. Although such halo orbits about artificial equilibrium points can be generated using a solar sail, there are points inside L ₁ and beyond L ₂ where a solar sail cannot be placed, so low-thrust, such as solar electric propulsion, is the only option to generate artificial halo orbits around points inaccessible to a solar sail. Analytical and numerical halo orbits for such low-thrust propulsion systems are obtained by using the Lindstedt Poincaré and differential corrector method respectively. Both the period and minimum amplitude of halo orbits about artificial equilibrium points inside L ₁ decreases with an increase in low-thrust acceleration. The halo orbits about artificial equilibrium points beyond L ₂ in contrast show an increase in period with an increase in low-thrust acceleration. However, the minimum amplitude first increases and then decreases after the thrust acceleration exceeds 0.415 mm/s². Using a continuation method, we also find stable artificial halo orbits which can be sustained for long integration times and require a reasonably small low-thrust acceleration 0.0593 mm/s².     

Reef-island topography and the vulnerability of atolls to sea-level rise

Low-lying reef islands on the rim of atolls are perceived as particularly vulnerable to the impacts of sea-level rise. Three effects are inferred: erosion of the shoreline, inundation of low-lying areas, and saline intrusion into the freshwater lens. Regional reconstruction of sea-level trends, supplementing the short observational instrumental record, indicates that monthly mean sea level is rising in the eastern Indian and western Pacific Oceans. This paper reviews the morphology and substrate characteristics of reef islands on Indo-Pacific atolls, and summarises their topography. On most atolls across this region, there is an oceanward ridge built by waves to a height of around 3 m above MSL; in a few cases these are topped by wind-blown dunes. The prominence of these ridges, together with radiocarbon dating and multi-temporal studies of shoreline position, indicate net accretion rather than long-term erosion on most of these oceanward shores. Less prominent lagoonward ridges occur, but their morphology and continuity are atoll-specific, being a function of the processes operating in each lagoon. Low-lying central areas are a feature of many islands, often locally excavated for production of taro. These lower-lying areas are already subject to inundation, which seems certain to increase as the sea rises. Tropical storms play an important role in the geomorphology of reef islands in those regions where they are experienced. Topographical differences, as well as features such as emergence of the reef flat and the stability of the substrate, mean that islands differ in terms of their susceptibility to sea-level rise. Further assessment of variations in shoreline vulnerability based on topography and substrate could form the basis for enhancing the natural resilience of these islands.