A word of caution on calculating market-based minimum capital risk requirements

This paper demonstrates that the use of GARCH-type models for the calculation of minimum capital risk requirements (MCRRs) may lead to the production of inaccurate and therefore inefficient capital requirements. We show that this inaccuracy stems from the fact that GARCH models typically overstate the degree of persistence in return volatility. A simple modification to the model is found to improve the accuracy of MCRR estimates in both back- and out-of-sample tests. Given that internal risk management models are currently in widespread usage in some parts of the world (most notably the USA), and will soon be permitted for EC banks and investment firms, we believe that our paper should serve as a valuable caution to risk management practitioners who are using, or intend to use this popular class of models.

Automated identification and tracking of polar-cap plasma patches at solar minimum

A method of automatically identifying and tracking polar-cap plasma patches, utilising data inversion and feature-tracking methods, is presented. A well-established and widely used 4-D ionospheric imaging algorithm, the Multi-Instrument Data Assimilation System (MIDAS), inverts slant total electron content (TEC) data from ground-based Global Navigation Satellite System (GNSS) receivers to produce images of the free electron distribution in the polar-cap ionosphere. These are integrated to form vertical TEC maps. A flexible feature-tracking algorithm, TRACK, previously used extensively in meteorological storm-tracking studies is used to identify and track maxima in the resulting 2-D data fields. Various criteria are used to discriminate between genuine patches and "false-positive" maxima such as the continuously moving day-side maximum, which results from the Earth's rotation rather than plasma motion. Results for a 12-month period at solar minimum, when extensive validation data are available, are presented. The method identifies 71 separate structures consistent with patch motion during this time. The limitations of solar minimum and the consequent small number of patches make climatological inferences difficult, but the feasibility of the method for patches larger than approximately 500 km in scale is demonstrated and a larger study incorporating other parts of the solar cycle is warranted. Possible further optimisation of discrimination criteria, particularly regarding the definition of a patch in terms of its plasma concentration enhancement over the surrounding background, may improve results.

September Arctic sea-ice minimum predicted by spring melt-pond fraction

The area of Arctic September sea ice has diminished from about 7 million km2 in the 1990s to less than 5 million km2 in five of the past seven years, with a record minimum of 3.6 million km2 in 2012 (ref. 1). The strength of this decrease is greater than expected by the scientific community, the reasons for this are not fully understood, and its simulation is an on-going challenge for existing climate models2, 3. With growing Arctic marine activity there is an urgent demand for forecasting Arctic summer sea ice4. Previous attempts at seasonal forecasts of ice extent were of limited skill5, 6, 7, 8, 9. However, here we show that the Arctic sea-ice minimum can be accurately forecasted from melt-pond area in spring. We find a strong correlation between the spring pond fraction and September sea-ice extent. This is explained by a positive feedback mechanism: more ponds reduce the albedo; a lower albedo causes more melting; more melting increases pond fraction. Our results help explain the acceleration of Arctic sea-ice decrease during the past decade. The inclusion of our new melt-pond model10 promises to improve the skill of future forecast and climate models in Arctic regions and beyond.

Sparse density estimation on the multinomial manifold

A new sparse kernel density estimator is introduced based on the minimum integrated square error criterion for the finite mixture model. Since the constraint on the mixing coefficients of the finite mixture model is on the multinomial manifold, we use the well-known Riemannian trust-region (RTR) algorithm for solving this problem. The first- and second-order Riemannian geometry of the multinomial manifold are derived and utilized in the RTR algorithm. Numerical examples are employed to demonstrate that the proposed approach is effective in constructing sparse kernel density estimators with an accuracy competitive with those of existing kernel density estimators.

On-line Gaussian mixture density estimator for adaptive minimum bit-error-rate beamforming receivers

We develop an on-line Gaussian mixture density estimator (OGMDE) in the complex-valued domain to facilitate adaptive minimum bit-error-rate (MBER) beamforming receiver for multiple antenna based space-division multiple access systems. Specifically, the novel OGMDE is proposed to adaptively model the probability density function of the beamformer’s output by tracking the incoming data sample by sample. With the aid of the proposed OGMDE, our adaptive beamformer is capable of updating the beamformer’s weights sample by sample to directly minimize the achievable bit error rate (BER). We show that this OGMDE based MBER beamformer outperforms the existing on-line MBER beamformer, known as the least BER beamformer, in terms of both the convergence speed and the achievable BER.

