Assessing cause and effect of multiple stressors on marine systems

An operational framework is developed to serve as a guideline for investigating causal relationships between environmental stressors and effects on marine biota. Because of the complexity and variability of many marine systems, multiple lines of evidence are needed to understand relationships between stressors and effects on marine resources. Within this framework, a weight of evidence approach based on multiple lines of evidence are developed and applied in a sequential manner by (1) characterizing the study system which involves determining if target biota are impaired, assessment of food and habitat availability, and measuring contaminant levels in the environment, (2) assessing direct effects of contaminant exposure on target biota using biomarkers and assessing indirect effects of exposure using suites of bioindicators, and (3) applying standard causal criteria based on epidemiological principles and diagnostic health profiling techniques to assess potential causes. Use of multiple lines of evidence should also reduce the risk of false positives (Type I error or falsely concluding that there is a causal relationship when there is none) and false negatives (Type II error or falsely concluding there is not a causal relationship when there actually is). Understanding causal relationships and the mechanistic processes between environmental stressors and effects on biota is important in the effective management and restoration of impaired marine ecosystems.     

Volcanic risk ranking for Auckland, New Zealand. I: Methodology and hazard investigation

Volcanic eruptions typically produce a number of hazards, and many regions are at risk from more than one volcano or volcanic field. So that detailed risk assessments can be carried out, it is necessary to rank potential volcanic hazards and events in terms of risk. As it is often difficult to make accurate predictions regarding the characteristics of future eruptions, a method for ranking hazards and events has been developed that does not rely on precise values. Risk is calculated individually for each hazard from each source as the product of likelihood, extent and effect, based on the parameters order of magnitude. So that multiple events and outcomes can be considered, risk is further multiplied by the relative probability of the event occurring (probability_e) and the relative importance of the outcome (importance_o). By adding the values obtained, total risk is calculated and a ranking can be carried out.This method was used to rank volcanic hazards and events that may impact the Auckland Region, New Zealand. Auckland is at risk from the Auckland volcanic field, Okataina volcanic centre, Taupo volcano, Tuhua volcano, Tongariro volcanic centre, and Mt. Taranaki volcano. Relative probabilities were determined for each event, with the highest given to Mt. Taranaki. Hazards considered were, for local events: tephra fall, scoria fall and ballistic impacts, lava flow, base surge and associated shock waves, tsunami, volcanic gases and acid rain, earthquakes and ground deformation, mudflows and mudfills, lightning and flooding; and for distal events: tephra fall, pyroclastic flows, poisonous gases and acid rain, mudflows and mudfills, climate variations and earthquakes. Hazards from each source were assigned values for likelihood, with the largest for tephra fall from all sources, earthquakes and ground deformation, lava flows, scoria fall and base surge for an Auckland eruption on land, and earthquakes and ground deformation from an Auckland eruption in the ocean. The largest values for extent were for tephra fall and climate variation from each of the distal centres. However, these parameters do not give a true indication of risk. In a companion paper the effect of each hazard is fully investigated and the risk ranking completed.     

Radiative Heat Transfer and Hydrostatic Stability in Nocturnal Fog

We have performed a one-dimensional and transient radiative heat transfer analysis in order to investigate interaction between atmospheric radiation and convective instability within a nocturnal fog. The radiation element method using the Ray Emission Model (REM²), which is a generalized numerical method, in conjunction with a line-by-line (LBL) method, is employed to attain high spectral resolution calculations for anisotropically scattering fog. The results show that the convective instability has a strong dependence on radiative properties of the fog. For the condition of a 20-m droplet diameter and liquid water content of 0.1 × 10^–3 kg m^–3;, the temperature profile within the fog becomes S shaped, and a convective instability layer forms in the middle or lower level of the fog. However, for the same water content and a 40-m diameter droplet, no strong convective instability layer forms, whereas for a 10-m diameter droplet a strong convective instability is observed.     

上甸子本底站气溶胶散射系数变化特征的初步分析

根据上甸子大气成分本底站从2004～2006年的连续观测数据,分析了区域大气粒子光散射系数的变化特点,以了解本底地区大气气溶胶粒子的基本特征.卜甸子地区3年数据比较表明,散射系数平均水平较低,但在2006年有所升高,主要表现在受春季沙尘天气的影响,污染次数增多.其次,散射系数日变化规律表现为白大低夜间高,午后出现最低值,此变化与大气层结的日变化趋势一致.从季节变化看,冬季和春季的散射系数相对较低,而夏季和秋季的值较高,这与气象条件及内外源的影响都有关系.天气对粒子变化的影响表现为,晴天无云时的大气状况有利于污染物的重直输送,散射系数远远低于阴天时的数值.另一方面,风向对本底站粒子浓度影响也较明显.来自东北东方向的空气较干净,散射系数值通常较低,而西南西风向通常会引起散射系数值的增大,说明位于西南方向上的北京等城镇的污染输送对上甸子本底站的大气状况有一定影响.最后,通过拟和散射系数小时平均值出现频率曲线,对上甸子地区本底浓度值作了初步估算,其范围在10~20 Mm~(-1).     

