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.

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.

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.