We present a measurement of the Hubble Constant based on Cepheid distances to 27 galaxies within 20 Mpc. We take the Cepheid data from published measurements by the Hubble Telescope Key Project on the Distance Scale (H0KP). We calibrate the Cepheid Period-Luminosity (PL) relation with data from over 700 Cepheids in the LMC obtained by the OGLE collaboration; we assume an LMC distance modulus of 18.50 mag (d=50.1 kpc). Using this PL calibration we obtain new distances to the H0KP galaxies. We correct the redshifts of these galaxies for peculiar velocities using two distinct velocity field models: the phenomenological model of Tonry et al. and a model based on the IRAS density field and linear gravitational instability theory. We combine the Cepheid distances with the corrected redshifts for the 27 galaxies to derive H_0, the Hubble constant. The results are H_0 = 85 +/- 5 km/s/Mpc (random error) at 95% confidence when the IRAS model is used, and 92 +/- 5 km/s/Mpc when the phenomenological model is used. The IRAS model is a better fit to the data and the Hubble constant it returns is more reliable. Systematic error stems mainly from LMC distance uncertainty which is not directly addressed by this paper. Our value of H_0 is significantly larger than that quoted by the H0KP, H_0 = 71 +/- 6 km/s/Mpc. Cepheid recalibration explains ~30% of this difference, velocity field analysis accounts for ~70%. We discuss in detail possible reasons for this discrepancy and future study needed to resolve it.

We study the reliability of dark matter (DM) profile determinations for rich clusters of galaxies based on weak gravitational lensing measurements. We consider three functional forms for the mass distribution: a softened isothermal sphere (SIS), the Navarro, Frenk, & White model (NFW) and the Klypin, Kravstov, Bullock, & Primack model (KKBP). A maximum likelihood parametric method is developed and applied to simulated images of background galaxies (z=2.0) lensed by clusters at z=0.2, 0.5, and 0.8. The images have field size (~10 arcmin) and angular resolution (~0.5 arcsec) characteristic of the new generation of 8-10 meter-class ground-based telescopes. The method recovers the correct virial mass of the DM halo, M_{200}, to within ~10% at all redshifts. The recovered values of \delta_c and r_s are only accurate to ~25%; however, their correlated errors are such that greater accuracy is obtained for the mass. Clusters with an SIS density profile are recognized by the method, while NFW and KKBP clusters are recognized as not being SIS. The method does not discriminate between NFW and KKBP profiles for clusters of a given mass, because they are too similar in their outer regions. Strongly lensed arcs near the cluster centers would be needed to distinguish NFW from KKBP.

We describe the scientific motivation behind, and the methodology of, the Stanford Cluster Search (StaCS), a program to compile a catalog of optically selected clusters of galaxies at intermediate and high (0.3 < z < 1) redshifts. The clusters are identified using an matched filter algorithm applied to deep CCD images covering approximately 60 square degrees of sky. These images are obtained from several data archives, principally that of the Berkeley Supernova Cosmology Project of Perlmutter et al. Potential clusters are confirmed with spectroscopic observations at the 9.2 m Hobby-Eberly Telescope. Follow-up observations at optical, sub-mm, and X-ray wavelengths are planned in order to estimate cluster masses. Our long-term scientific goal is to measure the cluster number density as a function of mass and redshift, which is sensitive to the cosmological density parameter, and the amplitude of density fluctuations on cluster scales. Our short-term goals are the detection of high-redshift cluster candidates over a broad mass range and the measurement of evolution in cluster scaling relations. The combined data set will contain clusters ranging over an order of magnitude in mass, and allow constraints on these parameters accurate to ~10%. We present our first spectroscopically confirmed cluster candidates and describe how to access them electronically.

We present the first results from the Shellflow program, an all-sky Tully-Fisher (TF) peculiar velocity survey of 276 Sb-Sc galaxies with redshifts between 4500 and 7000 km/s. Shellflow was designed to minimize systematic errors between observing runs and between telescopes, thereby removing the possibility of a spurious bulk flow caused by data inhomogeneity. A fit to the data yields a bulk flow amplitude V_bulk = 70{+100}{-70} km/s (1 sigma error) with respect to the Cosmic Microwave Background, i.e., consistent with being at rest. At the 95% confidence level, the flow amplitude is < 300 km/s. Our results are insensitive to which Galactic extinction maps we use, and to the parameterization of the TF relation. The larger bulk motion found in analyses of the Mark III peculiar velocity catalog are thus likely to be due to non-uniformities between the subsamples making up Mark III. The absence of bulk flow is consistent with the study of Giovanelli and collaborators and flow field predictions from the observed distribution of IRAS galaxies.