Optical properties of dust from Jupiter family comets

To try to define specific physical properties of the dust of Jupiter-family comets (JFCs), we compare the light scattered by them. Amongst the more than 1000 JFCs, less than 200 are numbered, 40 of them being rather bright. In the present work we use data from the latter. In situ observations of three nuclei show low albedo surfaces. The albedo of the dust particles in the coma is low, with generally a red colour. The A(α)fρ product is a measure of cometary activity and secular changes. Images of different regions (jets and fans) give indications on the nucleus rotation and position of the emitting areas, as compared to the position of the rotation axis. Differences in physical properties between the particles in different regions are pointed out by differences in the linear polarization of the scattered light and by spectral variations in brightness and polarization. Jupiter family comets are considered as dust-poor comets. Tails and trails’ studies give an estimation of the size distribution of the particles. However the dust production rates depend on the largest particles (up to centimetre size), which are mainly observed in the trails where large dark compact particles are found. These dark particles are also responsible for the high polarization in the inner most coma of some comets. The meaning, in terms of physical properties, of the linear polarization is discussed through different examples such as 2P/Encke, 9P/Tempel 1 or the fragments of 73P/Schwassmann-Wachmann 3. Cometary outbursts and splitting events show that the properties of the dust ejected from the interior of the nucleus are similar to the ones of more active comets (new or with larger semi-major axis).     

Annual cycles in the interstellar scintillation time-scales of PKS B1519–273 and PKS B1622–253

We have used the University of Tasmania's 30-m radio telescope at Ceduna in South Australia to regularly monitor the flux density of a number of southern blazars. We report the detection of an annual cycle in the variability time-scale of the centimetre radio emission of PKS B1622−253. Observations of PKS B1519−273 over a period of nearly 2 yr confirm the presence of an annual cycle in the variability time-scale in that source. These observations prove that interstellar scintillation is the principal cause of inter-day variability at radio wavelengths in these sources. The best-fitting annual cycle model for both sources implies a high degree of anisotropy in the scattering screen and that it has a large velocity offset with respect to the local standard of rest. This is consistent with a greater screen distance for these 'slow' intra-day variability (IDV) sources than for rapid scintillators such as PKS B0405−385 or J1819+3845.     

Formation of the radio profile components of the Crab pulsar

The induced Compton scattering of radio emission off the particles of the ultrarelativistic electron–positron plasma in the open field line tube of a pulsar is considered. We examine the scattering of a bright narrow radio beam into the background over a wide solid angle and specifically study the scattering in the transverse regime, which holds in a moderately strong magnetic field and gives rise to the scattered component nearly antiparallel to the streaming velocity of the scattering particles. Making use of the angular distribution of the scattered intensity and taking into account the effect of rotational aberration in the scattering region, we simulate the profiles of the backscattered components as applied to the Crab pulsar. It is suggested that the interpulse (IP), the high-frequency interpulse (IP') and the pair of so-called high-frequency components (HFC1 and HFC2) result from the backward scattering of the main pulse (MP), precursor (PR) and low-frequency component (LFC), respectively. The components of the high-frequency profiles, the IP' and HFCs, are interpreted for the first time. The HFC1 and HFC2 are argued to be a single component split by the rotational aberration close to the light cylinder. It is demonstrated that the observed spectral and polarization properties of the profile components of the Crab pulsar as well as the giant pulse phenomenon outside the MP can be explained in terms of our model.     

Multiple solutions in preliminary orbit determination from three observations

Charlier’s theory (1910) provides a geometric interpretation of the occurrence of multiple solutions in Laplace’s method of preliminary orbit determination, assuming geocentric observations. We introduce a generalization of this theory allowing to take into account topocentric observations, that is observations made from the surface of the rotating Earth. The generalized theory works for both Laplace’s and Gauss’ methods. We also provide a geometric definition of a curve that generalizes Charlier’s limiting curve, separating regions with a different number of solutions. The results are generically different from Charlier’s: they may change according to the value of a parameter that depends on the observations.