Long term storage of relativistic particles in the solar corona

A model is presented which shows that large numbers of energetic electrons (0.3-> 10 MeV) and protons (1–30 MeV) can be stored in the solar corona at altitudes around 3 × 10⁵ km for periods in excess of 5 days. Specific reference is made to the time period July 6–16 1968 as an excellent example of energetic solar particle storage. Time histories of interplanetary charged particle intensities observed by the IMP-4 and Pioneer 8 satellites are used to substantiate this contention. Detailed reference is also made to solar X-ray, optical and radio data obtained during the period in question, in addition to interplanetary magnetometer data. This model provides a unique solution to many hitherto unexplained solar particle events, and can also account for the lack of prompt particle emission from certain large solar flares recorded in the past.     

Brans-Dicke cosmology,II

We discuss the implications of the Brans-Dicke scalar-tensor theory for cosmology with particular emphasis on the primordial element abundances that would obtain. Two general classes of models are found. Models of one class expand through the nuclear burning stage slightly more rapidly than the general relativistic case: models of the other class may expand at any rate whatsoever. The first class of models yeilds primordial abundances of D, He³ and He⁴ in agreement with their general relativistic values if the present mass density is low. High-density cosmologies, however, would produce too much He⁴. The second class of models yields element abundances which are far too high unless the expansion rate was quite large: in this case no He⁴ at all is produced. Finally, we determine the rate of change of the constant of gravitationG at the present epoch. For all but a very small class of models is negative at the present epoch. Models with positive values of at the present epoch produce no primordial He⁴ whatsoever, and have ages significantly lower than the corresponding general relativistic ages.     

Soil–pile–structure interaction: experimental outcomes from shaking table tests

An effective way to study the complex seismic soil‐structure interaction phenomena is to investigate the response of physical scaled models in 1‐g or n‐g laboratory devices. The outcomes of an extensive experimental campaign carried out on scaled models by means of the shaking table of the Bristol Laboratory for Advanced Dynamics Engineering, University of Bristol, UK, are discussed in the present paper. The experimental model comprises an oscillator connected to a single or a group of piles embedded in a bi‐layer deposit. Different pile head conditions, that is free head and fixed head, several dynamic properties of the structure, including different masses at the top of the single degree of freedom system, excited by various input motions, e.g. white noise, sinedwells and natural earthquake strong motions recorded in Italy, have been tested. In the present work, the modal dynamic response of the soil–pile–structure system is assessed in terms of period elongation and system damping ratio. Furthermore, the effects of oscillator mass and pile head conditions on soil–pile response have been highlighted, when the harmonic input motions are considered. Copyright © 2015 John Wiley & Sons, Ltd.     

Statistical properties of aftershock rate decay: Implications for the assessment of continuing activity

Aftershock rates seem to follow a power law decay, but the assessment of the aftershock frequency immediately after an earthquake, as well as during the evolution of a seismic excitation remains a demand for the imminent seismic hazard. The purpose of this work is to study the temporal distribution of triggered earthquakes in short time scales following a strong event, and thus a multiple seismic sequence was chosen for this purpose. Statistical models are applied to the 1981 Corinth Gulf sequence, comprising three strong (M = 6.7, M = 6.5, and M = 6.3) events between 24 February and 4 March. The non-homogeneous Poisson process outperforms the simple Poisson process in order to model the aftershock sequence, whereas the Weibull process is more appropriate to capture the features of the short-term behavior, but not the most proper for describing the seismicity in long term. The aftershock data defines a smooth curve of the declining rate and a long-tail theoretical model is more appropriate to fit the data than a rapidly declining exponential function, as supported by the quantitative results derived from the survival function. An autoregressive model is also applied to the seismic sequence, shedding more light on the stationarity of the time series.     

One- and multi-component models of the upper photosphere based on molecular spectra

Non-LTE synthetic spectra derived from a detailed analysis of the formation of the CN (0, 0) λ13883 Å spectrum are compared with center-limb photoelectric spectra taken at Kitt Peak National Observatory. Kitt Peak National Observatory is operated by the Association of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation. Significant non-LTE effects are found and the Kurucz, Altrock-Cannon, Mount-Linsky II, and HSRA models are compared. We derive a solar carbon abundance of A _c =8.30±0.10 for the Mount-Linsky model and A _c =8.40±0.10 for the Altrock-Cannon model, compared to the HSRA value of A _c =8.55±0.10, assuming a nitrogen abundance of logA _N=7.93. In addition we specify the regions of formation for the CN(0, 0) 3883.35 Å bandhead at disc center and limb.     

EUV analysis of an active region

