The spreading and subduction of oceanic plates are key parts of the Wilson cycle. Starting with the birth at mid-ocean ridges, oceanic plates are controlled by thermal, rheological, and chemical conditions, followed by plate maturation as it ages and cools while crossing the seafloor, and finishing with the dynamics of plate destruction as it retires at the subduction zone to become a deeper part of Earth's convective system. In order to understand how oceanic plates form, evolve, and are destroyed, I employed 2D and 3D self-consistent magmatic-thermomechanical numerical models with both brittle/plastic strain weakening and grain size evolution to focus on the spatial and temporal evolution of mid-ocean ridges and subduction zones, systemically explore the effects of magmatism, strain-dependent friction coefficient, grain size evolution, and brittle-ductile damage on these processes, and then compare the results to available geological-geophysical observations. For the ocean spreading, the spatial and temporal evolution of the mid-ocean ridge is investigated in two separate directions – along the spreading ridges and along the spreading direction. Along the spreading ridges, two distinct sections, magmatic and amagmatic, are documented. They are caused by the dynamical boundary instability beneath spreading centers. The sustained along-ridge discontinuity starts with the limited melt supply and thick brittle layer, maintains through density and latent heating feedback and different spreading modes in magmatic compared to amagmatic sections, and vanishes with increasing spreading rate and mantle potential temperature and/or decreasing intensity of hydrothermal cooling. Along the spreading direction, a spectrum of tectonic patterns, from asymmetric long-lived detachment faults in rolling-hinge mode to symmetric conjugate faults in flip-flop mode, are observed. Fault strength reduction (induced by hydrothermal alteration and strain-dependent friction coefficient) and axial brittle layer thickness (controlled by spreading rates, hydrothermal circulation, and mantle potential temperature) are two pivotal factors in controlling the faulting patterns and spreading modes along the spreading direction. Strong strain weakening and thick brittle layer facilitate the formation of long-lived detachment faults in rolling-hinge mode, while weak strain weakening and thin brittle layer result in conjugated faults in flip-flop mode at spreading centers. Grain size reduction is observed at the root of detachment faults, but its effect was found negligible owing to the small fault strength reduction induced by grain damage. For oceanic plate subduction, the mechanism of ridge-inversed subduction initiation and the effect of grain size reduction on the lithospheric discontinuities at the ocean-continent subduction zone were explored. The forced near-ridge subduction initiation was investigated by modeling the inversion of three-dimensional inherited mature spreading patterns formed after long-term spreading. Forced compression predominantly reactivates and rotates inherited extensional faults at former spreading centers, shortening and thickening the weakest near-ridge region of the oceanic lithosphere, thereby producing ridge swellings. As a result, a new megathrust zone is developed, which accommodates further shortening and subduction initiation. Brittle/plastic strain weakening has a key impact on the collapse of the thickened ridge and the onset of near-ridge subduction initiation. In contrast, grain size evolution of the mantle only slightly enhances the localization of shear zones at the brittle-ductile transition and thus plays a subordinate role. For the origin of discontinuity layers at the ocean-continent subduction zone, it is related to the grain size reduction layer caused by subduction and formed within the brittle/ductile transition zone. This grain size reduction layer is a potential barrier zone to compact the rising volatiles from the deep mantle. The compacted volatiles below the barrier zone may lead to a sharp seismic velocity drop, which is the possible genesis of the discontinuity layer within the lithosphere. Furthermore, the enriched melts/volatiles at the Lithosphere-Asthenosphere Boundary, induced by small-scale convection, slab bending, or slab dehydration, are critical sources for compacted volatiles. My findings provide a better understanding of the formation, development, and destruction of oceanic plates. I clarify the mechanisms controlling variations in tectonic patterns and oceanic crustal thickness at mid-ocean ridges, elucidate the detailed evolution process of forced near-ridge subduction initiation, and propose a new mechanism for the origin of mid-lithosphere discontinuity. These new findings shed light on geological and seismic studies.