First results of operational ionospheric dynamics prediction for the Brazilian Space Weather program

It is shown the development and preliminary results of operational ionosphere dynamics prediction system for the Brazilian Space Weather program. The system is based on the Sheffield University Plasmasphere–Ionosphere Model (SUPIM), a physics-based model computer code describing the distribution of ionization within the Earth mid to equatorial latitude ionosphere and plasmasphere, during geomagnetically quiet periods. The model outputs are given in a 2-dimensional plane aligned with Earth magnetic field lines, with fixed magnetic longitude coordinate. The code was adapted to provide the output in geographical coordinates. It was made referring to the Earth’s magnetic field as an eccentric dipole, using the approximation based on International Geomagnetic Reference Field (IGRF-11). During the system operation, several simulation runs are performed at different longitudes. The original code would not be able to run all simulations serially in reasonable time. So, a parallel version for the code was developed for enhancing the performance. After preliminary tests, it was frequently observed code instability, when negative ion temperatures or concentrations prevented the code from continuing its processing. After a detailed analysis, it was verified that most of these problems occurred due to concentration estimation of simulation points located at high altitudes, typically over 4000 km of altitude. In order to force convergence, an artificial exponential decay for ion–neutral collisional frequency was used above mentioned altitudes. This approach shown no significant difference from original code output, but improved substantially the code stability. In order to make operational system even more stable, the initial altitude and initial ion concentration values used on exponential decay equation are changed when convergence is not achieved, within pre-defined values. When all code runs end, the longitude of every point is then compared with its original reference station longitude, and differences are compensated by changing the simulation point time slot, in a temporal adjustment optimization. Then, an approximate neighbor searching technique was developed to obtain the ion concentration values in a regularly spaced grid, using inverse distance weighting (IDW) interpolation. A 3D grid containing ion and electron concentrations is generated for every hour of simulated day. Its spatial resolution is 1° of latitude per 1° of longitude per 10 km of altitude. The vertical total electron content (VTEC) is calculated from the grid, and plotted in a geographic map. An important feature that was implemented in the system is the capacity of combining observational data and simulation outputs to obtain more appropriate initial conditions to the ionosphere prediction. Newtonian relaxation method was used for this data assimilation process, where ionosonde data from four different locations in South America was used to improve the system accuracy. The whole process runs every day and predicts the VTEC values for South America region with almost 24 h ahead.     

Fluid inclusion studies of chemosynthetic carbonates: strategy for seeking life on Mars

Fluid inclusions in minerals hold the potential to provide important data on the chemistry of the ambient fluids during mineral precipitation. Especially interesting to astrobiologists are inclusions in low-temperature minerals that may have been precipitated in the presence of microorganisms. We demonstrate that it is possible to obtain data from inclusions in chemosynthetic carbonates that precipitated by the oxidation of organic carbon around methane-bearing seepages. Chemosynthetic carbonates have been identified as a target rock for astrobiological exploration. Other surficial rock types identified as targets for astrobiological exploration include hydrothermal deposits, speleothems, stromatolites, tufas, and evaporites, each of which can contain fluid inclusions. Fracture systems below impact craters would also contain precipitates of minerals with fluid inclusions. As fluid inclusions are sealed microchambers, they preserve fluids in regions where water is now absent, such as regions of the martian surface. Although most inclusions are < 5 microns, the possibility to obtain data from the fluids, including biosignatures and physical remains of life, underscores the advantages of technological advances in the study of fluid inclusions. The crushing of bulk samples could release inclusion waters for analysis, which could be undertaken in situ on Mars.     

Near 13.5-day periodicity in Muon Detector data during late 2001 and early 2002

In this study we perform a continuous Morlet wavelet transform method in time series of secondary cosmic rays and 1 AU interplanetary medium parameters for the interval from October 2001 to October 2002. The near 13.5-day periodicity was obtained during late 2001, and it was remarkable for muon data. Even though some works have pointed out that the main activations of the 13.5 day recurrence in near-Earth solar wind are related, e.g., with the heliosheet crossings or to the occurrence at 1 AU of two high speed streams approximately 180° apart in solar longitude per solar rotation, we aim to show that the period of about half the solar rotation during the end months of 2001 present in muon time series was apparently due to the occurrence of non-recurrent interplanetary disturbances. The interconnections among successive Forbush decreases, recovery phases and gradual muon depressions (associated with corotating interaction regions) seem to play an important role in such 13.5-day periodicity.     

A method to test HF ray tracing algorithm in the ionosphere by means of the virtual time delay

It is well known that a 3D ray tracing algorithm furnishes the ray’s coordinates, the three components of the wave vector and the calculated group time delay of the wave along the path. The latter quantity can be compared with the measured group time delay to check the performance of the algorithm. Simulating a perfect reflector at an altitude equal to the virtual height of reflection, the virtual time delay is assumed as a real time delay. For a monotonic electronic density profile we find a very small relative difference between the calculated and the virtual delay for both analytic and numerical 3D electronic density models.     

MHD numerical simulations of proto-stellar jets

We will summarize in this paper the effects that the presence of the magnetic field can cause to proto-stellar jet dynamics, structure and emission line properties, and the differences between two- and three-dimensional numerical simulations will be emphasized. This revised version was published online in August 2006 with corrections to the Cover Date.