Numerical simulations of very large impacts on the Earth

Vertical impacts on the Earth of asteroids 500-3000 km in diameter at 15 km/s have been numerically modelled using the hydrodynamic SOVA code. This code has been modified for the spherical system of coordinates well suited for simulations of very large impacts when the entire Earth is involved in motion. The simulations include cratering process, upward motion of deep mantle layers, fall of ejecta on the Earth, escape of matter to space, and formation of rock vapour atmospheres. The calculations were made for the period preceding disappearance of rock vapour atmospheres caused by radiation several years after the largest impacts. For very large vertical impacts at 15 km/s, escaping masses proved to be negligibly small. Quantities of kinetic, internal, potential, and radiated away energies are obtained as functions of time and space. After the impacts, a global layer of condensed ejecta covers the whole of the Earth's surface and the ejecta energy is sufficient to vaporise an ocean 3 km deep. The mass of rock vapour atmosphere is 10-23% of the impactor mass. This atmosphere has a greater mass than the water atmosphere if impactor is 2000 km in diameter or larger.     

Low Visibility Formation and Forecasting on the Northern Coast of Brazil

Visibility analysis and forecast at the Maceio International Airport in the Brazilian Northeast (NEB) was the principal goal of this investigation. Surface meteorological data of the Maceio International Airport were used for low visibility frequency study. Low visibility in NEB was provoked more frequently by light fog (LF) formation (1,098 or 92 h month^?1 on average). Haze and fog were very rare (81 h and one event per year, respectively on average). Light fog with a visibility less than 2 km usually was detected together with rain or drizzle. Low visibility was observed more frequently at night and during the rainy season. Applications of the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model for light fog forecast were tested. Thermodynamic processes were studied by vertical profile, elaborated by: (1) National Centers for Environmental Prediction (NCEP) reanalysis data for Maceio (because of some radiosonde absence) and (2) forecast vertical temperature and humidity profiles were produced, using Air Parcels Trajectories of the HYSPLIT model at the pattern levels. The synoptic situations before and during low visibility phenomena were analyzed using different products of NCEP reanalysis, the high resolution (10 km) ETA model and infrared satellite images. Wave disturbance in the trade winds field, localized on the northwest periphery of the South Atlantic subtropical High, usually accompanied the phenomena. A humidity advection, weak ascendant movement and thermal inversion absence at the low levels were created by these waves. The middle level’s descendent movement provoked the humidity accumulation at levels below. Satisfactory results of the HYSPLIT model applications for light fog forecast were obtained with 12 h antecedence. In particular, stable level forecast by the ETA model was forecast satisfactorily with 12 h antecedence; vertical movements were predicted better with up to 48 h antecedence. The PSU/NCAR mesoscale model (MM5) and PAFOG models were tested for analysis and forecast of an intensive fog event. Intensive fog provoked a fatal accident of a small airplane near the Maceio Airport in 2007. These fog formation processes were studied by NCEP reanalysis data, the high resolution regional model MM5, and satellite and radar data. Fog formation was simulated by PAFOG model and satisfactory results were obtained with 10 h antecedence.     

Robotic optical telescopes global network MASTER II. Equipment, structure, algorithms

Presented paper describes the basic principles and features of the implementation of a robotic network of optical telescopes MASTER, designed to study the prompt (simultaneous with gamma radiation) optical emission of gamma-ray bursts and to perform the sky survey to detect unknown objects and transient phenomena. With joint efforts of Sternberg astronomical institute, High altitude astronomical station of the Pulkovo observatory, Ural state university, Irkutsk state university, Blagoveshchensk pedagogical university, the robotic telescopes MASTER?II near Kislovodsk, Yekaterinburg, Irkutsk and Blagoveshchensk were installed and tested. The network spread over the longitudes is greater than 6?h. A further expansion of the network is considered.     

Application of horizontal-to-vertical (H/V) Fourier spectral ratio for analysis of site effect on rock (NEHRP-class B) sites in Taiwan

