Modern Multimessenger Astronomy is a powerful instrument for performing cosmological crucial tests of the Standard Cosmological Model in the wide redshift interval up to $z \sim 10$. This is principally important for discussion related to discrepancies between local and global measurements of cosmological parameters. We present a review of multimessenger gamma-ray burst observations currently conducted and planed for THESEUS$-$BTA cooperative program. Such observations give a unique opportunity to test the fundamental foundations of cosmological models: gravitation theory; cosmological principle of homogeneity and isotropy of large-scale distribution of matter; and space expansion paradigm. Important role of various selection effects leading to systematic distortions of true cosmological relations is discussed.

We examine the prospects of the high redshift Long Gamma Ray Bursts (LGRB) Hubble Diagram as a test of the basic cosmological principles. Analysis of the Hubble Diagram allows us to test several fundamental cosmological principles using the directly observed flux-distance-redshift relation. Modern LGRB data, together with the correlation between the spectral peak energy and the isotropic-equivalent radiated energy (the so-called Amati relation) can be used for construction of the Hubble Diagram at the model-independent level. We emphasise observational selection effects, which inevitably exist and distort the theoretically predicted relations. An example is the weak and strong gravitational lensing bias effect for high redshift LGRB in the presence of limited observational sensitivity (Malmquist bias). After bias correction, there is a tendency to vacuum dominated models with $\Omega_\Lambda \rightarrow 0.9$, $\Omega_\text{m} \rightarrow 0.1$. Forthcoming gamma-ray observations by the space THESEUS mission together with ground and space based multimessenger facilities will allow us to improve essentially the restrictions on alternative basic principles of cosmological models.

A review of modern VLTI observations of the orbital motion of closest stars to the relativistic compact object Sgr A* and its ability to test gravitation theories in the conditions of Post-Newtonian approximation. The observed orbital parameters, second order Doppler effect and gravitational redshift, measured for several S-stars, are compared with theoretical PN predictions.

The spatial distribution of sources of gamma-ray bursts (GRB) with known red shifts is analyzed by the conditional density and pairwise distance methods. The sample of GRB is based on data from the Swift mission and contains fluxes, coordinates, and red shifts for 384 GRB sources. Selection effects that distort the true source distribution are taken into account by comparing the observed distribution with fractal and uniform model catalogs. The Malmqvist effect is modeled using an approximation for the visible luminosity function of the GRB. The case of absorption in the galactic plane is also examined.This approach makes it possible to study the spatial structure of the entire sample at one time without artificial truncations. The estimated fractal dimensionality is $D=2.55\pm0.06$ on scales of $2\div6$ Gpc.

We have had the chance to live through a fascinating revolution in measuring the fundamental empirical cosmological Hubble law. The key progress is analysed : 1) improvement of observational means (ground-based radio and optical observations, space missions) ; 2) understanding of the biases that affect both distant and local determinations of the Hubble constant; 3) new theoretical and observational results. These circumstances encourage us to take a critical look at some facts and ideas related to the cosmological red-shift. This is important because we are probably on the eve of a new understanding of our Universe, heralded by the need to interpret some cosmological key observations in terms of unknown processes and substances.

The Local Universe is the most detail studied part of the observable region of space with the radius R about 100 Mpc. There are two empirical fundamental cosmological laws directly established from observations in the Local Universe independently from cosmological theory: first, the Hubble-Humason-Sandage linear redshift-distance law and second, Carpenter- Karachentsev-deVaucouleurs density-radius power-law. Review of modern state of these empirical laws and their cosmological significance is given. Possible theoretical interpretations of the surprising coexistence of both laws at the spatial scales from 1 Mpc to 100 Mpc are discussed. Comparison of the standard space-expansion explanation of the cosmological redshift with possible global gravitational redshift model is given

We report the results of a statistical analysis of the space distribution of galaxies within distances about 300 Mpc using the 2MRS catalog, which contains redshifts of 43533 galaxies of the 2MASS all-sky IR survey. Because of the unique features of the 2MRS survey, such as its 90 percent sky coverage, galaxy selection in the IR, the complete incorporation of the old stellar population of galaxies, weakness of the dust extinction effects, and the smallness of the k- and e-corrections allowed us to determine the statistical properties of the global distribution of galaxies in the Local Universe. We took into account the main methodological factors that distort the theoretically expected relations compared to those actually observed. We construct the radial galaxy number counts N(R), SL(R, r) statistics, and the complete correlation function (conditional density) for volume-limited (VL) galaxy samples. The observed conditional density in the redshift space is independent of the luminosity of galaxies and has the form of a power-law function with slope ~ 1.0 over a large range scale-length spanning from 0.1 to 100 Mpc. We compare the statistical properties of the space distribution of galaxies of the 2MRS catalog with the corresponding properties of simulated catalogs: stochastic fractal distributions and galaxies of the Millennium catalog.

New arguments supporting the reality of large-scale fluctuations in the density of the visible matter in deep galaxy surveys are presented. A statistical analysis of the radial distributions of galaxies in the COSMOS and HDF-N deep fields is presented. Independent spectral and photometric surveys exist for each field, carried out in different wavelength ranges and using different observing methods. Catalogs of photometric redshifts in the optical (COSMOS-Zphot) and infrared (UltraVISTA) were used for the COSMOS field in the redshift interval $0.1 < z < 3.5$, as well as the zCOSMOS (10kZ) spectroscopic survey and the XMM-COSMOS and ALHAMBRA-F4 photometric redshift surveys. The HDFN-Zphot and ALHAMBRA-F5 catalogs of photometric redshifts were used for the HDF-N field. The Pearson correlation coefficient for the fluctuations in the numbers of galaxies obtained for independent surveys of the same deep field reaches $R = 0.70 \pm 0.16$. The presence of this positive correlation supports the reality of fluctuations in the density of visible matter with sizes of up to 1 000 Mpc and amplitudes of up to 20% at redshifts $z \sim 2$. The absence of correlations between the fluctuations in different fields (the correlation coefficient between COSMOS and HDF-N is $R = -0.20 \pm 0.31$) testifies to the independence of structures visible in different directions on the celestial sphere. This also indicates an absence of any influence from universal systematic errors (such as "spectral voids"), which could imitate the detection of correlated structures.

