Some factors affecting the growth and decay of plages

The Mount Wilson coarse array magnetograph data set is analyzed to examine the dependence of growth and decay rates on the tilt angles of the magnetic axes of the regions. It is found that there is a relationship between these quantities which is similar to that found earlier for sunspot groups. Regions near the average tilt angle show larger average (absolute) growth and decay rates. Thepercentage growth and decay rates show minima (in absolute values) at the average tilt angles because the average areas of regions are largest near this angle. This result is similar to that derived earlier for sunspot groups. As in the case of spot groups, this suggests that, for decay, the effect results from the fact that the average tilt angle may represent the simplest subsurface configuration of the flux loop or loops that make up the region. In the case of region growth, it was suggested that the more complicated loop configuration should result in increased magnetic tension in the flux loop, and thus in a slower ascent of the loop to the surface, and thus a slower growth rate. In order to examine this further, the growth and decay rates of plage regions were examined as functions of the magnetic complexity of the regions. In the case of decay, the result was as expected from the model suggested above - that is, the more complex regions decayed more slowly. But for growing regions the effect is the opposite to that expected (more complex regions grow faster, even in terms of percentage growth), so the explanation of the tilt angle effect for growing regions proposed earlier may not be valid.Operated by the Association of Universities for Research in Astronomy, Inc., under Cooperative Agreement with the National Science Foundation.     

Measurement of Kodaikanal white-light images

A program of digitization of the daily white-light solar images from the Kodaikanal station of the Indian Institute of Astrophysics is in progress. A similar set of white-light data from the Mount Wilson Observatory was digitized some years ago. In both cases, areas and positions of individual sunspot umbrae are measured. In this preliminary report, comparisons of these measurements from the two sites are made. It is shown that both area and position measurements are in quite good agreement. The agreement is sufficiently good that it is possible to measure motions and area changes of sunspots from one site to the next, involving time differences from about 12 hours to about 36 hours. This enables us to trace the motions of many more small sunspots than could be done from one site alone. Very small systematic differences in rotation rate between the two sites of about 0.4% are found. A portion of this discrepancy is apparently due to the difference in plate scales between the two sites. Another contributing factor in the difference is the latitude visibility of sunspots. In addition it is suggested that a small, systematic difference in the measured radii at the two sites may contribute a small amount to this discrepancy, but it has not been possible to confirm this hypothesis. It is concluded that in general, when dealing with high precision rotation results of this sort, one must be extremely careful about subtle systematic effects.Operated by the Association of Universities for Research in Astronomy, Inc., under Cooperative Agreement with the National Science Foundation.     

An Attempt to Detect Coronal Mass Ejections in Lyman-α Using SOHO Swan

In this study, the possibility that coronal mass ejections (CMEs) may be observed in neutral Lyman-α emission was investigated. An observing campaign was initiated for SWAN (Solar Wind ANisotropies), a Lyman-α scanning photometer on board the Solar and Heliospheric Observatory (SOHO) dedicated to monitoring the latitude distribution of the solar wind from its imprints on the interstellar sky background. This was part of SOHO Joint Observing Program (JOP) 159 and was an exploratory investigation as it was not known how, or even if, CMEs interact with the solar wind and interstellar neutral hydrogen at this distance (≈60 and 120 R _S). The study addresses the lack of methods for tracking CMEs beyond the field-of-view of current coronagraphs (30 R _S). In our first method we used LASCO, white-light coronagraphs on SOHO, and EIT, an extreme ultraviolet imaging telescope also on SOHO, to identify CME candidates which, subject to certain criteria, should have been observable in SWAN. The criteria included SWAN observation time and location, CME position angle, and extrapolated speed. None of the CME candidates that we discuss were identified in the SWAN data. For our second method we analyzed all of the SWAN data for 184 runs of the observing campaign, and this has yielded one candidate CME detection. The candidate CME appears as a dimming of the background Lyman-α intensity representing ≈10% of the original intensity, moving radially away from the Sun. Multiple candidate CMEs observed by LASCO and EIT were found which may have caused this dimming. Here we discuss the campaign, data analysis technique and statistics, and the results.     

Electrical properties of Titan's surface from Cassini RADAR scatterometer measurements

