Non-spherical source-surface model of the heliosphere: a scalar formulation

The source-surface method offers an alternative to full MHD simulation of the heliosphere. It entails specification of a surface from which the solar wind flows normally outward along straight lines. Compatibility with MHD results requires this (source) surface to be non-spherical in general and prolate (aligned with the solar dipole axis) in prototypical axisymmetric cases. Mid-latitude features on the source surface thus map to significantly lower latitudes in the heliosphere. The model is usually implemented by deriving the B field (in the region surrounded by the source surface) from a scalar potential formally expanded in spherical harmonics, with coefficients chosen so as to minimize the mean-square tangential component of B over this surface. In the simplified (scalar) version the quantity minimized is instead the variance of the scalar potential over the source surface. The scalar formulation greatly reduces the time required to compute required matrix elements, while imposing essentially the same physical boundary condition as the vector formulation (viz., that the coronal magnetic field be, as nearly as possible, normal to the source surface for continuity with the heliosphere). The source surface proposed for actual application is a surface of constant , where r is the heliocentric distance and B is the scalar magnitude of the B field produced by currents inside the Sun. Comparison with MHD simulations suggests that k 1.4 is a good choice for the adjustable exponent. This value has been shown to map the neutral line on the source surface during Carrington Rotation 1869 (May–June 1993) to a range of latitudes that would have just grazed the position of Ulysses during that month in which sector structure disappeared from Ulysses magnetometer observations.     

The orion-IRc2 outflow observed in COJ=4-3

A CO J=4-3 map of the Orion-IRc2 region (size 80, beam 15) is presented. The outflow has a bipolar structure observed nearly pole-on. We interpret our observations in terms of an hierarchy of tubes where the fastest gas proceeds through the narrowest, innermost cone. The gas has an inhomogeneous temperature distribution and must be highly clumped. Most of the properties derived comply with the picture developed for the few known extreme high velocity outflows; however, this outflow contains exceptionally dense gas ( > 10⁶ cm^–3).     

Airborne total eclipse observation of the extreme solar limb at 400 μm

The total solar eclipse of February 26, 1979 was monitored at far infrared wavelengths from the NASA Lear Jet Observatory flying at 12.9 km in the eclipse shadow. The resultant eclipse curve for radiation within a bandwidth of 20 cm^–1 centered upon 25 cm^–1 (400 m) was measured and analysed at an equivalent angular resolution of 1 arc sec over a 100 arc sec region adjacent to the limb to provide information on the intensity distribution of continuum radiation close to this limb. The curve has been compared to predictions derived from models of the solar atmosphere for the specific geometry of this eclipse, and is shown to match most closely that derived from a uniform distribution of radiation across the disk. This is in distinct contrast to the result obtained in the only other comparable experiment, carried out over Africa in 1973 from a supersonic Concorde aircraft, in which an intense but narrow spike of far infrared radiation at the extreme solar limb was inferred from the data. The absence also in the present observations of the significant limb brightening predicted by the HSRA model (in which homogeneity within the source region is assumed) is in substantial agreement with lower resolution results from mountain altitudes. This result is interpreted as further evidence for the presence in the Sun's lower chromosphere of significant inhomogeneity with a scale size of at least 1000 km at this depth.     

Narrow-line dynamics in Seyfert and starburst nuclei

Rotating disks and expanding supershells of ionized gas are supposed to be kinematical components essential for the formation of emission-like profiles from starburst nuclei and the narrow-line regions of Seyfert galaxies. Simple models predict a few typical features of line profiles by which disks may be recognized.A giant supershell is identified SE of the nucleus of the starburst galaxy NGC 253. We suggest that the variously reported outflow cones in Seyfert nuclei are also induced by superbubble dynamics. Observations of NGC 4151 are shown to be consistent with this scenario.Paper presented at the 11th European Regional Astronomical Meetings of the IAU on New Windows to the Universe, held 3–8 July, 1989, Tenerife, Canary Islands, Spain.     

Atmospheric nitrogen inputs into the North Sea: effect on productivity