We discuss how the redshift (Mattig) method in Friedmann cosmology relates to dynamical distance indicators based on Newton's gravity (Teerikorpi 2011). It belongs to the class of indicators where the relevant length inside the system is the distance itself (in this case the proper metric distance). As the Friedmann model has Newtonian analogy, its use to infer distances has instructive similarities to classical dynamical distance indicators. In view of the theoretical exact linear distance-velocity law, we emphasize that it is conceptually correct to derive the cosmological distance via the route: redshift (primarily observed) --> space expansion velocity (not directly observed) --> metric distance (physical length in "cm"). Important properties of the proper metric distance are summarized.

The idea of the global gravitational effect as the source of cosmological redshift was considered by de Sitter (1916, 1917), Eddington (1923), Tolman (1929) and Bondi (1947). Also Hubble (1929) called the discovered distance-redshift relation as "De Sitter effect". For homogeneous matter distribution cosmological gravitational redshift is proportional to square of distance: z_grav ~ r^2. However for a fractal matter distribution having the fractal dimension D=2 the global gravitational redshift is the linear function of distance: z_grav ~ r, which gives possibility for interpretation of the Hubble law without the space expansion. Here the field gravity fractal cosmological model (FGF) is presented, which based on two initial principles. The first assumption is that the Feynman's field gravity approach describes the gravitational interaction, which delivers a natural basis for the conceptual unity of all fundamental physical interactions within the framework of the relativistic and quantum fields in Minkowski space. The second hypothesis is that the spatial distribution of gravitating matter is a fractal at all scales up to the Hubble radius. The fractal dimension of matter distribution is assumed to be D = 2, which implies that the global gravitational redshift is the explanation of the observed linear Hubble law. In the frame of the FGF all three phenomena - the cosmic background radiation, the fractal large scale structure, and the Hubble law, - could be the consequence of a unique large scale structure evolution process of the initially homogeneous ordinary matter without nonbaryonic matter and dark energy.

One of the main problems of observational cosmology is to determine the range in which a reliable measurement of galaxy correlations is possible. This corresponds to determine the shape of the correlation function, its possible evolution with redshift and the size and amplitude of large scale structures. Different selection effects, inevitably entering in any observation, introduce important constraints in the measurement of correlations. In the context of galaxy redshift surveys selection effects can be caused by observational techniques and strategies and by implicit assumptions used in the data analysis. Generally all these effects are taken into account by using pair-counting algorithms to measure two-point correlations. We review these methods stressing that they are based on the a-priori assumption that galaxy distribution is spatially homogeneous inside a given sample. We show that, when this assumption is not satisfied by the data, results of the correlation analysis are affected by finite size effects.In order to quantify these effects, we introduce a new method based on the computation of the gradient of galaxy counts along tiny cylinders. We show, by using artificial homogeneous and inhomogeneous point distributions, that this method is to identify redshift dependent selection effects and to disentangle them from the presence of large scale density fluctuations. We then apply this new method to several redshift catalogs and we find evidences that galaxy distribution, in those samples where selection effects are small enough, is characterized by power-law correlations with exponent $\gamma=0.9$ up to $20$ Mpc/h followed by a change of slope that, in the range [20,100] Mpc/h, corresponds to a power-law exponent $\gamma=0.25$. Whether a crossover to spatial unformity occurs at $\sim 100$ Mpc/h cannot be clarified by the present data.

We have analyzed the spatial distribution of galaxies from the release of the Sloan Digital Sky Survey of galactic redshifts (SDSS DR7), applying the complete correlation function (conditional density), two-point conditional density (cylinder), and radial density methods. Our analysis demonstrates that the conditional density has a power-law form for scales lengths 0.5-30 Mpc/h, with the power-law corresponding to the fractal dimension D = 2.2+-0.2; for scale lengths in excess of 30 Mpc/h, it enters an essentially flat regime, as is expected for a uniform distribution of galaxies. However, in the analysis applying the cylinder method, the power-law character with D = 2.0+-0.3 persists to scale lengths of 70 Mpc/h. The radial density method reveals inhomogeneities in the spatial distribution of galaxies on scales of 200 Mpc/h with a density contrast of two, confirming that translation invariance is violated in the distribution of galaxies to 300 Mpc/h, with the sampling depth of the SDSS galaxies being 600 Mpc/h.

We present an analysis of BVRcIc observations of the field sized around 4' x 4' centered at the host galaxy of the gamma-ray burst GRB 021004 with the 6-m BTA telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences. We measured the magnitudes and constructed the color diagrams for 311 galaxies detected in the field (S/N > 3). The differential and integral counts of galaxies up to the limit, corresponding to 28.5 (B), 28.0 (V), 27.0 (Rc), 26.5 (Ic) were computed. We compiled the galaxy catalog, consisting of 183 objects, for which the photometric redshifts up to the limiting magnitudes 26.0 (B), 25.5 (V), 25.0 (Rc), 24.5 (Ic) were determined using the HyperZ code. We then examined the radial distribution of galaxies based on the z estimates. We have built the curves expected in the case of a uniform distribution of galaxies in space, and obtained the estimates for the size and contrast of the possible super-large-scale structures, which are accessible with the observations of this type.

Cosmological density fields are assumed to be translational and rotational invariant, avoiding any special point or direction, thus satisfying the Copernican Principle. A spatially inhomogeneous matter distribution can be compatible with the Copernican Principle but not with the stronger version of it, the Cosmological Principle which requires the additional hypothesis of spatial homogeneity. We establish criteria for testing that a given density field, in a finite sample at low redshifts, is statistically and/or spatially homogeneous. The basic question to be considered is whether a distribution is, at different spatial scales, self-averaging. This can be achieved by studying the probability density function of conditional fluctuations. We find that galaxy structures in the SDSS samples, the largest currently available, are spatially inhomogeneous but statistically homogeneous and isotropic up to ~ 100 Mpc/h. Evidences for the breaking of self-averaging are found up to the largest scales probed by the SDSS data. The comparison between the results obtained in volumes of different size allows us to unambiguously conclude that the lack of elf-averaging is induced by finite-size effects due to long-range correlated fluctuations. We finally discuss the relevance of these results from the point of view of cosmological modeling.

Initial principles of a method of analysis of the luminous matter spatial distribution with sizes about thousands Mpc are presented. The method is based on an analysis of the photometric redshift distribution N(z) in the deep fields with large redshift bins \Deltaz=0.1{\div}0.3. Number density fluctuations in the bins are conditioned by the Poisson's noise, the correlated structures and the systematic errors of the photo-z determination. The method includes covering of a sufficiently large region on the sky by a net of the deep multiband surveys with the sell size about 10^{\circ}x10^{\circ} where individual deep fields have angular size about 10'x10' and may be observed at telescopes having diameters 3-10 meters. The distributions of photo-z within each deep field will give information about the radial extension of the super large structures while a comparison of the individual radial distributions of the net of the deep fields will give information on the tangential extension of the super large structures. A necessary element of the method is an analysis of possible distortion effects related to the methodic of the photo-z determination.

