We develop a new and general method to prove the the existence of the random attractor (strong attractor) for the primitive equations (PEs) of large-scale ocean and atmosphere dynamics under non-periodic boundary conditions and driven by infinitedimensional additive fractional Wiener processes. In contrast to our new method, the common method, compact Sobolev embedding theorem, is to obtain the uniform a priori estimates in some Sobolev space whose regularity is high enough. But this is very complicated for the 3D stochastic PEs with the non-periodic boundary conditions. Therefore, the existence of universal attractor ( weak attractor) was established in previous work. The main idea of our method is that we first derive that P-almost surely the solution operator of stochastic PEs is compact. Then we construct a compact absorbing set by virtue of the compact property of the the solution operator and the existence of a absorbing set. We should point out that our method has some advantages over the common method of using compact Sobolev embedding theorem, i.e., if the random attractor in some Sobolev space do exist in view of the common method, our method would then further implies the existence of random attractor in this space. The present work provides a general way for proving the existence of random attractor for common classes of dissipative stochastic partial differential equations and improves the existing results concerning random attractor of stochastic PEs. In a forth coming paper, we use this new method to prove the existence of strong attractor for the stochastic moist primitive equations, improving the results, the existence of weak (universal) attractor of the deterministic model. M × (−h, 0) ⊂ R 3 , and consider the following 3D stochastic PEs of Geophysical Fluid Dynamics.The unknowns for the 3D stochastic viscous PEs are the fluid velocity field (υ, w) = (υ 1 , υ 2 , w) ∈ R 3 with υ = (υ 1 , υ 2 ) and υ ⊥ = (−υ 2 , υ 1 ) being horizontal, the temperature T and the pressure p. f = f 0 (β + y) is the given Coriolis parameter, Q is a given heat source. The viscosity and the heat diffusion operators L 1 and L 2 are given byHere the positive constants ν 1 , µ 1 are the horizontal and vertical Reynolds numbers, respectively, and ν 2 , µ 2 are positive constants which stand for the horizontal and vertical heat diffusivity, respectively. To simplify the nations, we assume ν i = µ i = 1, i = 1, 2. The results in this paper are still valid when we consider the general cases. We set ∇ = (∂ x , ∂ y ) to be the horizontal gradient operator and ∆ = ∂ 2x + ∂ 2 y to be the horizontal Laplacian. Here, we takeẆ H i (t, x, y, z), i = 1, 2, the informal derivative for the fractional Wiener process W H i given below. The boundary of ℧ is partitioned into three parts:Here h is a sufficiently smooth function. Without loss generality, we assume h = 1. We consider the following boundary conditions of the stochastic 3D viscous PEs.where η(x, y) is the wind stress on the surface of the ocean, α is a positive constant, τ is t...