Deeper model endgame analysis

A reference model of Fallible Endgame Play has been implemented and exercised with the chess-engine WILHELM. Past experiments have demonstrated the value of the model and the robustness of decisions based on it: experiments agree well with a Markov Model theory. Here, the reference model is exercised on the well-known endgame KBBKN.

Model endgame analysis

A reference model of Fallible Endgame Play has been implemented and exercised with the chess engine WILHELM. Various experiments have demonstrated the value of the model and the robustness of decisions based on it. Experimental results have also been compared with the theoretical predictions of a Markov model of the endgame and found to be in close agreement.

Testing the limits of quasi-geostrophic theory: application to observed laboratory flows outside the quasi-geostrophic regime

We compare laboratory observations of equilibrated baroclinic waves in the rotating two-layer annulus, with numerical simulations from a quasi-geostrophic model. The laboratory experiments lie well outside the quasi-geostrophic regime: the Rossby number reaches unity; the depth-to-width aspect ratio is large; and the fluid contains ageostrophic inertia–gravity waves. Despite being formally inapplicable, the quasi-geostrophic model captures the laboratory flows reasonably well. The model displays several systematic biases, which are consequences of its treatment of boundary layers and neglect of interfacial surface tension and which may be explained without invoking the dynamical effects of the moderate Rossby number, large aspect ratio or inertia–gravity waves. We conclude that quasi-geostrophic theory appears to continue to apply well outside its formal bounds.

Interfacial structure of wild-type and mutant forms of puroindoline-b bound to DPPG monolayers

The interaction of wild-type puroindoline-b (Pin-b+) and two mutant forms having single residue substitutions (G46S or W44R) with L-alpha-dipalmitoylphosphatidyl-dl-glycerol (DPPG) as a Langmuir monolayer at the air/water interface was investigated by neutron reflectivity (NR) and Brewster angle microscopy (BAM). NR profiles were fitted using a three-layer model to enable differences in penetration of protein between the lipid headgroup and acyl regions to be determined. The data showed similar surface excesses for each of the three proteins at the interface; however, it was revealed that the depth of penetration of protein into the lipid region differed for each protein with Pin-b+ penetrating further into the acyl region of the lipid compared to the mutant forms of the protein that interacted with the headgroup region only. BAM images revealed that the domain structure of the DPPG monolayers was disrupted when Pin-b+ adsorption had reached equilibrium, suggesting protein penetration had led to compression of the lipid region. In contrast, the domain structure was unaffected by the W44R mutant, suggesting no change in compression of the lipid region and hence little or no penetration of protein into the lipid layer.

Neural networks for identification of nonlinear systems under random piecewise polynomial disturbances

The problem of identification of a nonlinear dynamic system is considered. A two-layer neural network is used for the solution of the problem. Systems disturbed with unmeasurable noise are considered, although it is known that the disturbance is a random piecewise polynomial process. Absorption polynomials and nonquadratic loss functions are used to reduce the effect of this disturbance on the estimates of the optimal memory of the neural-network model.

Evaluation of satellite-based and model re-analysis rainfall estimates for Uganda

The dependence of much of Africa on rain fed agriculture leads to a high vulnerability to fluctuations in rainfall amount. Hence, accurate monitoring of near-real time rainfall is particularly useful, for example in forewarning possible crop shortfalls in drought-prone areas. Unfortunately, ground based observations are often inadequate. Rainfall estimates from satellite-based algorithms and numerical model outputs can fill this data gap, however rigorous assessment of such estimates is required. In this case, three satellite based products (NOAA-RFE 2.0, GPCP-1DD and TAMSAT) and two numerical model outputs (ERA-40 and ERA-Interim) have been evaluated for Uganda in East Africa using a network of 27 rain gauges. The study focuses on the years 2001 to 2005 and considers the main rainy season (February to June). All data sets were converted to the same temporal and spatial scales. Kriging was used for the spatial interpolation of the gauge data. All three satellite products showed similar characteristics and had a high level of skill that exceeded both model outputs. ERA-Interim had a tendency to overestimate whilst ERA-40 consistently underestimated the Ugandan rainfall.

Fast northward energy transfer in the Atlantic due to Agulhas rings

The adiabatic transit time of wave energy radiated by an Agulhas ring released in the South Atlantic Ocean to the North Atlantic Ocean is investigated in a two-layer ocean model. Of particular interest is the arrival time of baroclinic energy in the northern part of the Atlantic, because it is related to variations in the meridional overturning circulation. The influence of the Mid-Atlantic Ridge is also studied, because it allows for the conversion from barotropic to baroclinic wave energy and the generation of topographic waves. Barotropic energy from the ring is present in the northern part of the model basin within 10 days. From that time, the barotropic energy keeps rising to attain a maximum 500 days after initiation. This is independent of the presence or absence of a ridge in the model basin. Without a ridge in the model, the travel time of the baroclinic signal is 1300 days. This time is similar to the transit time of the ring from the eastern to the western coast of the model basin. In the presence of the ridge, the baroclinic signal arrives in the northern part of the model basin after approximately 10 days, which is the same time scale as that of the barotropic signal. It is apparent that the ridge can facilitate the energy conversion from barotropic to baroclinic waves and the slow baroclinic adjustment can be bypassed. The meridional overturning circulation, parameterized in two ways as either a purely barotropic or a purely baroclinic phenomenon, also responds after 1300 days. The ring temporarily increases the overturning strength. Th presence of the ridge does not alter the time scales.