Using the deep fields of COSMOS, FDF, HUDF, and HDF-N as an example, we discuss the prospects for and limitations on the method for searching for super large structures in the spatial distribution of galaxies proposed in the preceding article of this series. An analysis of the distribution N(z) of photometric redshifts in a grid of the deep fields of HUDF-FDF-COSMOS-HDFN reveals the possible existence of super large structures with a contrast dN/N~50% and tangential and radial dimensions of about 1000 Mpc. The reality of the detected candidate super large structures in the universe can be verified by further observations with a finer grid of deep fields. The influence of systematic errors can be reduced by observing the same deep fields with several 3-10 meter telescopes and utilizing different methods for determining the photometric redshifts.

Statistical analyses of finite sample distributions usually assume that fluctuations are self-averaging, i.e. that they are statistically similar in different regions of the given sample volume. By using the scale-length method, we test whether this assumption is satisfied in several samples of the Sloan Digital Sky Survey Data Release Six. We find that the probability density function (PDF) of conditional fluctuations, filtered on large enough spatial scales (i.e., r>30 Mpc/h), shows relevant systematic variations in different sub-volumes of the survey. Instead for scales r<30 Mpc/h the PDF is statistically stable, and its first moment presents scaling behavior with a negative exponent around one. Thus while up to 30 Mpc/h galaxy structures have well-defined power-law correlations, on larger scales it is not possible to consider whole sample average quantities as meaningful and useful statistical descriptors. This situation is due to the fact that galaxy structures correspond to density fluctuations which are too large in amplitude and too extended in space to be self-averaging on such large scales inside the sample volumes: galaxy distribution is inhomogeneous up to the largest scales, i.e. r ~ 100 Mpc/h, probed by the SDSS samples. We show that cosmological corrections, as K-corrections and standard evolutionary corrections, do not qualitatively change the relevant behaviors. Finally we show that the large amplitude galaxy fluctuations observed in the SDSS samples are at odds with the predictions of the standard LCDM model of structure formation.(Abridged version).

We consider the conditional galaxy density around each galaxy, and study its fluctuations in the newest samples of the Sloan Digital Sky Survey Data Release 7. Over a large range of scales, both the average conditional density and its variance show a nontrivial scaling behavior, which resembles to criticality. The density depends, for 10 < r < 80 Mpc/h, only weakly (logarithmically) on the system size. Correspondingly, we find that the density fluctuations follow the Gumbel distribution of extreme value statistics. This distribution is clearly distinguishable from a Gaussian distribution, which would arise for a homogeneous spatial galaxy configuration. We also point out similarities between the galaxy distribution and critical systems of statistical physics.

One of the most striking features predicted by standard models of galaxy formation is the presence of anti-correlations in the matter distribution at large enough scales (r>r_c). Simple arguments show that the location of the length-scale r_c, marking the transition from positive to negative correlations, is the same for any class of objects as for the full matter distribution, i.e. it is invariant under biasing. This scale is predicted by models to be at about the same distance of the scale signaling the baryonic acoustic oscillation scale r_{bao}. We test these predictions in the newest SDSS galaxy samples.We find that, in several MG samples, the correlation function remains positive at scales >250 Mpc/h, while in the concordance LCDM it should be negative beyond r_c\approx 120 Mpc/h. In other samples the correlation function becomes negative at scales <50 Mpc/h. To investigate the origin of these differences we consider in detail the propagation of errors on the sample density into the estimation of the correlation function. We conclude that these are important at large enough separations, and that they are responsible for the observed differences between different estimators and for the measured sample to sample variations of the correlation function. We conclude that, in the newest SDSS samples, the large scale behavior of the galaxy correlation function is affected by intrinsic errors andv olume-dependent systematic effects which make the detection of correlations to be only an estimate of a lower limit of their amplitude, spatial extension and statistical errors. We point out that these results represent an important challenge to LCDM models as they largely differ from its predictions.(Abridged version).

The properties of the galaxy distribution at large scales are usually studied using statistics which are assumed to be self-averaging inside a given sample. We present a new analysis able to quantitatively map galaxy large scale structures while testing for the stability of average statistical quantities in different sample regions. We find that the newest samples of the Sloan Digital Sky Survey provide unambiguous evidence that galaxy structures correspond to large amplitude density fluctuations at all scales limited only by sample sizes. The two-point correlations properties are self-averaging up to approximately 30 Mpc/h and are characterized by a fractal dimension D=2.1 +- 0.1. Then at all larger scales probed density fluctuations are too large in amplitude and too extended in space to be self-averaging inside the considered volumes. These inhomogeneities are compatible with a continuation of fractal correlations but incompatible with: (i) a homogeneity scale smaller than 100 Mpc/h, (ii) predictions of standard theoretical models, (iii) mock galaxy catalogs generated from cosmological Nbody simuations.

We study statistical properties of galaxy structures in several samples extracted from the 2dF Galaxy Redshift Survey. In particular, we measured conditional fluctuations by means of the scale-length method and determined their probability distribution. In this way we find that galaxy distribution in these samples is characterized by large amplitude fluctuations with a large spatial extension, whose size is only limited by the sample's boundaries. These fluctuations are quite typical and persistent in the sample's volumes, and they are detected in two independent regions in the northern and southern galactic caps. We discuss the relation of the scale-length method to several statistical quantities, such as counts of galaxies as a function of redshift and apparent magnitude. We confirm previous results, which have determined by magnitude and redshift counts that there are fluctuations of about 30% between the southern and the northern galactic caps and we relate explicitly these counts to structures in redshift space. We show that the estimation of fluctuation amplitude normalized to the sample density is biased by systematic effects, which we discuss in detail. We consider the type of fluctuations predicted by standard cosmological models of structure formation in the linear regime and, to study nonlinear clustering, we analyze several samples of mock-galaxy catalogs generated from the distribution of dark matter in cosmological N-body simulations. In this way we conclude that the galaxy fluctuations present in these samples are too large in amplitude and too extended in space to be compatible with the predictions of the standard models of structure formation.