Influence of inherent structure shear stress of supercooled liquids on their shear moduli

Configurations of supercooled liquids residing in their local potential minimum (i.e. in their inherent structure, IS) were found to support a non-zero shear stress. This IS stress was attributed to the constraint to the energy minimization imposed by boundary conditions, which keep size and shape of the simulation cell fixed. In this paper we further investigate the influence of these boundary conditions on the IS stress. We investigate its importance for the computation of the low frequency shear modulus of a glass obtaining a consistent picture for the low- and high frequency shear moduli over the full temperature range. Hence, we find that the IS stress corresponds to a non-thermal contribution to the fluctuation term in the Born-Green expression. This leads to an unphysical divergence of the moduli in the low temperature limit if no proper correction for this term is applied. Furthermore, we clarify the IS stress dependence on the system size and put its origin on a more formal basis.

Regional climate impacts of a possible future grand solar minimum

Any reduction in global mean near-surface temperature due to a future decline in solar activity is likely to be a small fraction of projected anthropogenic warming. However, variability in ultraviolet solar irradiance is linked to modulation of the Arctic and North Atlantic Oscillations, suggesting the potential for larger regional surface climate effects. Here, we explore possible impacts through two experiments designed to bracket uncertainty in ultraviolet irradiance in a scenario in which future solar activity decreases to Maunder Minimum-like conditions by 2050. Both experiments show regional structure in the wintertime response, resembling the North Atlantic Oscillation, with enhanced relative cooling over northern Eurasia and the eastern United States. For a high-end decline in solar ultraviolet irradiance, the impact on winter northern European surface temperatures over the late twenty-first century could be a significant fraction of the difference in climate change between plausible AR5 scenarios of greenhouse gas concentrations.

Possible impacts of a future grand solar minimum on climate: stratospheric and global circulation changes

It has been suggested that the Sun may evolve into a period of lower activity over the 21st century. This study examines the potential climate impacts of the onset of an extreme ‘Maunder Minimum like’ grand solar minimum using a comprehensive global climate model. Over the second half of the 21st century, the scenario assumes a decrease in total solar irradiance of 0.12% compared to a reference RCP8.5 experiment. The decrease in solar irradiance cools the stratopause (~1 hPa) in the annual and global mean by 1.4 K. The impact on global mean near-surface temperature is small (~−0.1 K), but larger changes in regional climate occur during the stratospheric dynamically active seasons. In Northern hemisphere (NH) winter-time, there is a weakening of the stratospheric westerly jet by up to ~3-4 m s1, with the largest changes occurring in January-February. This is accompanied by a deepening of the Aleutian low at the surface and an increase in blocking over northern Europe and the north Pacific. There is also an equatorward shift in the Southern hemisphere (SH) midlatitude eddy-driven jet in austral spring. The occurrence of an amplified regional response during winter and spring suggests a contribution from a top-down pathway for solar-climate coupling; this is tested using an experiment in which ultraviolet (200–320 nm) radiation is decreased in isolation of other changes. The results show that a large decline in solar activity over the 21st century could have important impacts on the stratosphere and regional surface climate.

Inferring the structure of the solar corona and inner heliosphere during the Maunder minimum using global thermodynamic magnetohydrodynamic simulations

Observations of the Sun’s corona during the space era have led to a picture of relatively constant, but cyclically varying solar output and structure. Longer-term, more indirect measurements, such as from 10Be, coupled by other albeit less reliable contemporaneous reports, however, suggest periods of significant departure from this standard. The Maunder Minimum was one such epoch where: (1) sunspots effectively disappeared for long intervals during a 70 yr period; (2) eclipse observations suggested the distinct lack of a visible K-corona but possible appearance of the F-corona; (3) reports of aurora were notably reduced; and (4) cosmic ray intensities at Earth were inferred to be substantially higher. Using a global thermodynamic MHD model, we have constructed a range of possible coronal configurations for the Maunder Minimum period and compared their predictions with these limited observational constraints. We conclude that the most likely state of the corona during—at least—the later portion of the Maunder Minimum was not merely that of the 2008/2009 solar minimum, as has been suggested recently, but rather a state devoid of any large-scale structure, driven by a photospheric field composed of only ephemeral regions, and likely substantially reduced in strength. Moreover, we suggest that the Sun evolved from a 2008/2009-like configuration at the start of the Maunder Minimum toward an ephemeral-only configuration by the end of it, supporting a prediction that we may be on the cusp of a new grand solar minimum.