I review the status of cosmic velocity analysis as of January 2000, with an emphasis on two key questions: (1) What is the scale of the largest bulk flows in the universe? and (2) What is the value of beta=Omega_m^{0.6}/b indicated by cosmic velocities, and what does this tell us about Omega_m itself? These are the most important issues for cosmic flow analysis, and each has been controversial in recent years. I argue that a preponderance of the evidence at present argues against very large scale (> 100 Mpc) bulk flows, and favors beta = 0.4-0.5, corresponding to a low-density (Omega_m=0.2-0.3) universe.

We present preliminary results from the Shellflow program, an all-sky Tully-Fisher (TF) survey of 297 Sb-Sc with redshifts between 4000 and 7500 km/s. The program was designed to ensure uniformity of the data between observing runs and between telescopes, thereby eliminating the possibility of a spurious bulk flow caused by data inhomogeneity. A simple bulk flow model suggests that the Shellflow galaxies are nearly at rest with respect to the cosmic microwave background (CMB). Taken at face value, this result suggests that most of the ~600 km/s motion of the Local Group with respect to the CMB is due to material within ~60 Mpc/h -- a striking confirmation of the gravitational instability picture and the notion of homogeneity in the Universe. Differences between our Shellflow ``null'' result and claims of larger bulk motions based on analyses of the MarkIII catalog are possibly due to small inhomogeneities between the sub-samples in MarkIII. Our preliminary Shellflow bulk velocity, V_bulk = 80+/-150 km/s in the CMB frame, must still be refined with Monte-Carlo simulations and tidal field analysis.

The closely related issues of (1) reconstructing density and velocity fields from redshift survey and distance indicator data and (2) determining cosmological parameters from such data sets are discussed. I focus on possible explanations for the discrepant values of $\beta\equiv\Omega^{0.6}/b$ resulting from analyses that employ various reconstruction methods. Although no firm conclusions are reached, possible routes towards resolution of these discrepancies are suggested by the discussion.

I describe a Tully-Fisher (TF) survey of galaxies in 15 Abell Clusters distributed around the sky in the redshift range 9000-13,000 km/sec. The observations were done between 1992 and 1995 at the Las Campanas (LCO) and Palomar Observatories. The program is known as the LP10K survey. The basic data set consists of R band CCD magnitudes and H-alpha rotation curves for over 200 spiral galaxies. The TF relation is found to have a mild but statistically significant surface brightness dependence, v_{TF} ~ L^{0.28} I_e^{0.13}, and the intrinsic TF scatter is found to be 0.28 +/- 0.07 mag. A maximum-likelihood method is applied to the LP10K data set to test for large-scale bulk streaming. A formal result of v_B=720 +/- 300 km/sec is found. However, this finding is in disagreement with several results on smaller scales, including the Shellflow result reported at the Cosmic Flows conference, which indicate convergence to the CMB frame by ~ 6000 km/sec redshift.

The implications of the massive, X-ray selected cluster of galaxies MS1054--03 at $z=0.83$ are discussed in light of the hypothesis that the primordial density fluctuations may be nongaussian. We generalize the Press-Schechter (PS) formalism to the nongaussian case, and calculate the likelihood that a cluster as massive as MS1054 would appear in the EMSS. The probability of finding an MS1054-like cluster depends only on $\omegam$ and the extent of primordial nongaussianity. We quantify the latter by adopting a specific functional form for the PDF, denoted $\psi_\lambda,$ which tends to Gaussianity for $\lambda\gg 1,$ and show how $\lambda$ is related to the more familiar statistic $T,$ the probability of $\ge 3\sigma$ fluctuations for a given PDF relative to a Gaussian. We find that Gaussian initial density fluctuations are consistent with the data on MS1054 only if $\omegam\simlt 0.2.$ For $\omegam\ge 0.25$ a significant degree of nongaussianity is required, unless the mass of MS1054 has been substantially overestimated by X-ray and weak lensing data. The required amount of nongaussianity is a rapidly increasing function of $\omegam$ for $0.25 \le \omegam \le 0.45,$ with $\lambda \le 1$ ($T \simgt 7$) at the upper end of this range. For a fiducial $\omegam=0.3,$ $\omegal=0.7$ universe, favored by several lines of evidence we obtain an upper limit $\lambda \le 10,$ corresponding to a $T\ge 3.$ This finding is consistent with the conclusions of Koyama, Soda, & Taruya (1999), who applied the generalized PS formalism to low ($z\simlt 0.1$) and intermediate ($z\simlt 0.6$) redshift cluster data sets.