The ANICE (Atmospheric Nitrogen Inputs into the Coastal Ecosystem) project addressed the atmospheric deposition of nitrogen to the North Sea, with emphasis on coastal effects. ANICE focused on quantifying the deposition of inorganic nitrogen compounds to the North Sea and the governing processes. An overview of the results from modelling and experimental efforts is presented. They serve to identify the role of the atmosphere as a source of biologically essential chemical species to the marine biota. Data from the Weybourne Atmospheric Observatory (UK) are used to evaluate the effect of short episodes with very high atmospheric nitrogen concentrations. One such episode resulted in an average deposition of 0.8 mmol N m⁻² day⁻¹, which has the potential to promote primary productivity of 5.3 mmol C m⁻² day⁻¹. This value is compared to long-term effects determined from model results. The total calculated atmospheric deposition to the North Sea in 1999 is 948 kg N km⁻¹, i.e. 0.19 mmol N m⁻² day⁻¹ which has the potential to promote primary productivity of 1.2 mmol C m⁻² day⁻¹. Detailed results for August 1999 show strong gradients across the North Sea due to adjacent areas where emissions of NO_x and NH₃ are among the highest in Europe. The average atmospheric deposition to the southern part of the North Sea in August 1999 could potentially promote primary production of 2.0 mmol C m⁻² day⁻¹, i.e. 5.5% of the total production at this time of the year in this area of the North Sea. For the entire study area the atmospheric contribution to the primary production per m² is about two-third of this value. Most of the deposition occurs during short periods with high atmospheric concentrations. This atmospheric nitrogen is almost entirely anthropogenic in origin and thus represents a human-induced perturbation of the ecosystem.     

Efficiency of base isolation systems in structural seismic protection and energetic assessment

This paper concerns the seismic response of structures isolated at the base by means of High Damping Rubber Bearings (HDRB). The analysis is performed by using a stochastic approach, and a Gaussian zero mean filtered non‐stationary stochastic process is used in order to model the seismic acceleration acting at the base of the structure. More precisely, the generalized Kanai–Tajimi model is adopted to describe the non‐stationary amplitude and frequency characteristics of the seismic motion. The hysteretic differential Bouc–Wen model (BWM) is adopted in order to take into account the non‐linear constitutive behaviour both of the base isolation device and of the structure. Moreover, the stochastic linearization method in the time domain is adopted to estimate the statistical moments of the non‐linear system response in the state space. The non‐linear differential equation of the response covariance matrix is then solved by using an iterative procedure which updates the coefficients of the equivalent linear system at each step and searches for the solution of the response covariance matrix equation. After the system response variance is estimated, a sensitivity analysis is carried out. The final aim of the research is to assess the real capacity of base isolation devices in order to protect the structures from seismic actions, by avoiding a non‐linear response, with associated large plastic displacements and, therefore, by limiting related damage phenomena in structural and non‐structural elements. In order to attain this objective the stochastic response of a non‐linear n‐dof shear‐type base‐isolated building is analysed; the constitutive law both of the structure and of the base devices is described, as previously reported, by adopting the BWM and by using appropriate parameters for this model, able to suitably characterize an ordinary building and the base isolators considered in the study. The protection level offered to the structure by the base isolators is then assessed by evaluating the reduction both of the displacement response and the hysteretic dissipated energy. Copyright © 2003 John Wiley & Sons, Ltd.     

Use of phytoplankton species indicators to track the origin of phytoplankton blooms in Willapa Bay, Washington

We focus on the question of whether high phytoplankton production events observed in a United States Pacific Northwest estuary consist of estuarine species blooms fueled by oceanic nutrient input or reflect offshore oceanic blooms that have advected into the estuary. Our approach is to use certain phytoplankton species as indicators associated with water mass origin, either estuarine or oceanic, to help resolve this question in Willapa Bay, Washington. We used species analysis and primary production data from 10 selected dates in May–September of 1998 and 1999, representing periods of high through low productivity. Out of 108 phytoplankton species identified from Willapa Bay, nine were selected and tested as indicators of oceanic species, six as estuarine, and two as surf zone. Our test results demonstrated the oceanic and estuarine species to be satifactory indicators of source waters. The prevalence of these species indicators in our samples revealed that the highest primary production and the appearance ofPseudo-nitzschia spp. were associated with oceanic intrusions of phytoplankton biomass into Willapa Bay. While the largest blooms were oceanic in origin, numerous medium-sized production events were from either oceanic, surf zone, or estuarine sources, indicating a complex situation.     

