Spread <Emphasis Type="Italic">F</Emphasis> in the Midlatitude Ionosphere According to DPS-4 Ionosonde Data

The results of studying spread F obtained from the DPS-4 ionosonde data at the observatory of the Pushkov Institute of Terrestrial Magnetism, Ionosphere, and Radio Wave Propagation (Moscow) are presented. The methodical questions that arise during the study of a spread F phenomenon in the ionosphere are considered; the current results of terrestrial observations are compared with previously published data and the results of sounding onboard an Earth-satellite vehicle. The automated algorithm for estimation of the intensity of frequency spread F, which was developed by the authors and was successfully verified via comparison of the data of the digisonde DPS-4 and the results of manual processing, is described. The algorithm makes it possible to quantify the intensity of spread F in megahertz (the dFs parameter) and in the number of points (0, 1, 2, 3). The strongest spread (3 points) is shown to be most likely around midnight, while the weakest spread (0 points) is highly likely to occur during the daytime. The diurnal distribution of a 1–2 point spread F in the winter indicates the presence of additional maxima at 0300–0600 UT and 1400–1700 UT, which may appear due to the terminator. Despite the large volume of processed data, we can not definitively state that the appearance of spread F depends on the magnetic activity indices Kp, Dst, and AL, although the values of the dFs frequency spread interval strongly increased both at day and night during the magnetic storm of March 17–22, 2015, especially in the phase of storm recovery on March 20–22.     

Study on uncertainty of inversion and resolution limit in the presence of noise

This study analyzed the uncertainty of inversion and the resolution limit in the presence of noise by means of statistical experiments. The exhaustive method is adopted to obtain the global optimal solution in each experiment. We found that even with small level of noise, solutions fluctuate in a large range for the thin bed. The distribution of solutions in the presence of noise is closely related to the spread of the cost function in the absence of noise. As a result, the area of a certain neighborhood around the true solution on the spread of the cost function in the absence of noise is used to evaluate the uncertainty of inversion and the resolution limit in the presence of noise. In the case that the SNR (signal-to-noise ratio) is 5 in this study, solutions focus around the true solution with a very small uncertainty only when the bed thickness is greater than the reciprocal of the double predominant frequency of the convoluting wavelet.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Evolution and observational signature of diffused antiworld

The existence of macroscopic regions with antibaryon excess in the baryon asymmetric Universe with general baryon excess is the possible consequence of practically all models of baryosynthesis. Diffusion of matter and antimatter to the border of antimatter domains defines the minimal scale of the antimatter domains surviving to the present time. A model of diffused antiworld is considered, in which the density within the surviving antimatter domains is too low to form gravitationally bound objects. The possibility to test this model by measurements of cosmic gamma ray fluxes is discussed. The expected gamma ray flux is found to be acceptable for modern cosmic gamma ray detectors and for those planned for the near future.     

Modeling of type III bursts dynamic spectra

The one-dimensional process of spatially limited electron stream propagation in the solar corona is simulated. It is shown that the beam instability development results either in strong relaxation in velocity space and inhibition of spatial diffusion (high-stream density) or in velocity space relaxation decrease and simultaneous growth of spatial stream length (low-stream density). Assuming a profile of background plasma density to be exponential, dynamic spectra of type III bursts are modeled, which shows that the emission source velocity is constant, and a duration of the burst emission at a given frequency reduces for high-stream densities.     

Tectonics of southwestern Mexico,isotopic evidence,nuclear Central America,Late Cretaceous break up

In the paleogeographic reconstruction of Mexico and northern Central America, an ever-increasing amount of evidence shows that the entire region is a collage of suspect terranes transported from abroad, whose timing and sense of motion are now beginning to be understood. Among these, the Chortis block (nuclear Central America) and the Baja California Peninsula have been proposed as pieces of continent separated from the Pacific coast of southwestern Mexico, that have moved either southeastward by the Farallon plate or northwestward by the Kula plate. Previous studies mainly confined to the northern margin of the Chortis block, confirmed a left-lateral displacement of 130 km in Neogene time. Further studies made northwestward along the Mexican coast provided a better understanding of magmatic and metamorphic processes in the area, and suggested times of detachment increased to 30 Ma, 40 Ma, and 66 Ma. The pre-detachment westernmost position of the block has changed, depending on the model chosen, from Puerto Vallarta and beyond, to the current position. Here we show that the isotopic mineral ages from coastal granites along the coast from Puerto Vallarta, Jalisco (80 Ma) to Puerto Angel, Oaxaca (11 Ma) record systematic decrease of cooling ages from NW to SE. This pattern is interpreted to result from the progressive uplift of rocks exposed at the present-day coast in that direction, such uplift occurred in response to the development of the Middle America Trench at the newly formed continental margin when the Chortis block was sliding at an average rate of 1.5 cm/year in a sinistral sense to its present position. Our results also constrain the position of the Kula-Farallon spreading axis north of Puerto Vallarta. These observations led us to conclude that several indicators point to this time and region for the onset of strike-slip drifting of the Chortis block toward its current position. Here, we also present several view points in terms of other possilble interpretations to different tectonic, geologic and isotopic data sets published recently by different authors.     