The LCO/Palomar 10,000 km/sec (LP10K) Tully-Fisher (TF) data set is used to test for bulk streaming motions on a ~150 Mpc scale. The sample consists of 172 cluster galaxies in the original target range of the survey, 9000-13,000 km/sec, plus an additional 72 galaxies with cz < 30,000 km/sec. A maximum-likelihood analysis that is insensitive to Malmquist and selection bias effects is used to constrain the bulk velocity parameters, and realistic Monte-Carlo simulations are carried out to correct residual biases and determine statistical errors. When the analysis is restricted to the original target range, the bias-corrected bulk flow is v_B=720 +/- 280 km/sec toward l=266, b=19. When all objects out to z=0.1 are included the result is virtually unchanged, v_B=700 +/- 250 km/sec toward l=272, b=10. The hypothesis that the Hubble flow has converged to the CMB frame at distances less than ~ 100 Mpc is ruled out at the 97% confidence level. The data are inconsistent with the flow vector found by Lauer & Postman. However, the LP10K bulk flow is consistent with that obtained from the SMAC survey of elliptical galaxies recently described by Hudson et al. If correct, the LP10K results indicate that the convergence depth for the Hubble flow is >~ 150 Mpc.

We present an improved POTENT method for reconstructing the velocity and mass density fields from radial peculiar velocities, test it with mock catalogs, and apply it to the Mark III Catalog. Method improvments: (a) inhomogeneous Malmquist bias is reduced by grouping and corrected in forward or inverse analyses of inferred distances, (b) the smoothing into a radial velocity field is optimized to reduce window and sampling biases, (c) the density is derived from the velocity using an improved nonlinear approximation, and (d) the computational errors are made negligible. The method is tested and optimized using mock catalogs based on an N-body simulation that mimics our cosmological neighborhood, and the remaining errors are evaluated quantitatively. The Mark III catalog, with ~3300 grouped galaxies, allows a reliable reconstruction with fixed Gaussian smoothing of 10-12 Mpc/h out to ~60 Mpc/h. We present maps of the 3D velocity and mass-density fields and the corresponding errors. The typical systematic and random errors in the density fluctuations inside 40 Mpc/h are \pm 0.13 and \pm 0.18. The recovered mass distribution resembles in its gross features the galaxy distribution in redshift surveys and the mass distribution in a similar POTENT analysis of a complementary velocity catalog (SFI), including the Great Attractor, Perseus-Pisces, and the void in between. The reconstruction inside ~40 Mpc/h is not affected much by a revised calibration of the distance indicators (VM2, tailored to match the velocities from the IRAS 1.2Jy redshift survey). The bulk velocity within the sphere of radius 50 Mpc/h about the Local Group is V_50=370 \pm 110 km/s (including systematic errors), and is shown to be mostly generated by external mass fluctuations. With the VM2 calibration, V_50 is reduced to 305 \pm 110 km/s.

The first results from a Tully-Fisher (TF) survey of cluster galaxies are presented. The galaxies are drawn from fifteen Abell clusters that lie in the redshift range 9000-12,000 km/sec and are distributed uniformly around the celestial sky. The data set consists of R-band CCD photometry and long- slit H-alpha spectroscopy. The rotation curves (RCs) are characterized by a turnover radius (r_t) and an asymptotic velocity v_a, while the surface brightness profiles are characterized in terms of an effective exponential surface brightness I_e and a scale length r_e. The TF scatter is minimized when the rotation velocity is measured at 2.0 +/- 0.2 r_e; a significantly larger scatter results when the rotation velocity is measured at > 3 or < 1.5 scale lengths. This effect demonstrates that RCs do not have a universal form, as has been suggested by Persic, Salucci, and Stel. In contrast to previous studies, a modest but statistically significant surface-brightness dependence of the TF relation is found, log v = const + 0.28*log L + 0.14*log I_e. This indicates a stronger parallel between the TF relation and the FP relations of elliptical galaxies than has previously been recognized. Future papers in this series will consider the implications of this cluster sample for deviations from Hubble flow on 100-200 Mpc scales.