The frequency-dependent amplification for rock (NEHRP-class B) sites was studied using earthquake ground-motion database collected in Taiwan during implementation of the Taiwan Strong Motion Instrumentation Program. The database used includes several hundred records from earthquakes of M_L 4.0–7.3 occurred between 1993 and 2004. The characteristics of amplification were evaluated using the well-known technique of horizontal-to-vertical Fourier spectral ratio (H/V) of the S-wave phase [Lermo J, Chavez-Garcia FJ. Site effect evaluation using spectral ratios with only one station. Bull Seism Soc Am 1993;83:1574–94]. The study allows us to analyze peculiarities of rock sites amplification in Northern and Eastern Taiwan. It was suggested to divide the NEHRP-class B site amplification into four types based on frequency of maximum amplification and the shape of amplification function. The applicability of the technique was also checked for a few stiff and soft soil sites (NEHRP-classes D and E).     

Determination of ecological significance based on geostatistical assessment: a case study from the Slovak Natura 2000 protected area

The Sitno Natura 2000 Site covers an area of 935,56 hectares. The Sitno region is significant due to the number of rare and endangered species of plants, and as a result is considered a location of great importance to the maintenance of floral gene pools. The study area suffers human impacts in the form of tourism. The main purpose of this study is to the measure landscape elements, determine the ecological significance of habitats within the Sitno area, and from this data, organize the study area into conservation zones. The results of this landscape quantification are numerical values that can be used to interpret the quality of ongoing ecological processes within individual landscape types. Interpretation of this quantified data can be used to determine the ecological significance of landscapes in other study areas. This research examines the habitats of Natura 2000 Sites by a set of landscape metrics for habitat area, size, density, and shape, such as Number of patches (NP), Patch density (PD), Mean patch size (MPS), Patch size standard deviation (PSSD) and Mean shape index (MSI). The classification of land cover patches is based on the Annex Code system.     

On the modeling of ground-motion field for assessment of multiple-location hazard,damage, and loss: example of estimation of electric network performance during scenario earthquake

Consideration of within-earthquake ground-motion correlation is essential for the estimation of seismic hazards, damage, and loss for spatially distributed systems. In many seismically active regions, the strong motion data of real engineering significance are completely unavailable or very scarce. The absence of necessary data does not allow developing regional empirical correlation models, and the models obtained for other regions should be used with correspondent sensitivity analysis. The level of within-earthquake correlation may vary in broad range; therefore, development of correspondent criteria for selection from available models is important. In this paper, we analyzed the performance of a system of critical elements of electric power network (substations) depending on variations in within-earthquake correlation. The performance is described as probability of different levels of non-functionality, i.e., percentage of area suffering power outage, and probability of expected number of customers without power. We have shown that the proper choice of the within-earthquake correlation model is critical in comprehensive estimations of functionality of substations in electrical power system. Skipping the ground-motion variability and within-earthquake correlation may lead to unreliable results. Relevance of geology-based within-earthquake correlation models has been demonstrated, and a scheme, which allows reducing uncertainty in the choice of realistic correlation, has been proposed.     

Non-Fickian mass transport in fractured porous media

The paper provides an introduction to fundamental concepts of mathematical modeling of mass transport in fractured porous heterogeneous rocks. Keeping aside many important factors that can affect mass transport in subsurface, our main concern is the multi-scale character of the rock formation, which is constituted by porous domains dissected by the network of fractures. Taking into account the well-documented fact that porous rocks can be considered as a fractal medium and assuming that sizes of pores vary significantly (i.e. have different characteristic scales), the fractional-order differential equations that model the anomalous diffusive mass transport in such type of domains are derived and justified analytically. Analytical solutions of some particular problems of anomalous diffusion in the fractal media of various geometries are obtained. Extending this approach to more complex situation when diffusion is accompanied by advection, solute transport in a fractured porous medium is modeled by the advection-dispersion equation with fractional time derivative. In the case of confined fractured porous aquifer, accounting for anomalous non-Fickian diffusion in the surrounding rock mass, the adopted approach leads to introduction of an additional fractional time derivative in the equation for solute transport. The closed-form solutions for concentrations in the aquifer and surrounding rocks are obtained for the arbitrary time-dependent source of contamination located in the inlet of the aquifer. Based on these solutions, different regimes of contamination of the aquifers with different physical properties can be readily modeled and analyzed.     

