amongst the highest quality data available for exploring galaxy evolution are the long exposures of the hubble deep field ( williams et al .
1996 ) , which register the faintest and sharpest images of galaxies ever detected .
the counts of faint galaxies are found to increase to the completeness limit ( @xmath10 ) , roughly doubling in number per magnitude , @xmath11 ( williams et al .
1996 ; see tables 9 - 10 ) , in all bands .
this trend is accompanied by a marked decline in angular size , so that most faint galaxies are barely resolved .
while much of this light must arise from the epoch when galaxies and stars were first forming , the photometric limit of the hdf far exceeds the practical limit of spectroscopy ( some @xmath12 of the hdf galaxies do not have measured redshifts ) , hence for the foreseeable future interpretation of these images rests principally on photometric information .
a major impediment in analysing the hdf images is our ignorance of the ultraviolet ( uv ) properties of ordinary galaxies , since at high redshift it is the restframe uv light that is detected in the optical bands . despite the rather modest uv performance of imaging satellites
it has become apparent that the appearance of even ordinary hubble - sequence galaxies is knotty and irregular in the uv , governed by the spatial distribution of high contrast star - forming regions ( oconnell & markum 1996 ; giavalisco et al .
hence , it is unclear to what extent the large numbers of blue and irregular looking objects reported in hst images ( glazebrook et al .
1995 ; driver , windhorst , & griffiths 1995 ; abraham et al .
1996b ; colley et al.1997 ) simply result from redshifted uv light rather than evolution .
there have been several efforts to obtain images of representative nearby galaxies .
the foca balloon experiment at 2000 @xmath5 has measured six uv images of local galaxies ( blecha et al .
1990 ) , the astro-1 ( astro-2 ) mission measured 27 ( 45 ) uv images of objects using the ultraviolet imaging telescope ( stecher et al .
1992 ) , and maoz et al .
( 1997 ) have recently compiled uv images of 110 galaxies using the foc on hst ( @xmath13 ) .
cognizant of the importance of these uv images for assessing galaxy evolution , researchers ( bohlin et al.1991 ; giavalisco et al .
1996a ; maoz 1997 ) have already used these data to comment on deep hst images .
while these studies have been very illustrative , the samples are either incomplete or register only the inner regions of galaxies so that a quantification of galaxy evolution in deep hst images has hitherto not been attempted .    due to the inadequacies of uv instrumentation , it is interesting to note that even the gross effect of the lyman - limit is more easily observed in the _ optical _ , by virtue of the high redshifts of some faint galaxies ( steidel et al .
the restframe uv morphologies of several tens of these high redshift galaxies are measured directly with hst , particularly in the hubble deep field ( steidel et al .  1996b ; giavalisco et al .
1996b ; lowenthal et al .
1997 ) down to the lyman limit .
these galaxies lie in the range @xmath14 and are reported to have a generally compact appearance , but many have faint plumes and some are up to 3  in length ( object c4 - 06 of steidel et al .
the uv continuum flux places some constraint on the luminosity density of star - forming galaxies at high - z ( madau et al.1996 ) at least for galaxies blue enough to be dominated by the light of o - stars .
comparison of the implied high mass star - formation rate with local estimates from from h@xmath15 emission ( gallego et al.1995 ) , provides an indirect means of estimating the evolution of the integrated uv - luminosity density ( madau et al .
1996 ) .    in this paper
, we outline a simple and practical method for measuring evolution in deep space images , by constructing empirical `` no - evolution '' fields of galaxies as a benchmark , matched in image properties to the data , in this case the hdf .
our method has the distinct advantages of being purely empirical and entirely free of parameterizations , so that our conclusions are liberated from the usual caveats .
we first establish the k - correction for each pixel belonging to a given galaxy image drawn from a complete sample of relatively bright galaxies with redshifts ( selected from the hdf ) .
we then scale the numbers of these k - corrected galaxy images out to high redshift in proportion to the product of the volume and the space density of each galaxy , the latter being set by its value of @xmath0 . finally , we correct the psf , add the appropriate noise and sampling of the hdf observations and compare these `` no - evolution '' fields with the hdf at fainter magnitudes ( @xmath6 ) . this approach is feasible due to the relatively high quality @xmath16 images of the brighter galaxies detected in the hdf ( @xmath17 ) , as well as the many subsequent redshift measurements of hdf galaxies with the keck telescope ( cohen et al .
1996 ; lowenthal et al .
1997 ) . together
these observations allow us to construct a statistically useful complete sample of 31 galaxies with sufficient uv - optical coverage for construction of an entirely empirical and parameter free `` no - evolution '' deep field for comparison with the hdf .
the universal geometry is the only unknown , dictating the volume and angle - redshift relations , and bearing somewhat on the degree of evolution derived in the @xmath18 band .    by following this procedure
, we make no recourse to the usual parameterizations of morphological or colour classes , galaxy luminosity functions , and light profiles .
implicitly , such parameterizations have tended to assume one - parameter relationships , particularly between morphological type and spectral type , reducing the inherent richness of the general galaxy population and potentially introducing artificial model dependencies . in light of this , it is not very surprising that simulations of two - dimensional images look artificial , appearing more `` hygienic '' than reality , and hence of limited value in reliably quantifying evolution . designating a morphological class
has in itself proven to be problematic for well - resolved galaxies .
naim et al .  ( 1995 ) , for example , have found from a comparison of the classification of six human ` experts ' that classification exhibits a variance of 2 morphological types relative to one another .
hence there is still no clear way to map an observed image onto a given morphological type or to produce a realistic image given a morphological type .
finally , we emphasize that our procedure is free of the many internal selection biases which have long haunted low - redshift to high - redshift comparisons , e.g. , where deep ccd data have lower surface brightness thresholds than photographic surveys used to construct the local luminosity functions .
we begin this paper by describing the manner in which we select our bright galaxy sample in  2 . in  3 , we discuss our empirical procedure for deriving pixel - by - pixel spectral energy distributions and number densities for each of these bright galaxies . in  4 , we assess both the fairness of our sample and the viability of the method by making some simple predictions chiefly with regard to the canada france redshift survey ( cfrs ; lilly et al .
1995 ) . in  5
, we describe the simulation procedure in detail , and in  6 we describe our analysis of the simulations and the hdf data using sextractor . in  7 we present our results , in  8 we discuss these results , and in  9 we present a summary .
we adopt @xmath19 and express all magnitudes in this paper in the ab ( oke 1974 ) ] magnitude system ( defined in terms of a flat spectrum in frequency ) .
also , to associate the hdf bands with their more familiar optical counterparts , we shall refer to the @xmath20 , @xmath21 , @xmath22 , and @xmath16 bands as @xmath23 , @xmath24 , @xmath25 , and @xmath26 , respectively , throughout this paper .
we select our sample of bright galaxies from the hdf .
first we measure the photometry of all galaxies detected in the hdf by applying the sextractor photometry package 1.2b5 ( bertin & arnouts 1996 ) to the publicly released version 2 images of the hdf . in using sextractor
we require objects to be at least 2 @xmath27 above the sky noise in the @xmath23 image over an area of at least 10 contiguous pixels , after first mildly smoothing the images by a gaussian of width 0.06-arcsec sigma ( the approximate width of the psf ) .
we take the deblending parameter ( @xmath28 ) to be 0.04 , which we found to be a good compromise between splitting too many faint objects apart and merging too many different objects together . for photometry
, we use isophotal magnitudes ( @xmath29 ) corrected with a gaussian extrapolation to approximate the light beyond the isophotes for all non - crowded ( almost all ) objects , where we take our isophotes to be equal to 2 @xmath27 times the sky noise ( @xmath30 ) .
objects with half - light radii less than 0.15 arcsec are excluded as stars ( 5 objects with @xmath31 ) .
only objects on the three wide field camera chips lying within 4 arcsecs of the edges are included , beyond which the image quality is diminished .
we match up our resulting photometric catalogues redshifts measured by various groups ( cohen et al .
1996 ; lowenthal et al .
1997 ) to construct a strictly magnitude - limited subset , bright enough to be highly complete in redshift .
we settle for a conservative limit of @xmath31 , which contains 32 objects in total  only one of which is without a redshift and for which we know of no attempt to measure a redshift nor any reason to suppose that its measurement would prove difficult .
consequently , we shall assume that the 31 objects with redshifts are representative , and simply scale volume densities by @xmath32 .
as more redshifts become available this analysis can be extended , although clearly the limited area of the hdf means that new fields are required in order to significantly enlarge the bright galaxy sample and to help average over the possible effects of clustering .
the resultant bright sample is listed in table 1 , along with our determination of its coordinates , @xmath23 magnitude , absolute magnitude , surface brightness ( taken to be @xmath33 ) , ( where @xmath34 is the rest - frame @xmath35 apparent magnitude and @xmath36 is the half - light radius ) , our eye - ball determination of the morphological type , and the known redshift . for reference
we also include our determination of the k - correction over the surface of each galaxy at @xmath37 .
the calculation of k - corrections , rest - frame absolute magnitudes , and rest - frame surface brightnesses are described below .
in addition , the redshift distribution of our bright subset is shown in figure 1 , for comparison with all of the redshifts measured in the compilation by cohen et al .
( 1996 ) , which includes the hdf and its flanking fields .
we display the images of the bright sample in figures 2a and figure 2b , as a set of colour images ( @xmath25 , @xmath24 , and @xmath23 bands ) as they appear if k - corrected to the mean bright sample redshift , @xmath38 , and also to @xmath37 , using the pixel - by - pixel k - corrections described in the next section .
we treat our bright galaxies as two - dimensional pixelated light - emitting surfaces for the purpose of k - correction , avoiding the unnecessary step of parameterizing their light profiles .
we do not attempt to alter the inclination of each galaxy in the simulations , with the attendant problems regarding extinction this would entail , since by default our sample already contains galaxies of all orientations .
the question arises as to the best method of k - correction when the redshift extends to wavelengths short enough that no uv information is available .
the spectrum of any given pixel is only sampled by 4 broad passbands and no information exists shortward of 2000 @xmath5 in the restframe of the typical bright galaxy ( for which @xmath39 ) , allowing a reliable k - correction to @xmath40 in the longest wavelength band , @xmath23 . in order to extend this work to higher redshift and in order to implement a smooth interpolation between the four passbands we adopt empirical template spectra . to indicate the robustness of the results to the choice of templates , we compare simulations for two independent sets of templates , one based on the small but large aperture uv sample of coleman , wu & weedman ( 1980 ) and the other on the larger but small aperture sample of template spectra carefully compiled and combined from the iue archive by kinney et al .
these data sets provide useful uv spectra for a range of optically selected galaxies nicely spanning the observed spectral range ( see connolly et al .
1996 ) , so that a reasonable empirical interpolation and extrapolation can be performed by matching these templates to the 4 passbands of the hdf and improved by interpolating between these spectra to provide a smooth result .
additionally , we include the average effects of the hi lyman continuum and series forest absorption at high redshift , as parameterized by madau ( 1995 ) for distant qso s , since they significantly affect the broad - band colours of high redshift galaxies .
formally , we represent each element @xmath41 in the two - dimensional template of a galaxy by @xmath42 where @xmath43 is the surface brightness template at element @xmath44 in the @xmath45 band and where @xmath46 is the spectral energy distribution shape at the element @xmath44 .
we take the surface brightness template to equal that observed in the hdf in band @xmath45 , the band with the highest integrated signal - to - noise ( either the @xmath23 or the @xmath24 band for the objects in our sample . ) for clarity we use the word `` element '' when referring to a pixel in the original bright galaxy image because these pixels are transformed in area and flux by redshift , and therefore distinct from the pixel scale of the simulations which is set by the hdf .
we take the sed shape , @xmath46 , for each element @xmath44 from two different compilations of spectral templates : the coleman , wu , & weedman ( 1980 ) ( hereinafter , cww ) set supplemented with the ngc4449 spectrum ( bruzual & ellis 1985 , unpublished ) and the kinney et al .
1996 ( k96 ) spectra .
cww includes seds intended to be representative of the dust - free e , sbc , scd , sdm , and starburst ( ngc4449 ) galaxies , and the k96 set includes seds intended to represent ellipticals , sa , sb , and sb1 ( starburst ; @xmath47 ) galaxies . we extend these observed sed templates below 1200 @xmath5 and 1400 @xmath5 respectively , by extrapolating the slope of the sed at these wavelengths down to the lyman break at 912 @xmath5 , below which we set the observed flux to zero ( see figure 3 ) .
we also interpolate linearly between spectra to form a smoothly continuous set for better fitting the observed 4 passband flux measurements .    integrating the seds over the transmission curve of the 4 passbands
, we then find the most - likely spectral template for each element @xmath44 , which we shall call @xmath48 , such that the sum of the squares of the differences between the sed template fluxes and the observed fluxes divided by the expected error in these fluxes is a minimum .
the expected error is taken to be the error in the value of the pixel fluxes added in quadrature with the typical cosmic variance in the model sed fluxes ( @xmath49 ) we find .
note that where the signal in a pixel is less than twice the sky noise , we have set the signal in that pixel equal to the value given by eq .
( 1 ) , adopting the surface brightness profile and the spectral energy template which best matches the mean value determined in a 1.0 arcsec diameter aperture .
we have not attempted to realign the hdf images in the different bands , nor have we tried to correct for the different shapes that the psf has in the different bands , since checks indicate that these effects are small and induce changes on scales smaller than the hdf pixel scale , i.e. , 0.04 arcsec .
this wavelength independence is due to the undersampling of images by wfpc , the effective hdf psf being more set by the wavelength independent sub - sampling grid pattern than by the wavelength dependent diffraction of the telescope .
clearly , we should take the spatial extent of each prototype image to be as large as possible within the limits of the signal . in practice , we take this to be the radius at which the mean pixel signal within an elliptical annulus is equal to the sky noise
. inspection of each image allows us to filter out the small number of obviously unrelated galaxies within this extended aperture , the pixels affected being replaced by their reflected counterparts . to be sure of including the `` whole '' object we extend its 2d profile by 50% , using the radial gradient . for this extended region
we produce the noise by adding poisson noise from modeled signal in quadrature with the general background noise . except for the very brightest `` cloned '' galaxies ,
this extrapolation is not relevant , since the vast majority of replicated images generated from a given template galaxy are much fainter than the original , lying at higher redshift where these outer regions of the profile lie further into the noise .
the question arises as to the utility of the hdf @xmath26 images for generating the higher redshift images in the redder bands , since the sensitivity of the redder bands of the hdf is greater than that of the @xmath26 band .
as it turns out , the signal in the hdf @xmath26 images is more than sufficient for our purposes , due to the strong @xmath50 cosmological surface brightness dimming . to see this ,
consider the s / n at a given @xmath27 in a pixel through the @xmath26 filter at the observed redshift .
this translates into a lower s / n at higher redshift in the other bands , given their relative sensitivities and exposure times . taking the worst case that the spectrum of a given pixel is flat in frequency ( though in practice the spectral indices of some starburst systems are somewhat steeper than this ) , the surface brightness dims by a factor @xmath51 , or @xmath52 , or a decrease of the surface brightness by 1.3 , 2.3 , and 3.3 @xmath53 , for the @xmath25 , @xmath24 and @xmath23 bands respectively .
however , the measured 1-@xmath27 @xmath26 noise level for the hdf is 32.4 @xmath54 , compared to the 32.2 , 34.0 , and 32.6 @xmath54 levels in the @xmath25 , @xmath24 , @xmath23 bands .
clearly , then , since few pixels are as blue or bluer than a flat spectrum in our template galaxy sample , the s / n of any pixel in a redder band at higher redshift will always be less than that of the observed @xmath26 image .    while it is apparent that a number of objects in our input sample have extremely low values of s / n in the @xmath26 band ( see the ellipticals in figure 2a )
, we emphasize that the s / n in @xmath26 is _ always _ sufficient to determine the appearance of these galaxies at higher redshifts in the redder bands of the hdf given the relative band sensitivities and exposure times used in these passbands .
greater inequities between bands would limit the simulations to depths less than the limiting magnitudes of these passbands .
we set the space density of each galaxy in our sample equal to @xmath55 , determined by the maximum redshift , @xmath56 , to which each galaxy could have potentially been selected given our chosen magnitude limit .
@xmath57 is given by @xmath58 where @xmath59 where @xmath60 is the solid angle of the surveyed area and where @xmath61 is the well - known luminosity distance at redshift @xmath56 .
@xmath56 is determined by : @xmath62 where @xmath63 is the observed redshift and @xmath64 is the @xmath65-correction in the @xmath23 band at redshift @xmath66 .
we determine the @xmath65-correction for the galaxy from the total spectral energy distribution , which incorporates the contribution of all the elements in our two - dimensional representation of each galaxy .
note this procedure avoids the usual steps of constructing a luminosity function , parameterizing it , and then assigning k - corrections .
we simply treat each galaxy as a class of its own with a space density set by its own value of @xmath0 , allowing a fully unbinned , individual treatment of each template galaxy .
astronomers have long used the @xmath3 distribution ( schmidt 1968 ) to assess the completeness and uniformity of a sample , where @xmath67 is the volume up to and including the redshift of the galaxy in question and @xmath57 is the volume in which this galaxy could have been observed , given the selection criteria . for a uniform distribution of galaxies in a volume one expects a uniform distribution of @xmath3 between 0 and 1 , in the absence of evolution and clustering , and for a sample with @xmath68 galaxies the average value of this quantity us given by : @xmath69 the distribution of @xmath3 for our template sample is shown in figure 4 , for @xmath70/@xmath71 , @xmath70 , and @xmath72 , where the k - correction is formed from the sed averaged over the whole object .
the average value for the @xmath3 distribution is 0.48 , 0.51 , and 0.53 for the @xmath70/@xmath73 , @xmath70 , and @xmath74 geometries , which is within the expected deviations for the @xmath3 statistic @xmath75 and not inconsistent with the rate of evolution detected in a larger redshift survey to a similar magnitude limit ( lilly et al .
we could have chosen to derive @xmath57 by redshifting each pixel of the object and performing the photometry on the redshifted image for greater self consistency , but this turns out to be virtually irrelevant because the magnitudes we recover from the redshifted images placed at @xmath56 are on average displaced by only 0.07 magnitudes faintward of our magnitude limit ( @xmath76 ) .
a slight faintward shift is expected for aperture magnitudes since the lower surface brightnesses at @xmath56 results in some small loss of the light on the wings . to illustrate this point
, we compare the recovered magnitudes at @xmath56 with the magnitude limit chosen for our bright sample in figure 5 .    to provide a basic context for understanding the no - evolution simulations described above we display the luminosity function obtained for this sample in figure 6 , for two different representative cosmologies .
we compare this luminosity function with the @xmath77-band luminosity functions determined by loveday et al .
( 1992 ) and zucca et al .
( 1997 ) in the apm and esp surveys , respectively correcting our sample to restframe @xmath77 directly using our pixel - by - pixel best - fit seds .
we see that our luminosity function is shifted to larger luminosity and/or space density than that of the local universe , broadly consistent with the findings of ellis et al .
( 1996 ) and lilly et al .
( 1995 ) where the luminosity function of blue objects is observed both to brighten and to steepen .
we can also compare our sample with the redshift distribution , n(z ) , of the canada - france redshift survey ( cfrs ) , a survey of similar depth , by simulating ( as described in detail below ) its survey parameters .
the cfrs sample covers @xmath78 over 112 @xmath79 .
scaling by these criteria and the 19% redshift incompleteness of the cfrs , we can construct a prediction for n(z ) using our method ( figure 7 ) .
we make a small correction , inferable from figure 5 of lilly et al .
( 1995 ) , to convert from isophotal to total magnitudes by a uniform 0.1 mag . a noticeable difference between our predicted n(z ) and that of the cfrs is found in the amplitude of the distributions , our redshift distribution being 28% higher , consistent with the differences in the respective luminosity functions ( figure 6 ) .
despite this , the shape of our predicted n(z ) is very similar to that of the cfrs , which is not surprising given the wide spread of @xmath3 for our sample ( figure 4 ) .
this basic agreement is quite satisfactory , showing that our bright galaxy sample fairly samples the range , if not quite the average density , of galaxies comprising the general field population .
according to the number density derived for each template galaxy , @xmath0 , and assuming homogeneity , we generate monte - carlo catalogues of the objects out to very large redshift ( @xmath80 ) over the solid angle of the hdf , each object being assigned a random redshift , position , and position angle .
we then generate mock images on the basis of these catalogues . each redshifted image
must be scaled in size and resampled with more noise and additional psf smoothing to account for its higher redshift and generally smaller size relative to its brighter counterpart .
since the two - dimensional surface brightness profiles of the bright galaxies already contain noise , we find it useful to simultaneously generate both a signal @xmath18 and noise @xmath68 image .
the noise image @xmath68 keeps track of how much noise has been implicitly added to each pixel in the signal image @xmath60 by virtue of each galaxy template implicitly containing noise .
such an accounting allows the proper amount of noise to be added to each pixel in the signal image after all the scaled galaxy templates have been laid down .    to generate an image at a chosen redshift
, we must calculate the change in size and surface brightness of each element @xmath44 in our two - dimensional galaxy templates . for a galaxy with redshift , @xmath66
, we take the angular size of each element in our two - dimensional galaxy template , @xmath81 , to be equal to @xmath82 where @xmath63 is the redshift of the object as measured in the hdf , @xmath83 is the angular distance ( peebles 1993 ) , and @xmath84 is the angular size of each element in our two - dimensional galaxy template image at @xmath63 , i.e. the pixel size of galaxies in the hdf ( 0.04 arcsec ) .
we calculate the surface brightness @xmath85 that each element @xmath44 of our two - dimensional galaxy template has at redshift @xmath66 for a given passband @xmath86 as @xmath87 where @xmath88 and @xmath89 are the @xmath65-corrections and @xmath90 @xmath91 colours calculated based on the spectral energy distribution @xmath92 of element @xmath44 .
we have included the line blanketing by the lyman - alpha forest as given by madau ( 1995 ) as well as lyman - limit absorption in our calculations of the @xmath65-corrections since these corrections become important at @xmath93 in the @xmath26 band and at @xmath94 in the @xmath25 band .    in calculating both the signal and the noise at each pixel on the images @xmath18 and @xmath68 , we break up each pixel into 25 smaller subpixels , and calculate the contribution of each redshifted element @xmath44 to each of these smaller subpixels to account for the generally much larger area covered by a data pixel than that of the redshifted image element , for the purpose of rebinning .
because the real psf has smoothed each galaxy template , to correctly calculate the appearance of each prototype galaxy in our sample at higher redshifts we must add more smoothing present in each object in the @xmath18 image , depending on the reduction in angular size for the object in question .
for simplicity , we simply smooth each object with a kernel derived from a relatively isolated , unsaturated star from the hubble deep field reduced in size so that its scale length is simply @xmath95 times that of the original scale length .
the above expression is exact for the case of a perfectly gaussian psf and is also close to exact in those cases where the angular size of the simulated galaxy laid down in the @xmath18 image is much smaller than the hdf .
of course , it is true that the real hdf psf differs from a gaussian in that it has much more extended wings , but for the most part the differences that this makes are small . to verify this
, we compared the angular sizes recovered using the present procedure and from assuming the psf to be exactly gaussian ( the @xmath27 s for which we determined by fitting to the same unsaturated star in the hdf ) , and we found little if any dependence on these differences .    for a pixel in which the contribution to the signal image @xmath18 from an element in the template image is @xmath96 , the noise contribution to the same pixel in the noise image @xmath68 is taken to be @xmath97 added in quadrature , where @xmath98 is the s / n ratio calculated for the element @xmath44 for the @xmath43 template drawn from the hdf , and account is made of both the background and poissonian noise associated with the measured pixel - by - pixel signal .
having generated the signal image , we want to calculate the appropriate noise for each pixel in our simulated hdf , which we take to be : @xmath99 where @xmath27 is the background noise level in the hdf , @xmath100 is the signal at pixel @xmath44 , and @xmath101 is the gain .
having already calculated the amount of noise which each pixel in the signal image had by virtue of implicit noise in the templates , i.e. the noise image @xmath68 , we can bring the noise in each pixel up to the appropriate value @xmath102 by adding pixel - by - pixel gaussian - distributed noise with standard deviation @xmath103 , smoothing this added noise with the noise kernel specified in table 4 of williams et al .  (
1996 ) so as to approximate the observed correlation properties of the noise in the hdf .
while the outlined procedure will add the appropriate amount of noise to pixels if noise is lacking in those pixels , it is possible that some pixels will already have more noise added to them than is present per pixel in the hdf .
in particular , this occurs when the redshift for an object is lower than that of the prototype , since of course we can not improve on the s / n of the prototype galaxy image . to artificially prevent this possibility
, we simply replace every galaxy in our mock catalogue whose redshift is lower than its observed redshift , with another galaxy from our bright sample placed at the same redshift whose observed redshift is lower than the mock catalogue redshift .
though this results in a bright galaxy population which is slightly unrepresentative ( and , in fact , this is the cause of the slight disagreement observed in the distribution of bright hdf galaxies and the simulations found later in the paper ) , it does not appreciably bias the measured properties of the cloned galaxies at the fainter magnitudes of interest .
the simulations cover a sky area four times that of the hdf in each of the four broadbands @xmath26 , @xmath25 , @xmath24 , and @xmath23 .
we perform these simulations self - consistently for @xmath104/@xmath71 , @xmath70 , and @xmath74 , by which we mean that the @xmath57 of the bright sample uses the same geometry as the volume used in constructing the simulation . in figure 8
we compare a no - evolution simulation assuming @xmath72 , @xmath104 , and @xmath105/@xmath106 with an area of the same size selected from the hdf .
a simple test of our procedure is to construct images of the 31 prototype galaxies at their observed redshifts ( @xmath63 ) .
good agreement is found between these images and their originals at @xmath63 , within the errors expected on the basis of our pixel - by - pixel sed fits . to illustrate this agreement more quantitatively
, we plot a comparison of the recovered values for the half - light radius , petrosian radius ( petrosian 1976 ) , and apparent @xmath23 magnitudes in figure 9 .
on both the simulations and the hdf itself , we perform the object identification and photometry using sextractor version 1.2b5 ( bertin & arnouts 1996 ) . after smoothing the @xmath23 images with a gaussian of 0.06-arcsec radius ( the approximate psf ) within sextractor
, we require that objects be at least 2 sigma above the noise over an area of at least 10 contiguous pixels , and we use a cleaning parameter of 1.0 .
we select our apparent magnitudes to be equal to sextractor s @xmath29 estimate , which for non - crowded objects is equivalent to an isophotal magnitude extrapolated beyond the isophotes .
we exclude objects with half - light radii less than 0.15 arcsec and @xmath107 as stars , fainter than which the contamination is relatively small . then , based upon the detected objects , we used sextractor to determine the apparent magnitudes in the other bands with the @xmath23 apertures .
we set the deblending parameter ( @xmath28 ) , which is important in determining the extent to which sextractor breaks up objects , to be equal to 0.04 .
we found some dependence of our results on the value chosen for this parameter , particularly the break - up rate , but for the most part this dependence was small .
our chosen value of the deblending parameter is very close to that ( 0.05 ) used by clements & couch ( 1997 ) ( couch , private communication ) . with this choice of parameter
, we obtained the reasonable result that no bright galaxy from our sample broke up into smaller pieces when placed at its @xmath56 , in agreement with our qualitative impressions in looking at these images .
we derive two different measures of the angular size for each object in our image , a petrosian radius ( petrosian 1976 ) and a half - light radius .
we take the petrosian radius to equal the smallest radius for which the surface brightness at that radius equals half the average surface brightness interior to that radius .
the half - light radius is equal to the radius of the aperture which contained half the light as determined by sextractor s best estimate of the total light .
we performed some simulations to test our method for recovering half - light radii and found that our recovered half - light radii are only slightly scattered about sizes 10% smaller than the input half - light radii at @xmath108 .
the scatter in this relationship partially derives from the uncertainty in the overall photometry performed by sextractor .
it is natural to begin our comparison of the faint galaxy population in the hdf by looking at the counts in @xmath23 , the longest wavelength band and hence least affected by uncertainty in the k - correction .
figure 10 compares the number counts recovered from both our simulations and the observations ( a description of the plotted @xmath109 estimates based on the size of our bright sample is given in appendix a. ) for all geometries , we find that our no - evolution predictions fall steadily short of the observations with increasing apparent magnitude . for @xmath105/@xmath106 ,
the shortfall is noticeably less pronounced due to the relatively larger volume available . at @xmath110 ,
for example , the counts fall short by a factor @xmath111 for @xmath105/@xmath106 compared to @xmath112 for @xmath113 and @xmath114 for @xmath106 .
of course , the counts recovered at bright magnitudes , @xmath115 , are in fair agreement with the observations , as one would expect given the definition of our sample up to approximately the same bright limit , @xmath76 .
it is reassuring to find that there little dependence of the number counts on the choice of spectral template , whether it be the cww set or the k96 set . finally ,
as discussed in appendix b , we note that the present no - evolution faint counts ( @xmath116 ) are compromised of a non - negligible number of the input prototype galaxies ( @xmath117 ) .
given our knowledge of the positions , redshifts , magnitudes , and types of galaxies laid down in each simulated image , it is simple to determine systematic uncertainties such as incompleteness in the number counts , by matching up the object catalogues from the simulations created by sextractor with our input monte - carlo catalogues .
figure 11 shows that the incompleteness becomes significant in the range @xmath118 for both our no - evolution simulations .
this incompleteness stems from the fact that at fainter magnitudes , we are sensitive to smaller and hence the higher surface brightness galaxies , so that large fraction of the population of galaxies at faint magnitudes subtend areas too large for detection given their redshifted surface brightnesses .    in a similar manner , we can determine the rate of overcounting of in our simulated fields .
as explained in colley et al .  ( 1997 ) , it is possible to overcount the number of faint galaxies by misidentifying individual parts of a galaxy as distinct galaxies , especially at higher redshifts , where the rest - frame ultraviolet light is accessed and hii regions have a higher contrast ( cf . the galaxy at @xmath119 in figure 2b ) .
for simplicity , rather than perform an angular correlation analysis like colley et al .
( 1997 ) , we have chosen to measure this number directly . in figure 12
, we display the rate of overcounting for all the simulations performed by comparison of the input random catalogue with those recovered .
clearly , in our no - evolution simulations , using the parameters chosen for the photometry , overcounting is never an important effect for any of the cosmologies examined .
the clustering seen in the hdf redshift data at bright magnitudes ( figure 1 ) is responsible for the overdensity at bright magnitudes @xmath120 in the hdf relative to the mean field counts measured in wider field surveys .
compared to the cfrs , for example , this amounts to a @xmath121 overdensity ( figure 7 ) . at fainter magnitudes ,
however , the count variance should be less affected by clustering , given the greater projected volume and the expectation of less well developed structure at earlier times .
we provide a quick test of the plausibility of this hypothesis using other deep hst @xmath23 data .
figure 13 shows the count variance between 6 pointings ( two sets of 3 hst pointings separated by 3 hours on the sky from another program ) .
these data show that the count variance in @xmath23 band decreases steadily to the poissonian limit ( @xmath122 ) over the area of wfpc2 , consistent with isotropy at faint magnitudes .
we compare the distributions of half - light radii recovered from the hdf with those of our `` no - evolution '' simulations in figure 14 .
the hatched area represents the @xmath123 uncertainty range based on the finite size of our bright sample for the no - evolution model using the cww sed templates .
the solid curve indicates the angular size distribution recovered from simulations using the k96 sed templates . at bright magnitudes ( @xmath124 ) ,
the angular sizes recovered from the simulations agree quite well with the observations as expected given the fact that we defined our bright sample in terms of many of these same galaxies .
in contrast , at fainter magnitudes , the half - light radii from the no - evolution simulations become significantly larger than those from the observations : larger by 43% , 44% , and 57% in the magnitude range @xmath125 and larger by 35% , 42% , and 53% in the magnitude range @xmath126 for the @xmath127/@xmath70 , @xmath70 , and @xmath74 geometries , respectively , where the median of the angular size distribution in each magnitude bin is used for comparison . as expected
, we see that galaxies in the @xmath71/ @xmath70 geometry tend to possess smaller angular sizes than galaxies in the @xmath70 geometry and especially galaxies in the @xmath72 geometry because of the somewhat larger angular - diameter distances .
we repeat the above comparisons using the petrosian radii instead of the half - light radii .
ideally , petrosian radii provide a more reliable surface brightness - independent estimator of the angular size than half - light radii , though one might question its meaningfulness given the increasingly lumpy appearance and small sizes of faint galaxies . in any case , an inspection of figure 15 shows that the same general trends and conclusions hold here as for the half - light radii .      in figure 16 - 17
, we plot a comparison between the colour distributions recovered from the observations and our ` no - evolution ' simulations in three different magnitude ranges for two different colours , @xmath128 and @xmath129 .
good agreement is found at the bright magnitudes ( @xmath124 ) as expected given that our bright sample is selected up to this magnitude limit . at fainter magnitudes , however , a clear excess of bluer galaxies is observed relative to our no - evolution simulations . in addition , in the faintest magnitude bin ( @xmath126 ) , there also appears to be an excess of red galaxies relative to that found in the no - evolution simulations .
this is somewhat unexpected because one would expect the real universe , for which our sample is representative at @xmath9 to have a younger and therefore bluer appearance than that of our extrapolated sample .
this excess may , therefore , reflect the presence of dust and lyman - series forest absorption at moderate to high - redshift . conceivably , an ad - hoc maximal dwarf - model , might also account for this red excess in terms of faded dwarfs at low redshift .      in figure 18 , we plot the redshift distributions of the galaxies recovered by sextractor by matching them up with our input catalogues .
it is apparent that in the absence of evolution very few galaxies lie beyond a redshift @xmath130 , even at the faintest magnitude , with little dependence on the choice of geometry .
a dramatic illustration of the relative unobservability of high redshift galaxies for @xmath70 is given in figure 19 , where we have broken up the simulation of a 52  x 72  hdf exposure into 4 redshift slices .
it is interesting to see how these distributions compare with the photometric redshift estimates derived from fits to the 4 passbands of the hdf , and from subsequent near - ir imaging of this field from the ground .
surprisingly enough , the redshift distribution recovered from our `` no - evolution '' simulations agree remarkably well with the approximate distribution of lanzetta , yahil , & fernandez - soto ( 1996 ) . given that similar coincidences have been found for redshift surveys at somewhat brighter magnitude limits ( broadhurst et al .
1988 ; colless et al .  1993
; cowie et al .  1996 ) , it is somewhat tempting to suppose that this result might hold to yet fainter magnitudes . despite such hopes ,
clearly this relationship begins to break down in the @xmath8 band at fainter magnitudes ( @xmath131 ) , at least , as cowie et al .
( 1996 ) have shown , and presumably a similar break down will occur in the redder bands at fainter magnitudes .      in figures 20 and 21 , we plot several colour - colour diagrams , and overplot the @xmath26-dropout and @xmath25-dropout selection criteria given by madau et al .
( 1996 ) , criteria , useful for selecting galaxies at redshifts @xmath132 and @xmath133 , respectively .
we count the number of recovered galaxies satisfying these criteria and provide a sumary in table 2 . as expected from the recovered redshift distributions plotted in figure 18 , it is clear that the no - evolution simulations contain manifestly fewer dropouts than the observations .
clearly then , the high - redshift universe contains a larger number of objects which of greater luminosity in the ultraviolet than the galaxy population in our redshift - complete sample .
we have compared the numbers of dropouts in the hdf found by madau et al .
( 1996 ) with our values listed in table 2 , using the same colour - magnitude window .
our dropout rate is higher for the hdf , a finding which might result from our use of sextractor for photometry rather than focas as employed by madau et al .
given that these programs likely differ with regard to the degree to which they successfully estimate the magnitudes of detected objects , our analyses probably probe slightly different depths , even though we choose for comparison the same nominal magnitude limit ( @xmath134 ) .
imposing a magnitude limit @xmath135 brighter , we can reproduce their dropout rate .
also , as stated by madau et al .  ( 1996 ) , their focas magnitudes are in general estimated to be too faint by about @xmath136 .      in the light of recent work attempting to quantify the morphology of faint galaxies , we compare the properties of the faint galaxies from the hdf against those from our simulations with the @xmath137 statistic proposed by abraham et al .
( 1996a,1996b ) , in order to provide an approximate estimate of the extent to which evolution may affect the asymmetry of galaxies , where @xmath137 is defined as the sum of the absolute value of the differences between a galaxy image and itself rotated 180 degrees about the center of the image .
the approximate contribution of the noise to the apparent asymmetry is calculated and subtracted .
we plot the distribution of this asymmetry parameter recovered both from our simulations and from the hdf in figure 22 , for two different magnitude bins .
we find a clear trend toward larger values of this parameter with increasing magnitude . at the simplest level
, this would seem to imply that the faint galaxies are less smooth than the redshifted field population at @xmath138 used in our simulations .
some caution must be exercised in interpreting the @xmath137 statistic since it appears to be extremely sensitive to the manner in which one determines the center about which to rotate for evaluating the statistic ( abraham & brinchmann , private communication ) and can exhibit a fair amount of scatter ( @xmath139 ) depending upon whether one takes the centroid or the maximum as the for the center . because of this , the statistic , or at least the present implementation , seems to become increasingly unreliable for the smallest galaxies and hence we restrict our comparison to the bright magnitudes ( @xmath140 ) where we are more confident in the reliability of the statistic from our own internal tests .    note that along similar lines , based upon pixel - by - pixel k - corrections of the frei , guhathakurta , & gunn ( 1996 ) sample and systematic applications of asymmetry ( @xmath137 ) & central - concentration ( @xmath141 ) statistics , abraham et al .
( 1996a,1996b ) have argued that the fainter galaxy population is quite unlike the local galaxy population broadly represented by the frei et al .
( 1996 ) sample , which is nevertheless somewhat ill - defined , lacking a firm magnitude limit or well - defined selection criteria .
since the novelty of the hdf is principally its angular size information obtained to unprecedentedly faint magnitudes , it is not too surprising that our most interesting finding relates to the sizes of the faint images .
the count excess is clearly composed of galaxies with smaller projected areas than expected on the basis of relatively low - redshift galaxies ( @xmath2 ) for any interesting geometry .
the evolution of angular - size has not been completely clear in the literature to date .
previously , in the range @xmath142 , roche et al .
( 1996 ) find an apparent excess of galaxies with small angular sizes relative to their no - evolution models which are designed to replicate the angular sizes of galaxies in the local universe . at brighter magnitudes ( @xmath143 ) , the study by i m et al .
( 1995 ) on the prerefurbished hst data indicates that the angular sizes of this population is also smaller than expected based upon more local , brighter measures .
however , there have also been some studies based on this parallel data indicating no changes in size ( mutz et al .
1994 ; casertano et al .  1995 ) .
we also find that the number density of dropouts discovered at faint magnitudes exceeds that recovered from our no - evolution simulation , a finding which clearly indicates that the real high - redshift universe contains a larger density of galaxies with high uv surface brightnesses .
note , however , that despite this finding , the estimated integrated star - formation rates in these redshift ranges are similar ( madau et al .
1996 ) . clearly , then
, as argued in steidel et al.(1996 ) and lowenthal et al .
( 1997 ) ( cf . figure 7 ) , more star - formation took place on relatively small scales at high - redshift .
it remains unclear whether this is the result of the intrinsically smaller baryonic aggregations or the increasing prominence of starbursts at high redshifts .
our finding that overcounting rates do not seem to be very important ( see also the companion paper : bouwens , broadhurst , & silk 1998 ) is somewhat contrary to the claims of colley et al .
( 1997 ) , even though we used a similar program ( sextractor ) with similar values of the deblending parameter as one catalogue ( clements & couch 1997 ) used in their study .
of course , the galaxy population we have used for our study is more evolved than expected in the real universe at higher redshifts , so we might expect our value to be an underestimate .
note that in any case , the small scale correlation claimed is a statistically marginal result , ( @xmath144 , colley et al 1997 ) .
less contentious are our findings regarding the colours and the number counts . in agreement with most authors
we find clear evidence that the predictions without evolution are redder than observed ( lilly et al .
1995 ) and that the number counts fall short of the data at faint magnitudes ( lilly 1993 ; pozzetti et al .
1996 ) .
as we have already argued , it is possible to see the present work as a demonstration that the trends found in comparing the local ( @xmath145 ) population with intermediate - redshift galaxies ( @xmath9 ) extend to yet fainter magnitudes and presumably higher redshifts . at these intermediate magnitudes a larger number of faint , blue , irregular , and presumably small galaxies ( ellis et al .  1996 ; lilly et al .
1995 ; brinchmann et al .  1998 ; guzman et al .
1997 ) have been reported at @xmath146 relative to local expectations .
the present study demonstrates that the these trends extend to yet fainter magnitudes , in that there continues to be both an increase in the number and apparent starburst activity of the faint galaxy population as well as a decrease in mean galaxy size .
uncertainty estimates for the present study are somewhat `` internal '' and , therefore , might underestimate a field - averaged variance in our measurements .
certainly , at bright magnitudes the hdf is @xmath121 overdense with respect to the cfrs .
however , since 90% of our bright galaxy sample ranges over a wide spread in distance , some 1500 mpc / h , and @xmath3 is evenly spread over the interval [ 0,1 ] , most bright galaxies should be considered as independent .
clearly , we would prefer a larger , more local sample of redshift - complete bright galaxies for our simulations , and other deep multicolour hdf imaging would be useful in the first two respects .
we have developed a technique for generating a model - independent faint galaxy population based on pixel - by - pixel k - corrected images of the bright galaxy population , for comparison with deep high - resolution images .
our technique has the virtue of being model - independent and completely empirical , except for the usual choice of geometry .
we have made use of the hdf for the purposes of defining a redshift - complete bright galaxy sample from which to construct empirical simulations and also for evaluating the evolution of the much fainter galaxies detected in the hdf .
we have made use of all four passbands and have been able to make concrete statements about the evolution of the image properties with estimated uncertainties based on the size of a volume - limited sample .
we find that , relative to our simulations based on the bright galaxy sample , the faint galaxies are smaller , more numerous , bluer , less regular , and contain more dropouts for any interesting geometry . of particular note is our finding with regard to angular sizes , for which a variety of nebulous and seemingly contradictory statements may be found in the literature .
our bright galaxy sample is small , and it will be important to see if our conclusions remain valid when a larger sample is available . nonetheless , we believe this work is an important foundation on which to build , in particular , as we look forward to independent deep fields and the very significant improvement in the quality of uv and optical imaging promised by the advanced camera ( ford et al .
1996 ) .
we would like to thank daniela calzetti , marc davis , mike fall , holland ford , andy fruchter , nick kaiser , and alex szalay for some very useful conversations , emmanuel bertin for answering several of our questions regarding sextractor , gordon squires and nick kaiser for the use of several routines from their software package imcat , harry ferguson for his help in producing colour images and answering a few of our questions , jarle brinchmann and roberto abraham for a number of useful tips and comments with regard to the routine we incorporated into sextractor for measuring the @xmath137 parameter , and finally harry ferguson , steve zepf , eric gawiser , and jonathan tan for some helpful comments on near - final drafts of this document .
rjb acknowledges support from an nsf graduate fellowship , tjb acknowledges the nasa grant go-05993.01 - 94a , and js acknowledges support from nsf and nasa grants .
here we quantify the empirical model uncertainties due to the finite size of our input sample .
we look at this uncertainty in terms of an arbitrary quantity describing the surface density of objects on the sky satisfying some observational criteria .    for any faint sample of objects satisfying a certain set of observable criteria ( @xmath147 ) ,
it is possible to express the probability distribution for the number of these objects as the sum of the probability distributions for this variable over all possible galaxy morphologies , luminosities , and sizes , i.e. , @xmath148 .
we suppose that this faint sample of galaxies is derived in a monte - carlo manner from our bright sample on the basis of no - evolution assumptions , so that it is possible to express this probability distribution for each galaxy type as the product of the probability distribution for the number of times this object will appear in our bright sample , call it @xmath149 , and the probability distribution for the number of times such a galaxy would be expected to appear in the faint sample in question ( @xmath147 ) , call it @xmath150 . in this way ,
the overall probability distribution for the number of galaxies observed in a given faint sample ( @xmath147 ) can be expressed as @xmath151 we shall suppose that each galaxy in our bright sample is sufficiently unique to be represented by its own term in the above equation .
furthermore , since each of these galaxies appeared in our bright sample once and only once , we suppose that @xmath149 is poisson - distributed with a mean of 1 , which we approximate as a normal distribution about this mean . for lack of something better
, we shall take the probability distributions @xmath149 for all other galaxy types to be delta functions at 0 .
, i.e. , with relatively larger contributions at faint magnitudes than at bright magnitudes , could result in large systematic errors in this distribution . ]    now , it is possible to estimate the uncertainty in the expectation value for the number of galaxies in our faint sample based on the finite number of galaxies we have in our bright sample . because our simulated area is , in principle , unlimited ,
the uncertainty in the expectation value of @xmath150 is vanishingly small , and in practice we take this uncertainty to be zero . as such , the above probability distribution is simply the sum of normal distributions with various weights , and therefore its 1-sigma uncertainty is simply equal to @xmath152 . we have used this expression throughout our paper in the estimation of errors for each observable calculated from our no - evolution simulations .
at fainter magnitudes , one expects a smaller number of the galaxy prototypes to make up an increasing fraction of our faint galaxy sample , as a result both of the differential k - corrections and of the relatively greater volume available to lower luminosity objects at fainter magnitudes . to quantify the importance of this trend , we compute a weighted measure of the effective number of galaxies which contribute to the counts in each magnitude bin , a quantity we take to equal @xmath153 where @xmath154 is the number of galaxies of a given prototype per magnitude per square degree recovered from our simulations and @xmath68 is the number of prototypes considered ( i.e. 31 ) .
basically , the above expression equals the square of the number counts divided by the estimated poissonian error based on the size of the contributing sample at the magnitude in question , the derivation of which is given in appendix a. in the case that all the @xmath154 s are the same , the expression reduces to @xmath68 as required .
we compute the above expression from our simulations , and we plot the results in figure 23 . at bright magnitudes ,
the bright galaxy catalogues are still dominated by shot noise so that the number of galaxies estimated from the above expression to contribute at any magnitude is lower than the real number ( @xmath155 ) which contribute . at fainter magnitudes ,
the effective of number of galaxies decreases again for the reasons stated above , i.e. , a relatively smaller differential volume for the more luminous galaxies and the k - correction . for the sake of clarity , we emphasize that at each magnitude , every galaxy prototype contributes some fraction to the number counts there .
abraham , r.g . ,
tanvir , n.r . ,
santiago , b.x . ,
ellis , r.s . ,
glazebrook , k. , van den bergh , s. , , 1996a , , 279 , l47 .
abraham , r.g .
, van den bergh , s. , glazebrook , k. , ellis , r.s . , santiago , b. , surma , p. , & griffiths , r.e . , 1996b , , 107 , 1 .
bertin , e. , & arnouts , s. 1996 , , 117 , 393 .
blecha , a. , golay , m. , huguenin , d. , reichen , d. , bersier , d. 1990 , , 233 , l9 - 12 .
bohlin , r.c .
, et al .
1991 , , 368 , 12 .
bouwens , r.j . , broadhurst , t.j .
, & silk , j. 1998 , , submitted .
brinchmann , j. , et al .
1998 , , submitted .
broadhurst , t.j .
, ellis , r.s , & shank , t. 1988 , , 235 , 827 .
casertano , s. , ratnatunga , k.u . , griffiths , r.e .
, i m . , m. , neuschaeffer , l.w . ,
ostrander , e.j . , windhorst , r.a .
1995 , , 453 , 699 .
clements , d. , & couch , w. 1997 , , 280 , l43 - 48 .
cohen , j.g . ,
cowie , l.l . ,
hogg , d.w .
, songaila , a. , blandford , r. , hu , e.m . , shopbell , p. 1996
, , 471 , l5 .
colless , m. , ellis , r.s . ,
taylor , k. , broadhurst , t.j . , peterson , b.a .
1993 , , 261 , 19 .
colley , w.n . ,
rhoads , j.e . ,
ostriker , j.p . , spergel , d.n .
1996 , , 473 , l63 .
connolly , a.j . ,
szalay , a.s .
, bershady , m.a . ,
kinney , a.l . , & calzetti , d. 1996 , , 110 , 1071 .
coleman , g.d . ,
wu , c .- c . , & weedman , d.w .
1980 , , 43 , 393 .
cowie , l.l . ,
songaila , a. , hu , e.m . , cohen , j.d .
1996 , , 112 , 839 .
driver , s.p . ,
windhorst , r.a . , & griffiths , r.e .
1995 , , 453 , 48 .
ellis , r.s . , colless , m. , broadhurst , t.j . , heyl , j.s . ,
glazebrook , k. 1996 , , 280 : 235 - 251 .
ford , h.c . , et al .
1996 . , 2807 , 184 .
frei , z. , guhathakurta , p. , gunn , j.e .
1996 , , 111 , 174 .
gallego , j. , zamorano , j. , aragon - salamanca , a. , rego , m. 1995 , , 445 , l1 .
giavalisco , m. , livio , m. , bohlin , r.c . ,
macchetto , f.d .
, stecher , t.p . 1996a , , 112 , 369 .
giavalisco , m. , steidel , c.c . ,
macchetto , f.d .
1996b , , 470 , 189 .
glazebrook , k. , ellis , r.s . ,
santiago , b. , griffiths , r. 1995b , , 275 , l19 .
guzman , r. , gallego , j. , koo , d.c . , phillips , a.c . , lowenthal , j.d . ,
faber , s.m . , illingworth , g.d . ,
vogt , n.p .
1997 , , 489 , 559 .
i m , m. , casertano , s. , griffiths , r.e .
, & ratnatunga , k.u .
1995 , , 441 , 494 .
kinney , a.l . ,
calzetti , d. , bohlin , r.c . , mcquade , k. , storchi - bergmann , t. , & schmitt , h.r .
1996 , , 467 , 38 .
lanzetta , k.m . ,
yahil , a. , fernandez - soto , a. 1996 , , 381 , 759 .
leitherer et al .
1996 , , 108 , 996 .
lilly , s.j .
1993 , , 411 , 502 .
lilly , s.j .
, tresse , l. , hammer , f. , crampton , d. , lefevre , o. 1995 , , 455 , 108 .
loveday , j. , peterson , b.a . , efstathiou , g. , & maddox , s.j .
1992 , , 390 , 338 .
lowenthal , j. , et al .
1997 , , 481 , 673 .
madau , p. 1995
, , 441 , 18 .
madau , p. , ferguson , h.c . ,
dickinson , m.e .
, giavalisco , m. , steidel , c.c . ,
& fruchter , a. 1996 , , 283 , 1388 .
maoz , d. 1997 , , 490 , l135 .
maoz , d. , filippenko , a.v .
, ho , l.c . , macchetto , f.d . ,
rix , h .- w . ,
schneider , d.p .
1997 , , 107 , 215 .
mutz , s.b .
, et al .
1994 , , 434 , l55 .
naim , a. , et al .
1995 , , 274 , 1107 .
oconnell , r.w . , & marcum , p. 1996 , in _ hst and the high redshift universe _
( 37th hertstmonceux conference ) , eds .
tanvir , a. aragon - salamanca , and j.v . wall , 1996 , astro - ph/9609101 .
oke , j.b .
1974 , , 27 , 21 .
peebles , p.j.e .
1993 , principles of physical cosmology , princeton university press , princeton .
petrosian , v. 1976 , , 209 , l1 .
pozzetti , l. , bruzual , g. , zamorani , g. 1996 , , 274 , 832 .
roche , n. , ratnatunga , k. , griffiths , r.e .
, i m , m. , & neuschaefer , l. 1996 , , 282 , 1247 .
schmidt , m. 1968 , , 151 , 393 .
stecher , t. , et al .
1992 , , 395 , l1 .
steidel , c.c . , dickinson , m. , persson , s.e .
1994 , , 437 , l75 .
steidel , c.c .
, giavalisco , m. , pettini , m. , dickinson , m. , adelburger , k.l .
1996a , , 462 , l17 .
steidel , c.c . ,
giavalisco , m. , dickinson , m. , adelburger , k.i .
1996b , , 112 , 352 .
williams , r.e . , et al .
1996 , , 112 , 1335 .
zucca , e. , et al .
1997 , , 326 , 477      12:36:49.351 & 62:13:47.934 & 0.089 & 18.2 & -19.2 & 4.4 & 19.7 & e + 12:36:50.971 & 62:13:21.738 & 0.199 & 19.6 & -20.2 & 1.0 & 21.9 & sp + 12:36:56.572 & 62:12:46.478 & 0.517 & 20.1 & -21.8 & 4.1 & 21.6 & sp + 12:36:47.992 & 62:13:10.107 & 0.475 & 20.5 & -21.2 & 4.3 & 19.9 & e + 12:36:50.138 & 62:12:40.811 & 0.474 & 20.7 & -21.4 & 1.1 & 20.7 & mrg + 12:36:42.827 & 62:12:17.400 & 0.454 & 20.7 & -21.1 & 1.2 & 20.5 & sp + 12:36:41.654 & 62:11:32.961 & 0.089 & 20.8 & -18.4 & 0.6 & 21.9 & sp + 12:36:43.713 & 62:11:43.907 & 0.764 & 20.9 & -22.5 & 3.8 & 19.2 & e + 12:36:53.820 & 62:12:55.070 & 0.642 & 20.9 & -21.8 & 2.0 & 20.4 & sp + 12:36:41.850 & 62:12:06.478 & 0.432 & 21.0 & -20.7 & 1.5 & 20.6 & sp + 12:36:46.952 & 62:12:37.901 & 0.318 & 21.0 & -19.9 & 0.9 & 20.7 & e + 12:36:51.708 & 62:13:54.832 & 0.557 & 21.1 & -21.2 & 2.1 & 20.5 & ir + 12:36:46.082 & 62:11:43.119 & 1.012 & 21.3 & -23.1 & 1.5 & 20.3 & sp + 12:36:58.690 & 62:12:53.457 & 0.319 & 21.3 & -19.7 & 1.0 & 22.6 & sp + 12:36:46.262 & 62:14:05.706 & 0.960 & 21.3 & -23.0 & 2.9 & 18.5 & e + 12:36:57.230 & 62:13:00.701 & 0.473 & 21.4 & -20.8 & 0.5 & 21.7 & sp + 12:37:00.485 & 62:12:35.720 & 0.562 & 21.4 & -20.8 & 3.3 & 20.1 & e + 12:36:50.193 & 62:12:46.807 & 0.677 & 21.5 & -21.5 & 4.4 & 20.1 & e + 12:36:49.634 & 62:13:14.095 & 0.475 & 21.5 & -20.4 & 2.3 & 21.7 & sp + 12:36:44.287 & 62:11:34.292 & 1.013 & 21.6 & -23.5 & 4.7 & 18.9 & e 12:36:51.640 & 62:12:21.247 & 0.299 & 21.6 & -19.0 & 2.5 & 21.1 & sp + 12:36:44.097 & 62:12:48.868 & 0.555 & 21.7 & -20.8 & 0.4 & 20.6 & sp + 12:36:43.066 & 62:12:43.258 & 0.847 & 21.8 & -22.2 & 4.8 & 18.9 & e + 12:36:49.416 & 62:14:07.778 & 0.752 & 21.8 & -21.4 & 1.7 & 19.5 & e + 12:36:55.493 & 62:12:46.551 & 0.790 & 22.0 & -21.5 & 1.6 & 20.4 & sp + 12:36:46.422 & 62:11:52.360 & 0.504 & 22.0 & -19.9 & 3.5 & 19.6 & e + 12:36:49.559 & 62:12:58.629 & 0.475 & 22.0 & -19.9 & 1.7 & 20.7 & sp + 12:36:49.288 & 62:13:12.300 & 0.478 & 22.2 & -19.9 & 0.9 & 19.3 & e + 12:36:43.936 & 62:12:50.496 & 0.556 & 22.2 & -21.3 & 2.0 & 20.5 & irr + 12:36:38.882 & 62:12:20.817 & 0.608 & 22.2 & -20.4 & 1.0 & 21.2 & sp + 12:36:39.920 & 62:12:08.424 & 1.015 & 22.3 & -22.5 & 4.5 & 18.3 & e +    ccc observations ( madau et al .
1996 ) & 58 & 14 + observations ( this work ) & 90 & 19 + ne ( @xmath70/@xmath71/cww ) & @xmath156 & @xmath157 + ne ( @xmath70/cww ) & @xmath158 & @xmath159 + ne ( @xmath74/cww ) & @xmath160 & @xmath161 + ne ( @xmath70/@xmath71/k96 ) & @xmath162 & @xmath157 + ne ( @xmath70/k96 ) & @xmath158 & @xmath159 + ne ( @xmath74/k96 ) & @xmath157 & @xmath162 +