Images of the Hubble Ultra Deep Field are analyzed to obtain a catalog of galaxies for which the angular sizes, surface brightness, photometric redshifts, and absolute magnitudes are found. The catalog contains a total of about 4000 galaxies identified at a high signal-to-noise ratio, which allows the cosmological relations angular size{redshift and surface brightness-redshift to be analyzed. The parameters of the evolution of linear sizes and surface brightness of distant galaxies in the redshift interval 0.5-6.5 are estimated in terms of a grid of cosmological models with different density parameters. The distribution of photometric redshifts of galaxies is analyzed and possible superlarge inhomogeneities in the radial distribution of galaxies are found with scale lengths as large as 2000 Mpc.

We apply the scale-length method to several three dimensional samples of the Two degree Field Galaxy Redshift Survey. This method allows us to map in a quantitative and powerful way large scale structures in the distribution of galaxies controlling systematic effects. By determining the probability density function of conditional fluctuations we show that large scale structures are quite typical and correspond to large fluctuations in the galaxy density field. We do not find a convergence to homogeneity up to the samples sizes, i.e. ~ 75 Mpc/h. We then measure, at scales r <~ 40 Mpc/h, a well defined and statistically stable power-law behavior of the average number of galaxies in spheres, with fractal dimension D=2.2 +- 0.2. We point out that standard models of structure formation are unable to explain the existence of the large fluctuations in the galaxy density field detected in these samples. This conclusion is reached in two ways: by considering the scale, determined by the linear perturbation analysis of a self-gravitating fluid, below which large fluctuations are expected in standard models and through the determination of statistical properties of mock galaxy catalogs generated from cosmological N-body simulations of the Millenium consortitum.

The space expansion physics contains several paradoxes which were clearly demonstrated by Edward Harrison (1981, 1995, 2000), who emphasized that the cooling of homogeneous hot gas (including photon gas of CBR) in the standard cosmological model based on the violation of energy conservation by the expanding space. In modern version of SCM the term "space expansion" actually means continuous creation of vacuum, something that leads to conceptual problems. Recent discussion by Francis, Barnes, James, and Lewis (2007) on the physical sense of the increasing distance to a receding galaxy without motion of the galaxy is just a particular consequence of the arising paradoxes. Here we present an analysis of the following conceptual problems of the SCM: the violation of energy conservation for local comoving volumes, the exact Newtonian form of the Friedmann equation, the absence of an upper limit on the receding velocity of galaxies which can be greater than the speed of light, and the presence of the linear Hubble law deeply inside inhomogeneous galaxy distribution. The common cause of these paradoxes is the geometrical description of gravity, where there is no a well defined concept of the energy-momentum tensor for the gravitational field, no energy quanta - gravitons, and no energy-momentum conservation for matter plus gravity because gravity is not a material field.

Gravitational lensing of background compact objects like active galactic nuclei and quasars, by extended intermediate mass lenses such as globular clusters and and dark matter clumps with masses 10^5 - 10^8 M_sun, is considered. It is shown that observational study of the galaxy-quasar's associations is a powerful direct observational test of the nature of massive galaxy haloes. Optical interferometric observations with VLTI and Keck instruments are able to constrain masses and number of substructure halo objects. Evidence of gravitational lensing by globular clusters in haloes of spiral and elliptical galaxies is presented.

Recent extensive, multi-color deep surveys of galaxies open a possibility to get observational estimation of sizes for the largest structures in the Universe. Photometric redshift accuracy (about 0.03(1+z)) allows directly study clustering at scales about 1000 Mpc. Thanks to large number of galaxies in each redshift bin one may detect super-large structures if they really exist. Here we show that the observed behavior of the redshift distribution of galaxies in deep surveys such as HUDF and FDF is consistent with existence of super-large structures of luminous matter with scales about 2000 Mpc. We detect a large underdense region in radial galaxy distribution at redshift interval z=1.2 - 2.2 which separate our "Local Hubble Volume" from the neighboring over-density region at z=2.2 - 3.5. This result can also explain the observed deficiency of gamma ray sources at redshift about 2. Observational test on the reality of the supper-large structures may be obtained by organizing sky covering net (cells about 10n x 10n degrees) of very deep narrow angle (1n x 1n arc-minutes) multi-band photometric surveys of galaxies which is achievable for large ground-based telescopes.

A history of the problem of mathematical and physical definition for the energy-momentum of the gravity field is reviewed. As it was noted 90 years ago by Hilbert (1917), Einstein (1918), Schrodinger (1918) and Bauer (1918) within Geometrical Gravity approach (General Relativity) there is no tensor characteristics of the energy-momentum for the gravity field. Landau & Lifshitz (1971) called this quantity pseudo-tensor of energy-momentum and noted that Einstein's equations does not express the energy conservation for matter plus gravity field. This has crucial consequences for gravity physics and cosmology, such as negative energy density for static gravity field and violation of energy conservation in expanding space. However there is alternative Field Gravity approach for description of gravitation as a symmetric tensor field in Minkowski space, which is similar to description of all other physical interactions and based on well-defined positive, localizable energy-momentum of the gravity field. This relativistic quantum Field Gravity approach was partially developed by Firz & Pauli (1939), Birkhoff (1944), Thirring (1961), Kalman (1961), Feynman (1963) and others. Here it is shown that existence of well-defined positive energy-momentum of the gravity field leads to radical changes in gravity physics and cosmology, including such new possibilities as two-component nature of gravity - attraction (spin 2) and repulsion (spin 0), absence of black holes and singularities, scalar gravitational radiation caused by spherically symmetric gravitational collapse.

Crucial observational tests of gravity physics are reviewed. Such tests are able to clarify the key question on the nature of gravitational interaction: is gravity the curvature of space? or is gravity a matter field in Minkowski flat space as other physical forces? Up to now all actually performed experiments do not allow to distinguish between these two alternatives in gravity physics. The existence of well-defined positive energy-momentum of the gravity field in Poincare-Feynman approach leads to radical changes in gravity physics and cosmology which may be tested by laboratory experiments and astrophysical observations. New possibilities for observational distinction between geometrical general relativity and field gravity theories are discussed. Among them: the contribution of the scalar repulsive force into Newtonian gravitational interaction, post-Newtonian translational motion of rotating bodies, gravitational deflection of light by small mass bodies, scalar gravitational radiation from spherically pulsating stars, existence of limiting radius, surface, magnetic field for massive bodies and absence of singularities and horizons for relativistic compact objects.

We present a summary of the International conference "Problems of practical cosmology", held at Russian Geographical Society, 23-27 June 2008, St.-Petersburg, Russia, where original reports were offered for discussion of new developments in modern cosmological physics, including the large scale structure of the Universe, the evolution of galaxies, cosmological effects in the local stellar systems, gravity physics for cosmology, cosmological models, and crucial observational tests of rival world models. The term "Practical Cosmology" was introduced by Allan Sandage in 1995 when he formulated "23 astronomical problems for the next three decades" at the conference on "Key Problems in Astronomy and Astrophysics" held at Canary Islands. Now when the first decade has passed, we can summarise the present situation in cosmological physics emphasizing interesting hot problems that have arisen during the last decade. Full texts of all reports are available at the website of the conference.

We report the detection of dark energy near the Milky Way made with precision observations of the local Hubble flow of expansion. We estimate the local density of dark energy and find that it is near, if not exactly equal to, the global dark energy density. The result is independent of, compatible with, and complementary to the horizon-scale observations in which dark energy was first discovered. Together with the cosmological concordance data, our result forms direct observational evidence for the Einstein antigravity as a universal phenomenon -- in the same sense as the Newtonian universal gravity.

We discuss the estimation of galaxy correlation properties in several volume limited samples, in different sky regions, obtained from the Fourth Data Release of the Sloan Digital Sky Survey. The small scale properties are characterized through the determination of the nearest neighbor probability distribution. By using a very conservative statistical analysis, in the range of scales [0.5,~30] Mpc/h we detect power-law correlations in the conditional density in redshift space, with an exponent \gamma=1.0 \pm 0.1. This behavior is stable in all different samples we considered thus it does not depend on galaxy luminosity. In the range of scales [~30,~100] Mpc/h we find evidences for systematic unaveraged fluctuations and we discuss in detail the problems induced by finite volume effects on the determination of the conditional density. We conclude that in such range of scales there is an evidence for a smaller power-law index of the conditional density. However we cannot distinguish between two possibilities: (i) that a crossover to homogeneity (corresponding to \gamma=0 in the conditional density) occurs before 100 Mpc/h, (ii) that correlations extend to scales of order 100 Mpc/h (with a smaller exponent 0 < \gamma <1). We emphasize that galaxy distributions in these samples present large fluctuations at the largest scales probed, corresponding to the presence of large scale structures extending up to the boundaries of the present survey. Finally we discuss several differences between the behavior of the conditional density in mock galaxy catalogs built from cosmological N-body simulations and real data. We discuss some theoretical implications of such a fact considering also the super-homogeneous features of primordial density fields.

The typical spectra of gamma-ray bursts (GRBs) are discussed in the context of the compactness problem for GRB sources and how it is resolved in the popular fireball model. In particular, observational (model-independent) constraints on the collimation of the gamma-rays and the dependence of the collimation angle on the photon energy are considered. The fact that the threshold for the creation of $e^{-}e^{+}$ pairs depends on the angle between the momenta of the annihilating photons in the GRB source provides an alternative solution to the compactness problem. A new approach to explaining GRBs, taking into account the angular dependence for pair creation, is proposed, and the main features of a scenario describing a GRB source with a total (photon) energy smaller or of the order of $10^{49}$ erg are laid out. Thus, we are dealing with an alternative to an ultra-relativistic fireball, if it turns out (as follows from observations) that all "long" GRBs are associated with normal (not peculiar) core-collapse supernovae. The effects of radiation pressure and the formation of jets as a consequence of even a small amount of anisotropy in the total radiation field in a (compact) GRB source are examined in this alternative model. Possible energy release mechanisms acting in regions smaller or of the order of $10^{8}$ cm in size (a compact model for a GRB) are discussed. New observational evidence for such compact energy release in the burst source is considered.

Typical observational gamma-ray burst (GRB) spectra are discussed and, in this connection, what is the origin of the compactness problem and how it was solved at first. If the threshold for $e^{-}e^{+}$ pair production depends on an angle between photon momenta, then another solution of the compactness problem is possible. We discuss a possibility of the $\gamma$-rays collimation and the dependence of photon beaming on photon energies. The list of basic assumptions of the scenario describing the GRB source with energy $< 10^{49}$ ergs are adduced: the matter is about an alternative to the ultrarelativistic fireball if {\it all} long-duration GRBs are related or physically connected with {\it normal/unpeculiar} core-collapse supernovae (SNe). Namely, we consider the questions about radiation pressure and how the jet arises on account of even small asymmetry of the radiation field in a compact GRB source. The possibility of a new approach to explanation of GRB phenomenon is shown. Possible mechanisms of their generation in regions of size $< 10^8$cm are discussed (a compact model of GRBs). Observational consequences of the compact GRB energy release are considered.

Basing on the increasing evidence for the cosmological relevance of the local Hubble flow, we consider a simple analytical cosmological model for the Local Universe. This is a non-Friedmann model with a non-uniform static space-time. The major dynamical factor controlling the local expansion is the antigravity produced by the omnipresent and permanent dark energy of the cosmic vacuum (or the cosmological constant). The antigravity dominates at distances larger than 1-2 Mpc from the center of the Local Group. The model gives a natural explanation of the two key quantitative characteristics of the local expansion flow, which are the local Hubble constant and the velocity dispersion of the flow. The observed kinematical similarity of the local and global flows of expansion is clarified by the model. We demonstrate analytically the efficiency of the vacuum cooling mechanism that allows one to see the Hubble flow so close to the Local Group. Special significance is argued for the 'universal Hubble constant' H_V, depending only on the vacuum density (H_V \approx 60 km/s/Mpc). The model makes a number of testable predictions. It also shows unexpectedly that it is the dwarf galaxies of the local flow with the shortest distances and lowest redshifts that may be the most sensitive indicators of dark energy in our neighborhood.

We study galaxy correlations from samples extracted from the 2dFGRS final release. Statistical properties are characterized by studying the nearest neighbor probability density, the conditional density and the reduced two-point correlation function. The result is that the conditional density has a power-law behavior in redshift space described by an exponent \gamma=0.8 \pm 0.2 in the interval from about 1 Mpc/h, the average distance between nearest galaxies, up to about 40 Mpc/h, corresponding to radius of the largest sphere contained in the samples. These results are consistent with other studies of the conditional density and are useful to clarify the subtle role of finite-size effects on the determination of the two-point correlation function in redshift and real space

The physics of the expansion of the universe is still a poorly studied subject of the standard cosmological model. This because the concept of expanding space can not be tested in the laboratory and because ``expansion'' means continuous creation of space, something that leads to several paradoxes. We re-consider and expand here the discussion of conceptual problems, already noted in the literature, linked to the expansion of space. In particular we discuss the problem of the violation of energy conservation for local comoving volumes, the exact Newtonian form of the Friedmann equations, the receding velocity of galaxies being greater than the speed of light, and the Hubble law inside inhomogeneous galaxy distribution. Recent discussion by Kiang, Davis \& Lineweaver, and Whiting of the non-Doppler nature of the Lemaitre cosmological redshift in the standard model is just a particular consequence of the paradoxes mentioned above. The common cause of these paradoxes is the geometrical description of gravity (general relativity), where there is not a well defined concept of the energy-momentum tensor for the gravitational field and hence no energy-momentum conservation for matter plus gravity.

We review the increasing evidence for the cosmological relevance of the cold local Hubble flow. New observations, N-body simulations and other theoretical arguments are discussed, supporting our previous suggestion that the cosmological vacuum or uniform dark energy can have locally observable consequencies, especially a lower velocity scatter in DE dominated regions. The apparent contradiction between the slight dependence of the growth factor on $\Omega_{\Lambda}$ and the significant influence of dark energy in realistic N-body calculations is clarified. An interesting new result is that in the standard $\Lambda$ cosmology, gravitation dominates around a typical matter fluctuation up to about the correlation length $r_0$, and we tentatively link this with the high pairwise velocity dispersion on scales up to several Mpc, as measured from galaxy redshift-space correlations. Locally, the smooth Hubble flow on similar scales is consistent with N-body simulations including $\Omega_{\Lambda} \approx 0.7$ and a low density contrast in the Local Volume, which make it generally vacuum-dominated beyond 1 - 2 Mpc from galaxies and groups. We introduce a useful way to view the Hubble flow in terms of "zero gravity" spheres aound galaxies: e.g. a set of non-intersecting spheres, observed to be expanding, actually participates in accelerating expansion. The observed insensitiveness of the local velocity dispersion to galaxy mass is explained as an effect of the vacuum, too.

We present a review of the history and the present state of the fractal approach to the large-scale distribution of galaxies. Angular correlation function was used as a general instrument for the structure analysis. It was realized later that a normalization condition for the reduced correlation function estimator results in distorted values for both R_{hom} and fractal dimension D. Moreover, according to a theorem on projections of fractals, galaxy angular catalogues can not be used for detecting a structure with the fractal dimension D>2. For this 3-d maps are required, and indeed modern extensive redshift-based 3-d maps have revealed the ``hidden'' fractal dimension of about 2, and have confirmed superclustering at scales even up to 500 Mpc (e.g. the Sloan Great Wall). On scales, where the fractal analysis is possible in completely embedded spheres, a power--law density field has been found. The fractal dimension D =2.2 +- 0.2 was directly obtained from 3-d maps and R_{hom} has expanded from 10 Mpc to scales approaching 100 Mpc. In concordance with the 3-d map results, modern all sky galaxy counts in the interval 10^m - 15^m give a 0.44m-law which corresponds to D=2.2 within a radius of 100h^{-1}_{100} Mpc. We emphasize that the fractal mass--radius law of galaxy clustering has become a key phenomenon in observational cosmology.

We suggest a new method for estimating the fractal dimension of the spatial distribution of galaxies, the method of selected cylinders. We show the capabilities of this method by constructing a two-point conditional column density for galaxies with known redshifts from the LEDA database. The fractal dimension of a sample of LEDA and EDR SDSS galaxies has been estimated to be D = 2.1 for cylinder lengths of 200 Mpc. Amajor advantage of the suggested method is that it allows scales comparable to the catalog depth to be analyzed for galaxy surveys in the form of conical sectors and small fields in the sky.

Gravitational wave (GW) detectors operating on a long time range can be used for the study of space distribution of sources of GW bursts or to put strong upper limits on the GW signal of a wide class of source candidates. For this purpose we propose here a sidereal time analysis to analyze the output signal of GW detectors. Using the characteristics of some existing detectors, we demonstrate the capability of the sidereal time analysis to give a clear signature of different localizations of GW sources: the Galactic Center, the Galactic Plane, the Supergalactic plane, the Great Attractor. On the contrary, a homogeneous 3D-distribution of GW sources gives a signal without features. In this paper we consider tensor gravitational waves with randomly oriented polarization. We consider GW detectors at fixed positions on the Earth, and a fixed orientation of the antenna.

Recent observations support the view that the universe is described by a FLRW model with $\Omega_m^0 \approx 0.3$, $\Omega_{\Lambda}^0 \approx 0.7$, and $w \leq -1/3$ at the present epoch. There are several theoretical suggestions for the cosmological $\Lambda$ component and for the particular form of the energy transfer between this dark energy and matter. This gives a strong motive for a systematic study of general properties of two-fluid FLRW models. We consider a combination of one perfect fluid, which is quintessence with negative pressure ($p_Q = w\epsilon_Q$), and another perfect fluid, which is a mixture of radiation and/or matter components with positive pressure ($p = \beta \epsilon_m$), which define the associated one-fluid model ($p = \gamma \epsilon$). We introduce a useful classification which contains 4 classes of models defined by the presence or absence of energy transfer and by the stationarity ($w = const.$ and $\beta = const.$) or/and non stationarity ($w$ or $\beta$ time dependent) of the equations of state. It is shown that, for given $w$ and $\beta$, the energy transfer defines $\gamma$ and, therefore, the total gravitating mass and dynamics of the model. We study important examples of two-fluid FLRW models within the new classification. The behaviour of the energy content, gravitating mass, pressure, and the energy transfer are given as functions of the scale factor. We point out three characteristic scales, $a_E$, $a_{\cal P}$ and $a_{\cal M}$, which separate periods of time in which quintessence energy, pressure and gravitating mass dominate. Each sequence of the scales defines one of 6 evolution types.

The probability of strong gravitational lensing of distant compact objects by transparent lenses with King mass distributions is calculated. Examples of such lenses are globular clusters, dwarf galaxies, and clusters of hidden mass. The lens masses thus range from 10^3 to 10^9 Solar mass, intermediate between stellar and galactic masses. Taking account of the conic caustic of King lenses and the fractal large-scale distribution of matter along the line of sight, it is possible to explain the correlation between quasars and galaxies as the effect of strong gravitational lensing of the active nuclei of distant galaxies by globular clusters in the halos of nearby galaxies. Observational tests for verifing this mesolensing hypothesis are described.

We use data on the local 3-dimensional galaxy distribution for studying the statistics of the detection rates of gravitational waves (GW) coming from supernova explosions. We consider both tensor and scalar gravitational waves which are possible in a wide range of relativistic and quantum gravity theories. We show that statistics of GW events as a function of sidereal time can be used for distinction between scalar and tensor gravitational waves because of the anisotropy of spatial galaxy distribution. For calculation of the expected amplitudes of GW signals we use the values of the released GW energy, frequency and duration of GW pulse which are consistent with existing scenarios of SN core collapse. The amplitudes of the signals produced by Virgo and the Great Attractor clusters of galaxies is expressed as a function of the sidereal time for resonant bar detectors operating now (IGEC) and for forthcoming laser interferometric detectors (VIRGO).Then, we calculate the expected number of GW events as a function of sidereal time produced by all the galaxies within 100 Mpc. In the case of axisymmetric rotational core collapse which radiates a GW energy of $10^{-9}M_{\odot}c^2$, only the closest explosions can be detected. However, in the case of nonaxisymmetric supernova explosion, due to such phenomena as centrifugal hangup, bar and lump formation, the GW radiation could be as strong as that from a coalescing neutron-star binary. For radiated GW energy higher than $10^{-6}M_{\odot}c^2$ and sensitivity of detectors at the level $h \approx 10^{-23}$ it is possible to detect Virgo cluster and Great Attractor, and hence to use the statistics of GW events for testing gravity theories.

Cepheid distances of local galaxies (<7 Mpc) are used to study the very nearby velocity field, as pioneered by Sandage (Sandage 1986) who also pointed out its remarkable properties: linearity and quietness. The new data show that the velocity dispersion in the distance range as seen from the barycentre of the Local Group 1-8 Mpc is as low as 38 km/sec. The local rate of expansion coincides with the global Hubble constant. Down to 1.5 Mpc we cannot detect a deviation from the linear Hubble flow. This puts an upper limit for the mass of the Local Group, for a wide class of Friedman models, including those with the cosmological constant.

The cosmological vacuum, which is perfectly uniform, dominates by density over all the forms of cosmic matter. It makes the Universe be actually more uniform than it could be seen from the visible picture of the highly non-uniform matter distribution, especially inside the observed cell of uniformity (100-150 Mpc). This uniformity reveals itself in the structure of the Hubble matter flow which extends over a giant range of cosmic space scales -- from few Mpc to a thousand Mpc, -- preserving its kinematical identity. According to Sandage (1999), this flow is mysteriously regular and quiet even deep inside the cell of uniformity. An answer we propose to the question in the title above is as follows: This is most probably because the flow is dynamically controlled by the cosmological vacuum. An additional conjecture of cosmological intermittency, that addresses a complex statistical structure of initial chaotic perturbations, is also suggested in this context.

Our local environment at $r<10$ Mpc expands linearly and smoothly, as if ruled by a uniform matter distribution, while observations show the very clumpy local universe. This is a long standing enigma in cosmology. We argue that the recently discovered vacuum or quintessence (dark energy (DE) component with the equation of state $p_Q = w \rho_Q c^2$, $w \in [-1,0)$) from observations of the high-redshift universe may also manifest itself in the properties of the very local Hubble flow. We introduce the concept of the critical distance $r_Q$ where the repulsive force of dark energy starts to dominate over the gravity of a mass concentration. For the Local Group $r_Q$ is about 1.5 Mpc. Intriguingly, at the same distance 1.5 Mpc the linear and very "cold" Hubble flow emerges, with about the global Hubble constant. We also consider the critical epoch $t_Q$, when the DE antigravity began to dominate over the local matter gravity for a galaxy which at the present epoch is in the local DE dominated region. Our main result is that the homogeneous dark energy component, revealed by SNIa observations, resolves the old confrontation between the local Hubble flow and local highly non-uniform, fractal matter distribution. It explains why the Hubble law starts on the outskirts of the Local Group, with the same Hubble constant as globally and with a remarkably small velocity dispersion.

This report continues recent Peebles-Turner debate "Is cosmology solved?" and considers the first results for Sandage's program for "Practical cosmology". A review of conceptual problems of modern cosmological models is given, among them: the nature of the space expansion; recession velocities of distant galaxies more than velocity of light; cosmological Friedmann force; continuous creation of gravitating mass in Friedmann's equation; cosmological pressure is not able to produce a work; cosmological gravitational frequency shift; Friedmann-Holtsmark paradox; the problem of the cosmological constant; Einstein's and Mandelbrot's Cosmological Principles; fractality of observed galaxy distribution; Sandage's 21st problem: Hubble - de Vaucouleurs paradox; quantum nature of gravity force.

A retrospective analysis of the field theory of gravitation, describing gravitational field in the same way as other fields of matter in the flat space-time, is done. The field approach could be called "quantum gravidynamics" to distinguish it from the "geometrodynamics" or general relativity. The basic propositions and main conclusions of the field approach are discussed with reference to classical works of Birkhoff, Moshinsky, Thirring, Kalman, Feynman, Weinberg, Deser. In the case of weak fields both "gravidynamics" and "geometrodynamics" give the same predictions for classical relativistic effects. However, in the case of strong field, and taking into account quantum nature of the gravitational interaction, they are profoundly different. Contents of the paper: 1) Introduction; 2) Two ways in gravity theory: 2.1.Hypotheses of Poincar\'e and Einstein, 2.2. Gravity as a geometry of space, 2.3. Gravitation as a material field in flat space-time; 3) Classical theory of tensor field: 3.1.Works of Birkhoff and Moshinsky, 3.2.Works of Thirring and Kalman, 3.3.Thirring and Deser about identity of GR and FTG; 4) Quantum theory of tensor field; 5) Modern problems in field theory of gravitation: 5.1.Multicomponent nature of tensor field, 5.2.Choice of energy-momentum tensor of gravitational field, 5.3.Absence of black holes in FTG, 5.4.Astrophysical tests of FTG; 6) Conclusions.

