We present a full forward-modeled $w\mathrm{CDM}$ analysis of the KiDS-1000 weak lensing maps using graph-convolutional neural networks (GCNN). Utilizing the cosmogrid, a novel massive simulation suite spanning six different cosmological parameters, we generate almost one million tomographic mock surveys on the sphere. Due to the large dataset size and survey area, we perform a spherical analysis while limiting our map resolution to HEALPix $n_{\mathrm{side}}=512$. We marginalize over systematics such as photometric redshift errors, multiplicative calibration and additive shear bias. Furthermore, we use a map-level implementation of the nonlinear intrinsic alignment model along with a novel treatment of baryonic feedback to incorporate additional astrophysical nuisance parameters. We also perform a spherical power spectrum analysis for comparison. The constraints of the cosmological parameters are generated using a likelihood-free inference method called Gaussian process approximate Bayesian computation (GPABC). Finally, we check that our pipeline is robust against choices of the simulation parameters. We find constraints on the degeneracy parameter of $S_8\equiv {\sigma }_8\sqrt{{\mathrm{\Omega }}_M/0.3}=0.78_{-0.06}^{+0.06}$ for our power spectrum analysis and $S_8=0.79_{-0.05}^{+0.05}$ for our GCNN analysis, improving the former by 16%. This is consistent with earlier analyses of the 2-point function, albeit slightly higher. Baryonic corrections generally broaden the constraints on the degeneracy parameter by about 10%. These results offer great prospects for full machine learning based analyses of ongoing and future weak lensing surveys.

• Received 20 January 2022
• Accepted 31 March 2022

DOI:https://doi.org/10.1103/PhysRevD.105.083518