Ocean-atmosphere CO2 exchange: An accessible lab simulation for considering biological effects

Phytoplankton is considered a key component mediating the ocean-atmospheric exchange of carbon dioxide and oxygen. Lab simulations which model biological responses to atmospheric change are difficult to translate into natural settings owing in part to the vertical migration of phytoplankton. In the sea this vertical migration acts to regulate actual carbon dioxide consumption. To capture some critical properties of this vertical material transfer, we monitored the effects of atmospheric CO₂ on dense suspensions of bioconvecting microorganisms. Bioconvection refers to the spontaneous patterns of circulation which arise among such upwardly swimming cells as alga, protozoa, zoospore and large bacteria. Gravity, phototaxis and chemotaxis have all been implicated as affecting pattern-forming ability. The ability of a biologically active suspension to detect atmospheric changes offers a unique method to quantify organism adjustment and vertical migration. With increasing CO₂, bioconvection patterns in alga (P. parva) and protozoa (T. pyriformis) lose their robustness, and surface cell populations retreat from the highest CO₂ regions. Cell movement (both percent motile and mean velocity) generally diminishes. A general program of image analysis yields statistically significant variations in macroscopic migration patterns; both fractal dimension and various crystallographic parameters correlate strongly with carbon dioxide content.     

An Efficient Non-iterative Bulk Parametrization of Surface Fluxes for Stable Atmospheric Conditions Over Polar Sea-Ice

In climate and weather prediction models the near-surface turbulent fluxes of heat and momentum and related transfer coefficients are usually parametrized on the basis of Monin–Obukhov similarity theory (MOST). To avoid iteration, required for the numerical solution of the MOST equations, many models apply parametrizations of the transfer coefficients based on an approach relating these coefficients to the bulk Richardson number \(Ri_{b}\). However, the parametrizations that are presently used in most climate models are valid only for weaker stability and larger surface roughnesses than those documented during the Surface Heat Budget of the Arctic Ocean campaign (SHEBA). The latter delivered a well-accepted set of turbulence data in the stable surface layer over polar sea-ice. Using stability functions based on the SHEBA data, we solve the MOST equations applying a new semi-analytic approach that results in transfer coefficients as a function of \(Ri_{b}\) and roughness lengths for momentum and heat. It is shown that the new coefficients reproduce the coefficients obtained by the numerical iterative method with a good accuracy in the most relevant range of stability and roughness lengths. For small \(Ri_{b}\), the new bulk transfer coefficients are similar to the traditional coefficients, but for large \(Ri_{b}\) they are much smaller than currently used coefficients. Finally, a possible adjustment of the latter and the implementation of the new proposed parametrizations in models are discussed.     

Bubble growth in visco-elastic magma: implications to magma fragmentation and bubble nucleation

We present a visco-elastic bubble growth model, accounting for viscous and elastic deformations and for volatile mass transfer between bubbles and melt. We define the borders between previous bubble growth models accounting for incompressible viscous melt, and our new model accounting also for elastic deformation; this is done by a set of end-member analytical solutions and numerical simulations. Elastic deformation is most prominent for magma of small vesicularity, where the growth regime depends on the shear modulus. For high shear modulus, bubble growth is slow and follows an exponential law in a viscous growth regime, while for low shear modulus bubbles quickly follow a square-root diffusive solution. Our model provides all the elastic components (stresses, strains and strain rates) required for defining criteria for failure and magma fragmentation. We suggest two failure criteria, a stress related one based on the internal friction and the Mohr-Coulomb failure theory, and a strain related one based on fibre elongation experiments. We argue that both criteria are equivalent if we consider their shear modulus dependency and its effect on magma rheology. Last, we apply our model to the process of bubble nucleation. In the incompressible case, following nucleation, growth is slow and leads to long incubation times during which bubbles may be dissolved back into the melt. The elastic response in magmas with low shear modulus results in a short incubation time, increasing the probability of survival. The above effects emphasize the significance of visco-elasticity for the dynamic processes occurring in magmas during volcanic activity.