It is shown that detecting or setting an upper limit on the scalar gravitational radiation is a good experimental test of relativistic gravity theories. The relativistic tensor-field theory of gravitation is revised and it is demonstrated that the scalar monopole gravitational radiation must be added to the usual quadrupole radiation. In the case of the binary pulsar PSR 1913+16 it is predicted the existence 0.735 % excess of the gravitational radiation due to the scalar gravitational waves.

The KLUN (Kinematics of the Local UNiverse) sample of 6600 spiral galaxies is used in studying dark matter in different scales: Type dependence of the zero-point of Tully-Fisher relation indicates M/L \approx 9 -- 16 in galactic scale. Preliminary results from a study of selection effect influencing double galaxies give a larger value M/L \approx 30 -- 50. Study of the Perseus-Pisces supercluster, using Malmquist bias corrected TF distances and Tolman-Bondi solutions indicates M/L \approx 200 -- 600 for the large clusters. Application of developed version of Sandage-Tammann-Hardy test of the linearity of Hubble law inside the observed hierarchical (fractal) galaxy distribution up to 200 Mpc suggests that either \Omega_0 is very small (0.01) or major part of the matter is uniformly distributed dark matter.

Late-time observations of GRB970508 with the SAO-RAS 6-m telescope in July-August 1998 show that the optical transient (OT) varies very little since November 1997. Here we report the $BVR_cI_c$ light-curves behaviour up to 470 days from the GRB occurence with independent $BVR_cI_c$ power-law fits. After ~=200 days the power-law decay has ceased and the OT contamination to the host galaxy flux is already less than the observational errors of the brightness estimates. Spectral energy distribution for the host galaxy alone (without OT) and comparison of its broad-band spectrum to the average continuum spectra of different Hubble types galaxies are shown. The best fits to the BVR_cI_c spectral energy distribution of the GRB host with local starburst galaxies are found for the blue compact S0 dwarf galaxy NGC1510, the Scd starburst HII galaxy NGC7793 and the blue compact galaxy Mrk 1267. Position of the host galaxy in the M_B vs. log(D_25) diagram for local late-type LEDA galaxies allows us to attribute GRB970508 host galaxy to the blue compact galaxies.