Selfishness and altruism can coexist when help is subject to diminishing returns

Altruism and selfishness are 30–50% heritable in man in both Western and non-Western populations. This genetically based variation in altruism and selfishness requires explanation. In non-human animals, altruism is generally directed towards relatives, and satisfies the condition known as Hamilton's rule. This nepotistic altruism evolves under natural selection only if the ratio of the benefit of receiving help to the cost of giving it exceeds a value that depends on the relatedness of the individuals involved. Standard analyses assume that the benefit provided by each individual is the same but it is plausible in some cases that as more individuals contribute, help is subject to diminishing returns. We analyse this situation using a single-locus two-allele model of selection in a diploid population with the altruistic allele dominant to the selfish allele. The analysis requires calculation of the relationship between the fitnesses of the genotypes and the frequencies of the genes. The fitnesses vary not only with the genotype of the individual but also with the distribution of phenotypes amongst the sibs of the individual and this depends on the genotypes of his parents. These calculations are not possible by direct fitness or ESS methods but are possible using population genetics. Our analysis shows that diminishing returns change the operation of natural selection and the outcome can now be a stable equilibrium between altruistic and selfish alleles rather than the elimination of one allele or the other. We thus provide a plausible genetic model of kin selection that leads to the stable coexistence in the same population of both altruistic and selfish individuals. This may explain reported genetic variation in altruism in man.

Nonlinear stability of multilayer quasi-geostrophic flow

New nonlinear stability theorems are derived for disturbances to steady basic flows in the context of the multilayer quasi-geostrophic equations. These theorems are analogues of Arnol’d's second stability theorem, the latter applying to the two-dimensional Euler equations. Explicit upper bounds are obtained on both the disturbance energy and disturbance potential enstrophy in terms of the initial disturbance fields. An important feature of the present analysis is that the disturbances are allowed to have non-zero circulation. While Arnol’d's stability method relies on the energy–Casimir invariant being sign-definite, the new criteria can be applied to cases where it is sign-indefinite because of the disturbance circulations. A version of Andrews’ theorem is established for this problem, and uniform potential vorticity flow is shown to be nonlinearly stable. The special case of two-layer flow is treated in detail, with particular attention paid to the Phillips model of baroclinic instability. It is found that the short-wave portion of the marginal stability curve found in linear theory is precisely captured by the new nonlinear stability criteria.

The coevolution of two phytoplankton species on a single resource: allelopathy as a pseudo-mixotrophy

Without the top-down effects and the external/physical forcing, a stable coexistence of two phytoplankton species under a single resource is impossible — a result well known from the principle of competitive exclusion. Here I demonstrate by analysis of a mathematical model that such a stable coexistence in a homogeneous media without any external factor would be possible, at least theoretically, provided (i) one of the two species is toxin producing thereby has an allelopathic effect on the other, and (ii) the allelopathic effect exceeds a critical level. The threshold level of allelopathy required for the coexistence has been derived analytically in terms of the parameters associated with the resource competition and the nutrient recycling. That the extra mortality of a competitor driven by allelopathy of a toxic species gives a positive feed back to the algal growth process through the recycling is explained. And that this positive feed back plays a pivotal role in reducing competition pressures and helping species succession in the two-species model is demonstrated. Based on these specific coexistence results, I introduce and explain theoretically the allelopathic effect of a toxic species as a ‘pseudo-mixotrophy’—a mechanism of ‘if you cannot beat them or eat them, just kill them by chemical weapons’. The impact of this mechanism of species succession by pseudo-mixotrophy in the form of alleopathy is discussed in the context of current understanding on straight mixotrophy and resource-species relationship among phytoplankton species.

Impact of atmospheric forcing on Antarctic continental shelf water masses

The Antarctic continental shelf seas feature a bimodal distribution of water mass temperature, with the Amundsen and Bellingshausen Seas flooded by Circumpolar Deep Water that is several degrees Celsius warmer than the cold shelf waters prevalent in the Weddell and Ross Seas. This bimodal distribution could be caused by differences in atmospheric forcing, ocean dynamics, ocean and ice feedbacks, or some combination of these factors. In this study, a highly simplified coupled sea ice–mixed layer model is developed to investigate the physical processes controlling this situation. Under regional atmospheric forcings and parameter choices the 10-yr simulations demonstrate a complete destratification of the Weddell Sea water column in winter, forming cold, relatively saline shelf waters, while the Amundsen Sea winter mixed layer remains shallower, allowing a layer of deep warm water to persist. Applying the Weddell atmospheric forcing to the Amundsen Sea model destratifies the water column after two years, and applying the Amundsen forcing to the Weddell Sea model results in a shallower steady-state winter mixed layer that no longer destratifies the water column. This suggests that the regional difference in atmospheric forcings alone is sufficient to account for the bimodal distribution in Antarctic shelf-sea temperatures. The model prediction of mixed layer depth is most sensitive to the air temperature forcing, but a switch in all forcings is required to prevent destratification of the Weddell Sea water column.