We present a new optical Tully-Fisher (TF) investigation for a complete, full-sky sample of 297 Sb-Sc spirals with redshifts between 4500 and 7000 km/s. The survey was specifically designed to provide uniform, well- calibrated data over both hemispheres. All previous TF surveys within the Supergalactic shell (cz < 6000 km/s) have relied on matching separate data sets in the Northern and Southern hemispheres and thus cannot attain full- sky homogeneity. Analyses of the cosmological dipole and peculiar velocities based on these studies have produced contradictory claims for the amplitude of the bulk flow and whether it is generated by internal or external mass fluctuations. With Shellflow, and further zero-point calibration of existing TF data sets, we expect a high-accuracy detection of the bulk flow amplitude and an unambiguous characterization of the tidal field at 6000 km/s. One of the key features of this study is not only its careful, all-sky sample selection but the independent duplication, by at least 2 of us, of all the data reductions. These reductions and a first flow analysis based on the Shellflow sample alone will be published in the Fall.

This is the second in a series of papers in which we compare Tully-Fisher (TF) data from the Mark III Catalog with predicted peculiar velocities based on the IRAS galaxy redshift survey and gravitational instability theory, using a rigorous maximum likelihood method called VELMOD. In Paper I (Willick et al. 1997b), we we applied the method to a $cz_{LG} \leq 3000$ km/sec, 838-galaxy TF sample and found $\beta_I=0.49\pm 0.07,$ where $\beta_I\equiv \Omega^{0.6}/b_I$ and $b_I$ is the linear biasing parameter for IRAS galaxies. In this paper we increase the redshift limit to $cz_{LG}=7500$ km/sec, thereby enlarging the sample to 1876 galaxies. The expanded sample now includes the W91PP and CF subsamples of the Mark III catalog, in addition to the A82 and MAT subsamples already considered in Paper I. We implement VELMOD using both the forward and inverse forms of the TF relation, and allow for a more general form of the quadrupole velocity residual detected in Paper I. We find $\beta_I=0.50\pm 0.04$ (1-sigma error) at 300 km/sec smoothing of the IRAS-predicted velocity field. The fit residuals are spatially incoherent for $\beta_I=0.5,$ indicating that the IRAS plus quadrupole velocity field is a good fit to the TF data. If we eliminate the quadrupole we obtain a worse fit, but a similar value for $\beta_I$ of $0.54\pm 0.04.$ Changing the IRAS smoothing scale to 500 km/sec has almost no effect on the best $\beta_I.$ We find evidence for a density-dependence of the small-scale velocity dispersion, $\sigma_v(\delta_g)\simeq (100 + 35 \delta_g)$ km/sec.

We compare Tully-Fisher (TF) data for 838 galaxies within cz=3000 km/sec from the Mark III catalog to the peculiar velocity and density fields predicted from the 1.2 Jy IRAS redshift survey. Our goal is to test the relation between the galaxy density and velocity fields predicted by gravitational instability theory and linear biasing, and thereby to estimate $\beta_I = \Omega^{0.6}/b_I,$ where $b_I$ is the linear bias parameter for IRAS galaxies. Adopting the IRAS velocity and density fields as a prior model, we maximize the likelihood of the raw TF observables, taking into account the full range of selection effects and properly treating triple-valued zones in the redshift-distance relation. Extensive tests with realistic simulated galaxy catalogs demonstrate that the method produces unbiased estimates of $\beta_I$ and its error. When we apply the method to the real data, we model the presence of a small but significant velocity quadrupole residual (~3.3% of Hubble flow), which we argue is due to density fluctuations incompletely sampled by IRAS. The method then yields a maximum likelihood estimate $\beta_I=0.49\pm 0.07$ (1-sigma error). We discuss the constraints on $\Omega$ and biasing that follow if we assume a COBE-normalized CDM power spectrum. Our model also yields the 1-D noise noise in the velocity field, including IRAS prediction errors, which we find to be be 125 +/- 20 km/sec.