We present direct imaging data of 22 ultra steep spectrum radio sources obtained at (or near) a subarcsecond seeing. The basic sample of 40 double radio sources was selected from the RATAN-600 catalogue. The FRII-structure has been confirmed with VLA and preliminary optical identifications which come from the 6 m-telescope. As the RATAN-600 flux limit at 3.9 GHz ($\approx$ 10mJy) is fainter than that of major surveys, the sample may have high-$z$ contents. This is also suggested by the faint magnitudes in the Hubble diagram. The final aim is to create a homogeneous sample of high-$z$ radio galaxies in a well defined strip around the sky, with faint radio limit and subarcsecond morphology down to $m_{R}=24$. We could confirm 16 identifications down to $m_{R}\sim$24. Most of the extended objects have multicomponent structures as expected from other surveys of high-redshift radio galaxies. We found five unresolved objects even with a subarcsecond seeing. Of the remaining six objects, three are extremely faint and the other three have such a complex environment that further observations are needed to confirm the optical identification.

Optical observations of the GRB 970508 optical remnant were continued with the 6-m telescope of SAO RAS in standard BVRcIc bands in Oct.-Dec.1997 and in Jan. 1998. The results of the photometry of GRB970508 remnant and of three nearby galaxies are presented. The BVRcIc light curves of the GRB970508 remnant may be described by a power law plus a constant F(t)=F_0(t-t_0)^{\alpha}+C. In In parameters determination of the faint host galaxy we used the results of our BVRcIc photometry of May-August, 1997, the data of recent observations Keck-II and WHT telescopes and also the data of Rc and B bands photometry obtained by other authors based on our secondary standards. The level-off from the initial power-law decline seen in the first months after the burst observed in all bands. The effect is the most strong in Ic band. The best \chi^2-fits for Fo, \alpha, C parameters of the data in each of 4 bands point to the presence of a faint constant source with Ic=24.13\pm0.28, Rc = 25.55\pm0.19, V=25.80\pm 0.14, B=26.68\pm 0.14. The average \alpha of light curve power law is -1.23\pm0.04. The optical remnant has the power law F_\nu \nu^\beta spectrum with spectral slope = -1.10\pm0.08 and does not change after optical curve maximum. BVRcIc spectrum together with an absolute magnitude of the constant component M_B_rest=-17.5\pm0.3 and the limit of the linear size of underlying host galaxy d<1.1 kpc, calculated for z=0.835,H_0=75km s^-1 Mpc^-1, q_0=0.5, correspond to a host galaxy, such as a starburst dwarf, red starburst dwarf, irregular dwarf, HII galaxy, or blue compact dwarf galaxy. All these types of dwarf galaxies show evidence of starburst or post starburst activity. Galaxy G2 has a spectrum similar to one of the host GRB galaxy and lies at the projected distance of ~20 kpc from GRB.