Impact of non-hydrostatic effects and trapped lee waves on mountain wave drag in directionally sheared flow

The orographic gravity wave drag produced in flow over an axisymmetric mountain when both vertical wind shear and non-hydrostatic effects are important was calculated using a semi-analytical two-layer linear model, including unidirectional or directional constant wind shear in a layer near the surface, above which the wind is constant. The drag behaviour is determined by partial wave reflection at the shear discontinuity, wave absorption at critical levels (both of which exist in hydrostatic flow), and total wave reflection at levels where the waves become evanescent (an intrinsically non-hydrostatic effect), which produces resonant trapped lee wave modes. As a result of constructive or destructive wave interference, the drag oscillates with the thickness of the constant-shear layer and the Richardson number within it (Ri), generally decreasing at low Ri and when the flow is strongly non-hydrostatic. Critical level absorption, which increases with the angle spanned by the wind velocity in the constant-shear layer, shields the surface from reflected waves, keeping the drag closer to its hydrostatic limit. While, for the parameter range considered here, the drag seldom exceeds this limit, a substantial drag fraction may be produced by trapped lee waves, particularly when the flow is strongly non-hydrostatic, the lower layer is thick and Ri is relatively high. In directionally sheared flows with Ri = O(1), the drag may be misaligned with the surface wind in a direction opposite to the shear, a behaviour which is totally due to non-trapped waves. The trapped lee wave drag, whose reaction force on the atmosphere is felt at low levels, may therefore have a distinctly different direction from the drag associated with vertically propagating waves, which acts on the atmosphere at higher levels.

A risk-based framework for assessing the effectiveness of stratospheric aerosol geoengineering

Geoengineering by stratospheric aerosol injection has been proposed as a policy response to warming from human emissions of greenhouse gases, but it may produce unequal regional impacts. We present a simple, intuitive risk-based framework for classifying these impacts according to whether geoengineering increases or decreases the risk of substantial climate change, with further classification by the level of existing risk from climate change from increasing carbon dioxide concentrations. This framework is applied to two climate model simulations of geoengineering counterbalancing the surface warming produced by a quadrupling of carbon dioxide concentrations, with one using a layer of sulphate aerosol in the lower stratosphere, and the other a reduction in total solar irradiance. The solar dimming model simulation shows less regional inequality of impacts compared with the aerosol geoengineering simulation. In the solar dimming simulation, 10% of the Earth’s surface area, containing 10% of its population and 11% of its gross domestic product, experiences greater risk of substantial precipitation changes under geoengineering than under enhanced carbon dioxide concentrations. In the aerosol geoengineering simulation the increased risk of substantial precipitation change is experienced by 42% of Earth’s surface area, containing 36% of its population and 60% of its gross domestic product.

The dynamics of Baroclinic Zonal jets

Multiple alternating zonal jets are a ubiquitous feature of planetary atmospheres and oceans. However, most studies to date have focused on the special case of barotropic jets. Here, the dynamics of freely evolving baroclinic jets are investigated using a two-layer quasigeostrophic annulus model with sloping topography. In a suite of 15 numerical simulations, the baroclinic Rossby radius and baroclinic Rhines scale are sampled by varying the stratification and root-mean-square eddy velocity, respectively. Small-scale eddies in the initial state evolve through geostrophic turbulence and accelerate zonally as they grow in horizontal scale, first isotropically and then anisotropically. This process leads ultimately to the formation of jets, which take about 2500 rotation periods to equilibrate. The kinetic energy spectrum of the equilibrated baroclinic zonal flow steepens from a −3 power law at small scales to a −5 power law near the jet scale. The conditions most favorable for producing multiple alternating baroclinic jets are large baroclinic Rossby radius (i.e., strong stratification) and small baroclinic Rhines scale (i.e., weak root-mean-square eddy velocity). The baroclinic jet width is diagnosed objectively and found to be 2.2–2.8 times larger than the baroclinic Rhines scale, with a best estimate of 2.5 times larger. This finding suggests that Rossby wave motions must be moving at speeds of approximately 6 times the turbulent eddy velocity in order to be capable of arresting the isotropic inverse energy cascade.

Recent lake-level and outflow variations at Lake Viljandi, Estonia: validation of a coupled lake-catchment modelling scheme for climate change studies

Changes in the depth of Lake Viljandi between 1940 and 1990 were simulated using a lake water and energy-balance model driven by standard monthly weather data. Catchment runoff was simulated using a one-dimensional hydrological model, with a two-layer soil, a single-layer snowpack, a simple representation of vegetation cover and similarly modest input requirements. Outflow was modelled as a function of lake level. The simulated record of lake level and outflow matched observations of lake-level variations (r = 0.78) and streamflow (r = 0.87) well. The ability of the model to capture both intra- and inter-annual variations in the behaviour of a specific lake, despite the relatively simple input requirements, makes it extremely suitable for investigations of the impacts of climate change on lake water balance.