This is the third in a series of papers in which we assemble and analyze a homogeneous catalog of peculiar velocity data. In Papers I and II, we described the Tully-Fisher (TF) redshift-distance samples that constitute the bulk of the catalog, and our methodology for obtaining mutually consistent TF calibrations for these samples. In this paper, we supply further technical details of the treatment of the data, and present a subset of the catalog in tabular form. The full catalog, known as the Mark III Catalog of Galaxy Peculiar Velocities, is available in accessible on-line databases, as described herein. The electronic catalog incorporates not only the TF samples discussed in Papers I and II, but also elliptical galaxy Dn-sigma samples originally presented elsewhere. The relative zeropointing of the elliptical and spiral data sets is discussed here. The basic elements of the Mark III Catalog are the observables for each object (redshift, magnitude, velocity width, etc.) and inferred distances derived from the TF or Dn-sigma relations. Malmquist bias-corrected distances are computed for each catalog object using density fields obtained from the IRAS 1.2 Jy redshift survey. Distances for both individual objects and groups are provided. A variety of auxiliary data, including distances and local densities predicted from the IRAS redshift survey reconstruction method, are tabulated as well. We study the distributions of TF residuals for three of our samples and conclude that they are well-approximated as Gaussian. For the Mathewson et al. sample we demonstrate a significant decrease in TF scatter with increasing velocity width. We test for, but find no evidence of, a correlation between TF residuals and galaxy morphology.

Six of the principal galaxy distance indicators are discussed: Cepheid variables, the Tully-Fisher relation, the Dn-sigma relation, Surface Brightness Fluctuations, Brightest Cluster Galaxies, and Type Ia Supernovae. The role they play in peculiar velocity surveys and Hubble constant determination is emphasized. Past, present, and future efforts at constructing catalogs of redshift-independent distances are described. The chapter concludes with a qualitative overview of Malmquist and related biases.

We consider a measure of the peculiar velocity field derived from the Mark III compilation of 2900 spiral galaxies (Willick \etal 1996b), using an analysis that is substantially free of bias (Nusser and Davis 1995). We expand the velocity field in a set of orthogonal, smooth modes, reducing the data to a set of 56 coefficients fitted to a maximum redshift of 6000 km/s, and maximum spherical harmonic of l=3. Equivalent mode coefficients can be computed for the gravity field derived from any whole-sky redshift catalog of galaxies, such as the \iras 1.2 Jy survey (Fisher \etal 1995). Detailed intercomparison shows the two independent fields to be remarkably aligned in general. There are, however, systematic discrepancies in the fields that lead to coherence in the residuals between them. These residuals take the form of a dipole field in the LG frame that grows with distance; it is not consistent with a bulk flow residual. A chi-squared analysis of the mode-mode comparison designed to determine which value of $\beta\equiv \Omega^{0.6}/b$ for the \iras gravity field best fits the Mark III velocity field shows instead that no value of $\beta$ is acceptable. This is in contrast to the simulated catalogs, and suggests that there must be some systematic discrepancy between the two fields.

We describe the construction of an N-body simulation that mimics the true velocity and mass-density fields in a box of side 256\hmpc about the Local Group , and the production of mock catalogs that mimic in detail current catalogs of redshifts and peculiar velocities. Our main purpose is to provide a tool for developing and testing reconstruction methods, but the different components of the method can be used on their own in other applications. The initial conditions are based on the IRAS 1.2Jy redshift survey, assuming that galaxies trace mass and \Omega=1. A density field smoothed is recovered from the redshift survey, using quasi-linear theory and a power-preserving filter. The corresponding potential field is traced back to the linear regime using the Zel'dovich-Bernoulli equation. Small-scale power is added by means of constrained realization to mimic fluctuations on galactic scales. The gravitating system is evolved forward in time with a PM code of 2\hmpc resolution. The result reproduces the real dynamical structures on large scales and the statistical properties of the structure down to galactic scales. "Galaxies" are identified via a linear biasing scheme (b=1.35) and they are divided into "spirals" and "ellipticals" to obey Dressler's morphology-density relation. The galaxies are assigned internal-velocity parameters and absolute magnitudes scattered about an assumed mean Tully-Fisher relation. They are then "observed" as magnitude-limited samples, trying to mimic the selection criteria of the data sets constituting the Mark III catalog of peculiar velocities. Artificial IRAS 1.2Jy redshift surveys are also compiled. The simulations and mock catalogs will be made available electronically as bench marks for testing reconstruction methods.