The Maunder minimum (1645-1715) was indeed a grand minimum: a reassessment of multiple datasets

Aims. Although the time of the Maunder minimum (1645–1715) is widely known as a period of extremely low solar activity, it is still being debated whether solar activity during that period might have been moderate or even higher than the current solar cycle (number 24). We have revisited all existing evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. Methods. We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of the sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum. Results. The level of solar activity during the Maunder minimum is reassessed on the basis of all available datasets. Conclusions. We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely lower than during the Dalton minimum of around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected with a high confidence level.

Conditioning of the weak-constraint variational data assimilation problem for numerical weather prediction

4-Dimensional Variational Data Assimilation (4DVAR) assimilates observations through the minimisation of a least-squares objective function, which is constrained by the model flow. We refer to 4DVAR as strong-constraint 4DVAR (sc4DVAR) in this thesis as it assumes the model is perfect. Relaxing this assumption gives rise to weak-constraint 4DVAR (wc4DVAR), leading to a different minimisation problem with more degrees of freedom. We consider two wc4DVAR formulations in this thesis, the model error formulation and state estimation formulation. The 4DVAR objective function is traditionally solved using gradient-based iterative methods. The principle method used in Numerical Weather Prediction today is the Gauss-Newton approach. This method introduces a linearised `inner-loop' objective function, which upon convergence, updates the solution of the non-linear `outer-loop' objective function. This requires many evaluations of the objective function and its gradient, which emphasises the importance of the Hessian. The eigenvalues and eigenvectors of the Hessian provide insight into the degree of convexity of the objective function, while also indicating the difficulty one may encounter while iterative solving 4DVAR. The condition number of the Hessian is an appropriate measure for the sensitivity of the problem to input data. The condition number can also indicate the rate of convergence and solution accuracy of the minimisation algorithm. This thesis investigates the sensitivity of the solution process minimising both wc4DVAR objective functions to the internal assimilation parameters composing the problem. We gain insight into these sensitivities by bounding the condition number of the Hessians of both objective functions. We also precondition the model error objective function and show improved convergence. We show that both formulations' sensitivities are related to error variance balance, assimilation window length and correlation length-scales using the bounds. We further demonstrate this through numerical experiments on the condition number and data assimilation experiments using linear and non-linear chaotic toy models.

Data Assimilation and Inverse Methods in Terms of a Probabilistic Formulation

The weak-constraint inverse for nonlinear dynamical models is discussed and derived in terms of a probabilistic formulation. The well-known result that for Gaussian error statistics the minimum of the weak-constraint inverse is equal to the maximum-likelihood estimate is rederived. Then several methods based on ensemble statistics that can be used to find the smoother (as opposed to the filter) solution are introduced and compared to traditional methods. A strong point of the new methods is that they avoid the integration of adjoint equations, which is a complex task for real oceanographic or atmospheric applications. they also avoid iterative searches in a Hilbert space, and error estimates can be obtained without much additional computational effort. the feasibility of the new methods is illustrated in a two-layer quasigeostrophic model.

Sparse density estimation on multinomial manifold combining local component analysis

A new sparse kernel density estimator is introduced based on the minimum integrated square error criterion combining local component analysis for the finite mixture model. We start with a Parzen window estimator which has the Gaussian kernels with a common covariance matrix, the local component analysis is initially applied to find the covariance matrix using expectation maximization algorithm. Since the constraint on the mixing coefficients of a finite mixture model is on the multinomial manifold, we then use the well-known Riemannian trust-region algorithm to find the set of sparse mixing coefficients. The first and second order Riemannian geometry of the multinomial manifold are utilized in the Riemannian trust-region algorithm. Numerical examples are employed to demonstrate that the proposed approach is effective in constructing sparse kernel density estimators with competitive accuracy to existing kernel density estimators.