In the early 1900s, Dickson introduced what he called generalized quaternion algebras over any field $K$ of characteristic not 2. These are exactly what we'd call quaternion algebras over $K$. His definition was in terms of a basis with rules for their products, and he gave a criterion for these to be division rings. In particular, these were the first noncommutative division rings besides the quaternions of Hamilton, aside from subrings of Hamilton's quaternions.
Three of Dickson's works where he introduces these algebras are
(1) Linear Algebras, Trans. AMS ${\bf 13}$ (1912) 59-73.
(2) Linear Associative Algebras and Abelian Equations, Trans. AMS ${\bf 15}$ (1914), 31-46.
(3) Algebras and Their Arithmetics, Univ. of Chicago Press, 1923.
In (1) he gives the defining equations for a generalized quaternion algebra and the norm criterion for it to be a division algebra (pp. 65-66), though without using the label "generalized quaternion algebra." He writes in a footnote that this work goes back to 1906.
In (2) he constructs cyclic algebras, without using that name, and calls the special case of dimension 4 a generalized quaternion algebra. (Archaic terminology alert: Dickson refers to equations defining cyclic Galois extensions as uniserial abelian equations.)
In (3) he defines cyclic algebras again without using that name (p. 65), writes them as $D$, refers to them as algebras of type $D$ (p. 68), and remarks in a footnote on p. 66 that Wedderburn calls them Dickson algebras. Near the end of the book he looks at generalized quaternion algebras over the rationals with $\mathbf Q(i)$ as a maximal subfield, and as a particular example he uses the Hamilton quaternions over $\mathbf Q$ to describe all rational and integral solutions of certain quadratic Diophantine equations in several variables. Think about how a sum of four squares factors over the quaternions to imagine how arithmetic properties of quaternions could be useful to analyze a Diophantine equation involving a sum of four squares; this is similar in spirit to the way the Gaussian integers are useful in studying a Diophantine equation involving a sum of two squares. The term "arithmetics" in the title of the book is, as far as I can tell, Dickson's label for what we'd call maximal orders, so the book would be called today "Algebras and Their Maximal Orders."
The motivation for Dickson's interest in quaternion algebras was the earlier development of integral Hamilton quaternions, due first to Lipschitz (all integral coefficients, which is clunky in the same way that $\mathbf Z[\sqrt{-3}]$ is compared to $\mathbf Z[(1+\sqrt{-3})/2]$) and then to Hurwitz (all integral or all half-integral coefficients). To study quadratic Diophantine equations in several variables going beyond a sum of four squares, Dickson was led to extend the original definition of quaternions. His main interest was developing the right theory of generalized integral quaternions, rather than just a theory over a field. A further reference in this direction is Dickson's paper
(4) On the theory of numbers and generalized quaternions, Amer. J. Math ${\bf 46}$ (1924), 1-16.
Concerning the connection to central simple algebras, Dickson was the first to show any division algebra that is 4-dimensional over its center is cyclic, at least outside characteristic 2 since he didn't have a good definition in characteristic 2. Taking into account Wedderburn's theorem that every finite-dimensional central simple algebra over a field is a matrix algebra over a division algebra, and that matrix algebras are cyclic, Dickson had shown that every 4-dimensional CSA over a field not having characteristic 2 is an algebra of "his" type. (Wedderburn proved the analogous theorem for dimension 9.)