We report regional-scale low-resolution backscatter images of Titan's surface acquired by the Cassini RADAR scatterometer at a wavelength of 2.18-cm. We find that the average angular dependence of the backscatter from large regions and from specific surface features is consistent with a model composed of a quasi-specular Hagfors term plus a diffuse cosine component. A Gaussian quasi-specular term also fits the data, but less well than the Hagfors term. We derive values for the mean dielectric constant and root-mean-square (rms) slope of the surface from the quasi-specular term, which we ascribe to scattering from the surface interface only. The diffuse term accommodates contributions from volume scattering, multiple scattering, or wavelength-scale near-surface structure. The Hagfors model results imply a surface with regional mean dielectric constants between 1.9 and 3.6 and regional surface roughness that varies between 5.3° and 13.4° in rms-slope. Dielectric constants between 2 and 3 are expected for a surface composed of solid simple hydrocarbons, water ice, or a mixture of both. Smaller dielectric constants, between 1.6 and 1.9, are consistent with liquid hydrocarbons, while larger dielectric constants, near 4.5, may indicate the presence of water-ammonia ice [Lorenz, R.D., 1998. Icarus 136, 344-348] or organic heteropolymers [Thompson, W.R., Squyres, S.W., 1990. Icarus 86, 336-354]. We present backscatter images corrected for angular effects using the model residuals, which show strong features that correspond roughly to those in 0.94-μm ISS images. We model the localized backscatter from specific features to estimate dielectric constant and rms slope when the angular coverage is within the quasi-specular part of the backscatter curve. Only two apparent surface features are scanned with angular coverage sufficient for accurate modeling. Data from the bright albedo feature Quivira suggests a dielectric constant near 2.8 and rms slope near 10.1°. The dark albedo feature Shangri-La is best fit by a Hagfors model with a dielectric constant close to 2.4 and an rms slope near 9.5°. From the modeled backscatter curves, we find the average radar albedo in the same linear (SL) polarization to be near 0.34. We constrain the total-power albedo in order to compare the measurements with available groundbased radar results, which are typically obtained in both senses of circular polarization. We estimate an upper limit of 0.4 on the total-power albedo, a value that is significantly higher than the 0.21 total albedo value measured at 13 cm [Campbell, D., Black, G., Carter, L., Ostro, S., 2003. Science 302, 431-434]. This is consistent with a surface that has more small-scale structure and is thus more reflective at 2-cm than 13-cm. We compare results across overlapping observations and observe that the reduction and analysis are repeatable and consistent. We also confirm the strong correlations between radar and near-infrared images.     

Using an inverse modelling approach with equifinality control to investigate the dominant controls on snowmelt nutrient export

There is great interest in modelling the export of nitrogen (N) and phosphorus (P) from agricultural fields because of ongoing challenges of eutrophication. However, the use of existing hydrochemistry models can be problematic in cold regions because models frequently employ incomplete or conceptually incorrect representations of the dominant cold regions hydrological processes and are overparameterized, often with insufficient data for validation. Here, a process‐based N model, WINTRA, which is coupled to a physically based cold regions hydrological model, was expanded to simulate P and account for overwinter soil nutrient biochemical cycling. An inverse modelling approach, using this model with consideration of parameter equifinality, was applied to an intensively monitored agricultural basin in Manitoba, Canada, to help identify the main climate, soil, and anthropogenic controls on nutrient export. Consistent with observations, the model results suggest that snow water equivalent, melt rate, snow cover depletion rate, and contributing area for run‐off generation determine the opportunity time and surface area for run‐off–soil interaction. These physical controls have not been addressed in existing models. Results also show that the time lag between the start of snowmelt and the arrival of peak nutrient concentration in run‐off increased with decreasing antecedent soil moisture content, highlighting potential implications of frozen soils on run‐off processes and hydrochemistry. The simulations showed TDP concentration peaks generally arriving earlier than NO₃ but also decreasing faster afterwards, which suggests a significant contribution of plant residue Total dissolved Phosphorus (TDP) to early snowmelt run‐off. Antecedent fall tillage and fertilizer application increased TDP concentrations in spring snowmelt run‐off but did not consistently affect NO₃ run‐off. In this case, the antecedent soil moisture content seemed to have had a dominant effect on overwinter soil N biogeochemical processes such as mineralization, which are often ignored in models. This work demonstrates both the need for better representation of cold regions processes in hydrochemical models and the model improvements that are possible if these are included.     