An estimate of hypocentre location accuracy in a large network: possible implications for tectonic studies in Italy

Data recorded by the Italian Telemetered Seismic Network (ITSN) of the Istituto Nazionale di Geofisica (ING) have been widely used in recent years to image slab structures and to find evidence for active processes along the Italian Peninsula. However, the use of seismic data for geostructural purposes may be affected by the well-known trade-off between earthquake location and seismic-velocity parameters. Furthermore, the confidence ellipse predicted by standard procedures may be inadequate for the representation of the probable error of a computed localization. This paper evaluates the probable errors on the hypocentre determinations of the seismic events recorded by the ITSN, using a Monte Carlo method.
We compute synthetic arrival times using a 1-D velocity model appropriate as an average for the Italian area. The hypocentres used are all those recorded by the ITSN during the period January 1992 to March 1994 (1972 events). Station locations are those of the current ITSN configuration. The synthetic arrival times are perturbed with a Gaussian distribution of errors and input to ING's standard hypocentral location procedure, but using crustal velocities differing by 10 per cent from those used to generate them. Each simulation is repeated at least 30 times. Average absolute shifts of hypocentres are assessed in grid cells of linear dimension 33 km covering the whole Italian region.
For regions within the ITSN, shifts are typically 5–10 km in location and up to 20 km in depth. However, for offshore and coastal regions, they are much greater: 50 km or more in both location and depth (far exceeding the equivalent uncertainties quoted by ING bulletins). Possible consequences of this are highlighted by producing a cross-section of subcrustal hypocentres from the Adriatic to the Tyrrhenian Sea, where the large uncertainty in depth precludes any confident interpretation of dipping tectonic features.     

Pore-water response on seasonal environmental changes in intertidal sediments of the Weser Estuary,Germany

In order to determine time-dependent changes in estuarine pore-water chemistry and flux variations across the sediment-water interface, sediment cores of an intertidal mud flat in the Weser Estuary were taken monthly over a one-year period. Sediment temperature, pH, Eh, Cl^–, O₂, NO ₃ ^– , and SO ₄ ^2– pore-water concentrations were measured and showed variations that relate to the changes of surface temperature and estuarine water composition. Fick's first law was applied to quantify diffusive fluxes from concentration gradients in the diffusive boundary layer and in the pore water. Total nitrate fluxes were calculated from flux chamber experiments. Diffusive oxygen fluxes increased from 5 mmol m^–2 d^–1 in winter to 18 mmol m^–2 d^–1 in early summer, while nitrate fluxes into the sediment increased from 3 mmol m^–2 d^–1 in winter to 60 mmol m^–2 d^–1 in early summer. Oxygen and nitrate fluxes into the sediment correlated linearly to sediment temperature. Sulfate fluxes increased from 0.5 mmol m^–2 d^–1 in winter to 10 mmol m^–2 d^–1 in August and September. Converted into carbon fluxes, the sum of these oxidants ranged from 10 mmol m^–2 d^–1 in winter to 80 mmol m^–2 d^–1 in summer. An estimation of the upper limit of the annual nitrate flux into the sediment showed that about 10% of the 250,000 t of nitrate discharged annually by the river may be decomposed within the inner Weser Estuary.     

Doubtful Events Reported in Current Italian Earthquake Catalogues: The Medieval Earthquakes of Subiaco (Central Italy)

Within the framework of a study of the seismicity of the Aniene Valley (Central Italy), we analysed the medieval earthquakes of Subiaco (1216, 1227, 1299), the largest events reported for the area. Our main goal was to investigate some doubtful events reported in earthquake catalogues and, as such, currently utilised for seismic hazard estimates. A careful screening of the oldest available sources and their filiation pattern up to the present pointed out the uncertainty on the date and nature of these phenomena. A multidisciplinary approach based on the joint analysis of archaeological, geomorphologic and historical evidence allowed us to propose new interpretations concerning these events and their significance for the assessment of seismic hazard in the Aniene Valley. The main conclusion is that the dates of the 1216 and 1227 events are fairly unsupported. In particular, the 1216 earthquake could be dated back to between AD 1159 and 1181.