Numerical Modeling of Ice Fog in Interior Alaska Using the Weather Research and Forecasting Model

An ice microphysics parameterization scheme has been modified to better describe and understand ice fog formation. The modeling effort is based on observations in the Sub-Arctic Region of Interior Alaska, where ice fog occurs frequently during the cold season due to abundant water vapor sources and strong inversions existing near the surface at extremely low air temperatures. The microphysical characteristics of ice fog are different from those of other ice clouds, implying that the microphysical processes of ice should be changed in order to generate ice fog particles. Ice fog microphysical characteristics were derived with the NCAR Video Ice Particle Sampler during strong ice fog cases in the vicinity of Fairbanks, Alaska, in January and February 2012. To improve the prediction of ice fog in the Weather Research and Forecasting model, observational data were used to change particle size distribution properties and gravitational settling rates, as well as to implement a homogeneous freezing process. The newly implemented homogeneous freezing process compliments the existing heterogeneous freezing scheme and generates a higher number concentration of ice crystals than the original Thompson scheme. The size distribution of ice crystals is changed into a Gamma distribution with the shape factor of 2.0, using the observed size distribution. Furthermore, gravitational settling rates are reduced for the ice crystals since the crystals in ice fog do not precipitate in a similar manner when compared to the ice crystals of cirrus clouds. The slow terminal velocity plays a role in increasing the time scale for the ice crystals to settle to the surface. Sensitivity tests contribute to understanding the effects of water vapor emissions as an anthropogenic source on the formation of ice fog.     

Decay of a highly eccentric satellite

The re-entry phase of a highly eccentric satellite is discussed. Numerical simulations allowing the prediction of the exact date of re-entry of a highly eccentric satellite are exposed.It is shown that under very particular circumstances the life of the satellite can be extended by a few days. The number of final revolutions of the rapidly contracting orbit depends critically on the air density between 70 km and 100 km.Re-entry of the European scientific satellite HEOS-1 predicted for 28 October, 1975 is near such a situation.     

Integration of reflection seismic and sediment grain-size data from Lake Khubsugul (Northern Mongolia): a reply to Prokopenko and Kendall

Prokopenko and Kendall (J Paleolimnol doi:, 2008) criticise the work presented in Fedotov et al. (J Paleolimnol 39:335–348, 2008), and instead propose an alternative interpretation for the grain-size evolution recorded in the KDP-01 core, retrieved from the central part of Lake Khubsugul. Their interpretation is based (i) on a seismic-stratigraphic re-interpretation of sparker seismic profile khub012 (which they copied from Fedotov et al. (EOS Trans 87:246–250, 2006)), (ii) on the presupposition that changes in lake level are the dominant control on facies distribution in Lake Khubsugul, and (iii) on the invalidation of our age-depth model. In this reply to their comment, we demonstrate that they interpreted seismic artefacts and geometries caused by changes in profile orientation as true stratigraphic features and that the lake-level reconstruction they derive from this interpretation is therefore incorrect. We also demonstrate that their grain-size predictions, which they consider to be predominantly driven by changes in lake level, are inconsistent with the measured sulphate concentration, which is a demonstrated proxy of lake level in Lake Khubsugul, and with the measured grain-size record. Finally, we point out that even if there would be a problem with the age-depth model, this problem would not affect the part of the sedimentary sequence discussed in Fedotov et al. (J Paleolimnol 39:335–348, 2008).     

On the instability of weakly radially anisotropic star clusters

We discuss contradictions existing in the literature in the problem on the stability of collisionless spherical stellar systems, which are the simplest anisotropic generalization of the well-known polytropic models. On the one hand, calculations of the growth rates within the framework of a linear stability theory and N-body simulations suggest that these systems should become stable when the parameter s characterizing the degree of anisotropy of the stellar velocity distribution becomes lower than some critical value s _crit > 0. On the other hand, according to Palmer and Papaloizou, the growth rate should be nonzero up to the isotropic limit s = 0. Using our method of determining the eigenmodes of stellar systems, we show that even though the mode growth rates in weakly radially anisotropic systems of this type are nonzero, they are exponentially small, i.e., decrease as γ ∝ exp(−a/s) when s → 0. For slightly radially anisotropic systems with a finite lifetime, this actually implies stability.