Simulating the development of Martian highland landscapes through the interaction of impact cratering, fluvial erosion, and variable hydrologic forcing

On the highlands of Mars early in the history of the planet precipitation-driven fluvial erosion competed with ongoing impact cratering. This disruption, and the multiple enclosed basins produced by impacts, is partially responsible for a long debate concerning the processes and effectiveness of fluvial erosion. The role of fluvial erosion in sculpting the early Martian landscape is explored here using a simulation model that incorporates formation of impact craters, erosion by fluvial and slope processes, deposition in basins, and flow routing through depressions. Under assumed arid hydrologic conditions, enclosed basins created by cratering do not overflow, drainage networks are short, and fluvial bajadas infill crater basins with sediment supplied from erosion of interior crater slopes and, occasionally from adjacent steep slopes. Even under arid conditions adjacent crater basins can become integrated into larger basins through lateral erosion of crater rims or by rim burial by sediment infilling. Fluvial erosion on early Mars was sufficient to infill craters of 10 km or more in diameter with 500–1500 m of sediment. When the amount of runoff relative to evaporation is assumed to be larger, enclosed basins overflow and deeply incised valleys interconnect basins. Examples of such overflow and interconnection on the Martian highlands suggest an active hydrological cycle on early Mars, at least episodically. When fluvial erosion and cratering occur together, the drainage network is often disrupted and fragmented, but it reintegrates quickly from smaller impacts. Even when rates of impact are high, a subtle fluvial signature is retained on the landscape as broad, smooth intercrater plains that feature craters with variable amounts of infilling and rim erosion, including nearly buried “ghost” craters. The widespread occurrence of such intercrater plains on Mars suggests a strong fluvial imprint on the landscape despite the absence of deep, integrated valley networks. Indisputable deltas and alluvial fans are rare in the crater basins on Mars, in part because of subsequent destruction of surficial fluvial features by impact gardening and eolian processes. Simulations, however, suggest that temporally-varying lake levels and a high percentage of suspended to bedload supplied to the basins could also result in poor definition of fan–delta complexes.     

The solar flares of February 1986 and the ensuing intense geomagnetic storm

A very intense geomagnetic storm, the largest observed in 26 years, was observed in early February 1986 having just been preceded by a series of six solar flares during the period 3–7 February. The storm and its antecedent flares are currently a subject of great interest because of the unusually large magnitude of the various geomagnetic effects that obtained. The fact that the flares were moderate to large in soft X-ray intensity, but much smaller than the largest that the Sun is capable of producing, coupled with the fact that these events occurred near the minimum of the current solar activity cycle, adds to the uniqueness of the overall episode.This paper describes the special circumstances surrounding these events and offers an interpretation of the cause and effect relationships through a numerical simulation of the dynamical evolutionary processes that may have occurred in interplanetary space.     

Observations of coronal disturbances from 1 to 9 R ⊙

Observations of a coronal disturbance on 1973 January 11 commencing at 18^h01^m UT are described. The event is homologous with an earlier disturbance from the same region of the corona. The observations suggest that a cloud of coronal gas containing 4 × 10³⁹ electrons propagated outwards to 5 R behind a piston-driven shock wave travelling at a velocity of 800 to 1200 km s^–1.On leave from Division of Radiophysics, CSIRO, Sydney, Australia.     

Center-to-limb profiles of the aluminum autoionization lines in the solar spectrum

The rotation of the solar corona has been studied using recurrence properties of the green coronal line (5303 Å) for the interval 1947–1970. Short-lived coronal activity is found to show the same differential rotation as short-lived photospheric magnetic field features. Long-lived recurrences show rigid rotation in the latitude interval ±57°.5. It is proposed that at least part of the variability of rotational properties of the solar atmosphere may be understood as a consequence of coexistence of differential and rigid solar rotation.On leave from Torino University, Italy, as an ESRO-NASA Fellow.