It's all summing dx/dt + dy/dt + dz/dt, for i paths between satellites and ground stations (or more receivers for differential or rtk or vrs style). [2]
Which reduces most of the time to summing DELTA-Xi + DELTA-Yi + DELTA-Zi + delta-t(timeerrors). For i paths between each sat and ground receiver.
Which you should recognize the transformation if you've ever taken calculus. Even if you don't integrate every time you get a fix.
Part of what I describe as math 'magic' is that you can cancel out most of the unknowns and most of the unsolved calculus if you add a second fixed receiver.
Google and Apple location services 'cheat' and do this via subbing a nearby wifi MAC with known coordinates, which for them is good-enough. But augmented gps from FAA or DOT or coastguard etc work the same way, but with real gps receivers on the ground in realtime. Obviously without having to substitute anything.
Either way- the extra known variable greatly simplifies the math via canceling-out terms.
Plus there are both closed and open form solutions developed since initial GPS deployment that allow solving without direct integration.
Chapter 12 of [0] Surveying gets into the math, including transformations, if you want to see the math details.
Or [1] GPS by van Sickle for a good overview of the various methods/ technologies. (Also survey-centric).
[0]https://books.google.com/books/about/Surveying_theory_and_pr...
[1]https://books.google.com/books?id=J0fLBQAAQBAJ&pg=PA63&sourc...
[2] despite wgs84 and lat/lon being associated as default 'GPS coordinates', the 'raw' gps system data is xyz Cartesian in feet, then transformed to lat lon or whatever else.