The paper is concerned with existence, multiplicity and asymptotic behavior of (weak) solutions for nonlocal systems involving critical nonlinearities. More precisely, we consider {M([u](s)(2)-mu integral V-R3(x)vertical bar u vertical bar(2)dx) [(-Delta)(s)u-mu V(x)u] - phi vertical bar u vertical bar(2)*(s,)t(-2)u -lambda h(x)vertical bar u vertical bar(p-2)u+vertical bar u vertical bar(2)*(s-2)u in R-3 (-Delta)(t)phi = vertical bar u vertical bar(2)*(s,t) in R-3, where (-Delta)(s) is the fractional Lapalcian, [u](s) is the Gagliardo seminorm of u, M : R-0(+) -> R-0(+) is a continuous function satisfying certain assumptions, V (x) = jxj is the Hardy potential function, 2*(s,t) = (3 + 2t)/(3 - 2s), s, t is an element of (0, 1), lambda, mu are two positive parameters, 1 < p < 2*(s) = 6=(3 - 2s) and h is an element of L-2*(s) =(2(s)* - p (R-3). By using topological methods and the Krasnoleskii's genus theory, we obtain the existence, multiplicity and asymptotic behaviour of solutions for above problem under suitable positive parameters lambda and mu. Moreover, we also consider the existence of nonnegative radial solutions and non-radial sign-changing solutions. The main novelties are that our results involve the possibly degenerate Kirchhoff function and the upper critical exponent in the sense of Hardy-Littlehood-Sobolev inequality. We emphasize that some of the results contained in the paper are also valid for nonlocal Schrodinger-Maxwell systems on Cartan-Hadamard manifolds.