We develop a straightforward and quantitative two-step method for spectroscopically classifying galaxies from the low signal-to-noise (S/N) optical spectra typical of galaxy redshift surveys. First, using \chi^2-fitting of characteristic templates to the object spectrum, we determine the relative contributions of the old stellar component, the young stellar component, and various emission line spectra. Then, we classify the galaxy by comparing the relative strengths of the components with those of galaxies of known morphological type. In particular, we use the ratios of (1) the emission line to absorption line contribution, (2) the young to old stellar contribution, and (3) the oxygen to hydrogen emission line contribution. We calibrate and test the method using published morphological types for 32 galaxies from the long-slit spectroscopic survey of Kennicutt (1992) and for 304 galaxies from a fiber spectroscopic survey of nearby galaxy clusters. From an analysis of a sample of long-slit spectra of spiral galaxies in two galaxy clusters, we conclude that the majority of the galaxies observed in the fiber survey are sufficiently distant that their spectral classification is unaffected by aperture bias. Our spectral classification is consistent with the morphological classification to within one type (e.g. E to S0 or Sa to Sb) for \gtsim 80% of the galaxies. Disagreements between the spectral and morphological classifications of the remaining galaxies reflect a divergence in the correspondence between spectral and morphological types, rather than a problem with the data or method.

We review the quantitative science that can be and has been done with redshift and peculiar velocity surveys of galaxies in the nearby universe. After a brief background setting the cosmological context for this work, the first part of this review focuses on redshift surveys. The practical issues of how redshift surveys are carried out, and how one turns a distribution of galaxies into a smoothed density field, are discussed. Then follows a description of major redshift surveys that have been done, and the local cosmography out to 8,000 km/s that they have mapped. We then discuss in some detail the various quantitative cosmological tests that can be carried out with redshift data. The second half of this review concentrates on peculiar velocity studies, beginning with a thorough review of existing techniques. After discussing the various biases which plague peculiar velocity work, we survey quantitative analyses done with peculiar velocity surveys alone, and finally with the combination of data from both redshift and peculiar velocity surveys. The data presented rule out the standard Cold Dark Matter model, although several variants of Cold Dark Matter with more power on large scales fare better. All the data are consistent with the hypothesis that the initial density field had a Gaussian distribution, although one cannot rule out broad classes of non-Gaussian models. Comparison of the peculiar velocity and density fields constrains the Cosmological Density Parameter. The results here are consistent with a flat universe with mild biasing of the galaxies relative to dark matter, although open universe models are by no means ruled out.

The mass density field in the local universe, recovered by the POTENT method from peculiar velocities of $\sim$3000 galaxies, is compared with the density field of optically-selected galaxies. Both density fields are smoothed with a Gaussian filter of radius 12 $h^{-1}$ Mpc. Under the assumptions of gravitational instability and a linear biasing parameter $b\sbo$ between optical galaxies and mass, we obtain $\beta\sbo \equiv \om^{0.6}/b\sbo = 0.74 \pm 0.13$. This result is obtained from a regression of POTENT mass density on optical density after correcting the mass density field for systematic biases in the velocity data and POTENT method. The error quoted is just the $1\sigma$ formal error estimated from the observed scatter in the density--density scatterplot; it does not include the uncertainty due to cosmic scatter in the mean density or in the biasing relation. We do not attempt a formal analysis of the goodness of fit, but the scatter about the fit is consistent with our estimates of the uncertainties.

We have combined five Tully-Fisher (TF) redshift-distance samples for peculiar velocity analysis: the cluster data of Han, Mould and coworkers (1991-93, HM) and Willick (1991, W91CL), and the field data of Aaronson et al. (1992), Willick (1991), Courteau & Faber (1992), and Mathewson et al. (1992), totaling over 3000 spiral galaxies. We treat the cluster data in this paper, which is the first of a series; in Paper II we treat the field TF samples. These data are to be combined with elliptical data (e.g., Faber et al. 1989) to form the MARK III CATALOG OF GALAXY PECULIAR VELOCITIES, which we will present in Paper III. The catalog will be used as input for POTENT reconstruction of velocity and density fields, described in later papers, as well as for alternative velocity analyses. Our main goal in Papers I & II is to place the TF data onto a self-consistent system by (i) applying a uniform set of corrections to the raw observables, (ii) determining the TF slopes and scatters separately for each sample, and (iii) adjusting the TF zeropoints to ensure mutually consistent distances. The global zeropoint is set by the HM sample, chosen because of its depth and uniformity on the sky and its substantial overlap with each of the other samples. In this paper, we calibrate the ``forward'' and ``inverse'' TF relations for HM and W91CL. We study the selection criteria for these samples and correct for the resultant statistical biases. The bias corrections are validated by comparing forward and inverse cluster distances. We find that many sample clusters are better modeled as ``expanding'' than relaxed, which significantly affects the TF calibrations. Proper corrections for internal extinction are derived self-consistently from the data.