Two fundamental empirical laws have been established in the analysis of galaxy space distribution. First, recent analyses have revealed that the three dimensional distribution of galaxies and clusters is characterized by large scale structures and huge voids: such a distribution shows fractal correlations up to the limits of the available samples. This has confirmed the earlier de Vaucouleurs power-law density - distance relation, now corresponding to a fractal structure with dimension $D \approx 2$, at least, in the range of scales $ \sim 1 \div 200 Mpc$ ($H_0 = 55 km/sec/Mpc$). An eventual cut-off towards homogenization has not been yet identified. Second, since Hubble's discovery, the linear redshift-distance law has been well established within $200 Mpc$ and also much deeper. The co-existence of these laws within the same scales is a challenge for the standard cosmology where the linear Hubble law is a strict consequence of homogeneity of the expanding universe. This puzzle is now sufficiently strong to raise doubts for the standard cosmology.

We have studied using the KLUN sample of 5171 spiral galaxies having Tully-Fisher distance moduli, the average radial space distribution of galaxies out to a distance of about 200 Mpc (for H_0=50 km/s/Mpc). One motivation came from the current debate on the fractal dimension D and maximum scale of fractality. To study this question, we used a new method based on photometric TF distances, independent of redshift, to construct the number density distribution. Our main results are: (1) While scattered below about 20 Mpc, at larger distances the radial distribution starts to follow, in terms of distance modulus mu_TF, the law log N = 0.46 mu + const., using diameter TF relation, and log N = 0.40 mu + const. for magnitudes. These are the predictions based on fractal dimensions 2.3 and 2.0, respectively. These radial density gradients are valid up to the limits of KLUN, or about 200 Mpc. (2) We have tried to understand the derived radial density behaviour as a result of some bias in KLUN or our analysis, however, without success. Numerical simulations have shown that the method itself works. If the density law is caused by incompleteness in the diameter limited KLUN sample, then the incompleteness should start at widely different angular diameters d_25 for different values of rotation parameter log V_M, which would be quite unexpected. (3) If correlation studies favoring long scale fractality (200 Mpc or more) and D ~ 2 are correct, the position of our Galaxy would be close to average in the Universe, with the galaxy density decreasing around us according to the expected law.