A sequence of extreme ultraviolet (EUV) spectroheliograms of McMath region No. 10283 were obtained by the Harvard College Observatory experiment on OSO-6. The lines Ovi λ1032 Mg × λ625, Si xii λ499 and Fe xvi λ 335 were used to determine coronal temperatures and densities above the active region. A comparison of theoretical and observed line ratios yielded coronal temperatures of 2.2 to 2.3 × 10⁶K above the active region and 2.0 to 2.1 × 10⁶K in the surrounding area. The temperatures derived from ratios involving the O vi intensities are systematically higher than the others. This is attributed to an error in the theoretical O vi intensities. The intensities observed above the limb are compared with intensities predicted with a simple model based on cylindrical geometry. The overall agreement shows that the assumption of an iso-thermal corona in hydrostatic equilibrium above the active region is a reasonable working hypothesis and that the adopted geometrical model for the electron density distribution is adequate.     

VHF double-peaked spectra at negative Doppler shifts in the morning sector of radio aurora

This paper presents a new Doppler spectral type of VHP (42 MHz) radio auroral backscatter. This spectrum, which has a double-peaked structure, was observed repeatedly during the morning sector of an exceptionally strong event (A_p = 48) and is due to irregularities moving northwards with quite different velocities. The stronger spectral component, which has a smaller Doppler shift, is centred in frequency at ～?130 Hz, corresponding to the ion-acoustic velocity range in the medium; the weak component, which has a greater frequency shift, usually is centred at about ?300 Hz (～ 1050 m s^?1). Evidence based on spectral analysis of sequential short time sequences shows that the spectral power alternates in time between the two distinct frequency bands where the peaks are located, suggesting that the double-peaked spectrum may result from two competing processes which cannot operate simultaneously. The possibility exists that the theoretical model proposed by Sato (1977), which predicts two different quasi-linear stabilization mechanisms for the two-stream instability, could explain the observed double-peaked spectral type.     

Tracing the Electron Density from the Corona to 1 au

We derive the electron density distribution in the ecliptic plane, from the corona to 1 AU, using observations from 13.8 MHz to a few kHz by the radio experiment WAVES aboard the spacecraft Wind. We concentrate on type III bursts whose trajectories intersect the spacecraft, as determined by the presence of burst-associated Langmuir waves, or by energetic electrons observed by the 3-D Plasma experiment. For these bursts we are able to determine the mode of emission, fundamental or harmonic, the electron density at 1 AU, the distance of emission regions along the spiral, and the time spent by the beams as they proceed from the low corona to 1 AU. For all of the bursts considered, the emission mode at burst onset was the fundamental; by contrast, in deriving many previous models, harmonic emission was assumed.By measuring the onset time of the burst at each frequency we are able to derive an electron density model all along the trajectory of the burst. Our density model, after normalizing the density at 1 AU to be ne(215 R0)=7.2 cm–3 (the average value at the minimum of solar activity when our measurements were made), is ne=3.3×105 r–2+4.1×106 r–4+8.0×107 r–6 cm–3, with r in units of R0. For other densities at 1 AU our result implies that the coefficients in the equation need to be multiplied by n _e (1 AU)/7.2.We compare this with existing models and those derived from direct, in-situ measurements (normalized to the same density at 1 AU) and find that it agrees very well with in-situ measurements and poorly with radio models based on apparent source positions or assumptions of the emission mode. One implication of our results is that isolated type III bursts do not usually propagate in dense regions of the corona and solar wind, as it is still sometimes assumed.     

Explicit integration of bounding surface model for the analysis of earthquake soil liquefaction

This paper presents a new plasticity model developed for the simulation of monotonic and cyclic loading of non‐cohesive soils and its implementation to the commercial finite‐difference code FLAC, using its User‐Defined‐Model (UDM) capability. The new model incorporates the framework of Critical State Soil Mechanics, while it relies upon bounding surface plasticity with a vanished elastic region to simulate the non‐linear soil response. Stress integration of constitutive relations is performed using a recently proposed explicit scheme with automatic error control and substepping, which so far has been employed in the literature only for constitutive models aiming at monotonic loading. The overall accuracy of this scheme is evaluated at element level by simulating cyclic loading along complex stress paths and by using iso‐error maps for paths involving change of the Lode angle. The performance of the new constitutive model and its stress integration scheme in complex boundary value problems involving earthquake‐induced liquefaction is evaluated, in terms of accuracy and computational cost, via a number of parametric analyses inspired by the successful simulation of the VELACS centrifuge Model Test No. 2 studying the lateral spreading response of a liquefied sand layer. Copyright © 2009 John Wiley & Sons, Ltd.     

On-line estimation of nonlinear physical systems

Recursive algorithms for estimating states of nonlinear physical systems are presented. Orthogonality properties are rediscovered and the associated polynomials are used to linearize state and observation models of the underlying random processes. This requires some key hypotheses regarding the structure of these processes, which may then take account of a wide range of applications. The latter include streamflow forecasting, flood estimation, environmental protection, earthquake engineering, and mine planning. The proposed estimation algorithm may be compared favorably to Taylor series-type filters, nonlinear filters which approximate the probability density by Edgeworth or Gram-Charlier series, as well as to conventional statistical linearization-type estimators. Moreover, the method has several advantages over nonrecursive estimators like disjunctive kriging. To link theory with practice, some numerical results for a simulated system are presented, in which responses from the proposed and extended Kalman algorithms are compared.