In the realm of astrophysics, black holes exist within nonvacuum cosmological backgrounds, making it crucial to investigate how these backgrounds influence the properties of black holes. In this work, we first introduce a novel static spherically-symmetric exact solution of Einstein field equations representing a surrounded hairy black hole. This solution represents a generalization of the hairy Schwarzschild solution recently derived using the extended gravitational decoupling method. Then, we discuss how the new induced modification terms attributed to the primary hairs and various background fields affect the geodesic motion in comparison to the conventional Schwarzschild case. Although these modifications may appear insignificant in most cases, we identify specific conditions where they can be comparable to the Schwarzschild case for some particular background fields.

We investigate here the final state of gravitational collapse of a non-spherical and non-marginally bound dust cloud as modeled by the Szekeres spacetime. We show that a directionally globally naked singularity can be formed in this case near the collapsing cloud boundary and not at its geometric center, as is typically the case for a spherical gravitational collapse. This singularity is a strong curvature naked singularity in the sense of Tipler criterion on gravitational strength. The null geodesics escaping from the singularity would be less scattered in this case in certain directions since the singularity is close to the boundary of the cloud, as is the case in the current scenario. The physical implications are pointed out.

In this paper, we apply the gravitational decoupling method for dynamical systems in order to obtain a new type of solution that can describe a hairy dynamical black hole. We consider three cases of decoupling. The first one is the simplest and most well known when the mass function is the function only of space coordinate r. The second case is a Vaidya spacetime case when the mass function depends on time v. Finally, the third case represents the generalization of these two cases: the mass function is the function of both r and v. We also calculate the apparent horizon and singularity locations for all three cases.