We report results from an intensive multiwavelength monitoring campaign on the TeV blazar Mrk 421 over the period of [2003][2004]. The source was observed simultaneously at TeV energies with the Whipple 10 m telescope and at X-ray energies with the Rossi X-Ray Timing Explorer (RXTE) during each clear night within the Whipple observing windows. Supporting observations were also frequently carried out at optical and radio wavelengths to provide simultaneous or contemporaneous coverages. The large amount of simultaneous data has allowed us to examine the variability of Mrk 421 in detail, including cross-band correlation and broadband spectral variability, over a wide range of flux. The variabilities are generally correlated between the X-ray and gamma-ray bands, although the correlation appears to be fairly loose. The light curves show the presence of flares with varying amplitudes on a wide range of timescales at both X-ray and TeV energies. Of particular interest is the presence of TeV flares that have no coincident counterparts at longer wavelengths, because the phenomenon seems difficult to understand in the context of the proposed emission models for TeV blazars. We have also found that the TeV flux reached its peak days before the X-ray flux did during a giant flare (or outburst) in 2004 (with the peak flux reaching $135 mcrab in X-rays, as seen by the RXTE ASM, and $3 crab in gamma rays). Such a difference in the development of the flare presents a further challenge to both the leptonic and hadronic emission models. Mrk 421 varied much less at optical and radio wavelengths. Surprisingly, the normalized variability amplitude in the optical seems to be comparable to that in the radio, perhaps suggesting the presence of different populations of emitting electrons in the jet. The spectral energy distribution of Mrk 421 is seen to vary with flux, with the two characteristic peaks moving toward higher energies at higher fluxes. We have failed to fit the measured spectral energy distributions (SEDs) with a one-zone synchrotron self-Compton model; introducing additional zones greatly improves the fits. We have derived constraints on the physical properties of the X-ray/gamma-ray flaring regions from the observed variability (and SED) of the source. The implications of the results are discussed.

[1] We assess the accuracy and resolution of topography data provided by the Shuttle Radar Topography Mission (SRTM) through spectral comparisons with the National Elevation Dataset (NED) and a high-resolution laser data set of the 1999 Hector Mine earthquake rupture. We find that SRTM and the NED are coherent for wavelengths greater than 200 m, however the spatial resolution of the NED data is superior to the SRTM data for wavelengths shorter than 350 m, likely due to the application of a boxcar filter applied during final SRTM processing stages. From these results, a low-pass filter/decimation algorithm can be designed in order to expedite large-area SRTM applications.

[1] Stress accumulation rates along the primary segments of the San Andreas Fault system are computed using a three-dimensional (3-D) elastic half-space model with realistic fault geometry. The model is developed in the Fourier domain by solving for the response of an elastic half-space due to a point vector body force and analytically integrating the force from a locking depth to infinite depth. This approach is then applied to the San Andreas Fault system using published slip rates along 18 major fault strands of the fault zone. GPS-derived horizontal velocity measurements spanning the entire 1700 Â 200 km region are then used to solve for apparent locking depth along each primary fault segment. This simple model fits remarkably well (2.43 mm/yr RMS misfit), although some discrepancies occur in the Eastern California Shear Zone. The model also predicts vertical uplift and subsidence rates that are in agreement with independent geologic and geodetic estimates. In addition, shear and normal stresses along the major fault strands are used to compute Coulomb stress accumulation rate. As a result, we find earthquake recurrence intervals along the San Andreas Fault system to be inversely proportional to Coulomb stress accumulation rate, in agreement with typical coseismic stress drops of 1-10 MPa. This 3-D deformation model can ultimately be extended to include both time-dependent forcing and viscoelastic response.

[1] Exploring the earthquake cycle for large, complex tectonic boundaries that deform over thousands of years requires the development of sophisticated and efficient models. In this paper we introduce a semianalytic three-dimensional (3-D) linear viscoelastic Maxwell model that is developed in the Fourier domain to exploit the computational advantages of the convolution theorem. A new aspect of this model is an analytic solution for the surface loading of an elastic plate overlying a viscoelastic half-space. When fully implemented, the model simulates (1) interseismic stress accumulation on the upper locked portion of faults, (2) repeated earthquakes on prescribed fault segments, and (3) the viscoelastic response of the asthenosphere beneath the plate following episodic ruptures. We verify both the analytic solution and computer code through a variety of 2-D and 3-D tests and examples. On the basis of the methodology presented here, it is now possible to explore thousands of years of the earthquake cycle along geometrically complex 3-D fault systems.

We report on observations of very high energy γ rays from the shell-type supernova remnant (SNR) Cassiopeia A with the Very Energetic Radiation Imaging Telescope Array System stereoscopic array of four imaging atmospheric Cherenkov telescopes in Arizona. The total exposure time for these observations is 22 hr, accumulated between September and November of 2007. The γ-ray source associated with the SNR Cassiopeia A was detected above 200 GeV with a statistical significance of 8.3σ. The estimated integral flux for this γ-ray source is about 3% of the Crab-Nebula flux. The photon spectrum is compatible with a power law dN/dE ∝ E −Γ with an index Γ = 2.61 ± 0.24 stat ± 0.2 sys. The data are consistent with a point-like source. We provide a detailed description of the analysis results and discuss physical mechanisms that may be responsible for the observed γ-ray emission.