On the construction of discrete approximations to linear differential expressions

Algorithm for generating discrete approximations in terms of ordinates for linear differential expression

Patents and Competition in the Automobile Industry

Methane (CH4) fluxes from world rivers are still poorly constrained, with measurements restricted mainly to temperate climates. Additional river flux measurements, including spatio-temporal studies, are important to refine extrapolations. Here we assess the spatio-temporal variability of CH4 fluxes from the Amazon and its main tributaries, the Negro, Solimoes, Madeira, Tapajos, Xingu, and Para Rivers, based on direct measurements using floating chambers. Sixteen of 34 sites were measured during low and high water seasons. Significant differences were observed within sites in the same river and among different rivers, types of rivers, and seasons. Ebullition contributed to more than 50% of total emissions for some rivers. Considering only river channels, our data indicate that large rivers in the Amazon Basin release between 0.40 and 0.58 Tg CH4 yr(-1). Thus, our estimates of CH4 flux from all tropical rivers and rivers globally were, respectively, 19-51% to 31-84% higher than previous estimates, with large rivers of the Amazon accounting for 22-28% of global river CH4 emissions.Funding Agencies|FAPESP [08/58089-9, 2011/06609-1, 2011/14502-2, 2012/17359-9]</p

Global shallow water magnetohydrodynamic waves in the solar tachocline

We derive analytical solutions and dispersion relations of global magnetic Poincar\'e (magneto-gravity) and magnetic Rossby waves in the approximation of shallow water magnetohydrodynamics. The solutions are obtained in a rotating spherical coordinate system for strongly and weakly stable stratification separately in the presence of toroidal magnetic field. In both cases magnetic Rossby waves split into fast and slow magnetic Rossby modes. In the case of strongly stable stratification (valid in the radiative part of the tachocline) all waves are slightly affected by the layer thickness and the toroidal magnetic field, while in the case of weakly stable stratification (valid in the upper overshoot layer of the tachocline) magnetic Poincar\'e and fast magnetic Rossby waves are found to be concentrated near the solar equator, leading to equatorially trapped waves. However, slow magnetic Rossby waves tend to concentrate near the poles, leading to polar trapped waves. The frequencies of all waves are smaller in the upper weakly stable stratification region than in the lower strongly stable stratification one

Long-term variation in the Sun's activity caused by magnetic Rossby waves in the tachocline

Long-term records of sunspot number and concentrations of cosmogenic radionuclides (10Be and 14C) on the Earth reveal the variation of the Sun's magnetic activity over hundreds and thousands of years. We identify several clear periods in sunspot, 10Be, and 14C data as 1000, 500, 350, 200 and 100 years. We found that the periods of the first five spherical harmonics of the slow magnetic Rossby mode in the presence of a steady toroidal magnetic field of 1200-1300 G in the lower tachocline are in perfect agreement with the time scales of observed variations. The steady toroidal magnetic field can be generated in the lower tachocline either due to the steady dynamo magnetic field for low magnetic diffusivity or due to the action of the latitudinal differential rotation on the weak poloidal primordial magnetic field, which penetrates from the radiative interior. The slow magnetic Rossby waves lead to variations of the steady toroidal magnetic field in the lower tachocline, which modulate the dynamo magnetic field and consequently the solar cycle strength. This result constitutes a key point for long-term prediction of the cycle strength. According to our model, the next deep minimum in solar activity is expected during the first half of this century.Comment: 4 pages, 4 figures, accepted in ApJ

Magnetic Rossby waves in the solar tachocline and Rieger-type periodicities

Apart from the 11-year solar cycle, another periodicity around 155-160 days was discovered during solar cycle 21 in high energy solar flares, and its presence in sunspot areas and strong magnetic flux has been also reported. This periodicity has an elusive and enigmatic character, since it usually appears only near the maxima of solar cycles, and seems to be related with a periodic emergence of strong magnetic flux at the solar surface. Therefore, it is probably connected with the tachocline, a thin layer located near the base of the solar convection zone, where strong dynamo magnetic field is stored. We study the dynamics of Rossby waves in the tachocline in the presence of a toroidal magnetic field and latitudinal differential rotation. Our analysis shows that the magnetic Rossby waves are generally unstable and that the growth rates are sensitive to the magnetic field strength and to the latitudinal differential rotation parameters. Variation of the differential rotation and the magnetic field strength throughout the solar cycle enhance the growth rate of a particular harmonic in the upper part of the tachocline around the maximum of the solar cycle. This harmonic is symmetric with respect to the equator and has a period of 155-160 days. A rapid increase of the wave amplitude could give place to a magnetic flux emergence leading to observed periodicities in solar activity indicators related with magnetic flux.Comment: 34 pages, 5 figures, accepted in Ap

Automated data reduction workflows for astronomy

Data from complex modern astronomical instruments often consist of a large number of different science and calibration files, and their reduction requires a variety of software tools. The execution chain of the tools represents a complex workflow that needs to be tuned and supervised, often by individual researchers that are not necessarily experts for any specific instrument. The efficiency of data reduction can be improved by using automatic workflows to organise data and execute the sequence of data reduction steps. To realize such efficiency gains, we designed a system that allows intuitive representation, execution and modification of the data reduction workflow, and has facilities for inspection and interaction with the data. The European Southern Observatory (ESO) has developed Reflex, an environment to automate data reduction workflows. Reflex is implemented as a package of customized components for the Kepler workflow engine. Kepler provides the graphical user interface to create an executable flowchart-like representation of the data reduction process. Key features of Reflex are a rule-based data organiser, infrastructure to re-use results, thorough book-keeping, data progeny tracking, interactive user interfaces, and a novel concept to exploit information created during data organisation for the workflow execution. Reflex includes novel concepts to increase the efficiency of astronomical data processing. While Reflex is a specific implementation of astronomical scientific workflows within the Kepler workflow engine, the overall design choices and methods can also be applied to other environments for running automated science workflows.Comment: 12 pages, 7 figure

Brachytherapy structural shielding calculations using Monte Carlo generated, monoenergetic data

To provide a method for calculating the transmission of any broad photon beam with a known energy spectrum in the range of 20 keV-1090 keV, through concrete and lead, based on the superposition of corresponding monoenergetic data obtained from Monte Carlo simulation

Collective and static properties of model two-component plasmas

Classical MD data on the charge-charge dynamic structure factor of two-component plasmas (TCP) modeled in Phys. Rev. A 23, 2041 (1981) are analyzed using the sum rules and other exact relations. The convergent power moments of the imaginary part of the model system dielectric function are expressed in terms of its partial static structure factors, which are computed by the method of hypernetted chains using the Deutsch effective potential. High-frequency asymptotic behavior of the dielectric function is specified to include the effects of inverse bremsstrahlung. The agreement with the MD data is improved, and important statistical characteristics of the model TCP, such as the probability to find both electron and ion at one point, are determined.Comment: 25 pages, 6 figures, 5 tables. Published in Physical Review E http://link.aps.org/abstract/PRE/v76/e02640